The genome of Dioscorea zingiberensis sheds light on the biosynthesis, origin and evolution of the medicinally important diosgenin saponins

Abstract

Diosgenin saponins isolated from Dioscorea species such as D. zingiberensis exhibit a broad spectrum of pharmacological activities. Diosgenin, the aglycone of diosgenin saponins, is an important starting material for the production of steroidal drugs. However, how plants produce diosgenin saponins and the origin and evolution of the diosgenin saponin biosynthetic pathway remain a mystery. Here we report a high-quality, 629-Mb genome of D. zingiberensis anchored on 10 chromosomes with 30 322 protein-coding genes. We reveal that diosgenin is synthesized in leaves (‘source’), then converted into diosgenin saponins, and finally transported to rhizomes (‘sink’) for storage in plants. By evaluating the distribution and evolutionary patterns of diosgenin saponins in Dioscorea species, we find that diosgenin saponin-containing may be an ancestral trait in Dioscorea and is selectively retained. The results of comparative genomic analysis indicate that tandem duplication coupled with a whole-genome duplication event provided key evolutionary resources for the diosgenin saponin biosynthetic pathway in the D. zingiberensis genome. Furthermore, comparative transcriptome and metabolite analysis among 13 Dioscorea species suggests that specific gene expression patterns of pathway genes promote the differential evolution of the diosgenin saponin biosynthetic pathway in Dioscorea species. Our study provides important insights and valuable resources for further understanding the biosynthesis, evolution, and utilization of plant specialized metabolites such as diosgenin saponins.

References

[1] Marker RE, Tsukamoto T, Turner DL . Sterols. C. Diosgenin1. J Am Chem Soc. 1940; 62: 2525-32.

[2] Li H, Huang W, Wen Y et al. Anti-thrombotic activity and chemical characterization of steroidal saponins from Dioscorea zingiberensis C.H. Wright. Fitoterapia. 2010; 81: 1147-56.

[3] Tan Y, Xiao X, Yao J et al. Syntheses and anti-cancer activities of glycosylated derivatives of diosgenin. Chem Res Chin Univ. 2017; 33: 80-6.

[4] Sethi G, Shanmugam M, Warrier S et al. Pro-apoptotic and anti-cancer properties of diosgenin: a comprehensive and critical review. Nutrients 2018; 10: 645.

[5] Yi T, Fan LL, Chen HL et al. Comparative analysis of diosgenin in Dioscorea species and related medicinal plants by UPLC-DAD-MS. BMC Biochem. 2014; 15: 19.

[6] Sautour M, Mitaine-Offer AC, Lacaille-Dubois MA . The Dioscorea genus: a review of bioactive steroid saponins. J Nat Med. 2007; 61: 91-101.

[7] Hsu KM, Tsai JL, Chen MY et al. Molecular phylogeny of Dioscorea (Dioscoreaceae) in east and Southeast Asia. Blumea. 2013; 58: 21-7.

[8] Viruel J, Segarra-Moragues JG, Raz L et al. Late cretaceous-early Eocene origin of yams (Dioscorea, Dioscoreaceae) in the Laurasian Palaearctic and their subsequent Oligocene-Miocene diversification . J Biogeogr. 2016; 43: 750-62.

[9] Shen L, Xu J, Luo L et al. Predicting the potential global distribution of diosgenin-contained Dioscorea species. Chin Med. 2018; 13: 58.

[10] Zhang X, Jin M, Tadesse N et al. Dioscorea zingiberensis C. H. Wright: an overview on its traditional use, phytochemistry, pharmacology, clinical applications, quality control, and toxicity . J Ethnopharmacol. 2018; 220: 283-93.

[11] Wang W, Zhao Y, Jing W et al. Ultrahigh-performance liquid chromatography-ion trap mass spectrometry characterization of the steroidal saponins of Dioscorea panthaica Prain et Burkill and its application for accelerating the isolation and structural elucidation of steroidal saponins . Steroids. 2015; 95: 51-65.

