泡桐维管形成层及木质部区的转录组分析

doi:10.7525/j.issn.1673-5102.2021.03.005

摘要/Abstract

摘要：

采用Illumina HiSeq^TM2500高通量测序技术，获得泡桐维管形成层及木质部区组织转录组的109 021 918条clean reads（16.35 Gb）。将clean reads从头组装得到104 432个单基因簇（Unigene），平均长度662?nt。将组装得到的Unigenes与公共数据库进行序列比对，分别有40 789（Nr：39.05%）、31 675（NT：30.33%）、15 539（COG：14.87%）、29 168（GO：27.93%）、16 316（KEGG：15.62%）、30 499（SwissProt：29.20%）以及28 828（Pfam：27.6%）个Unigenes获得功能注释。通过与GO数据库的比对分析，注释的29 168个Unigenes归于生物过程、细胞组分及分子功能三大类的55个功能组；15 539条Unigenes注释到COG数据库，被分为25个类别中；基于KEGG数据库可将16 316个Unigenes归于130个代谢途径。此外，在泡桐的维管形成层及木质部区转录组中共检测出16 118个简单序列重复（SSR）位点。为进一步挖掘泡桐重要功能基因提供了大量数据。

关键词: 泡桐, 维管形成层, 高通量测序, 转录组

Abstract:

We obtained about 16.35 Gb transcriptome data in total from vascular cambium and xylem zone of Paulownia tomentosa by Illumina Hiseq^TM 2500 sequencing. After assembled， 104 432 unigenes were obtained， the average length of unigene were 662 nt. Then， the unigenes were blastn and annotated by public databases， and 40 789（Nr： 39.05%）， 31 675（NT：30.33%）， 15 539（COG：14.87%）， 29 168（GO：27.93%）， 16 316（KEGG：15.62%）， 30 499（SwissProt：29.20%） and 28 828（Pfam：27.6%） unigenes were annotated respectively. By GO analysis， 29 168 unigenes were potentially involved in 55 functional categories of biological processes， cellular components， and molecular functions respectively. By COG analysis， 15 539 unigenes were classified into 25 functional groups； whereas， by KEGG analysis， 16 316 unigenes were assigned to 130 known pathways. In addition， 16 118 simple sequence repeat（SSR） loci were determined in the transcriptome of P.tomentosa. This study provided a basic data for further exploration of important functional genes in P. tomentosa.

Key words: Paulownia, vascular cambium, high-throughput sequencing, transcriptome

中图分类号:

S792.43

滕云, 邱宗波, 王稳利, 张雪如. 泡桐维管形成层及木质部区的转录组分析[J]. 植物研究, 2021, 41(3): 354-361.

Yun TENG, Zong-Bo QIU, Wen-Li WANG, Xue-Ru ZHANG. Transcriptome Analysis in Vascular Cambium and Xylem of Paulownia tomentosa[J]. Bulletin of Botanical Research, 2021, 41(3): 354-361.

图/表 6

表1

表2

图1

图2

表3

图3

参考文献 19

1	Qiu Z B，Wan L C，Chen T，et al.The regulation of cambial activity in Chinese fir（Cunninghamia lanceolata）involves extensive transcriptome remodeling［J］.New Phytologist，2013，199（3）：708-719.
2	许会敏，王莉，曹德昌，等.维管形成层活动周期调控研究进展［J］.科学通报，2015，60（7）：619-629.
	Xu H M，Wang L，Cao D C，et al.Research progress on the regulation of cambium activity periodicity［J］.Chinese Science Bulletin，2015，60（7）：619-629.
3	Hertzberg M，Aspeborg H，Schrader J，et al.A transcriptional roadmap to wood formation［J］.Proceedings of the National Academy of Sciences of the United States of America，2001，98（25）：14732-14737.
4	Gaete-Loyola J，Lagos C，Beltrán M F，et al.Transcriptome profiling of Eucalyptus nitens reveals deeper insight into the molecular mechanism of cold acclimation and deacclimation process［J］.Tree Genetics & Genomes，2017，13（2）：37.
5	曹玉婷，莫家兴，欧阳磊，等.柳杉维管形成层及木质部区转录组序列及基因表达分析［J］.分子植物育种，2019，17（3）：762-770.
	Cao Y T，Mo J X，Ouyang L，et al.Transcriptome sequence and gene expression in vascular cambium and xylem of Cryptomeria fortunei［J］.Molecular Plant Breeding，2019，17（3）：762-770.
6	杨光，梁坤南，黄桂华，等.基于高通量测序的柚木边材转录组分析［J］.分子植物育种，2020，18（13）：4274-4282.
	Yang G，Liang K N，Huang G H，et al.Transcriptome analysis of sapwood xylem in the Teak（Tectona grandis L.F.） based on high-throughput sequencing［J］.Molecular Plant Breeding，2020，18（13）：4274-4282.
