孝顺竹笋箨全长转录组测序分析

doi:10.3969/j.issn.1000-2006.202002025

国家林草科技领军期刊
中国精品科技期刊
中国高校百佳科技期刊
江苏省新闻出版政府奖期刊奖
RCCSE林学权威期刊（A+）
CSCD核心期刊
Scopus数据库收录期刊
中文核心期刊
SCD核心期刊

作者加群：102861116

微信公众号：南京林业大学学报

高级检索

南京林业大学学报（自然科学版） ›› 2020, Vol. 44 ›› Issue (6): 175-183.doi: 10.3969/j.issn.1000-2006.202002025

孝顺竹笋箨全长转录组测序分析

鞠烨¹^,²(), 江建平³, 尹增芳¹^,², 魏强¹^,^*()

1.南京林业大学,南方现代林业协同创新中心,江苏南京 210037
2.南京林业大学生物与环境学院,江苏南京 210037
3.万载县林业局,江西万载 336100

收稿日期:2020-02-15 修回日期:2020-03-22 出版日期:2020-11-30 发布日期:2020-12-07
通讯作者: 魏强
基金资助:
国家自然科学基金项目(31670602);江苏高校“青蓝工程”资助项目;江苏高校优势学科建设工程资助项目(PAPD)

Full-length transcriptome sequencing and annotation analyses of Bambusa multiplex sheath

JU Ye¹^,²(), JIANG Jianping³, YIN Zengfang¹^,², WEI Qiang¹^,^*()

1. Co-Innovation Center for the Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
2. College of Biology and the Environment, Nanjing Forestry University, Nanjing 210037, China
3. Wanzai Forestry Bureau, Wanzai 336100, China

Received:2020-02-15 Revised:2020-03-22 Online:2020-11-30 Published:2020-12-07
Contact: WEI Qiang

摘要/Abstract

摘要：

【目的】揭示竹子笋箨的生物学功能与其衰老的分子基础。【方法】利用PacBio Sequel 3代全长转录组测序技术结合生物信息学方法,对孝顺竹不同衰老阶段笋箨全长转录本进行分析。【结果】共获得106 148 条平均长度为3 615 bp的全长转录本。多数转录本长度分布在1.4~7.0 kb。NCBI等注释显示,97.34%的转录本具有注释结果。Mercator注释分析显示笋箨转录本覆盖了其全部34 个类别,其中与蛋白类别相关的序列最多,达到了18 224 条;与microRNA类别相关的最少,仅有9 条。在重要功能基因方面,共获得了2 489 条激素代谢及信号转导相关序列; 8 804 条转录因子序列中393 条与光合作用相关的转录本。其中,注释到的189 条共30 类NAC(NAM/ATAF/CUC)转录因子基因,93.33%的NAC基因,共计28类,其表达量与笋箨衰老呈相关性, 包括已被证实在叶片衰老中具有重要调控作用的NAC002、NAC016、NAC017、NAC029(NAP)、NAC042、NAC055与NAC083等7个NAC转录因子与NAC014等14个被报道在叶片衰老中具有潜在作用的NAC基因。NAC025、NAC028、NAC045、NAC061、NAC086、NAC103与NAC1L (NAC 1 Like) 7 个NAC转录因子表达与孝顺竹笋箨衰老呈正相关,为新发现的正调控笋箨衰老的潜在转录因子。此外,利用MISA程序在76 499 个序列中检测到简单序列重复(simple sequence repeats,SSR)位点,主要以单核苷酸重复(SSRs)为主,占据了全部(SSRs)的55.2%。利用PLEK等软件共获得2 769个长链非编码RNA(long non-coding RNA, lncRNA)。【结论】竹子笋箨基因表达多样化,并具有完整的光合系统基因,显示其具有潜在光合功能;NAC转录因子在孝顺竹笋箨衰老中具有潜在重要调控作用。本研究首次解析了竹子笋箨的全长转录本特征,为今后深入分析竹子笋箨功能及其衰老的分子机制奠定基础。

