Chromosome-scale assembly of the Dendrobium chrysotoxum genome enhances the understanding of orchid evolution

Abstract

As one of the largest families of angiosperms, the Orchidaceae family is diverse. Dendrobium represents the second largest genus of the Orchidaceae. However, an assembled high-quality genome of species in this genus is lacking. Here, we report a chromosome-scale reference genome of Dendrobium chrysotoxum, an important ornamental and medicinal orchid species. The assembled genome size of D. chrysotoxum was 1.37 Gb, with a contig N50 value of 1.54 Mb. Of the sequences, 95.75% were anchored to 19 pseudochromosomes. There were 30,044 genes predicted in the D. chrysotoxum genome. Two whole-genome polyploidization events occurred in D. chrysotoxum. In terms of the second event, whole-genome duplication (WGD) was also found to have occurred in other Orchidaceae members, which diverged mainly via gene loss immediately after the WGD event occurred; the first duplication was found to have occurred in most monocots (tau event). We identified sugar transporter (SWEET) gene family expansion, which might be related to the abundant medicinal compounds and fleshy stems of D. chrysotoxum. MADS-box genes were identified in D. chrysotoxum, as well as members of TPS and Hsp90 gene families, which are associated with resistance, which may contribute to the adaptive evolution of orchids. We also investigated the interplay among carotenoid, ABA, and ethylene biosynthesis in D. chrysotoxum to elucidate the regulatory mechanisms of the short flowering period of orchids with yellow flowers. The reference D. chrysotoxum genome will provide important insights for further research on medicinal active ingredients and breeding and enhances the understanding of orchid evolution.

References

[1] Christenhusz, M. J. M. & Byng, J. W. The number of known plants species in the world and its annual increase. Phytotaxa 261, 201-217 (2016).

[2] Zhang, G. et al. The Apostasia genome and the evolution of orchid. Nature 549, 379-383 (2017).

[3] Pridgeon, A. M. et al. Genera Orchidacearum: Epidendroideae (Part three) (Oxford University, 2014).

[4] Luo, J.-P., Deng, Y.-Y. & Zha, X.-Q. Mechanism of polysaccharides from Dendrobium huoshanense. on streptozotocin-induced diabetic cataract. Pharm. Biol. 46, 243-249 (2008).

[5] Leitch, I. J. et al. Genome size diversity in orchids: consequences and evolution. Ann. Bot. 104, 469-481 (2009).

[6] Ng, T. B. et al. Review of research on Dendrobium, a prized folk medicine. Appl. Microbiol. Biotechnol. 93, 1795-1803 (2012).

[7] Pan, L. H. et al. Comparison of hypoglycemic and antioxidative effects of polysaccharides from four different Dendrobium species. Int. J. Biol. Macromol. 64, 420-427 (2014).

[8] Tian, C.-C., Zha, X.-Q. & Luo, J.-P. A polysaccharide from Dendrobium huoshanense prevents hepatic inflammatory response caused by carbon tetrachloride. Biotechnol. Biotechnol. Equip. 29, 132-138 (2014).

[9] Huang, K. et al. Purification, characterization and biological activity of polysaccharides from Dendrobium officinale. Molecules 21, 701 (2016).

[10] Lu, J. et al. High-density genetic map construction and stem total polysaccharide content-related QTL exploration for Chinese endemic Dendrobium (Orchidaceae). Front. Plant Sci. 9, 398 (2018).

[11] Zhao, Y. et al. Antioxidant and anti-hyperglycemic activity of polysaccharide isolated from Dendrobium chrysotoxum Lindl. J. Biochem. Mol. Biol. 40, 670-677 (2007).

[12] Li, S. et al. Elution-extrusion counter-current chromatography separation of five bioactive compounds from Dendrobium chrysototxum Lindl. J. Chromatogr. A 1218, 3124-3128 (2011).

[13] Hu, J., Fan, W., Dong, F., Miao, Z. & Zhou, J. Chemical components of Dendrobium chrysotoxum. Chin. J. Chem. 30, 1327-1330 (2012).