7	Qiu Z B，He Y Y，Zhang Y M，et al.Genome-wide identification and profiling of microRNAs in Paulownia tomentosa cambial tissues in response to seasonal changes［J］.Gene，2018，677：32-40.
8	Dong Y P，Fan G Q，Deng M J，et al.Genome-wide expression profiling of the transcriptomes of four Paulownia tomentosa accessions in response to drought［J］.Genomics，2014，104（4）：295-305.
9	贾新平，孙晓波，邓衍明，等.鸟巢蕨转录组高通量测序及分析［J］.园艺学报，2014，41（11）：2329-2341.
	Jia X P，Sun X B，Deng Y M，et al.Sequencing and analysis of the transcriptome of Asplenium nidus［J］.Acta Horticulturae Sinica，2014，41（11）：2329-2341.
10	Grabherr M G，Haas B J，Yassour M，et al.Full-length transcriptome assembly from RNA-Seq data without a reference genome［J］.Nature Biotechnology，2011，29（7）：644-652.
11	Altschul S F，Gish W，Miller W，et al.Basic local alignment search tool［J］.Journal of Molecular Biology，1990，215（3）：403-410.
12	许蔷薇，楼雄珍，杨彬，等.云锦杜鹃转录组分析［J］.浙江农林大学学报，2019，36（6）：1190-1198.
	Xu Q W，Lou X Z，Yang B，et al.Transcriptome sequencing and analysis of Rhododendron fortunei［J］.Journal of Zhejiang A&F University，2019，36（6）：1190-1198.
13	马婧，邓楠，褚建民，等.泡泡刺高通量转录组鉴定及其黄酮类代谢途径初步分析［J］.林业科学研究，2016，29（1）：61-66.
	Ma J，Deng N，Chu J M，et al.High-throughput transcriptome identification and flavonoids metabolic pathways in Nitraria sphaerocarpa［J］.Forest Research，2016，29（1）：61-66.
14	蔡年辉，邓丽丽，许玉兰，等.基于高通量测序的云南松转录组分析［J］.植物研究，2016，36（1）：75-83.
	Cai N H，Deng L L，Xu Y L，et al.Transcriptome analysis for Pinus yunnanensis based on high throughput sequencing［J］.Bulletin of Botanical Research，2016，36（1）：75-83.
15	王帅，邵芬娟，李论，等.穗花杉的转录组测序及其转录组特性分析［J］.林业科学研究，2017，30（5）：759-764.
	Wang S，Shao F J，Li L，et al.Transcriptome sequencing and analysis of Amentotaxus argotaenia（Hance） Pilger［J］.Forest Research，2017，30（5）：759-764.
16	李太强，刘雄芳，万友名，等.基于高通量测序的极小种群野生植物长梗杜鹃转录组分析［J］.植物研究，2017，37（6）：825-834.
	Li T Q，Liu X F，Wan Y M，et al.Transcriptome analysis for Rhododendron longipedicellatum（plant species with extremely small populations） based on high throughput sequencing［J］.Bulletin of Botanical Research，2017，37（6）：825-834.
17	王兴春，谭河林，陈钊，等.基于RNA-Seq技术的连翘转录组组装与分析及SSR分子标记的开发［J］.中国科学：生命科学，2015，45（3）：301-310.
	Wang X C，Tan H L，Chen Z，et al.Assembly and characterization of the transcriptome and development of SSR markers in Forsythia suspensa based on RNA-Seq technology［J］.Scientia Sinica Vitae，2015，45（3）：301-310.
18	林开勤，李岩，赵德刚.杜仲转录组SSR发掘及标记开发［J］.分子植物育种，2016，14（6）：1548-1558.
	Lin K Q，Li Y，Zhao D G.SSR mining and marker development in Eucommia ulmoides oliver transcriptome［J］.Molecular Plant Breeding，2016，14（6）：1548-1558.
19	魏秀清，许玲，章希娟，等.莲雾转录组SSR信息分析及其分子标记开发［J］.园艺学报，2018，45（3）：541-551.
	Wei X Q，Xu L，Zhang X J，et al.Analysis on SSR information in transcriptome and development of molecular markers in wax apple［J］.Acta Horticulturae Sinica，2018，45（3）：541-551.

长度区间 Length interval	转录本 Transcript	单基因簇 Unigene
200~500	77 452	68 865
500~1 000	26 725	17 742
1 000~2 000	24 155	10 848
>2 000	17 549	6 977
总数量Total number	145 881	104 432
总长度Total length	129 567 733	69 112 491
N₅₀长度N₅₀ length	1 605	1 105
平均长度Mean length	888	662

数据库 Database	注释Unigenes数 Number of annotated Unigenes	占注释Unigenes百分数 Percentage of annotated Unigenes (%)
注释到Nr Annotated in Nr	40 789	39.05
注释到NT Annotated in NT	31 675	30.33
注释到SwissProt Annotated in SwissProt	30 499	29.20
注释到GO Annotated in GO	29 168	27.93
注释到Pfam Annotated in Pfam	28 828	27.6
注释到COG Annotated in COG	15 539	14.87
注释到KEGG Annotated in KEGG	16 316	15.62
注释到至少一个数据库 Annotated in at least one database	48 214	46.16
注释到所有数据库 Annotated in all database	7 386	7.07
Unigenes总量 Total Unigenes	104 432	100

