Haplotype-resolved T2T genome and population resequencing provide insights into the domestication and mogroside biosynthesis of Siraitia grosvenorii (Cucurbitaceae)

Abstract

Monk fruit (Siraitia grosvenorii, Cucurbitaceae) is globally renowned for its triterpenoid glycoside mogroside V, a high-intensity, non-caloric natural sweetener. However, its domestication and mogroside biosynthesis remain largely unknown. Here, we report a haplotype-resolved telomere-to-telomere (T2T) gapless genome for monk fruit, consisting of 14 chromosomes with genome sizes of 316.21 Mb (Hap1) and 316.07 Mb (Hap2). Comparative genomic analyses of the haplotypes revealed that structural variations and transposable elements have significantly contributed to genomic variation and architecture in monk fruit. Population genomic analyses based on 173 re-sequenced genomes indicated that cultivated monk fruit was mainly domesticated in situ from local wild populations in northern Guangxi of China, and that it likely experienced a mild domestication bottleneck, while exhibiting low genetic diversity. Demographic inference further revealed that the low genetic diversity is largely attributed to demographic changes driven by historical climate shifts. Selective sweeps were identified across all chromosomes of cultivated monk fruit, among which are genes exhibiting diverse putative functions and involved in various biosynthetic processes and secondary metabolism. This pattern of selective sweeps demonstrates the joint role of artificial selection and demographic changes in shaping the genomic landscape of cultivated monk fruit. Furthermore, comparative transcriptome analyses showed a pronounced temporally specific expression pattern among mogroside biosynthesis genes during fruit development and delineated additional candidate genes potentially involved in mogroside biosynthesis. This study not only provides insights into the domestication and mogroside biosynthesis of monk fruit but also lays a valuable genomic foundation for its molecular breeding and mogroside-targeted synthetic biology.

References

[1] Guo J, Xu W, Hu Y, et al. Phylotranscriptomics in Cucurbitaceae reveal multiple whole-genome duplications and key morphological and molecular innovations. Mol Plant. 2020; 13: 1117-33

[2] Lu AM, Zhang ZY. The genus Siraitia Merr. China Guihaia. 1984; 4: 27-33

[3] Itkin M, Davidovich-Rikanati R, Cohen S, et al. The biosynthetic pathway of the nonsugar, high-intensity sweetener mogroside V from Siraitia grosvenorii. Proc Natl Acad Sci USA. 2016; 113: E7619-28

[4] Cui S, Zhang S, Wang N, et al. Structural insights into the catalytic selectivity of glycosyltransferase SgUGT94-289-3 towards mogrosides. Nat Commun. 2024; 15: 6423

[5] Shivani TBK, Mallikarjun CP, et al. Introduction, adaptation and characterization of monk fruit (Siraitia grosvenorii): a non-caloric new natural sweetener. Sci Rep. 2021; 11: 6205

[6] Li H, Li R, Jiang W, et al. Research progress of pharmacological effects of Siraitia grosvenorii extract. J Pharm Pharmacol. 2022; 74: 953-60

[7] Guo Y, Chen X, Gong P, et al. Siraitia grosvenorii as a homologue of food and medicine: a review of biological activity, mechanisms of action, synthetic biology, and applications in future food. J Agric Food Chem. 2024; 72: 6850-70

[8] Huang H, Peng Z, Zhan S, et al. A comprehensive review of Siraitia grosvenorii (Swingle) C. Jeffrey: chemical composition, pharmacology, toxicology, status of resources development, and applications. Front Pharmacol. 2024;15:1388747

[9] Gong X, Chen N, Ren K, et al. The fruits of Siraitia grosvenorii: a review of a Chinese food-medicine. Front Pharmacol. 2019; 10: 1400

[10] Tang SQ, Bin XY, Peng YT, et al. Assessment of genetic diversity in cultivars and wild accessions of Luohanguo (Siraitia grosvenorii [Swingle] A. M. Lu et Z. Y. Zhang), a species with edible and medicinal sweet fruits endemic to southern China, using RAPD and AFLP markers. Genet Resour Crop Evol. 2007; 54: 1053-61

[11] Xie B, Lai B, Chen L, et al. Phylogeographic analysis of Siraitia grosvenorii in subtropical China provides insights into the origin of cultivated monk fruit and conservation of genetic resources. Ecol Evol. 2023; 13: e10181

[12] Tang S, Li Y, Geng Y, et al. Clonal and spatial genetic structure in natural populations of Luohanguo (Siraitia grosvenorii), an economic species endemic to South China, as revealed by RAPD markers. Biochem Syst Ecol. 2007; 35: 557-65

