Generation and classification of transcriptomes in two Croomia species and molecular evolution of CYC/TB1 genes in Stemonaceae

doi:10.1016/j.pld.2018.11.006

Abstract

Abstract: The genus Croomia (Stemonaceae) is an excellent model for studying the evolution of the Eastern Asia (EA) eEastern North America (ENA) floristic disjunction and the genetic mechanisms of floral zygomorphy formation. In addition to the presence of both actinomorphic and zygomorphic flowers within the genus, species are disjunctively distributed between EA and ENA. However, due to the limited availability of genomic resources, few studies of Croomia have examined these questions. In this study, we sequenced the floral and leaf transcriptomes of the zygomorphic flowered Croomia heterosepala and the actinomorphic flowered Croomia japonica, and used comparative genomic approaches to investigate the transcriptome evolution of the two closely related species. The sequencing and de novo assembly of transcriptomes from flowers of C. heterosepala (ChFlower), flowers of C. japonica (CjFlower), and leaves of C. japonica (CjLeaf) yielded 57,193, 62,131 and 64,448 unigenes, respectively. In addition, estimation of Ka/Ks ratios for 11,566 potential orthologous groups between ChFlower and CjFlower revealed that only six pairs had Ka/Ks ratios significantly greater than 1 and are likely under positive selection. A total of 429 single copy nuclear genes (SCNGs) and 21,460 expression sequence tags-simple sequence repeats (ESTSSRs) were identified in this study. Specifically, we identified seven CYC/TB1-like genes from Stemonaceae. Phylogenetic and molecular evolution analyses indicated that these CYC/TB1-like genes formed a monophyletic clade (SteTBL1) and were subject to strong purifying selection. The shifts of floral symmetry in Stemonaceae do not appear to be correlated with TBL copy number.

Key words: Croomia, Transcriptome, SCNGs, EST-SSRs, Flower symmetry, CYC/TB1

Ruisen Lu, Wuqin Xu, Qixiang Lu, Pan Li, Jocelyn Losh, Faiza Hina, Enxiang Li, Yingxiong Qiu. Generation and classification of transcriptomes in two Croomia species and molecular evolution of CYC/TB1 genes in Stemonaceae[J]. Plant Diversity, 2018, 40(06): 253-264.

References