关键词: 孝顺竹, PacBio, 简单序列重复(SSR), lncRNA, NAC, 衰老

Abstract:

【Objective】Currently, few reports have described the biological function of the bamboo sheath. This work aimed to elucidate the possible functions and the underlying biological mechanisms of the bamboo sheath from the perspective of molecular biology.【Method】Using PacBio Sequel full-length sequencing technology, we analyzed the transcriptomic profile of the culm sheath at different senescence stages of Bambusa multiplex.【Result】A total of 106 148 transcripts with an average length of 3 615 base pairs (bp) were obtained. Of those transcripts, 97.34% had annotated results. Mercator annotation further revealed that the transcripts covered all 34 functional subgroups associated with the bamboo sheath. The protein subgroup contained the most transcripts, with 18 224 transcripts. While the microRNA subgroup contained the least transcripts, with only nine being present in this subgroup. Additionally, 2 489 hormone-related transcripts and 8 804 transcription factors were detected. Interestingly, 393 transcripts related to photosynthesis were discovered in the full-length transcriptome of the B. multiplex sheath. A total of 189 transcripts encoding 30 NAC transcription factor genes were also identified. The expression of 93.33% of NAC transcription factor genes was found to be positively or negatively correlated with sheath senescence. Of them, seven NACs, including NAC002, NAC016, NAC017, NAC029(NAP), NAC042, NAC055 and NAC083 have been reported to play important roles in regulating leaf senescence; 14 NACs, including NAC014, have been reported to have potential roles in leaf senescence; and seven bamboo sheath senescence-associated NACs, including NAC025, NAC028, NAC045, NAC061, NAC086, NAC103 and NAC1L were newly discovered. The MISA software analysis revealed that 76 499 transcripts contain simple sequence repeats (SSRs). Most are single nucleotide repeat SSRs, accounting for 55.2% of all discovered SSRs. Additionally, 2 769 long non-coding RNAs (lncRNAs) were detected by using software such as PLRK et al.【Conclusion】These results revealed that various genes are actively expressed in the bamboo sheath. A number of transcripts related to all aspects of photosynthesis were detected in the bamboo sheath, indicating that the bamboo sheath may have the photosynthetic capacity. The majority of the identified NACs was associated with leaf senescence, suggesting a potential role of NACs in bamboo sheath senescence. Taken together, our work provides the first insights into the transcriptomic features of the bamboo sheath and provides a foundation for future studies to reveal the biological function of the bamboo sheath and the underlying molecular mechanisms of bamboo sheath senescence.

Key words: Bambusa multiplex, PacBio, simple sequence repeats(SSR), lncRNA, NAC, senescence

中图分类号:

S795
Q781

鞠烨,江建平,尹增芳,等. 孝顺竹笋箨全长转录组测序分析[J]. 南京林业大学学报（自然科学版）, 2020, 44(6): 175-183.

JU Ye, JIANG Jianping, YIN Zengfang, WEI Qiang. Full-length transcriptome sequencing and annotation analyses of Bambusa multiplex sheath[J].Journal of Nanjing Forestry University (Natural Science Edition), 2020, 44(6): 175-183.DOI: 10.3969/j.issn.1000-2006.202002025.

图/表 7

图1

图2

表1

光合作用相关转录本"

