High-quality de novo assembly of the Eucommia ulmoides haploid genome provides new insights into evolution and rubber biosynthesis

Abstract

We report the acquisition of a high-quality haploid chromosome-scale genome assembly for the first time in a tree species, Eucommia ulmoides, which is known for its rubber biosynthesis and medicinal applications. The assembly was obtained by applying PacBio and Hi–C technologies to a haploid that we specifically generated. Compared to the initial genome release, this one has significantly improved assembly quality. The scaffold N50 (53.15 MB) increased 28-fold, and the repetitive sequence content (520 Mb) increased by 158.24 Mb, whereas the number of gaps decreased from 104,772 to 128. A total of 92.87% of the 26,001 predicted protein-coding genes identified with multiple strategies were anchored to the 17 chromosomes. A new whole-genome duplication event was superimposed on the earlier γ paleohexaploidization event, and the expansion of long terminal repeats contributed greatly to the evolution of the genome. The more primitive rubber biosynthesis of this species, as opposed to that in Hevea brasiliensis, relies on the methylerythritol-phosphate pathway rather than the mevalonate pathway to synthesize isoprenyl diphosphate, as the MEP pathway operates predominantly in trans-polyisoprene-containing leaves and central peels. Chlorogenic acid biosynthesis pathway enzymes were preferentially expressed in leaves rather than in bark. This assembly with higher sequence contiguity can foster not only studies on genome structure and evolution, gene mapping, epigenetic analysis and functional genomics but also efforts to improve E. ulmoides for industrial and medical uses through genetic engineering.

References

[1] Baboo, M., Dixit, M., Sharma, K. & Saxena, N. S. Mechanical and thermal characterization of cis-polyisoprene and trans-polyisoprene blends. Polym. Bull. 66, 661-672 (2011).

[2] Mooibroek, H. & Cornish, K. Alternative sources of natural rubber. Appl. Microbiol. Biotechnol. 53, 355-365 (2000).

[3] Van Beilen, J. B. & Poirier, Y. Establishment of new crops for the production of natural rubber. Trends Biotechnol. 25, 522-529 (2007).

[4] Lieberei, R. South American leaf blight of the rubber tree (Hevea spp.): new steps in plant domestication using physiological features and molecular markers. Ann. Bot. 100, 1125-1142 (2007).

[5] Tang, C. et al. The rubber tree genome reveals new insights into rubber production and species adaptation. Nat. Plants 2, 16073 (2016).

[6] Bamba, T., Fukusaki, E., Nakazawa, Y. & Kobayashi, A. In-situ chemical analyses of trans-polyisoprene by histochemical staining and Fourier transform infrared microspectroscopy in a rubber-producing plant, Eucommia ulmoides Oliver. Planta 215, 934-939 (2002).

[7] Du, H. Y., Hu, W. Z. & Yu, R. The Report on Development of China’s Eucommia Rubber Resources and Industry (2014-2015) (Social Sciences Academic Press, 2015).

[8] Kent, E. G. & Swinney, F. B. Properties and applications of trans-1,4-polyisoprene. Ind. Eng. Chem. Prod. Res. Dev. 5, 134-138 (1966).

[9] Enoki, M., Doi, Y. & Iwata, T. Oxidative degradation of cis- and trans-1,4-polyisoprenes and vulcanized natural rubber with enzyme-mediator systems. Biomacromolecules 4, 314-320 (2003).

[10] Rose, K. & Steinbuchel, A. Biodegradation of natural rubber and related compounds: recent insights into a hardly understood catabolic capability of microorganisms. Appl. Environ. Microbiol. 71, 2803-2812 (2005).

[11] Schlesinger, W. & Leeper, H. M. Chicle - cis- and trans-polyisoprenes from a single plant species. Ind. Eng. Chem. 43, 398-403 (1951).

[12] Roth, W. B., Carr, M. E., Davis, E. A. & Bagby, M. O. New sources of gutta-percha in Garrya flavescens and G. wrightii. Phytochemistry 24, 183-184 (1985).

[13] Yan, R. F. Prospects and research progress on Eucommia ulmoides gum. Prog. Chem. 7, 65-71 (1995).

[14] Nakazawa, Y. et al. Production of Eucommia-rubber from Eucommia ulmoides Oliv (Hardy Rubber Tree). Plant Biotechnol. 26, 71-79 (2009).

[15] Wuyun, T. N. et al. The hardy rubber tree genome provides insights into the evolution of polyisoprene biosynthesis. Mol. Plant 11, 429-442 (2018).