[12] Lin S, Wang D, Yang D et al. Characterization of steroidal saponins in crude extract from Dioscorea nipponica Makino by liquid chromatography tandem multi-stage mass spectrometry. Anal Chim Acta. 2007; 599: 98-106.

[13] Lee HJ, Watanabe B, Nakayasu M et al. Novel steroidal saponins from Dioscorea esculenta (Togedokoro) . Biosci Biotechnol Biochem. 2017; 81: 2253-60.

[14] Tapondjou LA, Jenett-Siems K, Böttger S et al. Steroidal saponins from the flowers of Dioscorea bulbifera var. sativa . Phytochemistry. 2013; 95: 341-50.

[15] Ye T, Song W, Zhang JJ et al. Identification and functional characterization of DzS3GT, a cytoplasmic glycosyltransferase catalyzing biosynthesis of diosgenin 3-O-glucoside in Dioscorea zingiberensis . Plant Cell Tissue Organ Cult. 2017; 129: 399-410.

[16] Ding ZZ, Tang SR, Qin HZ et al. Resource Plant of Steroid Hormone. Chapter 1 Dioscorea species. In: Zhai RK. Beijing: Publishing House of Science; 1983: 14-113.

[17] Song W, Yan S, Li Y et al. Functional characterization of squalene epoxidase and NADPH-cytochrome P450 reductase in Dioscorea zingiberensis . Biochem Biophys Res Commun. 2019; 509: 822-7.

[18] Li X, Wang Y, Sun J et al. Chemotaxonomic studies of 12 Dioscorea species from China by UHPLC-QTOF-MS/MS analysis. Phytochem Anal. 2020; 31: 164-82.

[19] Moses T, Papadopoulou KK, Osbourn A . Metabolic and functional diversity of saponins, biosynthetic intermediates and semi-synthetic derivatives. Crit Rev Biochem Mol Biol. 2014; 49: 439-62.

[20] Nakayasu M, Kawasaki T, Lee HJ et al. Identification of furostanol glycoside 26-O- β-glucosidase involved in steroidal saponin biosynthesis from Dioscorea esculenta . Plant Biotechnol (Tokyo). 2015; 32: 299-308.

[21] Sonawane PD, Pollier J, Panda S et al. Plant cholesterol biosynthetic pathway overlaps with phytosterol metabolism. Nat Plants. 2017; 3: 16205.

[22] Upadhyay S, Jeena GS, Shukla RK . Recent advances in steroidal saponins biosynthesis and in vitro production. Planta. 2018; 248: 519-44.

[23] Mohammadi M, Mashayekh T, Rashidi-Monfared S et al. New insights into diosgenin biosynthesis pathway and its regulation in Trigonella foenum-graecum L. Phytochem Anal. 2020; 31: 229-41.

[24] Christ B, Xu C, Xu M et al. Repeated evolution of cytochrome P450-mediated spiroketal steroid biosynthesis in plants. Nat Commun. 2019; 10: 1-11.

[25] Cheng J, Chen J, Liu X et al. The origin and evolution of the diosgenin biosynthetic pathway in yam. Plant Commun. 2021; 2: 100079.

[26] Simão FA, Waterhouse RM, Ioannidis P et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015; 31: 3210-2.

[27] Parra G, Bradnam K, Korf I . CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007; 23: 1061-7.

[28] Gottlieb A, Müller HG, Massa AN et al. Insular organization of gene space in grass genomes. PLoS One. 2013; 8: e54101.

[29] Luo MC, Gu YQ, Puiu D et al. Genome sequence of the progenitor of the wheat D genome Aegilops tauschii . Nature. 2017; 551: 498-502.

[30] Li L, Stoeckert CJJ, Roos DS . OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003; 13: 2178-89.