光合作用过程 photosynthesis part	MapMan类别 MapMan Bin	转录本数目 transcript number	光合作用过程 photosynthesis part	MapMan类别 MapMan Bin	转录本数目 transcript number
卡尔文循环 calvin cycle	醛缩酶aldolase	11	光反应 light reaction	转酮醇酶transketolase	3
卡尔文循环 calvin cycle	醛缩酶aldolase	11		ATP合成酶ATPsynthase	12
	D-果糖1,6-二磷酸酯酶 D-fructose-1,6-bisphosphate 1-phosphohydrolase	2		ATP合成酶ATPsynthase	12
				循环电子流-叶绿体呼吸 cyclic electron flow-chlororespiration	26
				细胞色素b6/f cytochrome b6/f	4
	甘油醛-3-磷酸脱氢酶 glyceraldehyde-3-phosphate dehydrogenase	1		细胞色素b6/f cytochrome b6/f	4
				NADH脱氢酶NADH DH	12
				铁氧化还原蛋白ferredoxin	25
				质体醌plastocyanin	1
				光系统Ⅰ. LHC-Ⅰ photosystemI. LHC-Ⅰ	12
	磷酸甘油酸激酶 phosphoglycerate kinase	10		光系统Ⅰ. LHC-Ⅰ photosystemI. LHC-Ⅰ	12
	磷酸甘油酸激酶 phosphoglycerate kinase	10		光系统Ⅰ多肽亚基 photosystem Ⅰ. PSI polypeptide subunits	28
	磷酸戊糖激酶 phosphopentokinase	4		光系统Ⅰ多肽亚基 photosystem Ⅰ. PSI polypeptide subunits	28
	磷酸戊糖激酶 phosphopentokinase	4		光系统Ⅱ.生物发生 photosystem Ⅱ. Biogenesis	1
	核糖-5-磷酸异构酶 ribose 5-phosphate isomerase	1		光系统Ⅱ.生物发生 photosystem Ⅱ. Biogenesis	1
	核糖-5-磷酸异构酶 ribose 5-phosphate isomerase	1		光系统Ⅱ. LHC-Ⅱ photosystem Ⅱ. LHC-Ⅱ	20
	核酮糖磷酸3-表异构酶 ribulose-phosphate 3-epimerase	7		光系统Ⅱ. LHC-Ⅱ photosystem Ⅱ. LHC-Ⅱ	20
				光系统Ⅱ多肽亚基 photosystem Ⅱ. PSII polypeptide subunits	42
				状态转换state transition	24
	rubisco互作相关 rubisco interacting	30		状态转换state transition	24
	rubisco互作相关 rubisco interacting	30	光呼吸 photorespiration	过氧化物酶体转氨酶 aminotransferases peroxisomal	7
	rubisco大亚基 rubisco large subunit	6	光呼吸 photorespiration	过氧化物酶体转氨酶 aminotransferases peroxisomal	7
				甘油酸激酶glycerate kinase	9
				甘油酸清除. H蛋白 glycine cleavage. H protein	12
	rubisco小亚基 rubisco small subunit	4		甘油酸清除. H蛋白 glycine cleavage. H protein	12
				乙醇酸氧化酶glycolateoxidase	31
				羟基丙酮酸还原酶 hydroxy pyruvate reductase	17
	景天庚酮糖二磷酸酶 sedoheptulose bisphosphatase	2		羟基丙酮酸还原酶 hydroxy pyruvate reductase	17
	景天庚酮糖二磷酸酶 sedoheptulose bisphosphatase	2		磷酸甘油磷酸酯酶 phosphoglycolate phosphatase	5
	磷酸丙糖异构酶 triosephosphate isomerase	11		磷酸甘油磷酸酯酶 phosphoglycolate phosphatase	5
	磷酸丙糖异构酶 triosephosphate isomerase	11		丝氨酸羟甲基转移酶 serinehydroxy methyltransferase	1

表1

图3

图4

表2

图5

参考文献 38