[13] Zhou Y, Massonnet M, Sanjak JS, et al. Evolutionary genomics of grape (Vitis vinifera ssp. vinifera) domestication. Proc Natl Acad Sci USA. 2017; 114: 11715-20

[14] Lin J, Zhang W, Zhang X, et al. Signatures of selection in recently domesticated macadamia. Nat Commun. 2022; 13: 242

[15] Tang Q, Ma X, Mo C, et al. An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011; 12: 343

[16] Xia M, Han X, He H, et al. Improved de novo genome assembly and analysis of the Chinese cucurbit Siraitia grosvenorii, also known as monk fruit or luo-han-guo. Gigascience. 2018; 7: giy067

[17] Van Dijk EL, Naquin D, Gorrichon K, et al. Genomics in the long-read sequencing era. Trends Genet. 2023; 39: 649-71

[18] Cheng H, Concepcion GT, Feng X, et al. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021; 18: 170-5

[19] Cheng H, Jarvis ED, Fedrigo O, et al. Haplotype-resolved assembly of diploid genomes without parental data. Nat Biotechnol. 2022; 40: 1332-5

[20] Li C, Yuan Y, Nie Z, et al. The haplotype-resolved telomere-to-telomere genome and OMICS analyses reveal genetic responses to tapping in rubber tree. Nat Commun. 2025; 16: 6255

[21] Liu Y, Dong Y, Fang X, et al. Pangenomics and single-cell transcriptomics uncover the genetic basis of continuous bearing trait in grapevine. Hortic Res. 2025; 12: uhaf228

[22] Zhong H, Shi X, Zhang F, et al. Haplotype-resolved assemblies provide insights into genomic makeup of the oldest grapevine cultivar (Munage) in China. Hortic Res. 2025; 13: uhaf274

[23] Zhang Z, Liu W, Zhang T, et al. Population genomics of Vitis pseudoreticulata reveals the genetic basis of fungal resistance in grapevine. Hortic Plant J. 2025

[24] Yisilam G, Zheng E, Li C, et al. The chromosome-scale genome of black wolfberry (Lycium ruthenicum) provides useful genomic resources for identifying genes related to anthocyanin biosynthesis and disease resistance. Plant Diversity. 2025; 47: 201-13

[25] Wang J, Sun P, Li Y, et al. An overlooked paleotetraploidization in Cucurbitaceae. Mol Biol Evol. 2018; 35: 16-26

[26] Xie D, Xu Y, Wang J, et al. The wax gourd genomes offer insights into the genetic diversity and ancestral cucurbit karyotype. Nat Commun. 2019; 10: 5158

[27] Montero-Pau J, Blanca J, Bombarely A, 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. 2018; 16: 1161-71

[28] Zhang X, Zhao Y, Kou Y, et al. Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae). Hortic Res. 2022; 10: uhac231

[29] Zheng B, Xu Q, Shen Y . The relationship between climate change and Quaternary glacial cycles on the Qinghai-Tibetan Plateau: review and speculation. Quatern Int. 2002; 97-98: 93-101

[30] Colleoni F, Wekerle C, Näslund J, et al. Constraint on the penultimate glacial maximum Northern Hemisphere ice topography (≈140 kyrs BP). Quaternary Sci Rev. 2016; 137: 97-112

[31] Peng Y, Wang Y, Liu Y, et al. The genomic and epigenomic landscapes of hemizygous genes across crops with contrasting reproductive systems. Proc Natl Acad Sci USA. 2025; 122: e2422487122

[32] Wang C, Tan L, Zhang Z, et al. Haplotype-resolved genome reveals haplotypic variation and the biosynthesis of medicinal ingredients in Areca catechu L. Mol Hortic. 2025; 5: 24

[33] Shi X, Cao S, Wang X, et al. The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding. Hortic Res. 2023; 10: uhad061

[34] Chen W, Wang X, Sun J, et al. Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis. Nat Commun. 2024; 15: 4295

[35] Falcone Ferreyra ML, Rius SP, Casati P . Flavonoids: biosynthesis, biological functions, and biotechnological applications. Front Plant Sci. 2012; 3: 222

[36] Dong NQ, Lin HX . Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions. J Integr Plant Biol. 2021; 63: 180-209

[37] Waadt R, Seller CA, Hsu PK, et al. Plant hormone regulation of abiotic stress responses. Nat Rev Mol Cell Biol. 2022; 23: 680-94

[38] Yang Z, Li X, Yang L, et al. Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax. Plant Commun. 2023; 4: 100591

[39] Miller AJ, Gross BL . From forest to field: perennial fruit crop domestication. Am J Bot. 2011; 98: 1389-414

