杜仲雄花芽2个发育时期转录组分析

doi:10.7525/j.issn.1673-5102.2020.02.016

摘要/Abstract

摘要： 杜仲（Eucommia ulmoides）雄花富含多种活性成分和营养成分，具有重要的药用和营养价值。为了揭示杜仲雄蕊原基发育相关基因的表达情况，为杜仲雄花芽发育分子调控机制研究提供理论参考。本文以杜仲良种"华仲11号"（ "Huazhong No.11" ）为材料，采用lllumina高通量测序技术，分别对苞叶原基分化期和雄蕊原基分化期的花芽进行转录组测序，通过生物信息学对2个发育时期的转录组进行比较分析，筛选出与雄花芽形态发育相关的差异基因。结果显示，转录组测序共获得40.48 Gb过滤数据，各样品的clean reads与杜仲基因组进行序列比对，比对效率为90.56%~93.01%。在2个发育时期筛选出583个差异表达基因，其中在雄蕊原基发育期上调基因315个，下调基因267个。差异基因GO和KEGG功能分析显示，差异基因富集在与生长发育、光周期途径、激素合成和信号传导、碳代谢等相关的生物过程和代谢通路。结果显示光周期途径是杜仲成花诱导的重要途径，同时雄花芽在形态分化过程中受碳水化合物、植物激素和其他代谢物质调控。此外，MADS-box家族成员FLC、SOC1、AGL3和AGL8参与杜仲雄蕊器官发育。本研究为杜仲花发育基因调控提供了基础数据，也为雄花用杜仲的分子育种提供了参考。

关键词: 杜仲, 雄花芽, 不同发育时期, 转录组测序, 差异表达基因

Abstract: The male floral bud Eucommia ulmoides is rich in healthily kinds of nutrition compositions and active matter, which with a high nutritional and medicinal value. The transcriptome of male floral bud during two differentiation stages was sequenced, aimed to comprehend the expression of genes involved in the stamen development of E.ulcommia. lllumina high-throughput sequencing technology was applied to detect the transcriptome sequencing data of the male flora bud of ‘Huazhong No.11’ during bract differentiation and stamen differentiation stages. Obtained clean data was analyzed including gene function annotation and differentially expressed genes(DEGs) screening was analyzed by bioinformatics methods. A total of 40.48 Gb clean data were obtained, approximately 90.00% of the clean reads were mapped to the E.ulmoides reference genome. It was found that 315 genes were significantly up-regulated and 269 genes were significantly down-regulated in the male flower buds at stamen differentiation stage. GO and KEGG analyses of DEGs showed the DEGs were enriched in development process, photoperiodism, hormone biosynthetic and signal transduction process, carbon metabolism and other processes and pathways associated with flower induction. The photoperiodic was perhaps the most important pathway for the floral induction of E.ulcommia, and the carbohydrates, plant hormone and other metabolic substances was required in the process of male floral bud differentiation. Importantly, the members of MADs-box including FLC, SOC1, AGL3 and AGL8 regulated the stamen formation of E.ulcommia. It provides a significant reference for studying the genes involved in the floral organs and molecular breeding of E.ulcommia.

Key words: Eucommia ulmoides, male floral bud, different developmental stages, transcriptome, differentially expressed genes

中图分类号:

S722.3

朱利利, 杜庆鑫, 何凤, 庆军, 杜红岩. 杜仲雄花芽2个发育时期转录组分析[J]. 植物研究, 2020, 40(2): 284-292.

ZHU Li-Li, DU Qing-Xin, HE Feng, QING Jun, DU Hong-Yan. Sequencing Analysis of Transcriptome of Male Floral Bud at Two Development Stages in Eucommia ulmoides[J]. Bulletin of Botanical Research, 2020, 40(2): 284-292.

参考文献 38