Abascal, F., Zardoya, R., Posada, D., 2005. ProtTest:selection of best-fit models of protein evolution. Bioinformatics 21, 2104-2105.
Ai, B., Gao, Y., Zhang, X., Tao, J., Kang, M., Huang, H., 2015. Comparative transcriptome resources of eleven Primulina species, a group of ‘stone plants’ from a biodiversity hot spot. Mol. Ecol. Resour. 15, 619-632.
Akashi, H.D., Cádiz, D.A., Shigenobu, S., Makino, T., Kawata, M., 2016. Differentially expressed genes associated with adaptation to different thermal environments in three sympatric Cuban Anolis lizards. Mol. Ecol. 25, 2273-2285.
Altschul, S.F., Altschul, S.F., Gish, W., Gish, W., Miller, W., Miller, W., et al., 1990. Basic local alignment search tool. J. Mol. Biol. 215, 403-410.
Andersen, J.R., Lübberstedt, T., 2003. Functional markers in plants. Trends Plant Sci. 8, 554-560.
Ashburner, M., Ball, C.A., Blake, J.A., Botstein, D., Butler, H., Cherry, J.M., et al., 2000. Gene ontology:tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25, 25.
Bartlett, M.E., Specht, C.D., 2011. Changes in expression pattern of the TEOSINTE BRANCHED1-like genes in the Zingiberales provide a mechanism for evolutionary shifts in symmetry across the order. Am. J. Bot. 98, 227-243.
Bustamante, C.D., Fledel-Alon, A., Williamson, S., Nielsen, R., Hubisz, M.T., Glanowski, S., et al., 2005. Natural selection on protein-coding genes in the human genome. Nature 437, 1153-1157.
Caddick, L.R., Rudall, P.J., Wilkin, P., Hedderson, T.A.J., Chase, M.W., 2002. Phylogenetics of Dioscoreales based on combined analyses of morphological and molecular data. Bot. J. Linn. Soc. 138, 123-144.
Chapuis, M., Estoup, A., 2007. Microsatellite null alleles and estimation of population differentiation. Mol. Biol. Evol. 24, 621-631.
Chase, M.W., Christenhusz, M.J.M., Fay, M.F., Byng, J.W., Judd, W.S., Soltis, D.E., et al., 2016. An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants:APG IV. Bot. J. Linn. Soc. 181, 1-20.
Chen, G., Brecker, L., Felsinger, S., Cai, X., Kongkiatpaiboon, S., Schinnerl, J., et al., 2017. Morphological and chemical variation of Stemona tuberosa from southern China e evidence for heterogeneity of this medicinal plant species. Plant Biol. 19, 835-842.
Chen, L., Cao, Y., Yuan, N., Nakamura, K., Wang, G., Qiu, Y., 2015. Characterization of transcriptome and development of novel EST-SSR makers based on nextgeneration sequencing technology in Neolitsea sericea (Lauraceae) endemic to East Asian land-bridge islands. Mol. Breed. 35, 1-15.
Christenhusz, M., Byng, J., 2016. The number of known plants species in the world and its annual increase. Phytotaxa 261, 201-217.
Citerne, H., Reyes, E., Le, G.M., Delannoy, E., Simonnet, F., Sauquet, H., et al., 2017. Characterization of CYCLOIDEA-like genes in Proteaceae, a basal eudicot family with multiple shifts in floral symmetry. Ann. Bot. 119, 367-378.
Clark, M.S., Thorne, M.A.S., Vieira, F.A., Cardoso, J.C.R., Power, D.M., Peck, L.S., 2010. Insights into shell deposition in the Antarctic bivalve Laternula elliptica:gene discovery in the mantle transcriptome using 454 pyrosequencing. BMC Genomics 11, 362.
Conesa, A., Götz, S., García, J.M., Teról, J., Tal on, M., Robles, M., 2005. Blast2GO:a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674-3676.
Dempster, A.P., Laird, N.M., Rubin, D.B., 1977. Maximum likelihood from incomplete data via the EM algorithm. J. Roy. Stat. Soc. B 39, 1-38.
Dieringer, D., Schlötterer, C., 2003. MICROSATELLITE ANALYSER (MSA):a platform independent analysis tool for large microsatellite data sets. Mol. Ecol. Notes 3, 167-169.
Dijk, E.L.V., Auger, H., Yan, J., Thermes, C., 2014. Ten years of next-generation sequencing technology. Trends Genet. 30, 418-426.
Dolman, G., Joseph, L., 2016. Multi-locus sequence data illuminate demographic drivers of Pleistocene speciation in semi-arid southern Australian birds (Cinclosoma spp.). BMC Evol. Biol. 16, 226.
Duarte, J.M., Wall, P.K., Edger, P.P., Landherr, L.L., Ma, H., Pires, J.C., et al., 2010. Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels. BMC Evol. Biol. 10, 61.
Duran, C., Appleby, N., Edwards, D., Batley, J., 2009. Molecular genetic markers:discovery, applications, data storage and visualization. Curr. Bioinf. 4, 16-27.
Estill, J.C., Cruzan, M.B., 2001. Phytogeography of rare plant species endemic to the Southeastern United States. Castanea 66, 3-23.
Fang, M., Fu, C., Fu, C., Zhu, Y., Naiki, A., Li, E., 2013. Development of microsatellite markers for Croomia japonica and cross-amplification in its congener. Sci. Hortic. 161, 228-232.
Greger, H., 2006. Structural relationships, distribution and biological activities of Stemona alkaloids. Planta Med. 72, 99-113.
Guindon, S., Dufayard, J., Lefort, V., Anisimova, M., Hordijk, W., Gascuel, O., 2010. New algorithms and methods to estimate maximum-likelihood phylogenies:assessing the performance of PhyML 3.0. Syst. Biol. 59, 307-321.
Guo, R., Mao, Y., Cai, J., Wang, J., Wu, J., Qiu, Y., 2014. Characterization and crossspecies transferability of ESTeSSR markers developed from the transcriptome of Dysosma versipellis (Berberidaceae) and their application to population genetic studies. Mol. Breed. 34, 1733-1746.
Haas, B.J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P.D., Bowden, J., et al., 2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494-1512.
Horn, S., Pabón-Mora, N., Theuß, V.S., Busch, A., Zachgo, S., 2015. Analysis of the CYC/TB1 class of TCP transcription factors in basal angiosperms and magnoliids. Plant J. 81, 559-571.
Hoshino, Y., Igarashi, T., Ohshima, M., Shinoda, K., Murata, N., Kanno, A., et al., 2014. Characterization of CYCLOIDEA-like genes in controlling floral zygomorphy in the monocotyledon Alstroemeria. Sci. Hortic. 169, 6-13.
Howarth, D.G., Donoghue, M.J., 2005. Duplications in CYC-like genes from dipsacales correlate with floral form. Int. J. Plant Sci. 166, 357-370.
Howarth, D.G., Martins, T., Chimney, E., Donoghue, M.J., 2011. Diversification of CYCLOIDEA expression in the evolution of bilateral flower symmetry in Caprifoliaceae and Lonicera (Dipsacales). Ann. Bot. 107, 1521-1532.
Huang, D., Zhang, Y., Jin, M., Li, H., Song, Z., Wang, Y., et al., 2014. Characterization and high cross-species transferability of microsatellite markers from the floral transcriptome of Aspidistra saxicola (Asparagaceae). Mol. Ecol. Resour. 14, 569-577.
Huang, J., Gao, Y., Jia, H., Zhang, Z., 2016. Characterization of the teosinte transcriptome reveals adaptive sequence divergence during maize domestication. Mol. Ecol. Resour. 16, 1465-1477.
Huang, L.K., Yan, H.D., Zhao, X.X., Zhang, X.Q., Wang, J., Frazier, T., et al., 2015. Identifying differentially expressed genes under heat stress and developing molecular markers in orchardgrass (Dactylis glomerata L.) through transcriptome analysis. Mol. Ecol. Resour. 15, 1497-1509.
Iseli, C., Jongeneel, C.V., Bucher, P., 1999. ESTScan:a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 138-148.
Jabbour, F., Cossard, G., Le, G.M., Sannier, J., Nadot, S., Damerval, C., 2014. Specific duplication and dorsoventrally asymmetric expression patterns of cycloidealike genes in zygomorphic species of Ranunculaceae. PLoS One 9, e95727.
Kalinowski, S.T., Taper, M.L., Marshall, T.C., 2007. Revising how the computer program cervus accommodates genotyping error increases success in paternity assignment. Mol. Ecol. 16, 1099-1106.
Kanehisa, M., Araki, M., Goto, S., Hattori, M., Hirakawa, M., Itoh, M., et al., 2008. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, 480-484.
Kawakami, T., Backström, N., Burri, R., Husby, A., Olason, P., Rice, A.M., et al., 2014. Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single-nucleotide polymorphism array. Mol. Ecol. Resour. 14, 1248-1260.
Koilkonda, P., Sato, S., Tabata, S., Shirasawa, K., Hirakawa, H., Sakai, H., et al., 2012. Large-scale development of expressed sequence tag-derived simple sequence repeat markers and diversity analysis in Arachis spp. Mol. Breed. 30, 125-138.
Kongkiatpaiboon, S., Schinnerl, J., Felsinger, S., Keeratinijakal, V., Vajrodaya, S., Gritsanapan, W., et al., 2011. Structural relationships of stemona alkaloids:assessment of species-specific accumulation trends for exploiting their biological activities. J. Nat. Prod. 74, 1931-1938.
Kumar, S., Stecher, G., Tamura, K., 2016. MEGA7:molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870-1874.
Kumpatla, S.P., Mukhopadhyay, S., 2005. Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species. Genome 48, 985-998.
Leavitt, S.D., Esslinger, T.L., Divakar, P.K., Lumbsch, H.T., 2012. Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea(Parmeliaceae, Ascomycota) and Pleistocene population expansions. BMC Evol. Biol. 12, 176.
Lee, S., Xiao, C., Pei, S., 2008. Ethnobotanical survey of medicinal plants at periodic markets of Honghe Prefecture in Yunnan Province, SW China. J. Ethnopharmacol. 117, 362-377.
Li, E., Yi, S., Qiu, Y., Guo, J., Comes, H.P., Fu, C., 2008. Phylogeography of two East Asian species in Croomia (Stemonaceae) inferred from chloroplast DNA and ISSR fingerprinting variation. Mol. Phylogenet. Evol. 49, 702-714.
Li, L., Stoeckert, C.J., Roos, D.S., 2003. OrthoMCL:identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178-2189.
Li, Y., Korol, A.B., Fahima, T., Nevo, E., 2004. Microsatellites within genes:structure, function, and evolution. Mol. Biol. Evol. 21, 991-1007.
Li, Z., Zou, J., Mao, K., Lin, K., Li, H., Liu, J., et al., 2012. Population genetic evidence for complex evolutionary histories of four high altitude juniper species in the Qinghai-Tibetan plateau. Evolution 66, 831-845.
Lin, G., Lin, L., Pak-Ho, L.H., Zhu, J., Tang, C., Ke, C., et al., 2008. Croomine-and tuberostemonine-type alkaloids from roots of Stemona tuberosa and their antitussive activity. Tetrahedron 64, 10155-10161.
Lin, L., Zhong, Q., Cheng, T., Tang, C., Ke, C., Lin, G., et al., 2006. Stemoninines from the roots of Stemona tuberosa. J. Nat. Prod. 69, 1051-1054.
Liu, Z., Ma, L., Nan, Z., Wang, Y., 2013. Comparative transcriptional profiling provides insights into the evolution and development of the zygomorphic flower of Vicia sativa (Papilionoideae). PLoS One 8, e57338.
Logacheva, M.D., Kasianov, A.S., Vinogradov, D.V., Samigullin, T.H., Gelfand, M.S., Makeev, V.J., et al., 2011. De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum). BMC Genomics 12, 30.
Lukens, L., Doebley, J., 2001. Molecular evolution of the teosinte branched gene among maize and related grasses. Mol. Biol. Evol. 18, 627-638.
Ma, C., Xin, M., Feldmann, K.A., Wang, X., 2014. Machine learning-based differential network analysis:a study of stress-responsive transcriptomes in Arabidopsis. Plant Cell 26, 520-537.
Mao, Y., Zhang, Y., Xu, C., Qiu, Y., 2016. Comparative transcriptome resources of two Dysosma species (Berberidaceae) and molecular evolution of the CYP719A gene in Podophylloideae. Mol. Ecol. Resour. 16, 228-241.
McCormack, J.E., Heled, J., Delaney, K.S., Peterson, A.T., Knowles, L.L., 2011. Calibrating divergence times on species trees versus gene trees:implication for speciation history of Aphelocoma Jays. Evolution 65, 184-202.
Morgante, M., Hanafey, M., Powell, W., 2002. Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat. Genet. 30, 194-200.
Okuyama, S., 1944. On the Japanese species of Croomia. J. Jpn. Bot. 20, 31-32.
Ohwi, J., 1965. Croomia. Flora of Japan. Smithsonian Institution, Washington, p. 279.
Parchman, T.L., Geist, K.S., Grahnen, J.A., Benkman, C.W., Buerkle, C.A., 2010. Transcriptome sequencing in an ecologically important tree species:assembly, annotation, and marker discovery. BMC Genomics 11, 180.
Pertea, G., Huang, X., Liang, F., Antonescu, V., Sultana, R., Karamycheva, S., et al., 2003. TIGR gene indices clustering tools (TGICL):a software system for fast clustering of large EST datasets. Bioinformatics 19, 651-652.
Poncet, V., Rondeau, M., Tranchant, C., Cayrel, A., Hamon, S., Kochko, A., et al., 2006. SSR mining in coffee tree EST databases:potential use of ESTeSSRs as markers for the Coffea genus. Mol. Genet. Genom. 276, 436-449.
Preston, J.C., Hileman, L.C., 2009. Developmental genetics of floral symmetry evolution. Trends Plant Sci. 14, 147-154.
Preston, J.C., Hileman, L.C., 2012. Parallel evolution of TCP and B-class genes in Commelinaceae flower bilateral symmetry. EvoDevo 3, 6.
Qiu, L., Yang, C., Tian, B., Yang, J., Liu, A., 2010. Exploiting EST databases for the development and characterization of EST-SSR markers in castor bean (Ricinus communis L.). BMC Plant Biol. 10, 278.
Reardon, W., Fitzpatrick, D.A., Fares, M.A., Nugent, J.M., 2009. Evolution of flower shape in Plantago lanceolata. Plant Mol. Biol. 71, 241-250.
Reyes, E., Sauquet, H., Nadot, S., 2016. Perianth symmetry changed at least 199 times in angiosperm evolution. Taxon 65, 945-964.
Rogers, A.R., Harpending, H., 1992. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552-569.
Ronquist, F., Huelsenbeck, J.P., 2003. MrBayes 3:Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 1572-1574.
Rousset, F., 2008. Genepop'007:a complete re-implementation of the genepop software for windows and linux. Mol. Ecol. Resour. 8, 103-106.
Rudall, P.J., Cunniff, J., Wilkin, P., Caddick, L.R., 2005. Evolution of dimery, pentamery and the monocarpellary condition in the monocot family Stemonaceae (Pandanales). Taxon 54, 701-711.
Schluepen, C., Malito, E., Marongiu, A., Schirle, M., McWhinnie, E., Lo Surdo, P., et al., 2013. Mining the bacterial unknown proteome:identification and characterization of a novel family of highly conserved protective antigens in Staphylococcus aureus. Biochem. J. 455, 273-284.
Swanson, W.J., Wong, A., Wolfner, M.F., Aquadro, C.F., 2004. Evolutionary expressed sequence tag analysis of Drosophila female reproductive tracts identifies genes subjected to positive selection. Genetics 168, 1457-1465.
Teasdale, L.C., Köhler, F., Murray, K.D., O'Hara, T., Moussalli, A., 2016. Identi fication and qualification of 500 nuclear, single-copy, orthologous genes for the Eupulmonata (Gastropoda) using transcriptome sequencing and exon capture. Mol. Ecol. Resour. 16, 1107-1123.
Thompson, J.D., Higgins, D.G., Gibson, T.J., 1994. ClustalW:improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673-4680.
Varshney, R.K., Graner, A., Sorrells, M.E., 2005. Genic microsatellite markers in plants:features and applications. Trends Biotechnol. 23, 48-55.
Wang, S., Xie, Y., 2004. Red List. In:China Species Red List, vol. 1. Higher Education Press, Beijing.
Wen, J., Egan, A.N., Dikow, R.B., Zimmer, E.A., 2015. Utility of transcriptome sequencing for phylogenetic inference and character evolution. In:Elvira, H., Marc, S.A. (Eds.), Next-generation Sequencing in Plant Systematics, pp. 1-42.
Xiang, Y., Huang, C., Hu, Y., Wen, J., Li, S., Yi, T., et al., 2016. Evolution of Rosaceae fruit types based on nuclear phylogeny in the context of geological times and genome duplication. Mol. Biol. Evol. 34, 262-281.
Yang, Z., 2007. PAML 4:phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591.
Yang, Z., Nielsen, R., 2000. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evol. 17, 32-43.
Ye, J., Fang, L., Zheng, H., Zhang, Y., Chen, J., Zhang, Z., et al., 2006. WEGO:a web tool for plotting GO annotations. Nucleic Acids Res. 34, W293-W297.
Yuan, Z., Gao, S., Xue, D., Luo, D., Li, L., Ding, S., et al., 2009. RETARDED PALEA1 controls palea development and floral zygomorphy in rice. Plant Physiol. 149, 235.
Zeng, L., Zhang, N., Zhang, Q., Endress, P.K., Huang, J., Ma, H., 2017. Resolution of deep eudicot phylogeny and their temporal diversification using nuclear genes from transcriptomic and genomic datasets. New Phytol. 214, 1338-1354.
Zhang, J., Franks, R., Liu, X., Kang, M., Keebler, J., Schaff, J., et al., 2013a. De novo sequencing, characterization, and comparison of inflorescence transcriptomes of Cornus canadensis and C. florida (Cornaceae). PLoS One 8, e82674.
Zhang, L., Yan, H., Wu, W., Yu, H., Ge, X., 2013b. Comparative transcriptome analysis and marker development of two closely related Primrose species (Primula poissonii and Primula wilsonii). BMC Genomics 14, 329.
Zhang, W., Kramer, E.M., Davis, C.C., Donoghue, M.J., 2010. Floral symmetry genes and the origin and maintenance of zygomorphy in a plant-pollinator mutualism. P. Natl. Acad. Sci. USA. 107, 6388-6393.
Zhang, Z., Li, J., Zhao, X., Wang, J., Wong, G.K.S., Yu, J., 2006. KaKs_Calculator:calculating Ka and Ks through model selection and model averaging. Dev. Reprod. Biol. 4, 259-263.
Zhou, S., Yan, B., Li, F., Zhang, J., Ma, H., Liu, W., et al., 2017. RNA-seq analysis provides the first insights into the phylogenetic relationship and interspecific variation between Agropyron cristatum and Wheat. Front. Plant Sci. 8, 1644.
Zhu, S., Ding, Y., Yap, Z., Qiu, Y., 2016. De novo assembly and characterization of the floral transcriptome of an economically important tree species, Lindera glauca(Lauraceae), including the development of EST-SSR markers for population genetics. Mol. Biol. Rep. 43, 1243-1250.