[16] Ong, K. W., Hsu, A. & Tan, B. K. Chlorogenic acid stimulates glucose transport in skeletal muscle via AMPK activation: a contributor to the beneficial effects of coffee on diabetes. PLoS ONE 7, e32718 (2012).

[17] Onakpoya, I. J., Spencer, E. A., Thompson, M. J. & Heneghan, C. J. The effect of chlorogenic acid on blood pressure: a systematic review and meta-analysis of randomized clinical trials. J. Hum. Hypertens. 29, 77-81 (2015).

[18] Li, F. D. & Du, H. Y. Eucommia ulmoides (China Press of Traditional Chinese Medicine, 2001).

[19] Kwan, C. Y., Chen, C. X., Deyama, T. & Nishibe, S. Endothelium-dependent vasorelaxant effects of the aqueous extracts of the Eucommia ulmoides Oliv. leaf and bark: implications on their antihypertensive action. Vasc. Pharmacol. 40, 229-235 (2003).

[20] Wuyun, T. N., Du, H. Y., Liu, H. M., Wang, L. & Bao, W. Q. Important Research Progress in Genome Sequencing of Eucommia ulmoides (2014-2015) (Social Sciences Academic Press, 2015).

[21] Chin, C. S. et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 10, 563-569 (2013).

[22] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[23] Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078-2079 (2009).

[24] Burton, J. N. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol. 31, 1119-1125 (2013).

[25] Parra, G., Bradnam, K. & Korf, I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23, 1061-1067 (2007).

[26] Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[27] Keller, O., Kollmar, M., Stanke, M. & Waack, S. A novel hybrid gene prediction method employing protein multiple sequence alignments. Bioinformatics 27, 757-763 (2011).

[28] Majoros, W. H., Pertea, M. & Salzberg, S. L. TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics 20, 2878-2879 (2004).

[29] Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403-410 (1990).

[30] Birney, E., Clamp, M. & Durbin, R. GeneWise and genomewise. Genome Res. 14, 988-995 (2004).

[31] Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644-652 (2011).

[32] Haas, B. J. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 31, 5654-5666 (2003).

[33] Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639-1645 (2009).

[34] Li, L., Stoeckert, C. J. Jr. & Roos, D. S. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178-2189 (2003).

[35] De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-1271 (2006).

[36] Liu, Z. H. Floral Organogenesis and Phylogenetic Analysis in Eucommia ulmoides Oilv (Peking University, 2010).

[37] Zhang, N., Zeng, L., Shan, H. & Ma, H. Highly conserved low-copy nuclear genes as effective markers for phylogenetic analyses in angiosperms. N. Phytol. 195, 923-937 (2012).

[38] Van de Peer, Y., Mizrachi, E. & Marchal, K. The evolutionary significance of polyploidy. Nat. Rev. Genet. 18, 411-424 (2017).

[39] Iorizzo, M. et al. A high-quality carrot genome assembly provides new insights into carotenoid accumulation and asterid genome evolution. Nat. Genet. 48, 657-666 (2016).

[40] Tang, H., Krishnakumar, V. & Li, J. jcvi: JCVI utility libraries (Version v0.5.7). Zenodo (2015).

[41] Chae, L., Kim, T., Nilo-Poyanco, R. & Rhee, S. Y. Genomic signatures of specialized metabolism in plants. Science 344, 510-513 (2014).

[42] Nutzmann, H. W., Huang, A. & Osbourn, A. Plant metabolic clusters-from genetics to genomics. N. Phytol. 211, 771-789 (2016).

[43] Kautsar, S. A., Suarez Duran, H. G., Blin, K., Osbourn, A. & Medema, M. H. plantiSMASH: automated identification, annotation and expression analysis of plant biosynthetic gene clusters. Nucleic Acids Res. 45, W55-W63 (2017).

[44] Nutzmann, H. W. & Osbourn, A. Gene clustering in plant specialized metabolism. Curr. Opin. Biotechnol. 26, 91-99 (2014).

[45] Bamba, T. et al. Contribution of mevalonate and methylerythritol phosphate pathways to polyisoprenoid biosynthesis in the rubber-producing plant Eucommia ulmoides Oliver. Z. Naturforsch. C. J. Biosci. 65, 363-372 (2010).

[46] Suzuki, N. et al. Construction and analysis of EST libraries of the trans-polyisoprene producing plant, Eucommia ulmoides Oliver. Planta 236, 1405-1417 (2012).