[14] Yu, Z. et al. Dendrobium chrysotoxum Lindl. alleviates diabetic retinopathy by preventing retinal inflammation and tight junction protein decrease. J. Diabetes Res. 2015, 518317 (2015).

[15] Yu, Z. et al. Erianin inhibits high glucose-induced retinal angiogenesis via blocking ERK1/2-regulated HIF-1alpha-VEGF/VEGFR2 signaling pathway. Sci. Rep. 6, 34306 (2016).

[16] Wu, J. et al. Gigantol from Dendrobium chrysotoxum Lindl. binds and inhibits aldose reductase gene to exert its anti-cataract activity: an in vitro mechanistic study. J. Ethnopharmacol. 198, 255-261 (2017).

[17] Lim, V., Schneider, E., Wu, H. & Pang, I. H. Cataract preventive role of isolated phytoconstituents: findings from adecade of research. Nutrients 10, 1580 (2018).

[18] Ren, Z. Y. et al. Functional analysis of a novel C-glycosyltransferase in the orchid Dendrobium catenatum. Hortic. Res. 7, 111 (2020).

[19] Robles-Rivera, R. R. et al. Adjuvant therapies in diabetic retinopathy as an early approach to delay its progression: the importance of oxidative stress and inflammation. Oxid. Med. Cell Longev. 2020, 1-23 (2020).

[20] Gong, Y. et al. Erianin induces a jnk/sapk-dependent metabolic inhibition in human umbilical vein endothelial cells. In Vivo 18, 223-228 (2004).

[21] Cushman, M. et al. Synthesis and evaluation of stilbene and dihydrostilbene derivatives as potential anticancer agents that inhibit tubulin polymerization. J. Med. Chem. 34, 2579-2588 (1991).

[22] Ma, G. et al. Iuhibitory effects of Dendrobium chrysotoxutm and its constiuents on the mouse HePA and ESC. J. China Pharm. Univ. 25, 188-189 (1994).

[23] Gong, Y. Q. et al. In vivo and in vitro evaluation of erianin, a novel antiangiogenic agent. Eur. J. Cancer 40, 1554-1565 (2004).

[24] Zhang, G. et al. The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution. Sci. Rep. 6, 19029 (2016).

[25] Hu, M. J. et al. Chromosome-scale assembly of the Kandelia obovata genome. Hortic. Res. 7, 75 (2020).

[26] Kang, M. et al. A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine: an Isatis genome. Hortic. Res. 7, 18 (2020).

[27] Chen, F. et al. The sequenced angiosperm genomes and genome databases. Front. Plant Sci. 9, 418 (2018).

[28] Chen, F. et al. Genome sequences of horticultural plants: past, present, and future. Hortic. Res. 6, 112 (2019).

[29] Yan, L. et al. The genome of Dendrobium officinale illuminates the biology of the important traditional Chinese orchid herb. Mol. Plant 8, 922-934 (2015).

[30] Si, C. et al. DoRWA3 from Dendrobium officinale plays an essential role in acetylation of polysaccharides. Int. J. Mol. Sci. 21, 6250 (2020).

[31] Zheng, S. G. et al. Genome-wide researches and applications on Dendrobium. Planta 248, 769-784 (2018).

[32] Simão, F. A. et al. BUSCO online supplementary information: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[33] Cai, J. et al. The genome sequence of the orchid Phalaenopsis equestris. Nat. Genet. 47, 65-72 (2015).

[34] Yuan, Y. et al. The Gastrodia elata genome provides insights into plant adaptation to heterotrophy. Nat. Commun. 9, 1615 (2018).

[35] Sena, J. S. et al. Evolution of gene structure in the conifer Picea glauca: a comparative analysis of the impact of intron size. BMC Plant Biol. 14, 95 (2014).

[36] De La Torre, A. R. et al. Insights into conifer giga-genomes. Plant Physiol. 166, 1724-1732 (2014).

[37] Castillo-Davis, C. I. et al. Selection for short introns in highly expressed genes. Nat. Genet. 31, 415 (2002).