[1]	Yu-Wen Zhang, Yu-Cen Shi, Shi-Bao Zhang. Metabolic and transcriptomic analyses elucidate a novel insight into the network for biosynthesis of carbohydrate and secondary metabolites in the stems of a medicinal orchid Dendrobium nobile [J]. Plant Diversity, 2023, 45(03): 326-336.
[2]	Ting Tang, Faqing Tao, Weiqi Li. Characterisation of manganese toxicity tolerance in Arabis paniculata [J]. Plant Diversity, 2021, 43(02): 163-172.
[3]	Jun-Chu Peng, Xiang-Guang Ma, Yue-Hua Wang, Hang Sun. New insights into the evolutionary history of Megacodon: Evidence from a newly discovered species [J]. Plant Diversity, 2020, 42(03): 198-208.
[4]	Yonglu Wei, Jianpeng Jin, Xiani Yao, Chuqiao Lu, Genfa Zhu, Fengxi Yang. Transcriptome Analysis Reveals Clues into leaf-like flower mutant in Chinese orchid Cymbidium ensifolium [J]. Plant Diversity, 2020, 42(02): 92-101.
[5]	Yun Jia, Ji-Qing Bai, Mi-Li Liu, Zhen-Fang Jiang, Yan Wu, Min-Feng Fang, Zhong-Hu Li. Transcriptome analysis of the endangered Notopterygium incisum: Cold-tolerance gene discovery and identification of EST-SSR and SNP markers [J]. Plant Diversity, 2019, 41(01): 1-6.
[6]	Tao Liu a, Xiaoxian Li b, Shiqing Xie a, Ling Wang a, Shengchao Yang a, . RNA-seq analysis of Paris polyphylla var. yunnanensis* roots identified candidate genes for saponin synthesis [J]. Plant Diversity, 2016, 38(03): 163-170.
[7]	ZhAO Lei, Zachary LARSONRABIN, CHEN Si-Yun, GUO Zhen-Hua. Comparing De Novo Transcriptome Assemblers Using Illumina RNASeq Reads [J]. Plant Diversity, 2012, 34(5): 487-501.