[47] Rohmer, M. The discovery of a mevalonate-independent pathway for isoprenoid biosynthesis in bacteria, algae and higher plants. Nat. Prod. Rep. 16, 565-574 (1999).

[48] Kuzuyama, T. Mevalonate and nonmevalonate pathways for the biosynthesis of isoprene units. Biosci. Biotechnol. Biochem. 66, 1619-1627 (2002).

[49] Cui, Y., Wang, M. & Sun, K. Morphological study of gutta containing cells in Eucommia ulmoides Oliv. Chin. Bull. Bot. 16, 439-443 (1999).

[50] Du, H. Y., Du, L. Y. & Li, F. D. Dynamic of gutta-percha formation and accumulation in samara of Eucommia ulmoides. For. Res. 17, 185-191 (2004).

[51] Oh, S. K. et al. Isolation, characterization, and functional analysis of a novel cDNA clone encoding a small rubber particle protein from Hevea brasiliensis. J. Biol. Chem. 274, 17132-17138 (1999).

[52] Lin, T. et al. Genome analysis of Taraxacum kok-saghyz Rodin provides new insights into rubber biosynthesis. Nat. Sci. Rev. 5, 78-87 (2017).

[53] Dennis, M. S. & Light, D. R. Rubber elongation factor from Hevea brasiliensis. Identification, characterization, and role in rubber biosynthesis. J. Biol. Chem. 264, 18608-18617 (1989).

[54] Priya, P., Venkatachalam, P. & Thulaseedharan, A. Differential expression pattern of rubber elongation factor (REF) mRNA transcripts from high and low yielding clones of rubber tree (Hevea brasiliensis Muell. Arg.). Plant Cell Rep. 26, 1833-1838 (2007).

[55] Unver, T. et al. Genome of wild Olive and the evolution of oil biosynthesis. Proc. Natl Acad. Sci. USA 114, E9413-E9422 (2017).

[56] Wu, L. Q., Zhu, W. X., Zhang, Y. X., Fu, D. D. & Zhang, Z. X. Analysis of chlorogenic acid content and extraction detection method in Eucommia ulmoides Oliv. Food Sci. 26, 187-192 (2005).

[57] Huang, M. T., Smart, R. C., Wong, C. Q. & Conney, A. H. Inhibitory effect of curcumin, chlorogenic acid, caffeic acid, and ferulic acid on tumor promotion in mouse skin by 12-O-tetradecanoylphorbol-13-acetate. Cancer Res. 48, 5941-5946 (1988).

[58] Tajik, N., Tajik, M., Mack, I. & Enck, P. The potential effects of chlorogenic acid, the main phenolic components in coffee, on health: a comprehensive review of the literature. Eur. J. Nutr. 56, 2215-2244 (2017).

[59] Franke, R. et al. The Arabidopsis REF8 gene encodes the 3-hydroxylase of phenylpropanoid metabolism. Plant J. 30, 33-45 (2002).

[60] Niggeweg, R., Michael, A. J. & Martin, C. Engineering plants with increased levels of the antioxidant chlorogenic acid. Nat. Biotechnol. 22, 746-754 (2004).

[61] Bartley, G. E. et al. Transcriptional regulation of chlorogenic acid biosynthesis in carrot root slices exposed to UV-B light. Plant Gene 7, 1-10 (2016).

[62] Wang, H. et al. Mutation of WRKY transcription factors initiates pith secondary wall formation and increases stem biomass in dicotyledonous plants. Proc. Natl Acad. Sci. USA 107, 22338-22343 (2010).

[63] Zha, L. et al. DNA methylation influences chlorogenic acid biosynthesis in Lonicera japonica by mediating LjbZIP8 to regulate phenylalanine ammonia-lyase 2 expression. Front. Plant. Sci. 8, 1178 (2017).

[64] Zhang, J. et al. Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus. N. Phytol. 220, 502-516 (2018).

[65] Hiruma, K. Roles of plant-derived secondary metabolites during interactions with pathogenic and beneficial microbes under conditions of environmental stress. Microorganisms 7, 362 (2019).

[66] Wink, M. Plant secondary metabolites modulate insect behavior-steps toward addiction? Front. Physiol. 9, 364 (2018).

[67] McLean, S. & Duncan, A. J. Pharmacological perspectives on the detoxification of plant secondary metabolites: implications for ingestive behavior of herbivores. J. Chem. Ecol. 32, 1213-1228 (2006).