[1]	WEI Q, JIAO C, GUO L, et al. Exploring key cellular processes and candidate genes regulatingthe primary thickening growth of Moso underground shoots[J]. New Phytologist, 2017,214:81-96. DOI: 10.1111/nph.14284. doi: 10.1111/nph.14284 pmid: 27859288
[2]	WEI Q, JIAO C, DING Y L, et al. Cellular and molecular characterizations of a slow-growth variant provide insights into the fast growth of bamboo[J]. Tree Physiology, 2018,38(4):641-654. DOI: 10.1093/treephys/ tpx129. pmid: 29077967
[3]	WEI Q, GUO L, JIAO C, et al. Characterization of the developmental dynamics of the elongation of a bamboo internode during the fast growth stage[J]. Tree Physiology, 2019,39(7):1201-1214. DOI: 10.1093/treephys/tpz063. pmid: 31135922
[4]	GUO L, SUN X, LI Z, et al. Morphological dissection and cellular and transcriptome characterizations of bamboo pith cavity formation reveal a pivotal role of genes related to programmed cell death[J]. Plant Biotechnology Journal, 2019,17(5):982-997. DOI: 10.1111/pbi.13033.
[5]	魏强, 丁雨龙. 矢竹地下茎转录组测序及节间生长相关基因表达分析[J]. 南京林业大学学报(自然科学版), 2017,41(5):42-48.
	WEI Q, DING Y L. Transcriptome sequencing, de novo assembly and expression analysis of several genes related to internode development in Pseudosasa japonica[J]. J Nanjing For Univ (Nat Sci Ed), 2017,41(5):42-48. DOI: 10.3969/j.issn.1000-2006.201601006.
[6]	JIAO Y, HU Q, ZHU Y, et al. Comparative transcriptomic analysis of the flower induction and development of the Lei bamboo (Phyllostachys violascens)[J]. BMC bioinformatics, 2019,20(25):1-12. DOI: 10.1186/s12859-019-3261-z.
[7]	GAMUYAO R, NAGAI K, AYANO M, et al. Hormone distribution and transcriptome profiles in bamboo shoots provide insights on bamboo stem emergence and growth[J]. Plant and Cell Physiology, 2017,58(4):702-716. DOI: 10.1093/pcp/pcx023.
[8]	LIU Y, WU C, HU X, et al. Transcriptome profiling reveals the crucial biological pathways involved in cold response in Moso bamboo (Phyllostachys edulis)[J]. Tree Physiology, 2019. DOI: 10.1093/treephys/tpz133.
[9]	LI X, XIE L, ZHENG H, et al. Transcriptome profiling of postharvest shoots identifies PheNAP2-and PheNAP3-promoted shoot senescence[J]. Tree Physiology, 2019,39(12):2027-2044. DOI: 10.1093/treephys/tpz100.
[10]	柳延虎, 王璐, 于黎. 单分子实时测序技术的原理与应用[J]. 遗传, 2015,37(3):259-268.
	LIU Y H, WANG L, YU L. The principle and application of the single-molecule real-time sequencing technology[J]. HEREDITAS, 2015,37(3):259-268. DOI: 10.16288/j.yczz.14-323.