[31] Abbas FA, Ateya AM . Estradiol, esteriol, estrone and novel flavonoids from date palm pollen. Aust J Basic Appl Sci. 2011; 5: 606-14.

[32] Wang L, Wang X, Yuan X et al. Simultaneous analysis of diosgenin and sarsasapogenin in Asparagus officinalis byproduct by thin-layer chromatography. Phytochem Anal. 2011; 22: 14-7.

[33] Jiao Y, Wickett NJ, Ayyampalayam S et al. Ancestral polyploidy in seed plants and angiosperms. Nature. 2011; 473: 97-100.

[34] Clark JW, Donoghue PCJ . Whole-genome duplication and plant macroevolution. Trends Plant Sci. 2018; 23: 933-45.

[35] Ren R, Wang H, Guo C et al. Widespread whole genome duplications contribute to genome complexity and species diversity in angiosperms. Mol Plant. 2018; 11: 414-28.

[36] Li J, Yang D, Yu K et al. Determination of diosgenin content in medicinal plants with enzyme-linked immunosorbent assay. Planta Med. 2010; 76: 1915-20.

[37] Wang X, Chen D, Wang Y et al. De novo transcriptome assembly and the putative biosynthetic pathway of steroidal sapogenins of Dioscorea composita . PLoS One. 2015; 10: e0124560.

[38] Gas-Pascual E, Berna A, Bach TJ et al. Plant oxidosqualene metabolism: cycloartenol synthase-dependent sterol biosynthesis in Nicotiana benthamiana . PLoS One. 2014; 9: e109156.

[39] Gang DR . Evolution of flavors and scents. Annu Rev Plant Biol. 2005; 56: 301-25.

[40] Chae L, Kim T, Nilo-Poyanco R et al. Genomic signatures of specialized metabolism in plants. Science. 2014; 344: 510-3.

[41] Singh P, Singh G, Bhandawat A et al. Spatial transcriptome analysis provides insights of key gene(s) involved in steroidal saponin biosynthesis in medicinally important herb Trillium govanianum . Sci Rep. 2017; 7: 45295.

[42] Ciura J, Szeliga M, Grzesik M et al. Changes in fenugreek transcriptome induced by methyl jasmonate and steroid precursors revealed by RNA-Seq. Genomics. 2018; 110: 267-76.

[43] Li J, Liang Q, Li C et al. Comparative transcriptome analysis identifies putative genes involved in dioscin biosynthesis in Dioscorea zingiberensis . Molecules. 2018; 23: 1-13.

[44] Chen S, Xu J, Liu C et al. Genome sequence of the model medicinal mushroom Ganoderma lucidum . Nat Commun. 2012; 3: 913.

[45] Guo YL . Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes. Plant J. 2013; 73: 941-51.

[46] Zhou C, Yang Y, Tian J et al. 22R-but not 22S-hydroxycholesterol is recruited for diosgenin biosynthesis. Plant J. 2022; 109: 940-51.

[47] Hou L, Yuan X, Li S et al. Genome-wide identification of CYP72A gene family and expression patterns related to jasmonic acid treatment and steroidal saponin accumulation in Dioscorea zingiberensis . Int J Mol Sci 2021; 22: 10953.

[48] Prall W, Hendy O, Thornton LE . Utility of a phylogenetic perspective in structural analysis of CYP72A enzymes from flowering plants. PLoS One. 2016; 11: e0163024.

[49] He J, Chen Q, Xin P et al. CYP72A enzymes catalyse 13-hydrolyzation of gibberellins. Nat Plants. 2019; 5: 1057-65.

[50] Zhao Q, Yang J, Cui MY et al. The reference genome sequence of Scutellaria baicalensis provides insights into the evolution of wogonin biosynthesis. Mol Plant. 2019; 12: 935-50.

[51] Cheng J, Wang X, Liu X et al. Chromosome-level genome of Himalayan yew provides insights into the origin and evolution of the paclitaxel biosynthetic pathway. Mol Plant. 2021; 14: 1199-209.

