文章信息
- 唐海霞, 杜淑辉, 邢世岩, 桑亚林, 李际红, 刘晓静, 孙立民
- Tang Haixia, Du Shuhui, Xing Shiyan, Sang Yalin, Li Jihong, Liu Xiaojing, Sun Limin
- 银杏性别决定相关基因的筛选
- Screening of Sex Determination Related Genes in Ginkgo biloba
- 林业科学, 2017, 53(2): 76-82.
- Scientia Silvae Sinicae, 2017, 53(2): 76-82.
- DOI: 10.11707/j.1001-7488.20170209
-
文章历史
- 收稿日期:2016-01-29
- 修回日期:2016-04-12
-
作者相关文章
2. 山西农业大学林学院 太谷 030800
2. College of Forestry, Shanxi Agricultural University Taigu 030800
银杏 (Ginkgo biloba) 是我国特有的珍稀树种,有“活化石”、植物中的“熊猫”之称,在我国种植范围广泛。作为裸子植物中独特的一类,银杏现仅存一科一属一种 (邢世岩等,2013)。银杏是集材用、果用、药用、观赏于一体的高效多功能经济林树种,具有较高的经济价值、生态效益和研究价值。随着国内外市场对银杏食用、材用、保健、药用等领域的不断拓展,银杏资源已供不应求 (Zhang et al., 2015)。银杏为典型的雌雄异株植物,雌雄株在生长习性、形态特征与生产应用等方面存在显著差异。雄株生长速度快、树形美观、抗性强、根系发达,常用作绿化观赏或用材树使用;银杏种子含有银杏酸等生理药理活性物质,具有重要的食用、药用价值,因此雌株主要以生产种子为主 (黄茜等,2013)。性别决定机制是雌雄性别差异产生的根本原因,而性别决定基因是性别决定的直接执行者 (Charlesworth,2013),但是目前国内外关于银杏性别决定基因的研究还未见报道。对银杏性别决定基因的研究不仅可为破解银杏性别决定机制打下坚实基础,同时也能够加深对植物性别决定机制的理解。
早期的性别决定基因研究大多以一些与性别有关的发育性突变体为材料,研究性别决定基因的作用机制。随着分子生物学技术,特别是分离目的基因的方法的发展,mRNA差异显示PCR、消减杂交法、基因芯片技术、基因表达的系统分析等方法已广泛应用到大规模基因分离与表达分析研究中,但是以上方法普遍存在需要已知目的基因序列、灵敏度低、难以检测低丰度的目标、无法捕捉到目的基因表达水平的微小变化等缺点 (Wang et al., 2009)。随着新一代测序技术 (next-generation sequencing,NGS) 的迅猛发展,RNA-Seq (RNA-sequencing) 技术的应用弥补了上述研究方法存在的各种不足。Wu等 (2010)利用RNA-Seq技术对黄瓜 (Cucumis sativus) 雌雄花芽的转录组进行了测序与分析,发现了大量在雌雄花芽中差异表达或特异表达的转录因子和基因,为进一步深入探明黄瓜的性别决定机制奠定了坚实基础。Liu等 (2013)利用RNA-Seq技术对簸箕柳 (Salix suchowensis) 雌雄花芽差异表达基因进行了分析,发现了在雌雄花芽中差异表达的806个基因,其中有33个可能定位于性染色体上,有12个基因可能直接参与了簸箕柳的性别决定,为最终发现杨柳科 (Salicaceae) 性别决定基因和未来的功能基因组学研究奠定了基础。本研究利用高通量测序技术平台Illumina Solexa对银杏雌雄花芽和大、小孢子叶球进行转录组测序、拼接组装,再利用生物信息学方法对得到的unigene进行功能注释与功能分类,以期为理解银杏性别决定机制奠定基础。
1 材料与方法 1.1 试验材料根据相关文献 (史继孔等,1998;张万萍等,2001;韩晨静,2011) 报道,在山东农业大学银杏种质资源圃内挑选来源于同一家系的25年生银杏雌雄株各5株,分别于2015年3月4日和4月10日采集花芽和大、小孢子叶球样本各10个 (每株树采集2个样本,共采集40个样本),每5个样本为一组 (每组样品设置2个生物学重复) 用于构建cDNA文库。所采样本用液氮速冻后存于-80 ℃冰箱备用。
1.2 RNA提取使用改良的CTAB-LiCl法 (Chang et al., 1993) 提取总RNA,然后使用RNeasy Plant Mini Kit (QIAGEN) 对RNA样品进一步纯化。利用Agilent的Nano 6000 Assay Kit和琼脂糖凝胶电泳检测RNA的完整性和质量,保证RNA样品浓度大于400 ng·μL-1,28S/18S大于1.8。
1.3 cDNA文库构建与测序利用NEB Next®UltraTM RNA Library Prep Kit for Illumina®(NEB, USA) 构建cDNA文库,分别为雄花芽 (XY)、雌花芽 (CY)、小孢子叶球 (XH) 和大孢子叶球 (CH),然后利用Illumina Hiseq 2500平台对构建的cDNA文库进行DNA双末端测序,得到2×125 bp的原始数据。
1.4 数据拼接和功能注释去掉原始数据中的低质量序列和接头序列后,将所有高质量序列混合后利用软件Trinity (Grabherr et al., 2011) 进行组装得到unigene。利用BLAST系列软件将银杏的unigene序列与Nr (NCBI非冗余蛋白数据库)、Nt (NCBI核苷酸数据库)、Swiss-Prot蛋白质数据库、Pfam蛋白质数据库、COG/KOG、GO、KEGG数据库进行比对,获得unigene的注释信息。
1.5 qRT-PCR验证随机挑选26个差异表达基因,利用软件Primer 5进行引物设计。根据SYBR Green Dye (Takara, Dalian, China) 试剂盒进行荧光定量PCR反应,每个反应重复3次。用GAPDH(Zhang et al., 2015) 作为内参基因,采用2-ΔΔCt方法 (Livak et al., 2001) 计算每个基因表达量的变化。
1.6 差异表达基因的识别与富集分析利用Trapnell等 (2010)的方法计算每个unigene的FPKM (每百万片段中比对到某一基因每千碱基长度的片段数目,fragments per transcript kilobase per million fragments mapped) 值,作为衡量unigene表达量的标准。利用DESeq R软件包 (Anders,2010) 分析XY/CY和XH/CH中unigene表达量的差异,显著性检验采用校正后的P值 (Storey et al., 2003)。当检验结果的P值小于0.05时,就认为这个unigene是差异表达基因。利用软件GOseq (Young et al., 2010) 和KOBAS (Mao et al., 2005) 对差异基因进行功能富集分析。
2 结果与分析 2.1 测序结果及数据拼接利用Illumina Hiseq 2500平台对8个cDNA文库进行了高通量测序,去掉接头序列及低质量序列后,分别得到了119 494 172(雄花芽,XY)、115 958 434(雌花芽,CY)、121 730 470(小孢子叶球,XH)、116 504 860(大孢子叶球,CH) 条高质量序列。所有原始测序数据均已上传到NCBI Sequence Read Archive (SRA) 数据库,登录号分别为SRR2147720(XH),SRR2147715(CH),SRR2147717(CY),SRR2147721(XY)。将所有高质量序列混合后采用软件Trinity进行数据拼接,最终共获得108 307条unigene。unigene的平均长度为796 bp,N50为1 648 bp,其中长度大于1 000 bp的unigene有23 624条,占全部unigene的21.81%。
2.2 qRT-PCR验证利用qRT-PCR对随机挑选的26个差异表达基因在XY,CY,XH和CH中的表达量进行验证,从图 1可以看出,RNA-Seq结果与qRT-PCR结果具有很高的相关性 (R2> 0.9), 表明RNA-Seq方法准确反映了基因在不同器官中表达水平的变化。
使用BLAST系列软件将unigene序列与Nr,Nt,Swiss-Prot,Pfam,COG/KOG,GO,KEGG数据库进行比对,获得unigene的注释信息。从表 1可以看出,6 347条unigene (5.86%) 在所有数据库中都获得了注释信息,有51 953条 (47.97%) 在至少1个数据库中获得了注释信息,在Nr和SwissProt数据库中得到注释信息的unigene数量最多,而在KEGG和COG/KOG数据库中得到注释信息的unigene数量相对偏少。unigene在GO和KEGG数据库中注释结果见图 2和图 3。
为了能够确定参与银杏性别决定过程的基因,利用软件DESeq对XY/CY和XH/CH中unigene的表达量进行了比较。结果表明,有4 709条unigene在XY/CY中差异表达,其中2 235条unigene在XY中上调表达而有2 474条unigene下调表达;有9 802条unigene在XH/CH中差异表达,其中5 849条unigene在XH中上调表达而3 953条unigene下调表达。
对差异表达的unigene进行GO富集分析,结果表明XY/CY中的4 709条差异表达基因可以划分为36个类别,其中代谢过程 (GO:0008152,2 154个unigenes) 和细胞代谢过程 (GO:0044237,1 678个unigenes) 所占比例最大;XH/CH中的9 802条差异表达基因可以划分为18个类别,其中细胞组分合成 (GO:0044085,427个unigenes) 和核糖核蛋白复合体合成 (GO:0022613,272个unigenes) 所占比例最大。差异表达基因的KEGG代谢通路富集分析结果表明在XH/CH中,核糖体 (143个unigenes)、淀粉和蔗糖代谢 (97个unigenes)、植物激素信号转导 (76个unigenes) 所占比例最大;在XY/CY中,核糖体 (190个unigenes) 和苯丙素生物合成 (33条unigenes) 所占比例最大。在76个参与植物激素信号转导的差异表达基因中,有11个基因在XY/CY和XH/CH中共同差异表达 (表 2)。这些基因可能在银杏性别决定中发挥着重要的作用。
植物的性别决定涉及到复杂的调控网络以及在不同发育阶段差异或特异表达的功能基因 (Aryal et al., 2014)。有研究表明参与花发育的功能基因 (如ABC类基因) 并不是性别决定基因 (Ainsworth et al., 1995;Hardenack et al., 1994),而参与植物激素合成与信号转导的基因在植物性别决定中发挥着重要的作用 (Aryal et al., 2014;Guo et al., 2010)。从表 2中可以看到,2个负调控ABA信号转导的基因PP2C(protein phosphatase 2C) 和PYL(abscisic acid receptor 2) 在雄芽 (XY) 和小孢子叶球 (XH) 中下调表达,而正调控ABA信号转导的基因SNRK2 (SNF1-related protein kinase subfamily 2) 在XY和XH中上调表达。在大麦 (Hordeum vulgare) 中,反义表达SNRK2基因会导致花粉败育,从而引起雄性不育 (Zhanget al.,2001)。由此可以推测,参与ABA信号转导的基因可能在银杏雄性性别决定中发挥作用 (Fujii et al., 2009;Park et al., 2009)。生长素诱导基因SAUR(small auxin up RNA) 和生长素应答基因GH3(Gretchen Hagen 3) 在XY和XH中分别下调和上调表达,说明参与生长素信号转导的基因在银杏性别决定中发挥的作用有一定的差异,这与Kong等 (2013)对拟南芥 (Arabidopsis thaliana) 的研究结果一致。在XY和XH中,参与茉莉酸信号转导的基因COI1(coronatine insensitive 1) 上调表达。拟南芥COI1突变体呈现出不同程度的雄性不育 (Ellis et al., 2002;Turner et al., 2002),这说明茉莉酸在银杏雄性性别决定和育性发育方面也有重要作用。参与水杨酸信号转导的基因NPR1(nonexpressor of PR) 在XY和XH中上调表达,但是这个基因的上调表达可能只与银杏雄株的抗性相关,而与性别决定的关系不大 (Mou et al., 2003)。参与赤霉素信号转导的转录因子DELLA蛋白和参与乙烯信号转导的转录因子EIN3(ethylene-insensitive protein 3) 在XY和XH中上调表达。有研究表明DELLA蛋白D8在玉米 (Zea mays) 的性别决定中发挥重要的作用 (Peng et al., 1999);黄瓜中参与乙烯生物合成与信号转导,以及乙烯诱导表达的大量基因都在性别决定中发挥重要的作用 (Boualem,2009;Wu et al., 2010);拟南芥EIN3突变体中乙烯信号转导受阻并且雄性不育 (Solano et al., 1998)。因此可以推测乙烯在银杏性别决定中同样发挥一定的作用。
银杏的性染色体为ZW型,雄性性染色体为ZZ而雌性为ZW (Lan et al., 2008)。杨树的性染色体组成同样为ZW型,第19号染色体可能是杨树的性染色体 (Yin et al., 2008),Song等 (2013)研究发现2个参与DNA甲基化的基因定位到毛白杨性染色体上的性别决定区域。DNA甲基化对植物的性别决定起到重要的调控作用 (Hultquist et al., 2008;Janoušek et al., 1996;Martin et al., 2009;Parkinson et al., 2007)。本次研究中在雌芽 (CY) 和大孢子叶球 (CH) 中发现了参与DNA甲基化的特异表达基因,其中包括已经报道的参与植物性别决定的基因,如MET1和COMT1(Do et al., 2007;Schmidt et al., 2013)。另外,差异基因的GO富集分析表明,编码DNA甲基化转移酶的基因在CY和CH中均上调表达。这一结果与杨树中研究结果 (Song et al., 2012) 类似,CY和CH中DNA甲基化水平的上升可能有助于银杏雌性的分化与发育。
4 结论本研究通过使用RNA-Seq技术对银杏雌、雄花芽和大、小孢子叶球进行转录组测序与生物信息学分析,发现11个参与植物激素信号转导的基因在雄芽/雌芽和小孢子叶球/大孢子叶球中共同差异表达,参与DNA甲基化的基因在雌芽和大孢子叶球中显著上调表达,这些基因的差异表达可能促进了银杏雌性或雄性的分化,从而在银杏的性别决定中发挥一定的作用。本研究表明银杏的性别决定涉及多种调控途径和功能基因,性别决定相关基因的筛选为全面理解银杏的性别决定机制奠定了基础。
[] |
韩晨静. 2011. 银杏与叶籽银杏胚珠发育的组织化学比较研究. 泰安: 山东农业大学硕士学位论文. ( Han C J. 2011. Comparative studies on histochemistry of the developmental ovule in Ginkgo biloba L. and Ginkgo biloba var. epiphylla Mak. Tai'an:MS thesis of Shandong Agricultural University.[in Chinese][in Chinese]) http://cdmd.cnki.com.cn/Article/CDMD-10434-1011097847.htm |
[] |
黄茜, 刘霁瑶, 曹敏, 等. 2013. 银杏性别特征表现与鉴别研究进展. 果树学报, 30(6): 1065–1071.
( Huang Q, Liu Q Y, Cao M, et al. 2013. Advance of gender characteristics expression and identification of Ginkgo biloba. Journal of Fruit Science, 30(6): 1065–1071. [in Chinese] ) |
[] |
史继孔, 樊卫国, 文晓鹏. 1998. 银杏雌花芽形态分化的研究. 园艺学报, 25(1): 33–36.
( Shi J K, Fan W G, Wen X P. 1998. Female flower bud differentiation of Ginkgo biloba. Acta Horticulture Sinica, 25(1): 33–36. [in Chinese] ) |
[] |
邢世岩, 张倩, 付兆军, 等. 2013. 银杏垂乳个体发生及系统学意义. 林业科学, 49(8): 108–116.
( Xing S Y, Zhang Q, Fu Z J, et al. 2013. Ontogenesis and systematics implications of Ginkgo biloba chichi. Scientia Silvae Sinicae, 49(8): 108–116. [in Chinese] ) |
[] |
张万萍, 史继孔, 樊卫国, 等. 2001. 银杏雄花芽的形态分化. 园艺学报, 28(3): 255–258.
( Zhang W P, Shi J K, Fan W G, et al. 2001. Morphological differentiation of male flower of Ginkgo biloba. Acta Horticulture Sinica, 28(3): 255–258. [in Chinese] ) |
[] | Ainsworth C, Crossley S, Buchanan-Wollaston V, et al. 1995. Male and female flowers of the dioecious plant sorrel show different patterns of MADS Box gene expression. Plant Cell, 7(10): 1583–1598. DOI:10.1105/tpc.7.10.1583 |
[] | Anders S. 2010. Differential expression analysis for sequence count data. Genome Biology, 11(10): 79–82. |
[] | Aryal R, Ming R. 2014. Sex determination in flowering plants:papaya as a model system. Plant Science, 217: 56–62. |
[] | Boualem A. 2009. A conserved ethylene biosynthesis enzyme leads to andromonoecy in two Cucumis species. PloS One, 4(7): e6144. DOI:10.1371/journal.pone.0006144 |
[] | Chang S J, Puryear J, Caimey J. 1993. A simple and efficient method for isolation RNA from pine trees. Plant Molecular Biology Reporter, 11(2): 113–116. DOI:10.1007/BF02670468 |
[] | Charlesworth D. 2013. Plant sex chromosome evolution. Journal of Experimental Botany, 64(2): 405–420. DOI:10.1093/jxb/ers322 |
[] | Do C T, Pollet B, Thévenin J, et al. 2007. Both caffeoyl coenzyme A 3-O-methyltransferase 1 and caffeic acid O-methyltransferase 1 are involved in redundant functions for lignin, flavonoids and sinapoyl malate biosynthesis in Arabidopsis. Planta, 226(5): 1117–1129. DOI:10.1007/s00425-007-0558-3 |
[] | Ellis C, Turner J G. 2002. A conditionally fertile COI1 allele indicates cross-talk between plant hormone signalling pathways in Arabidopsis thaliana seeds and young seedlings. Planta, 215(4): 549–556. DOI:10.1007/s00425-002-0787-4 |
[] | Fujii H, Zhu J K. 2009. Arabidopsis mutant deficient in 3 abscisic acid-activated protein kinases reveals critical roles in growth, reproduction, and stress. Proceedings of the National Academy of Sciences, 106(20): 8380–8385. DOI:10.1073/pnas.0903144106 |
[] | Grabherr M G, Haas B J, Yassour M, et al. 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology, 29(7): 644–652. DOI:10.1038/nbt.1883 |
[] | Guo S, Zheng Y, Joung J G, et al. 2010. Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types. BMC Genomics, 11(1): 384–395. DOI:10.1186/1471-2164-11-384 |
[] | Hardenack S, Ye D, Saedler H, et al. 1994. Comparison of MADS box gene expression in developing male and female flowers of the dioecious plant white campion. Plant Cell, 6(12): 1775–1787. DOI:10.1105/tpc.6.12.1775 |
[] | Hultquist J, Dorweiler J. 2008. Feminized tassels of maize mop1 and ts1 mutants exhibit altered levels of miR156 and specific SBP-box genes. Planta, 229(1): 99–113. DOI:10.1007/s00425-008-0813-2 |
[] | Janoušek B, Široký J, Vyskot B. 1996. Epigenetic control of sexual phenotype in a dioecious plant, Melandrium album. Molecular Genetics and Genomics, 250(4): 483–490. |
[] | Kong Y, Zhu Y, Gao C, et al. 2013. Tissue-specific expression of SMALL AUXIN UP RNA41 differentially regulates cell expansion and root meristem patterning in Arabidopsis. Plant& Cell Physiology, 54(4): 609–621. |
[] | Lan T, Chen R, Li X, et al. 2008. Microdissection and painting of the W chromosome in Ginkgo biloba showed different labelling patterns. Bot Studies, 49: 33–37. |
[] | Liu J, Yin T, Ye N, et al. 2013. Transcriptome analysis of the differentially expressed genes in the male and female shrub willows (Salix suchowensis). PloS One, 8(4): e60181. DOI:10.1371/journal.pone.0060181 |
[] | Livak K J, Schmittgen T D. 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods, 25(4): 402–408. DOI:10.1006/meth.2001.1262 |
[] | Mao X, Tao C, Wei L. 2005. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics, 21(19): 3787–3793. DOI:10.1093/bioinformatics/bti430 |
[] | Martin A, Troadec C, Boualem A, et al. 2009. A transposon-induced epigenetic change leads to sex determination in melon. Nature, 461(7267): 1135–1138. DOI:10.1038/nature08498 |
[] | Mou Z, Fan W, Dong X. 2003. Inducers of plant systemic acquired resistance regulate NPR1 function through redox changes. Cell, 113(7): 935–944. DOI:10.1016/S0092-8674(03)00429-X |
[] | Park S Y, Fung P, Nishimura N, et al. 2009. Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science, 324(5930): 1068–1071. |
[] | Parkinson S E, Gross S M, Hollick J B. 2007. Maize sex determination and abaxial leaf fates are canalized by a factor that maintains repressed epigenetic states. Developmental Biology, 308(2): 462–473. DOI:10.1016/j.ydbio.2007.06.004 |
[] | Peng J, Richards D E, Hartley N M, et al. 1999. 'Green revolution' genes encode mutant gibberellin response modulators. Nature, 400(6741): 256–261. DOI:10.1038/22307 |
[] | Schmidt A, Wohrmann H J P, Raissig M T, et al. 2013. The polycomb group protein MEDEA and the DNA methyltransferase MET1 interact to repress autonomous endosperm development in Arabidopsis. Plant Journal, 73(5): 776–787. DOI:10.1111/tpj.2013.73.issue-5 |
[] | Solano R, Stepanova A, Chao Q, et al. 1998. Nuclear events in ethylene signaling:a transcriptional cascade mediated by ETHYLENE-INSENSITIVE3 and ETHYLENE-RESPONSE-FACTOR1. Genes& Development, 12(23): 3703–3714. |
[] | Song Y, Ma K, Bo W, et al. 2012. Sex-specific DNA methylation and gene expression in andromonoecious poplar. Plant Cell Reports, 31(8): 1393–1405. DOI:10.1007/s00299-012-1255-7 |
[] | Song Y, Ma K, Ci D, et al. 2013. Sexual dimorphic floral development in dioecious plants revealed by transcriptome, phytohormone, and DNA methylation analysis in Populus tomentosa. Plant Molecular Biology, 83(6): 559–576. DOI:10.1007/s11103-013-0108-2 |
[] | Storey J D, Tibshirani R. 2003. Statistical significance for genomewide studies. Proceedings of the National Academy of Sciences, 100(16): 9440–9445. DOI:10.1073/pnas.1530509100 |
[] | Trapnell C, Williams B A, Pertea G, et al. 2010. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology, 28(5): 511–515. DOI:10.1038/nbt.1621 |
[] | Turner J G, Ellis C, Devoto A. 2002. The jasmonate signal pathway. Plant Cell, 14(1): 153–164. |
[] | Wang Z, Gerstein M, Snyder M. 2009. RNA-Seq:a revolutionary tool for transcriptomics. Nature Reviews Genetics, 10(1): 57–63. DOI:10.1038/nrg2484 |
[] | Wu T, Qin Z, Zhou X, et al. 2010. Transcriptome profile analysis of floral sex determination in cucumber. Journal of Plant Physiology, 167(11): 905–913. DOI:10.1016/j.jplph.2010.02.004 |
[] | Yin T, DiFazio S P, Gunter L E, et al. 2008. Genome structure and emerging evidence of an incipient sex chromosome in Populus. Genome Research, 18(3): 422–430. DOI:10.1101/gr.7076308 |
[] | Young M D, Wakefield M J, Smyth G K, et al. 2010. Gene ontology analysis for RNA-seq:accounting for selection bias. Genome Biology, 11(2): R14. DOI:10.1186/gb-2010-11-2-r14 |
[] | Zhang Q, Li J, Sang Y, et al. 2015. Identification and characterization of microRNAs in Ginkgo biloba var. epiphylla Mak. PloS One, 10(5): e0127184. DOI:10.1371/journal.pone.0127184 |
[] | Zhang Y, Shewry P R, Jones H, et al. 2001. Expression of antisense SnRK1 protein kinase sequence causes abnormal pollen development and male sterility in transgenic barley. Plant Journal, 28(4): 431–441. |