[11]	ZHAO H, GAO Z, WANG L, et al. Chromosome-level reference genome and alternative splicing atlas of moso bamboo (Phyllostachys edulis) [J]. Giga Science, 2018, 7(10): giy115. DOI: 10.1093/gigascience/giy115. doi: 10.1093/gigascience/giy062 pmid: 29860514
[12]	WANG T, WANG H, CAI D, et al. Comprehensive profiling of rhizome-associated alternative splicing and alternative polyadenylation in moso bamboo (Phyllostachys edulis)[J]. The Plant Journal, 2017,91(4):684-699. DOI: 10.1111/tpj.13597. doi: 10.1111/tpj.13597 pmid: 28493303
[13]	FU L, NIU B, ZHU Z, et al. CD-HIT: accelerated for clustering the next-generation sequencing data[J]. Bioinformatics, 2012,28(23):3150-3152. DOI: 10.1093/bioinformatics/bts565.
[14]	THIMM O, BLäSING O, GIBON Y, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes[J]. The Plant Journal, 2004,37(6):914-939. DOI: 10.1111/j.1365-313X.2004.02016.x. doi: 10.1111/j.1365-313x.2004.02016.x pmid: 14996223
[15]	LOHSE M, NAGEL A, HERTER T, et al. Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data[J]. Plant Cell Environment, 2014,37(5):1250-1258. DOI: 10.1111/pce.12231.
[16]	BEIER S, THIEL T, MÜNCH T, et al. MISA-web: a web server for microsatellite prediction[J]. Bioinformatics, 2017,33(16):2583-2585. DOI: 10.1093/bioinformatics/btx198. doi: 10.1093/bioinformatics/btx198 pmid: 28398459
[17]	SUN L, LUO H, BU D, et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts[J]. Nucleic Acids Research, 2013, 41(17): e166-. DOI: 10.1093/nar/gkt646.
[18]	KONG L, ZHANG Y, YE Z Q, et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine[J]. Nucleic Acids Research, 2007, 35(S2): 345-349. 10.1093/nar/gkm391.
[19]	FINN R D, COGGILL P, EBERHARDT R Y, et al. The Pfam protein families database: towards a more sustainable future[J]. Nucleic Acids Research, 2016,44(1):279-285. DOI: 10.1093/nar/gkv1344.
[20]	LI A, ZHANG J, ZHOU Z. PLEK: a tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme[J]. BMC Bioinformatics, 2014,15(1):311. DOI: 10.1186/1471-2105-15-311.
[21]	KIM H J, NAM H G, LIM P O. Regulatory network of NAC transcription factors in leaf senescence[J]. Current Opinion in Plant Biology, 2016,33:48-56. DOI: 10.1016/j.pbi.2016.06. 002. doi: 10.1016/j.pbi.2016.06.002 pmid: 27314623
[22]	LIU M, QIAO G, JIANG J, et al. Transcriptome sequencing and de novo analysis for Ma bamboo (Dendrocalamus latiflorus Munro) using the Illumina platform[J]. PLoS ONE, 2012,7(10):e46766. DOI: 10.1371/journal.pone.0046766. doi: 10.1371/journal.pone.0046766
[23]	CUI K, WANG H, LIAO S, et al. Transcriptome sequencing and analysis for culm elongation of the world’s largest bamboo (Dendrocalamus sinicus)[J]. PLoS ONE, 2016,11(6):e0157362. DOI: 10.1371/ journal.pone.0157362. doi: 10.1371/journal.pone.0157362 pmid: 27304219
[24]	王身昌, 胡尚连, 曹颖, 等. 梁山慈竹高通量转录组测序及差异表达基因分析[J]. 华北农学报, 2016,31(3):65-71. doi: 10.7668/hbnxb.2016.03.010
	WANG S C, HU S L, CAO Y, et al. High-throughput RNA-seq and analysis on differential expressed gene from Dendrocalamus farinosus[J]. Acta Agriculturae Boreali-Sinica, 2016,31(3):65-71. DOI: 10.7668 /hbnxb.2016.03.010.
[25]	THIEL T, MICHALEK W, VARSHNEY R, et al. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.)[J]. Theoretical and Applied Genetics, 2003,106(3):411-422. DOI: 10.1007/s00122-002-1031-0. doi: 10.1007/s00122-002-1031-0 pmid: 12589540
[26]	TRANBARGER T J, KLUABMONGKOL W, SANGSRAKRU D, et al. SSR markers in transcripts of genes linked to post-transcriptional and transcriptional regulatory functions during vegetative and reproductive development of Elaeis guineensis[J]. BMC Plant Biology, 2012,12(1):1. DOI: 10.1186/1471-2229-12-1.
[27]	TAHERI S, LEE ABDULLAH T, YUSOP M R, et al. Mining and development of novel SSR markers using next generation sequencing (NGS) data in plants[J]. Molecules, 2018,23(2):399.DOI: 10.3390/ molecules 23020399. doi: 10.3390/molecules23020399
[28]	CHEKANOVA J A. Long non-coding RNAs and their functions in plants[J]. Current Opinion in Plant Biology, 2015,27:207-216. DOI: 10.1016/j.pbi.2015.08.003. doi: 10.1016/j.pbi.2015.08.003 pmid: 26342908
[29]	MENG X, ZHANG P, CHEN Q, et al. Identification and characterization of ncRNA-associated ceRNA networks in Arabidopsis leaf development[J]. BMC genomics, 2018,19(1):607. DOI: 10.1186/s12864-018-4993-2. doi: 10.1186/s12864-018-4993-2 pmid: 30103673
[30]	郭兆武, 李合松, 王若仲, 等. 高产杂交稻 ‘两优培九’剑叶及其叶鞘的光合作用[J]. 植物生理与分子生物学学报, 2007,33(6):531-537.
	GUO Z W, LI H S, WANG R Z, et al. Photosynjournal of the flag leaf blade and its sheath in high-yielding hybrid rice ‘Liangyoupeijiu’[J]. Journal of Plant Physiology and Molecular Biology, 2007,33(6):531-537.