KEGG通路 KEGG pathway	注释基因 Annotated genes	占总基因数的比例 Percentage of total genes(%)	通路ID Pathway ID
核糖体Ribosome	1 216	7.45	ko03010
植物病原体相互作用Plant-bacterium warfare	646	3.96	ko04626
碳代谢Carbon metabolism	638	3.91	ko01200
氨基酸的生物合成Biosynthesis of amino acids	559	3.43	ko01230
内质网蛋白加工Protein processing in endoplasmic reticulum	533	3.27	ko04141
剪接Spliceosome	396	2.43	ko03040
淀粉和蔗糖代谢Starch and sucrose metabolism	384	2.35	ko00500
RNA运输 RNA transport	371	2.27	ko03013
氧化磷酸化Oxidative phosphorylation	370	2.27	ko00190
植物激素信号转导Plant hormone signal transduction	361	2.21	ko04075
胞吞作用Endocytosis	325	1.99	ko04144
糖酵解/糖异生Glycolysis/gluconeogenesis	298	1.83	ko00010
嘌呤代谢Purine metabolism	275	1.69	ko00230
氨基糖和核苷酸糖代谢Amino sugar and nucleotide sugar metabolism	268	1.64	ko00520
苯丙素生物合成Phenylpropanoid biosynthesis	266	1.63	ko00940
泛素介导的蛋白水解Ubiquitin mediated proteolysis	257	1.58	ko04120
吞噬Phagocytosis	248	1.52	ko04145
丙酮酸代谢Pyruvate metabolism	220	1.35	ko00620
RNA降解 RNA degradation	220	1.35	ko03018
mRNA监视通路 mRNA surveillance pathway	217	1.33	ko03015

[1]	倪馨宇, 贺俊英, 燕孟娇, 杜超. RNA-Seq技术在珍稀濒危植物研究中的应用进展[J]. 植物研究, 2023, 43(4): 481-492.
[2]	郑占敏, 商玉冰, 周广波, 肖迪, 刘轶, 由香玲. PsnHB13与PsnHB15在小黑杨中的遗传转化与功能分析[J]. 植物研究, 2023, 43(3): 340-350.
[3]	黄东梅, 陈颖, 白露, 倪迪安, 徐奕扬, 张志国, 秦巧平. 萱草叶片响应低温胁迫的转录组分析[J]. 植物研究, 2022, 42(3): 424-436.
[4]	余静雅, 夏铭泽, 徐浩, 张发起. 青藏高原地区3种蒿属植物转录组比较分析[J]. 植物研究, 2022, 42(2): 200-210.
[5]	赵雪艳, 李思锋, 柏国清, 李为民. 百合属5种植物的比较转录组分析[J]. 植物研究, 2021, 41(4): 614-625.
[6]	周晨晨, 李义良, 王哲, 赵奋成, 钟岁英, 林昌明, 谭志强, 李福明, 吴惠姗, 郭文冰, 廖仿炎, 王建革. 基于RNA-Seq技术的湿加松特定时期产脂差异基因筛选及表达分析[J]. 植物研究, 2021, 41(3): 419-428.
[7]	鲁琳, 杨尚谕, 刘维东, 鲁黎明. 基于转录组测序的花烟草低温胁迫响应转录因子挖掘[J]. 植物研究, 2021, 41(1): 119-129.
[8]	赵历强, 单春苗, 张声祥, 施圆圆, 马克龙, 吴家文. 基于转录组测序的细风轮花青素合成途径及关键酶基因分析[J]. 植物研究, 2020, 40(6): 886-896.
[9]	杨斌, 孟庆瑶, 张凯, 段义忠. 孑遗濒危植物矮扁桃叶绿体全基因组特征分析及亲缘关系鉴定[J]. 植物研究, 2020, 40(5): 686-695.
[10]	刘丹, 曾钦朦, 刘斌, 李煜, 陈世品. 闽楠转录组分析及基因功能注释[J]. 植物研究, 2020, 40(4): 613-622.
[11]	张雨, 夏铭泽, 张发起. 药用资源植物山莨菪的转录组信息分析[J]. 植物研究, 2020, 40(3): 458-467.
[12]	朱利利, 杜庆鑫, 何凤, 庆军, 杜红岩. 杜仲雄花芽2个发育时期转录组分析[J]. 植物研究, 2020, 40(2): 284-292.
[13]	熊阳阳, 罗琳莉, 赵士豪, 周安龙, 王锦. 大花铁线莲单重被花的转录组差异分析[J]. 植物研究, 2020, 40(1): 93-105.
[14]	段义忠, 杜忠毓, 王海涛. 西北干旱区孑遗濒危植物四合木(Tetraena mongolica)叶绿体基因组特征研究及比较分析[J]. 植物研究, 2019, 39(5): 653-663.
[15]	叶兴状, 刘丹, 罗佳佳, 范辉华, 张国防, 刘宝, 陈世品. 濒危珍稀植物半枫荷的转录组分析[J]. 植物研究, 2019, 39(2): 276-286.