[38] Keane, P. A. & Seoighe, C. Intron length coevolution across mammalian genomes. Mol. Biol. Evol. 33, 2682-2691 (2016).

[39] Swinburne, I. A. & Silver, P. A. Intron delays and transcriptional timing during development. Dev. Cell 14, 324-330 (2008).

[40] Artieri, C. G. & Fraser, H. B. Transcript length mediates developmental timing of gene expression across Drosophila. Mol. Biol. Evol. 31, 2879-2889 (2014).

[41] Lin, I. W. et al. Nectar secretion requires sucrose phosphate synthases and the sugar transporter SWEET9. Nature 508, 546-549 (2014).

[42] Chen, L. Q. et al. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science 335, 207-211 (2012).

[43] Sun, M. X. et al. Arabidopsis RPG1 is important for primexine deposition and functions redundantly with RPG2 for plant fertility at the late reproductive stage. Plant Reprod. 26, 83-91 (2013).

[44] Chen, L. Q. et al. A cascade of sequentially expressed sucrose transporters in the seed coat and endosperm provides nutrition for the Arabidopsis embryo. Plant Cell 27, 607-619 (2015).

[45] Jiao, Y., Li, J., Tang, H. & Paterson, A. H. Integrated syntenic and phylogenomic analyses reveal an ancient genome duplication in monocots. Plant Cell 26, 2792-2802 (2014).

[46] Tang, H. et al. jcvi: JCVI Utility Libraries (Zenodo, 2015). https://doi.org/10.5281/zenodo.31631.

[47] Sang, X. et al. CHIMERIC FLORAL ORGANS1, encoding a monocot-specific MADS box protein, regulates floral organ identity in rice. Plant Physiol. 160, 788-807 (2012).

[48] Tapia-López, R. et al. An AGAMOUS-related MADS-box gene, XAL1 (AGL12), regulates root meristem cell proliferation and flowering transition in Arabidopsis. Plant Physiol. 146, 1182-1192 (2008).

[49] Parenicová, L. et al. Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: new openings to the MADS world. Plant Cell 15, 1538-1551 (2003).

[50] Masiero, S. et al. The emerging importance of type I MADS box transcription factors for plant reproduction. Plant Cell 23, 865-872 (2011).

[51] Niu, S. C. et al. Morphological type identification of self-incompatibility in Dendrobium and its phylogenetic evolution pattern. Int. J. Mol. Sci. 19, 2595 (2018).

[52] Wu, R. H. Studies on the Germplasm Resources of Dendrobium in China and their Genetic Relationships. Doctoral dissertation, Chin. Acad. For. (2007).

[53] Watkins, J. L. & Pogson B. J. Prospects for carotenoid biofortification targeting retention and catabolism. Trends Plant Sci. 25, 501-512 (2020).

[54] Johnson, P. R. & Ecker, J. R. The ethylene gas signal transduction pathway: a molecular perspective. Annu. Rev. Genet. 32, 227 (1998).

[55] Yin, C. C. et al. Ethylene responses in rice roots and coleoptiles are differentially regulated by a carotenoid isomerase-mediated abscisic acid pathway. Plant Cell 27, 1061-1081 (2015).

[56] Yu, Z. et al. Genome-wide identification and expression profile of TPS gene family in Dendrobium officinale and the role of DoTPS10 in linalool biosynthesis. Int. J. Mol. Sci. 21, 5419 (2020).

[57] Jiang, S. Y., Jin, J. J., Sarojam, R. & Ramachandran, S. A comprehensive survey on the terpene synthase gene family provides new insight into its evolutionary patterns. Genome Biol. Evol. 11, 2078-2098 (2019).

[58] Chen, F. The family of terpene synthases in plants: A mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 66, 212-229 (2011).

[59] Tholl, D. Biosynthesis and biological functions of terpenoids in plants. Adv. Biochem. Eng. Biotechnol. 148, 63-106 (2015).