[68] Germanà, M. A. Anther culture for haploid and doubled haploid production. Plant Cell Tissue Organ Cult. 104, 283-300 (2011).

[69] Kang, X. Y. Heating treatment device for non-separated branches and buds of trees. CN Patent No. 200610113448.X (2006).

[70] Feuillet, C., Leach, J. E., Rogers, J., Schnable, P. S. & Eversole, K. Crop genome sequencing: lessons and rationales. Trends Plant Sci. 16, 77-88 (2011).

[71] Yahav, T. & Privman, E. A comparative analysis of methods for de novo assembly of hymenopteran genomes using either haploid or diploid samples. Sci. Rep. 9, 6480 (2019).

[72] Simpson, J. T. & Pop, M. The theory and practice of genome sequence assembly. Annu. Rev. Genomics Hum. Genet. 16, 153-172 (2015).

[73] Bourque, G. et al. Ten things you should know about transposable elements. Genome Biol. 19, 199 (2018).

[74] Chuong, E. B., Elde, N. C. & Feschotte, C. Regulatory activities of transposable elements: from conflicts to benefits. Nat. Rev. Genet. 18, 71-86 (2017).

[75] Momohara, A. Floral and paleoenvironmental history from the late Pliocene to middle Pleistocene in and around central Japan. Palaeogeogr. Palaeoclimatol. Palaeoecol. 108, 281-293 (1994).

[76] Cerling, T. E. in C4 Plant Biology (eds Sage, R. F. & Monson, R. K.), 445-469 (Academic Press, 1999).

[77] Grandbastien, M. A. Activation of plant retrotransposons under stress conditions. Trends Plant Sci. 3, 181-187 (1998).

[78] Lisch, D. How important are transposons for plant evolution? Nat. Rev. Genet. 14, 49-61 (2013).

[79] Panchy, N., Lehti-Shiu, M. & Shiu, S. H. Evolution of gene duplication in plants. Plant Physiol. 171, 2294-2316 (2016).

[80] Vranova, E., Coman, D. & Gruissem, W. Structure and dynamics of the isoprenoid pathway network. Mol. Plant 5, 318-333 (2012).

[81] Putter, K. M., Van Deenen, N., Unland, K., Prufer, D. & Schulze Gronover, C. Isoprenoid biosynthesis in dandelion latex is enhanced by the overexpression of three key enzymes involved in the mevalonate pathway. BMC Plant Biol. 17, 88 (2017).

[82] Koyama, T. Molecular analysis of prenyl chain elongating enzymes. Biosci. Biotechnol. Biochem. 63, 1671-1676 (1999).

[83] Epping, J. et al. A rubber transferase activator is necessary for natural rubber biosynthesis in dandelion. Nat. Plants 1, 15048 (2015).

[84] Takahashi, S. & Koyama, T. Structure and function of cis-prenyl chain elongating enzymes. Chem. Rec. 6, 194-205 (2006).

[85] Schmidt, T. et al. Characterization of rubber particles and rubber chain elongation in Taraxacum kok-saghyz. BMC Biochem. 11, 11 (2010).

[86] Hillebrand, A. et al. Down-regulation of small rubber particle protein expression affects integrity of rubber particles and rubber content in Taraxacum brevicorniculatum. PLoS ONE 7, e41874 (2012).

[87] Yamashita, S. et al. Identification and reconstitution of the rubber biosynthetic machinery on rubber particles from Hevea brasiliensis. Elife 5, e19022 (2016).

[88] Tangpakdee, J. et al. Structure and biosynthesis of trans-polyisoprene from Eucommia ulmoides. Phytochemistry 45, 75-80 (1997).

[89] Liang, P. H., Ko, T. P. & Wang, A. H. J. Structure, mechanism and function of prenyltransferases. Eur. J. Biochem. 269, 3339-3354 (2002).

[90] Liu, J. et al. The Chromosome-based rubber tree genome provides new insights into spurge genome evolution and rubber biosynthesis. Mol. Plant 13, 336-350 (2020).

[91] Pendleton, M. et al. Assembly and diploid architecture of an individual human genome via single-molecule technologies. Nat. Methods 12, 780-786 (2015).

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

[93] Ghurye, J., Pop, M., Koren, S., Bickhart, D. & Chin, C. S. Scaffolding of long read assemblies using long range contact information. BMC Genom. 18, 527 (2017).