[1]	杜红岩.中国杜仲图志[M].北京:中国林业出版社,2014:76-79.Du H Y.China Eucommia pictorial[M].Beijing:China Forestry Publishing House,2014:76-79.
[2]	王丙武,王雅清,莫华,等.杜仲雌雄株细胞学、顶芽及叶含胶量的比较[J].植物学报,1999,41(1):11-15.Wang B W,Wang Y Q,Mo H,et al.Comparison of cytology,apical buds and gutta content between staminate and pistillate of Eucommia ulmoides trees[J].Acta Botanica Sinica,1999,41(1):11-15.
[3]	杜庆鑫,刘攀峰,庆军,等.杜仲雄花形态性状的遗传多样性分析[J].林业科学研究,2016,29(5):670-675.Du Q X,Liu P F,Qing J,et al.Genetic diversity analysis on morphology traits of Eucommia ulmoides male flowers[J].Forest Research,2016,29(5):670-675.
[4]	杜庆鑫,刘攀峰,魏艳秀,等.基于主成分与聚类分析的杜仲雄花品质综合评价[J].植物研究,2016,36(6):846-852.Du Q X,Liu P F,Wei Y X,et al.Comprehensive evaluation of Eucommia ulmoides male flowers quality by principal component and cluster analysis[J].Bulletin of Botanical Research,2016,36(6):846-852.
[5]	杜庆鑫,魏艳秀,刘攀峰,等.杜仲雄花主要活性成分含量的多样性[J].林业科学,2017,53(2):35-43.Du Q X,Wei Y X,Liu P F,et al.Diversity of the content of main active components in Eucommia ulmoides male flowers[J].Scientia Silvae Sinicae,2017,53(2):35-43.
[6]	李欣,乔家驹,冯汉青,等.杜仲雄花乙酸乙酯提取物的抑菌活性研究[J].食品工业科技,2015,36(11):62-65.Li X,Qiao J J,Feng H Q,et al.Study on antibacterial activities of the ethyl acetate extract in male Eucommia flowers[J].Science and Technology of Food Industry,2015,36(11):62-65.
[7]	杜红岩,胡文臻,俞锐,等.杜仲产业绿皮书:中国杜仲橡胶资源与产业发展报告[M].北京:社会科学文献出版社,2013:188-195.Du H Y,Hu W Z,Yu R,et al.Green book of Eucommia industry[M].Beijing:Social Sciences Academic Press,2013:188-195.
[8]	Zik M,Irish V F.Flower development:initiation,differentiation,and diversification[J].Annual Review of Cell and Developmental Biology,2003,19:119-140.
[9]	刘朝晖.杜仲(杜仲科)花器官的发生和发育及其系统学意义[D].北京:北京大学,2010.Liu Z H.Floral organogenesis and phylogenetic analysis in Eucommia ulmoides Oilv[D].Beijing:Peking University,2010.
[10]	Kim D,Langmead B,Salzberg S L.HISAT:a fast spliced aligner with low memory requirements[J].Nature Methods,2015,12(4):357-360.
[11]	Pertea M,Pertea G M,Antonescu C M,et al.StringTie enables improved reconstruction of a transcriptome from RNA-seq reads[J].Nature Biotechnology,2015,33(3):290-295.
[12]	Altschul S F,Madden T L,Schäffer A A,et al.Gapped BLAST and PSI-BLAST:a new generation of protein database search programs[J].Nucleic Acids Research,1997,25(17):3389-3402.
[13]	邓泱泱,荔建琦,吴松锋,等.nr数据库分析及其本地化[J].计算机工程,2006,32(5):71-73,76.Deng Y Y,Li J Q,Wu S F,et al.Integrated nr database in protein annotation system and its localization[J].Computer Engineering,2006,32(5):71-73,76.
[14]	Apweiler R,Bairoch A,Wu C H,et al.UniProt:the universal protein knowledgebase[J].Nucleic Acids Research,2004,32(S1):D115-D119.
[15]	Ashburner M,Ball C A,Blake J A,et al.Gene ontology:tool for the unification of biology[J].Nature Genetics,2000,25(1):25-29.
[16]	Tatusov R L,Galperin M Y,Natale D A,et al.The COG database:a tool for genome-scale analysis of protein functions and evolution[J].Nucleic Acids Research,2000,28(1):33-36.
[17]	Koonin E V,Fedorova N D,Jackson J D,et al.A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes[J].Genome Biology,2004,5(2):R7.
[18]	Finn R D,Bateman A,Clements J,et al.Pfam:the protein families database[J].Nucleic Acids Research,2014,42(D1):D222-D230.
[19]	Kanehisa M,Goto S,Kawashima S,et al.The KEGG resource for deciphering the genome[J].Nucleic Acids Research,2004,32(S1):D277-D280.
[20]	Xie C,Mao X Z,Huang J J,et al.KOBAS 2.0:a web server for annotation and identification of enriched pathways and diseases[J].Nucleic Acids Research,2011,39(S2):W316-W322.
[21]	Eddy S R.Profile hidden Markov models[J].Bioinformatics,1998,14(9):755-763.
[22]	Pertea M,Kim D,Pertea G M,et al.Transcript-level expression analysis of RNA-seq experiments with HISAT,StringTie and Ballgown[J].Nature Protocols,2016,11(9):1650-1667.
[23]	彭向永.雌、雄蒿柳花芽分化机制及性别决定基因挖掘[D].北京:中国林业科学研究院,2017.Peng X Y.Mechanism of floral bud differentiation and discovery of sex determination genes in male and female Salix viminalis[D].Beijing:Chinese Academy of Forestry,2017.
[24]	常莉,薛建平.生长素极性运输研究进展[J].生物学杂志,2008,25(6):9-13.Chang L,Xue J P.Advances in study of polar auxin transport[J].Journal of Biology,2008,25(6):9-13.
[25]	Shi Q L,Dong Y B,Qiao D H,et al.Isolation and characterization of ZmERF1 encoding ethylene responsive factor-like protein 1 in popcorn(Zea mays L.)[J].Plant Cell,Tissue and Organ Culture,2015,120(2):747-756.
[26]	程华,李琳玲,袁红慧,等.内源激素含量变化与板栗花芽分化关系研究[J].北方园艺,2013(22):5-9.Cheng H,Li L L,Yuan H H,et al.Study on the relationship of changes of the endogenous hormone content in Castanea mollissima and bud differentiation[J].Northern Horticulture,2013(22):5-9.
[27]	Ledger S,Strayer C,Ashton F,et al.Analysis of the function of two circadian-regulated CONSTANS-LIKE genes[J].Plant Journal,2010,26(1):15-22.
[28]	苏亚中.小麦CONSTANS-like基因的克隆与分析[D].开封:河南大学,2012.Su Y Z.Cloning and analysis of CONSTANS-like Genes in wheat(Triticum aestivum L.)[D].Kaifeng:Henan University,2012.
[29]	Liu H,Gu F W,Dong S Y,et al.CONSTANS-like 9(COL9) delays the flowering time in Oryza sativa by repressing the Ehd1 pathway[J].Biochemical and Biophysical Research Communications,2016,479(2):173-178.
[30]	Valverde F,Mouradov A,Soppe W,et al.Photoreceptor regulation of CONSTANS Protein in photoperiodic flowering[J].Science,2004,303(5660):1002-1006.
[31]	Fujiwara S,Oda A,Yoshida R,et al.Circadian clock proteins LHY and CCA1 regulate SVP protein accumulation to control flowering in Arabidopsis[J].The Plant Cell,2008,20(11):2960-2971.
[32]	Hicks K A,Albertson T M,Wagner D R.EARLY FLOWERING3 encodes a novel protein that regulates circadian clock function and flowering in Arabidopsis[J].The Plant Cell,2001,13(6):1281-1292.
[33]	李元元,王鲁,苏振刚,等.MADS-box基因控制植物成花的分子机理[J].基因组学与应用生物学,2010,29(6):1122-1132.Li Y Y,Wang L,Su Z G,et al.The molecular mechanism of MADS-box genes regulates floral formation and flowering in plant[J].Genomics and Applied Biology,2010,29(6):1122-1132.
[34]	Michaels S D,Amasino R M.FLOWERING LOCUS C encodes a novel MADS domain protein that acts as a repressor of flowering[J].The Plant Cell,1999,11(5):949-956.
[35]	Samach A,Onouchi H,Gold S E,et al.Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis[J].Science,2000,288(5471):1613-1616.
[36]	Moon J,Suh S S,Lee H,et al.The SOC1 MADS-box gene integrates vernalization and gibberellin signals for flowering in Arabidopsis[J].The Plant Journal,2003,35(5):613-623.
[37]	Ditta G,Pinyopich A,Robles P,et al.The SEP4 gene of Arabidopsis thaliana functions in floral organ and meristem identity[J].Current Biology,2004,14(21):1935-1940.
[38]	Mandel M A,Yanofsky M F.The Arabidopsis AGL8 MADS box gene is expressed in inflorescence meristems and is negatively regulated by APETALA1[J].The Plant Cell,1995,7(11):1763-1771.