[40] Hufford MB, Xu X, van Heerwaarden J, et al. Comparative population genomics of maize domestication and improvement. Nat Genet. 2012; 44: 808-11

[41] Jing CY, Zhang FM, Wang XH, et al. Multiple domestications of Asian rice. Nat Plants. 2023; 9: 1221-35

[42] Meyer RS, Purugganan MD . Evolution of crop species: genetics of domestication and diversification. Nat Rev Genet. 2013; 14: 840-52

[43] Gaut BS, Seymour DK, Liu Q, et al. Demography and its effects on genomic variation in crop domestication. Nat Plants. 2018; 4: 512-20

[44] An Z, Kutzbach JE, Prell WL, et al. Evolution of Asian monsoons and phased uplift of the Himalaya-Tibetan plateau since Late Miocene times. Nature. 2001; 411: 62-6

[45] An Z, Clemens SC, Shen J, et al. Glacial-interglacial Indian summer monsoon dynamics. Science. 2011; 333: 719-23

[46] Van der Biezen EA, Freddie CT, Kahn K, et al. Arabidopsis RPP4 is a member of the RPP5 multigene family of TIR-NB-LRR genes and confers downy mildew resistance through multiple signalling components. Plant J. 2002; 29: 439-51

[47] Jeong HJ, Kim YJ, Kim SH, et al. Nonsense-mediated mRNA decay factors, UPF1 and UPF3, contribute to plant defense. Plant Cell Physiol. 2011; 52: 2147-56

[48] Cambiagno DA, Nota F, Zavallo D, et al. Immune receptor genes and pericentromeric transposons as targets of common epigenetic regulatory elements. Plant J. 2018; 96: 1178-90

[49] Karasov TL, Shirsekar G, Schwab R, et al. What natural variation can teach us about resistance durability. Curr Opin Plant Biol. 2020; 56: 89-98

[50] Long Q, Cao S, Huang G, et al. Population comparative genomics discovers gene gain and loss during grapevine domestication. Plant Phsysiol. 2024; 195: 1401-13

[51] Zhou Y, Minio A, Massonnet M, et al. The population genetics of structural variants in grapevine domestication. Nat Plants. 2019; 5: 965-79

[52] Zhang T, Peng W, Xiao H, et al. Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape. J Integr Plant Biol. 2024; 66: 1408-26

[53] Kou Y, Liao Y, Toivainen T, et al. Evolutionary genomics of structural variation in Asian rice (Oryza sativa) domestication. Mol Biol Evol. 2020; 37: 3507-24

[54] Marçais G, Kingsford C . A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics. 2011; 27: 764-70

[55] Ranallo-Benavidez TR, Jaron KS, Schatz MC . GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat Commun. 2020; 11: 1432

[56] Servant N, Varoquaux N, Lajoie BR, et al. HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015; 16: 259

[57] Dudchenko O, Batra SS, Omer AD, et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science. 2017; 356: 92-5

[58] Durand NC, Shamim MS, Machol I, et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 2016; 3: 95-8

[59] Hu J, Wang Z, Sun Z, et al. NextDenovo: an efficient error correction and accurate assembly tool for noisy long reads. Genome Biol. 2024; 25: 107

[60] Xu M, Guo L, Gu S, et al. TGS-GapCloser: a fast and accurate gap closer for large genomes with low coverage of error-prone long reads. Gigascience. 2020; 9: giaa094

[61] Jain C, Rhie A, Hansen NF, et al. Long-read mapping to repetitive reference sequences using Winnowmap2. Nat Methods. 2022; 19: 705-10

[62] Hu J, Wang Z, Liang F, et al. NextPolish2: a repeat-aware polishing tool for genomes assembled using HiFi long reads. Genom Proteom Bioinform. 2024; 22: qzad009

[63] Manni M, Berkeley MR, Seppey M, et al. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol Biol Evol. 2021; 38: 4647-54

[64] Rhie A, Walenz BP, Koren S, et al. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol. 2020; 21: 245

[65] Brown MR, Gonzalez M, de La Rosa P, et al. Tidk: a toolkit to rapidly identify telomeric repeats from genomic datasets. Bioinformatics. 2025; 41: btaf049

[66] Lin Y, Ye C, Li X, et al. quarTeT: a telomere-to-telomere toolkit for gap-free genome assembly and centromeric repeat identification. Hortic Res. 2023; 10: uhad127

[67] Bao Z, Eddy SR . Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 2002; 12: 1269-76

[68] Price AL, Jones NC, Pevzner PA . De novo identification of repeat families in large genomes. Bioinformatics. 2005; 21: i351-8