[31]	SAKURABA Y, JEONG J, KANG M Y, et al. Phytochrome-interacting transcription factors PIF4 and PIF5 induce leaf senescence in Arabidopsis[J]. Nature Communications, 2014,5(1):1-13. DOI: 10.1038/ncomms5636.
[32]	KIM Y S, SAKURABA Y, HAN S H, et al. Mutation of the Arabidopsis NAC016 transcription factor delays leaf senescence[J]. Plant and Cell Physiology, 2013,54(10):1660-1672. DOI: 10.1093/pcp/pct113. doi: 10.1093/pcp/pct113 pmid: 23926065
[33]	LIANG M, LI H, ZHOU F, et al. Subcellular distribution of NTL transcription factors in Arabidopsis thaliana[J]. Traffic, 2015,16(10):1062-1074. DOI: 10.1111/tra.12311. doi: 10.1111/tra.12311 pmid: 26201836
[34]	NG S, IVANOVA A, DUNCAN O, et al. A membrane-bound NAC transcription factor, ANAC017, mediates mitochondrial retrograde signaling in Arabidopsis[J]. The Plant Cell, 2013,25(9):3450-3471. DOI: 10.1105/tpc.113.113985. doi: 10.1105/tpc.113.113985 pmid: 24045017
[35]	GUO Y, GAN S. AtNAP, a NAC family transcription factor, has an important role in leaf senescence[J]. The Plant Journal, 2006,46(4):601-612. DOI: 10.1111/j.1365-313X.2006.02723.x. doi: 10.1111/j.1365-313X.2006.02723.x pmid: 16640597
[36]	WU A, ALLU A D, GARAPATI P, et al. JUNGBRUNNEN1, a reactive oxygen species-responsive NAC transcription factor, regulates longevity in Arabidopsis[J]. The Plant Cell, 2012,24(2):482-506. DOI: 10.1105/tpc.111.090894. doi: 10.1105/tpc.111.090894 pmid: 22345491
[37]	HICKMAN R, HILL C, PENFOLD C A, et al. A local regulatory network around three NAC transcription factors in stress responses and senescence in Arabidopsis leaves[J]. The Plant Journal, 2013,75(1):26-39. DOI: 10.1111/tpj.12194. pmid: 23578292
[38]	YANG S D, SEO P J, YOON H K, et al. The Arabidopsis NAC transcription factor VNI2 integrates abscisic acid signals into leaf senescence via the COR/RD genes[J]. The Plant Cell, 2011,23(6):2155-2168. DOI: 10.1105/tpc.111.084913. doi: 10.1105/tpc.111.084913 pmid: 21673078

重复次数 replicates	重复基元长度及其数目 length of repeat motif and number
重复次数 replicates	单核苷酸 mononucleotide	二核苷酸 dinucleotide	三核苷酸 trinucleotide	四核苷酸 tetranucleotide	五核苷酸 pentanucleotide	六核苷酸 hexanucleotide
≥ 5~8	0	9 735	16 985	918	337	130
≥9~12	27 294	3 025	6 56	50	6	2
≥13~16	8 658	1 034	132	12	0	0
≥17~20	3 826	716	30	35	0	2
≥21~24	2 450	433	23	10	0	0
总数 total number	42 228	14 943	17 826	1 025	343	134
占比/% percentage	55.20	19.53	23.30	1.34	0.45	0.18

