Whole-genome sequencing provides insights into the genetic diversity and domestication of bitter gourd (Momordica spp.)

Abstract

Bitter gourd (Momordica charantia) is a popular cultivated vegetable in Asian and African countries. To reveal the characteristics of the genomic structure, evolutionary trajectory, and genetic basis underlying the domestication of bitter gourd, we performed whole-genome sequencing of the cultivar Dali-11 and the wild small-fruited line TR and resequencing of 187 bitter gourd germplasms from 16 countries. The major gene clusters (Bi clusters) for the biosynthesis of cucurbitane triterpenoids, which confer a bitter taste, are highly conserved in cucumber, melon, and watermelon. Comparative analysis among cucurbit genomes revealed that the Bi cluster involved in cucurbitane triterpenoid biosynthesis is absent in bitter gourd. Phylogenetic analysis revealed that the TR group, including 21 bitter gourd germplasms, may belong to a new species or subspecies independent from M. charantia. Furthermore, we found that the remaining 166 M. charantia germplasms are geographically differentiated, and we identified 710, 412, and 290 candidate domestication genes in the South Asia, Southeast Asia, and China populations, respectively. This study provides new insights into bitter gourd genetic diversity and domestication and will facilitate the future genomics-enabled improvement of bitter gourd.

References

[1] Schaefer, H. & Renner, S. S. A three-genome phylogeny of Momordica (Cucurbitaceae) suggests seven returns from dioecy to monoecy and recent long-distance dispersal to Asia. Mol. Phylogenet. Evol. 54, 553-560 (2010).

[2] Decker-Walters, D. S. Cucurbits, Sanskrit, and the Indo-Aryas. Econ. Bot. 53, 98-112 (1999).

[3] Willdenow, C. L. Linnaei Species Plantarum http://www.botanicus.org/title/b1206998x (1805).

[4] Rheede, H. v. Hortus Malabaricus http://hortusmalabaricus.net/ (1688).

[5] John, K. J. & Antony, V. T. Collection and preliminary evaluation of small bitter gourds (Momordica charantia L.) a relict vegetable of Southern Peninsular India. Genet. Resour. Crop Evol. 56, 99-104 (2009).

[6] John, K. J. & Antony, V. T. A taxonomic revision of the genus Momordica L. (Cucurbitaceae) in India. Indian. J. Plant Genet. Resour. 23, 172-184 (2010).

[7] Chakravarty, H. Fascicles of Flora of India, Fascicle 11:Cucurbitaceae. 87-95 (Botanical Survey of India, Calcutta, 1982).

[8] Marr, K. L., Mei, X. Y. & Bhattarai, N. K. Allozyme, morphological and nutritional analysis bearing on the domestication of Momordica charantia L. (Cucurbitaceae). Econ. Bot. 58, 435-455 (2004).

[9] Behera, T. K. et al. Bitter gourd: botany, horticulture, breeding. Hortic. Rev. 37, 101-141 (2010).

[10] Chen, J. C., Chiu, M. H., Nie, R. L., Cordell, G. A. & Qiu, S. X. Cucurbitacins and cucurbitane glycosides: structures and biological activities. Nat. Prod. Rep. 22, 386-399 (2005).

[11] Tan, M. J. et al. Antidiabetic activities of triterpenoids isolated from bitter melon associated with activation of the AMPK pathway. Chem. Biol. 15, 263-273 (2008).

[12] Amirthaveni, M., Premakumari, S., Gomathi, K. & Yang, R. Hypoglycemic effect of bitter gourd (Momordica charantia L.) among pre diabetics in India: a randomized placebo controlled cross over study. Indian. J. Nutr. Diet. 55, 44-63 (2018).

[13] Krawinkel, M. B. et al. Bitter gourd reduces elevated fasting plasma glucose levels in an intervention study among prediabteics. J. Ethnopharmacol. 216, 1-7 (2018).

[14] Dhillon, N. P. S., Sanguansil, S., Schafleitner, R., Wang, Y. W. & Mccreight, J. D. Diversity among a wide Asian collection of bitter gourd landraces and their genetic relationships with commercial hybrid cultivars. J. Am. Soc. Hort. Sci. 141, 475-484 (2016).