[52] Xia EH, Zhang HB, Sheng J et al. The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis. Mol Plant. 2017; 10: 866-77.

[53] Belton JM, McCord RP, Gibcus JH et al. Hi-C: a comprehensive technique to capture the conformation of genomes. Methods. 2012; 58: 268-76.

[54] Kokot M, Dlugosz M, Deorowicz S . KMC 3: counting and manipulating k-mer statistics. Bioinformatics. 2017; 33: 2759-61.

[55] Vurture GW, Sedlazeck FJ, Nattestad M et al. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics. 2017; 33: 2202-4.

[56] Li H . Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences. Bioinformatics. 2016; 32: 2103-10.

[57] Hu J, Fan J, Sun Z et al. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics. 2020; 36: 2253-5.

[58] Coombe L, Zhang J, Vandervalk BP et al. ARKS: chromosome-scale scaffolding of human genome drafts with linked read kmers. BMC Bioinformatics. 2018; 19: 234.

[59] Warren RL, Yang C, Vandervalk BP et al. LINKS: scalable, alignment-free scaffolding of draft genomes with long reads. GigaScience. 2015; 4: 35.

[60] Burton JN, Adey A, Patwardhan RP et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31: 1119-25.

[61] Xu Z, Wang H . LTR-FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Symp Ser. 2007; 35: W265-8.

[62] Han Y, Wessler SR . MITE-hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences. Nucleic Acids Symp Ser. 2010; 38: e199-8.

[63] Stanke M, Schöffmann O, Morgenstern B et al. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinformatics. 2006; 7: 1-11.

[64] Keilwagen J, Wenk M, Erickson JL et al. Using intron position conservation for homology-based gene prediction. Nucleic Acids Res. 2016; 44: e89-11.

[65] Haas BJ, Salzberg SL, Zhu W et al. Automated eukaryotic gene structure annotation using evidence modeler and the program to assemble spliced alignments. Genome Biol. 2008; 9: R7-22.

[66] Stamatakis A . RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014; 30: 1312-3.

[67] Yang Z . PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24: 1586-91.

[68] De Bie T, Cristianini N, Demuth JP et al. CAFE: a computational tool for the study of gene family evolution. Bioinformatics. 2006; 22: 1269-71.

[69] Edgar RC . MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004; 5: 113-9.

[70] Suyama M, Torrents D, Bork P . PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006; 34: W609-12.

[71] Scrucca L, Fop M, Murphy TB et al. Mclust 5: clustering, classification and density estimation using Gaussian finite mixture models. R J. 2016; 8: 289-317.

[72] Dobin A, Davis CA, Schlesinger F et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013; 29: 15-21.

[73] Anders S, Pyl PT, Huber W . HTSeq - a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015; 31: 166-9.

[74] Shannon P, Markiel A, Ozier O et al. A software environment for integrated models. Genome Res. 1971; 13: 2498-504.

[75] Johnson LS, Eddy SR, Portugaly E . Hidden Markov model speed heuristic and iterative HMM search procedure. BMC Bioinformatics. 2010; 11: 431.

[76] Nelson DR . The cytochrome P450 homepage. Hum Genomics. 2009; 4: 59-65.

[77] Grabherr MG, Haas BJ, Yassour M et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29: 644-52.

[78] Li B, Dewey CN . RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011; 12: 323.

[79] Zanne AE, Tank DC, Cornwell WK et al. Three keys to the radiation of angiosperms into freezing environments. Nature. 2014; 506: 89-92.

[80] Paradis E, Claude J, Strimmer K . APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004; 20: 289-90.

[81] Revell LJ . Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012; 3: 217-23.

[82] Pagel M . Inferring evolutionary processes from phylogenies. Zool Scr. 1997; 26: 331-48.

[83] Blomberg SP, Garland T, Ives AR . Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003; 57: 717-45.