孝顺竹笋箨全长转录组测序分析

Full-length transcriptome sequencing and annotation analyses of Bambusa multiplex sheath

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 7

参考文献 38

相关文章 15

编辑推荐

Metrics

本文评价

作者

读者

审稿

[1]	却枫, 刘庆楠, 查若飞, 魏强. 孝顺竹中笋箨衰老相关WRKY转录因子的鉴定与分析[J]. 南京林业大学学报（自然科学版）, 2023, 47(6): 113-123.
[2]	王晓静, 王涛, 杨凯, 李潞滨. PEG和NaCl胁迫下毛竹萌发种子中环状RNA特征及其表达研究[J]. 南京林业大学学报（自然科学版）, 2023, 47(6): 17-24.
[3]	杨甲定, 刘雨节, 冯建元, 张远兰. 树木叶片衰老中的氮素再吸收机制研究进展[J]. 南京林业大学学报（自然科学版）, 2023, 47(5): 1-8.
[4]	张强, 周鹏, 刘昌来, 余永帆, 张敏, 杨甲定. NaCl处理下全缘冬青和红果冬青根系的转录组活性比较[J]. 南京林业大学学报（自然科学版）, 2022, 46(3): 99-108.
[5]	田雪瑶, 周洁, 王保松, 何开跃, 何旭东. 柳树NAC基因的克隆与表达模式分析[J]. 南京林业大学学报（自然科学版）, 2020, 44(1): 119-124.
[6]	刘晓威,杨秀艳,武海雯,刘肖艳,朱建峰,张华新. NaCl胁迫下红砂种子萌动期差异表达基因的转录组分析[J]. 南京林业大学学报（自然科学版）, 2019, 43(03): 28-36.
[7]	吕东林,林琳,郭译文,韩锐,姜静. 紫雨桦耐盐性及花色苷合成相关基因的表达特性[J]. 南京林业大学学报（自然科学版）, 2018, 42(02): 25-32.
[8]	朱泓,王小敏,黄涛,吴文龙,李维林. NaCl胁迫对滨梅根际细菌群落多样性及优势菌群的影响[J]. 南京林业大学学报（自然科学版）, 2017, 41(04): 49-54.
[9]	王浩然,李爽爽,乐丽娜,匡华琳,黄敏仁,陈英. miR164a及其靶基因PeNAC1相互作用研究[J]. 南京林业大学学报（自然科学版）, 2016, 40(05): 29-33.
[10]	康桂娟,黎瑜,曾日中. 巴西橡胶树HbNAM基因克隆和表达分析[J]. 南京林业大学学报（自然科学版）, 2016, 40(01): 59-64.
[11]	林树燕,李洁,赵荣,董晓波,丁雨龙. 孝顺竹花芽分化及小孢子发生与雄配子体发育[J]. 南京林业大学学报（自然科学版）, 2015, 39(04): 51-56.
[12]	王浩然,邵志龙,朱燕宇,匡华琳,黄敏仁,陈英. ‘南林895’杨PdNAC1基因克隆及蛋白的亚细胞定位[J]. 南京林业大学学报（自然科学版）, 2015, 39(03): 50-54.
[13]	林树燕,李洁,赵荣,徐强,丁雨龙. 南京地区孝顺竹的开花生物学特性研究[J]. 南京林业大学学报（自然科学版）, 2015, 39(02): 52-56.
[14]	贺春玲,杨东烁,崔家俊,李雪萍,杨莉. 长木蜂筑巢对孝顺竹竹皮厚度和竹节直径的选择[J]. 南京林业大学学报（自然科学版）, 2014, 38(01): 179-182.
[15]	张凯敏,王玉成,杨桂燕,高彩球. 柽柳ThPR1基因的克隆与表达分析[J]. 南京林业大学学报（自然科学版）, 2013, 37(02): 45-49.