[60] Bohlmann, J., Meyer-Gauen, G. & Croteau, R . Plant terpenoid synthases: molecular biology and phylogenetic analysis. Proc. Natl Acad. Sci. USA 95, 4126-4133 (1998).

[61] Jiang, S. Y. et al. Comprehensive survey on the terpene synthase gene family provides new insight into its evolutionary patterns. Genome Biol. Evol. 11, 2078-2098 (2019).

[62] Huang, X. L. Research on molecular regulation mechanism of formation of floral color and floral fragrance of Dendrobium chrysotoxum based on transcriptome sequencing. Doctoral dissertation, Chin. Acad. For. (2019).

[63] Ranallo-Benavidez, T. R., Jaron, K. S. & Schatz, M. C. GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat. Commun. 11, 1-10 (2020).

[64] Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722-736 (2017).

[65] Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).

[66] Chen, Y. et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. GigaScience 7, 120 (2017).

[67] Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 3, 95-98 (2016).

[68] Dudchenko, O. et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356, 92-95 (2017).

[69] Jurka, J. et al. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 110, 462-467 (2005).

[70] Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. Bioinformatics 21, i351-i358 (2005).

[71] Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 35, W265-W268 (2007).

[72] Benson, G. et al. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573-580 (1999).

[73] Slater, G. S. & Birney, E. Automated generation of heuristics for biological sequence comparison. BMC Bioinform. 6, 31 (2005).

[74] Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34, W435-W439 (2006).

[75] Johnson, A. D. et al. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 24, 2938-2939 (2008).

[76] Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290-295 (2015).

[77] Kim, D. et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907-915 (2019).

[78] Holt, C. & Yandell, M. MAKER2: an annotation pipeline and genome database management tool for second generation genome projects. BMC Bioinform. 12, 491 (2011).

[79] Kent, W. J. BLAT-the BLAST-like alignment tool. Genome Res. 12, 656-664 (2002).

[80] Boeckmann, B. et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 31, 365-370 (2003).

[81] Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236-1240 (2014).

[82] Ashburner, M. et al. Gene Ontology: tool for the unification of biology. Nat. Genet. 25, 25-29 (2000).

[83] Lowe, T. M. & Eddy, S. R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955-964 (1997).

[84] Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933-2935 (2013).

[85] Wang, Y. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49 (2012).

[86] Yang, Z. PAML: a program package for phylogenetic analysis by maximum likelihood. Bioinformatics 13, 555-556 (1997).

[87] De Bie, T. et al. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-1271 (2006).

[88] El-Gebali, S. et al. The Pfam protein families database in 2019. Nucleic Acids Res. 47, D427-D432 (2019).

[89] Eddy, S. R. Accelerated profile HMM searches. PLoS Comput. Biol. 7, e1002195 (2011).

[90] Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772-780 (2013).

[91] Guindon, S., Lethiec, F., Duroux, P. & Gascuel, O. PHYML Online-a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 33, W557-W559 (2005).

[92] Wu, S. S. et al. The genome sequence of star fruit (Averrhoa carambola). Hortic. Res. 7, 95 (2020).

[93] Chen, S. P. et al. The Phoebe genome sheds light on the evolution of magnoliids. Hortic. Res. 7, 146 (2020).

[94] Letunic, I., Doerks, T. & Bork, P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 43, D257-D260 (2015).

[95] Kumar, S. et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547-1549 (2018).

[96] Finkelstein, R. Abscisic Acid synthesis and response. Arabidopsis Book 11, e0166- e0166 (2013).

[97] Lin, Z., Zhong, S. & Grierson, D. Recent advances in ethylene research. J. Exp. Bot. 60, 3311-3336 (2009).

[98] Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307-321 (2010).

[99] Jones, D. T., Taylor, W. R. & Thornton, J. M. The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275-282 (1992).

[100] Patro, R. et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417-419 (2017).

[101] Nisar, N. et al. Carotenoid metabolism in plants. Mol. Plant 8, 68-82 (2015).

[102] Sun, T. et al. Carotenoid metabolism in plants: the role of plastids. Mol. Plant 11, 58-74 (2017).