Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant

Abstract

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of β-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

References

[1] Moses, T., Papadopoulou, K. K. & Osbourn, A. Metabolic and functional diversity of saponins, biosynthetic intermediates and semi-synthetic derivatives. Crit. Rev. Biochem. Mol. Biol. 49, 439-462 (2014).

[2] Thimmappa, R., Geisler, K., Louveau, T., O’Maille, P. & Osbourn, A. Triterpene biosynthesis in plants. Annu. Rev. Plant Biol. 65, 225-257 (2014).

[3] Glauert, A. M., Dingle, J. & Lucy, J. Action of saponin on biological cell membranes. Nature 196, 953-955 (1962).

[4] Denoeud, F. et al. The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. Science 345, 1181-1184 (2014).

[5] Zhang, D. et al. The medicinal herb Panax notoginseng genome provides insights into ginsenoside biosynthesis and genome evolution. Mol. Plant 10, 903-907 (2017).

[6] Seki, H., Tamura, K. & Muranaka, T. P450s and UGTs: Key players in the structural diversity of triterpenoid saponins. Plant Cell Physiol. 56, 1463-1471 (2015).

[7] Miettinen, K. et al. The ancient CYP716 family is a major contributor to the diversification of eudicot triterpenoid biosynthesis. Nat. Commun. 8, 14153 (2017).

[8] Tamura, K. et al. Cytochrome P450 monooxygenase CYP716A141 is a unique beta-amyrin C-16beta oxidase Involved in triterpenoid saponin biosynthesis in Platycodon grandiflorus. Plant Cell Physiol. 58, 874-884 (2017).

[9] Nyakudya, E., Jeong, J. H., Lee, N. K. & Jeong, Y. S. Platycosides from the roots of Platycodon grandiflorum and their health benefits. Prev. Nutr. Food Sci. 19, 59-68 (2014).

[10] Zhang, L. et al. Platycodon grandiflorus - an ethnopharmacological, phytochemical and pharmacological review. J. Ethnopharmacol. 164, 147-161 (2015).

[11] Ma, C. H. et al. Candidate genes involved in the biosynthesis of triterpenoid saponins in Platycodon grandiflorum identified by transcriptome analysis. Front. Plant Sci. 7, 673 (2016).

[12] Xu, J. et al. Panax ginseng genome examination for ginsenoside biosynthesis. Gigascience 6, 1-15 (2017).

[13] Kim, N. H. et al. Genome and evolution of the shade-requiring medicinal herb Panax ginseng. Plant Biotechnol. J. 16, 1904-1917 (2018).

[14] Mapleson, D., Garcia Accinelli, G., Kettleborough, G., Wright, J. & Clavijo, B. KAT: a K-mer analysis toolkit to quality control NGS datasets and genome assemblies. Bioinformatics 33, 574-576 (2017).

[15] Mertens, J. et al. The bHLH transcription factors TSAR1 and TSAR2 regulate triterpene saponin biosynthesis in Medicago truncatula. Plant Physiol. 170, 194-210 (2016).

[16] Jarvis, D. E. et al. The genome of Chenopodium quinoa. Nature 542, 307-312 (2017).

[17] 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).

[18] Badouin, H. et al. The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution. Nature 546, 148-152 (2017).

[19] Hong, C. P. et al. accD nuclear transfer of Platycodon grandiflorum and the plastid of early Campanulaceae. BMC Genomics 18, 607 (2017).

[20] Ghosh, S. Triterpene structural diversification by plant cytochrome P450 enzymes. Front. Plant Sci. 8, 1886 (2017).

[21] Zhang, L., Liu, Z.-H. & Tian, J.-K. Cytotoxic triterpenoid saponins from the roots of Platycodon grandiflorum. Molecules 12, 832-841 (2007).

[22] Kwon, J. et al. Effect of processing method on platycodin D content in Platycodon grandiflorum roots. Arch. Pharm. Res. 40, 1087-1093 (2017).

[23] Seki, H. et al. Triterpene functional genomics in licorice for identification of CYP72A154 involved in the biosynthesis of glycyrrhizin. Plant Cell 23, 4112-4123 (2011).

[24] Wang, J., Gao, W., Zuo, B., Zhang, L. & Huang, L. Effect of methyl jasmonate on the ginsenoside content of Panax ginseng adventitious root cultures and on the genes involved in triterpene biosynthesis. Res. Chem. Intermed. 39, 1973-1980 (2013).

[25] Wang, W., Zhao, Z. J., Xu, Y., Qian, X. & Zhong, J. J. Efficient induction of ginsenoside biosynthesis and alteration of ginsenoside heterogeneity in cell cultures of Panax notoginseng by using chemically synthesized 2-hydroxyethyl jasmonate. Appl. Microbiol. Biotechnol. 70, 298-307 (2006).

[26] Chen, W. et al. Whole-genome sequencing and analysis of the Chinese herbal plant Panax notoginseng. Mol. Plant 10, 899-902 (2017).

[27] Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562-578 (2012).

[28] Luo, R. et al. Erratum: SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience 4, 30 (2015).

[29] 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).

[30] Nadalin, F., Vezzi, F. & Policriti, A. GapFiller: a de novo assembly approach to fill the gap within paired reads. BMC Bioinformatics 13, S8 (2012).

[31] Denisov, G. et al. Consensus generation and variant detection by Celera Assembler. Bioinformatics 24, 1035-1040 (2008).

[32] Boetzer, M., Henkel, C. V., Jansen, H. J., Butler, D. & Pirovano, W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics 27, 578-579 (2011).

[33] 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).

[34] Lomsadze, A., Burns, P. D. & Borodovsky, M. Integration of mapped RNA-Seq reads into automatic training of eukaryotic gene finding algorithm. Nucleic Acids Res. 42, e119 (2014).

[35] Stanke, M., Diekhans, M., Baertsch, R. & Haussler, D. Using native and syntentically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24, 637-644 (2008).

[36] Trapnell, C., Pachter, L. & Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105-1111 (2009).

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

[38] Wu, T. D. & Watanabe, C. K. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21, 1859-1875 (2005).

[39] Jaillon, O. et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 449, 463-467 (2007).

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

[41] Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312-1313 (2014).

[42] Kumar, S., Stecher, G., Suleski, M. & Hedges, S. B. TimeTree: a resource for timelines, timetrees, and divergence times. Mol. Biol. Evol. 34, 1812-1819 (2017).

[43] Han, M. V., Thomas, G. W., Lugo-Martinez, J. & Hahn, M. W. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol. Biol. Evol. 30, 1987-1997 (2013).

[44] Urich, M. A., Nery, J. R., Lister, R., Schmitz, R. J. & Ecker, J. R. MethylC-seq library preparation for base-resolution whole-genome bisulfite sequencing. Nat. Protoc. 10, 475-483 (2015).

[45] Krueger, F. & Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27, 1571-1572 (2011).

[46] Hofmeister, B. T., Lee, K., Rohr, N. A., Hall, D. W. & Schmitz, R. J. Stable inheritance of DNA methylation allows creation of epigenotype maps and the study of epiallele inheritance patterns in the absence of genetic variation. Genome Biol. 18, 155 (2017).