[69] Flynn JM, Hubley R, Goubert C, et al. RepeatModeler2 for automated genomic discovery of transposable element families. Proc Natl Acad Sci USA. 2020; 117: 9451-7

[70] Bao W, Kojima KK, Kohany O . Repbase update, a database of repetitive elements in eukaryotic genomes. Mob DNA. 2015; 6: 11

[71] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; 25: 4.10.1-4.10.14

[72] Stanke M, Keller O, Gunduz I, et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006; 34: W435-9

[73] Korf I. Gene finding in novel genomes. BMC Bioinformatics. 2004; 5: 59

[74] Majoros WH, Pertea M, Salzberg SL . TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics. 2004; 20: 2878-9

[75] Brůna T, Lomsadze A, Borodovsky M . GeneMark-ETP significantly improves the accuracy of automatic annotation of large eukaryotic genomes. Genome Res. 2024; 34: 757-68

[76] Keilwagen J, Hartung F, Grau J . GeMoMa: homology-based gene prediction utilizing intron position conservation and RNA-seq data. Methods Mol Biol. 2019; 1962: 161-77

[77] Slater GS, Birney E . Automated generation of heuristics for biological sequence comparison. BMC bioinformatics. 2005; 6: 31

[78] Haas BJ, Delcher AL, Mount SM, et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003; 31: 5654-66

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

[80] Camacho C, Coulouris G, Avagyan V, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10: 421

[81] Moriya Y, Itoh M, Okuda S, et al. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007; 35: W182-5

[82] Quevillon E, Silventoinen V, Pillai S, et al. InterProScan: protein domains identifier. Nucleic Acids Res. 2005; 33: W116-20

[83] Sun P, Jiao B, Yang Y, et al. WGDI: a user-friendly toolkit for evolutionary analyses of whole-genome duplications and ancestral karyotypes. Mol Plant. 2022; 15: 1841-51

[84] Emms DM, Kelly S . OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019; 20: 238

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

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

[87] Mendes FK, Vanderpool D, Fulton B, et al. CAFE 5 models variation in evolutionary rates among gene families. Bioinformatics. 2021; 36: 5516-8

[88] Yu G, Wang LG, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7

[89] Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34: 3094-100

[90] Goel M, Sun H, Jiao WB, et al. SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies. Genome Biol. 2019;20:277

[91] Goel M, Schneeberger K . plotsr: visualizing structural similarities and rearrangements between multiple genomes. Bioinformatics. 2022; 38: 2922-6

[92] Thorvaldsdóttir H, Robinson JT, Mesirov JP . Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013; 14: 178-92

[93] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20

[94] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60

[95] Danecek P, Bonfield JK, Liddle J, et al. Twelve years of SAMtools and BCFtools. Gigascience. 2021; 10: giab008

[96] Danecek P, Auton A, Abecasis G, et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8

[97] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64

[98] Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81: 559-75

[99] Zhang C, Dong SS, Xu JY, et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019; 35: 1786-8

[100] Li H, Durbin R . Inference of human population history from individual whole-genome sequences. Nature. 2011; 475: 493-6

[101] Terhorst J, Kamm JA, Song YS . Robust and scalable inference of population history from hundreds of unphased whole genomes. Nat Genet. 2017; 49: 303-9

[102] Pickrell JK, Pritchard JK . Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012; 8: e1002967

[103] Fitak RR . OptM: estimating the optimal number of migration edges on population trees using Treemix. Biol Methods Protoc. 2021; 6: bpab017

[104] Malinsky M, Matschiner M, Svardal H . Dsuite-Fast D-statistics and related admixture evidence from VCF files. Mol Ecol Resour. 2021; 21: 584-95

[105] Fick SE, Hijmans RJ . WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017; 37: 4302-15

[106] Phillips SJ, Anderson RP, Schapire RE . Maximum entropy modeling of species geographic distributions. Ecol Model. 2006; 190: 231-59

[107] Kou YX, Liu ML, López-Pujol J, et al. Contrasting demographic history and mutational load in three threatened whitebark pines (Pinus subsect. Gerardianae): implications for conservation. Plant J. 2024; 119: 2967-81

[108] Fielding AH, Bell JF . A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997; 24: 38-49

[109] Pavlidis P, Živkovic D, Stamatakis A, et al. SweeD: likelihood-based detection of selective sweeps in thousands of genomes. Mol Biol Evol. 2013; 30: 2224-34

[110] Kim D, Paggi JM, Park C, et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37: 907-15

[111] Liao Y, Smyth GK, Shi W . featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014; 30: 923- 30

[112] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550