[15] Urasaki, N. et al. Draft genome sequence of bitter gourd (Momordica charantia), a vegetable and medicinal plant in tropical and subtropical regions. DNA Res. 24, 51-58 (2017).

[16] Huang, X. et al. A map of rice genome variation reveals the origin of cultivated rice. Nature 490, 497-501 (2012).

[17] Hufford, M. B. et al. Comparative population genomics of maize domestication and improvement. Nat. Genet. 44, 808-811 (2012).

[18] Zhou, Z. et al. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat. Biotechnol. 33, 408-414 (2015).

[19] Qi, J. et al. A genomic variation map provides insights into the genetic basis of cucumber domestication and diversity. Nat. Genet. 45, 1510-1515 (2013).

[20] Lin, T. et al. Genomic analyses provide insights into the history of tomato breeding. Nat. Genet. 46, 1220-1226 (2014).

[21] Li, R. et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 20, 265-272 (2010).

[22] Cui, J. et al. A RAD-based genetic map for anchoring scaffold sequences and identifying QTLs in bitter gourd (Momordica charantia). Front. Plant Sci. 9, 477 (2018).

[23] Li, R. et al. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966-1967 (2009).

[24] Simão, 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).

[25] Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).

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

[27] Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 12, 323 (2011).

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

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

[30] Smit, A., Hubley, R. & Green, P. RepeatMasker Open-3.0 http://www.repeatmasker.org (1996).

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

[32] Edgar, R. C. & Myers, E. W. PILER: identification and classification of genomic repeats. Bioinformatics 21(Suppl 1), i152-i158 (2005).

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

[34] Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402 (1997).

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

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

[37] Elsik, C. G. et al. Creating a honey bee consensus gene set. Genome Biol. 8, R13 (2007).

[38] Bairoch, A. & Apweiler, R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 28, 45-48 (2000).

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

[40] Kanehisa, M. & Goto, S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 27, 29-34 (2000).

[41] Zdobnov, E. M. & Apweiler, R. InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17, 847-848 (2001).

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

[43] Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591 (2007).

[44] Wikstrom, N., Savolainen, V. & Chase, M. W. Evolution of the angiosperms: calibrating the family tree. Proc. Biol. Sci. 268, 2211-2220 (2001).

[45] Sun, H. et al. Karyotype stability and unbiased fractionation in the paleo-allotetraploid cucurbita genomes. Mol. Plant 10, 1293-1306 (2017).

[46] Sebastian, P., Schaefer, H., Telford, I. R. & Renner, S. S. Cucumber (Cucumis sativus) and melon (C. melo) have numerous wild relatives in Asia and Australia, and the sister species of melon is from Australia. Proc. Natl Acad. Sci. USA 107, 14269-14273 (2010).

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

[48] Hasegawa, M., Kishino, H. & Yano, T. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22, 160-174 (1985).

[49] Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078-2079 (2009).

[50] Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589-595 (2010).

[51] Mckenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297-1303 (2010).

[52] Price, A. L. et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904-909 (2006).

[53] Tang, H., Peng, J., Wang, P. & Risch, N. J. Estimation of individual admixture: analytical and study design considerations. Genet. Epidemiol. 28, 289-301 (2005).

[54] Barrett, J. C., Fry, B., Maller, J. & Daly, M. J. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263-265 (2005).

[55] Schiffels, S. & Durbin, R. Inferring human population size and separation history from multiple genome sequences. Nat. Genet. 46, 919-925 (2014).

[56] Ossowski, S. et al. The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana. Science 327, 92-94 (2010).

[57] Watterson, G. A. On the number of segregating sites in genetical models without recombination. Theor. Popul. Biol. 7, 256-276 (1975).

[58] Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123, 585-595 (1989).

[59] Akey, J. M., Zhang, G., Zhang, K., Jin, L. & Shriver, M. D. Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 12, 1805-1814 (2002).

[60] Chen, H., Patterson, N. & Reich, D. Population differentiation as a test for selective sweeps. Genome Res. 20, 393-402 (2010).

[61] Huang, S. et al. The genome of the cucumber, Cucumis sativus L. Nat. Genet. 41, 1275-1281 (2009).

