摘要
为了获得甘薯耐盐转录组序列信息,挖掘差异表达基因及其相关代谢途径,本文以盐胁迫处理0 d、3 d和6d的耐盐甘薯品种‘榕薯819’以及不耐盐甘薯品种‘榕薯910’的叶片为材料,借助高通量测序技术进行转录组测序分析。结果表明, 2个品种共获得157,252条Unigenes,平均组装长度为576 bp。其中有83,264条Unigenes在七大数据库中得到注释,占总数的52.95%。NR注释分类结果显示,在牵牛花(Ipomoea nil)中比对到同源序列的Unigenes最多,共43,620条,占总数的57.05%。Unigenes在KOG数据库中的注释主要富集在普通功能预测(8752个)、信号转导机制(5067个)以及翻译后修饰、蛋白转换、分子伴侣(4471个)中。差异表达分析显示,在‘榕薯819’中,盐处理3 d和6d的样品差异表达基因数分别为323个和3752个,共参与了33个GO功能分类项和302条KEGG代谢通路。在‘榕薯910’中,差异表达基因数则分别为5554个和7395个,共参与了50个GO功能分类项,涉及了329条KEGG代谢通路。以部分差异表达基因的转录组数据为基础,进行了表达量热图绘制。结果显示,注释到淀粉和蔗糖代谢途径的7个差异表达β-葡萄糖苷酶基因均表现为在耐盐品种中上调,在不耐盐品种中下调;注释到Ca2+信号途径的7个类钙调素基因中,有2个在耐盐品种中特异上调,5个在不耐盐品种中特异下调。差异表达转录因子表达量热图显示,IbERF1仅在耐盐品种中特异表达,IbNAC3、IbNAC11、IbERF3和IbERF4仅在不耐盐品种中特异表达,而IbNAC29在耐盐品种中下调表达,在不耐盐品种中上调表达。综上分析,甘薯耐盐转录组获得的Unigenes数量较大,序列信息丰富,盐胁迫下获得的差异表达基因及高丰度转录因子可能在甘薯抵御盐胁迫过程中发挥着重要作用。
The objective of this study is to obtain the salt tolerant transcriptome sequence information of sweet potato, mine and identify the differentially expressed genes and their related metabolic pathways. Rongshu 819(salt-tolerant variety) and Rongshu 910(salt-sensitive variety) under salt treatment for 0, 3, and 6 days were conducted for the transcriptome by high-throughput sequencing. The results showed that a total of 157,252 Unigenes with an average length of 576 bp were obtained from the two varieties. The 83,264 Unigenes, accounting for 52.95% of the total Unigenes, were annotated in the 7 functional databases. NR annotation revealed that sweet potato Unigenes had the most homologous sequences in Ipomoea nil with a total of 43,620 accounting for 57.05%. The annotation of Unigenes in the KOG database is mainly concentrated in general functional prediction only(8752),signal transformation mechanisms(5067), and post translation modification, protein transformation, molecular chaperones(4471) .The qRT-PCR indicated that, in Rongshu 819, the number of differentially expressed genes in 3 days and 6 days of salt treatment was 3752 and 807, respectively, which classified into 33 GO functional categories and 302 KEGG metabolic pathways. In Rongshu 910, the number of differentially expressed genes in 3 days and 6 days of salt treatment was 5554 and 7395, respectively,which classified into 50 GO functional categories and 329 KEGG metabolic pathways. The heat map was drawn based on the transcriptome data of partial differentially expressed genes. The result showed that 7 differentially expressed β-Glucosidase genes annotated into starch and sucrose metabolic pathway were up-regulated in salt-tolerant variety while down-regulated in salt-sensitive variety. Seven calmodulin-like genes annotated into Ca2+ signaling pathway were specifically expressed in two varieties, among which 2 were specifically up-regulated in salt-tolerant variety and 5 were specifically down-regulated in salt-sensitive variety. The expression heat map of differentially expressed transcription factors indicated that IbERF1 was only specifically expressed in salt-tolerant varieties, IbNAC3, IbNAC11, IbERF3, and IbERF4 were only specifically expressed in salt-sensitive varieties, while IbNAC29 was down-regulated in salt-tolerant varieties and up-regulated in salt-sensitive varieties. In conclusion, the number of Unigenes obtained from transcriptome of sweet potato under salt stress was large and the sequence information was rich. The differentially expressed genes and high abundance transcription factors may play an important role in sweet potato response to salt stress.
作者
张小红
彭琼
鄢铮
ZHANG Xiao-Hong;PENG Qiong;YAN Zheng(Fuzhou Institute of Agricultural Sciences,Fuzhou 350018,Fujian,China)
出处
《作物学报》
CAS
CSCD
北大核心
2023年第5期1432-1444,共13页
Acta Agronomica Sinica
基金
福建省科技厅引导性项目(2019N0037)资助。
关键词
甘薯
盐胁迫
转录组
高通量测序
sweet potato
salt stress
transcriptome
high-throughput sequencing