[62] Garcia-Mas, J. et al. The genome of melon (Cucumis melo L.). Proc. Natl Acad. Sci. USA 109, 11872-11877 (2012).

[63] Guo, S. et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions. Nat. Genet. 45, 51-58 (2013).

[64] Monteropau, J. et al. De novo assembly of the zucchini genome reveals a whole‐genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol. J. 16, 1161-1171 (2018).

[65] Wu, S. et al. The bottle gourd genome provides insights into Cucurbitaceae evolution and facilitates mapping of a Papaya ring-spot virus resistance locus. Plant J. 92, 963-975 (2017).

[66] Martínez-García, P. J. et al. The walnut (Juglans regia) genome sequence reveals diversity in genes coding for the biosynthesis of nonstructural polyphenols. Plant J. 87, 507-532 (2016).

[67] Schaefer, H., Heibl, C. & Renner, S. S. Gourds afloat: a dated phylogeny reveals an Asian origin of the gourd family (Cucurbitaceae) and numerous oversea dispersal events. Proc. R. Soc. B 276, 843-851 (2009).

[68] Cui, J. et al. Genome-wide analysis of simple sequence repeats in bitter gourd (Momordica charantia). Front. Plant Sci. 8, 1103 (2017).

[69] Doebley, J. F., Gaut, B. S. & Smith, B. D. The molecular genetics of crop domestication. Cell 127, 1309-1321 (2006).

[70] Przeworski, M., Hudson, R. R. & Di, R. A. Adjusting the focus on human variation. Trends Genet. 16, 296-302 (2000).

[71] Meyer, R. S., Duval, A. E. & Jensen, H. R. Patterns and processes in crop domestication: an historical review and quantitative analysis of 203 global food crops. N. Phytol. 196, 29-48 (2012).

[72] Phillips, D. R., Rasbery, J. M., Bartel, B. & Matsuda, S. P. Biosynthetic diversity in plant triterpene cyclization. Curr. Opin. Plant Bio. 9, 305-314 (2006).

[73] Shibuya, M., Adachi, S. & Ebizuka, Y. Cucurbitadienol synthase, the first committed enzyme for cucurbitacin biosynthesis, is a distinct enzyme from cycloartenol synthase for phytosterol biosynthesis. Tetrahedron 60, 6995-7003 (2004).

[74] Shang, Y. et al. Biosynthesis, regulation, and domestication of bitterness in cucumber. Science 346, 1084-1088 (2014).

[75] Davidovich-Rikanati, R. et al. Recombinant yeast as a functional tool for understanding bitterness and cucurbitacin biosynthesis in watermelon (Citrullus spp.). Yeast 32, 103-114 (2015).

[76] Takase, S. et al. Identification of triterpene biosynthetic genes from Momordica charantia using RNA-seq analysis. Biosci. Biotech. Bioch. 83, 251-261 (2019).

[77] Zhou, Y. et al. Convergence and divergence of bitterness biosynthesis and regulation in Cucurbitaceae. Nat. Plants 2, 16183 (2016).

[78] Itkin, M. et al. The biosynthetic pathway of the nonsugar, high-intensity sweetener mogroside V from Siraitiagrosvenorii. Proc. Natl Acad. Sci. USA 113, 7619-7628 (2016).

[79] Xie, D. et al. The wax gourd genomes offer insights into the genetic diversity and ancestral cucurbit karyotype. Nat. Commun. 10, 5158 (2019).

[80] Enoch, A. D., Sognigb, N., Adam, A., Jean, G. & Frank, B. Phenetic analysis of wild populations of Momordica charantia L. (Cucurbitaceae) in West Africa and inference of the definition of the new subspecies macroloba Achigan-Dako & Blattner. Candollea 63, 153-167 (2008).

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

[82] Osbourn, A. Secondary metabolic gene clusters: evolutionary toolkits for chemical innovation. Trends Genet. 26, 449-457 (2010).

[83] Kliebenstein, D. J. & Osbourn, A. Making new molecules-evolution of pathways for novel metabolites in plants. Curr. Opin. Plant Biol. 15, 415-423 (2012).