A High-Quality Melon Genome Assembly Provides Insights into Genetic Basis of Fruit Trait Improvement

The β-galactosidase gene AtrBGAL2 regulates Akebia trifoliata fruit cracking

Cracking of Akebia trifoliata fruit at maturity is problematic for the cultivation of the horticultural crop, shortening shelf-life quality and compromising commercial value. However, the molecular mechanisms underlying this feature of A. trifoliata are not known. Genes involved in cell wall metabolism were identified by genome and transcriptome sequencing, which may play important roles in fruit cracking. One of the galactose metabolism related genes, β-galactosidase (AtrBGAL2), was identified in A. trifoliata, and overexpression (OE) of AtrBGAL2 resulted in early fruit cracking, higher water-soluble pectin contents, and lower acid-soluble pectin, cellulose, and hemicellulose content compared to the wild type. Whereas silencing of AtrBGAL2 in trifoliata by virus induced gene silencing showed opposite trends. The levels of AtrBGAL2 transcripts were 24.6 and 66.0-fold higher in OE A. trifoliata and tomato fruits, respectively, and the cell wall-related genes were also gradually greater than in control plants during fruit ripening. Whereas the expression levels of AtrBGAL2 was significantly down-regulated by 54.1 % and 73.7 % in gene silenced A. trifoliata and CRISPR/Cas9 tomato mutant plants, respectively, and cell wall-related genes were also significantly reduced. These results demonstrate that AtrBGAL2 plays important roles in regulating fruit cracking during fruit ripening.

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

A wild melon reference genome provides novel insights into the domestication of a key gene responsible for melon fruit acidity

KEY MESSAGE

A wild melon reference genome elucidates the genomic basis of fruit acidity domestication. Structural variants (SVs) have been reported to impose major effects on agronomic traits, representing a significant contributor to crop domestication. However, the landscape of SVs between wild and cultivated melons is elusive and how SVs have contributed to melon domestication remains largely unexplored. Here, we report a 379-Mb chromosome-scale genome of a wild progenitor melon accession "P84", with a contig N50 of 14.9 Mb. Genome comparison identifies 10,589 SVs between P84 and four cultivated melons with 6937 not characterized in previously analysis of 25 melon genome sequences. Furthermore, the population-scale genotyping of these SVs was determined in 1175 accessions, and 18 GWAS signals including fruit acidity, fruit length, fruit weight, fruit color and sex determination were detected. Based on these genotyped SVs, we identified 3317 highly diverged SVs between wild and cultivated melons, which could be the potential SVs associated with domestication-related traits. Furthermore, we identify novel SVs affecting fruit acidity and proposed the diverged evolutionary trajectories of CmPH, a key regulator of melon fruit acidity, during domestication and selection of different populations. These results will offer valuable resources for genomic studies and genetic improvement in melon.

Crop Landraces and Indigenous Varieties: A Valuable Source of Genes for Plant Breeding

Landraces and indigenous varieties comprise valuable sources of crop species diversity. Their utilization in plant breeding may lead to increased yield and enhanced quality traits, as well as resilience to various abiotic and biotic stresses. Recently, new approaches based on the rapid advancement of genomic technologies such as deciphering of pangenomes, multi-omics tools, marker-assisted selection (MAS), genome-wide association studies (GWAS), and CRISPR/Cas9 gene editing greatly facilitated the exploitation of landraces in modern plant breeding. In this paper, we present a comprehensive overview of the implementation of new genomic technologies and highlight their importance in pinpointing the genetic basis of desirable traits in landraces and indigenous varieties of annual, perennial herbaceous, and woody crop species cultivated in the Mediterranean region. The need for further employment of advanced -omic technologies to unravel the full potential of landraces and indigenous varieties underutilized genetic diversity is also indicated. Ultimately, the large amount of genomic data emerging from the investigation of landraces and indigenous varieties reveals their potential as a source of valuable genes and traits for breeding. The role of landraces and indigenous varieties in mitigating the ongoing risks posed by climate change in agriculture and food security is also highlighted.

Telomere-to-telomere genome assembly of melon (Cucumis melo L. var. inodorus) provides a high-quality reference for meta-QTL analysis of important traits

Melon is an important horticultural crop with extensive diversity in many horticultural groups. To explore its genomic diversity, it is necessary to assemble more high-quality complete genomes from different melon accessions. Meanwhile, a large number of QTLs have been mapped in several studies. Integration of the published QTLs onto a complete genome can provide more accurate information for candidate gene cloning. To address these problems, a telomere-to-telomere (T2T) genome of the elite melon landrace Kuizilikjiz (Cucumis melo L. var. inodorus) was de novo assembled and all the published QTLs were projected onto it in this study. The results showed that a high-quality Kuizilikjiz genome with the size of 379.2 Mb and N50 of 31.7 Mb was de novo assembled using the combination of short reads, PacBio high-fidelity long reads, Hi-C data, and a high-density genetic map. Each chromosome contained the centromere and telomeres at both ends. A large number of structural variations were observed between Kuizilikjiz and the other published genomes. A total of 1294 QTLs published in 67 studies were collected and projected onto the T2T genome. Several clustered, co-localized, and overlapped QTLs were determined. Furthermore, 20 stable meta-QTLs were identified, which significantly reduced the mapping intervals of the initial QTLs and greatly facilitated identification of the candidate genes. Collectively, the T2T genome assembly together with the numerous projected QTLs will not only broaden the high-quality genome resources but also provide valuable and abundant QTL information for cloning the genes controlling important traits in melon.

The haplotype-resolved T2T reference genome highlights structural variation underlying agronomic traits of melon

Melon (Cucumis melo L.) is an important vegetable crop that has an extensive history of cultivation. However, the genome of wild and semi-wild melon types that can be used for the analysis of agronomic traits is not yet available. Here we report a chromosome-level T2T genome assembly for 821 (C. melo ssp. agrestis var. acidulus), a semi-wild melon with two haplotypes of ~373 Mb and ~364 Mb, respectively. Comparative genome analysis discovered a significant number of structural variants (SVs) between melo (C. melo ssp. melo) and agrestis (C. melo ssp. agrestis) genomes, including a copy number variation located in the ToLCNDV resistance locus on chromosome 11. Genome-wide association studies detected a significant signal associated with climacteric ripening and identified one candidate gene CM_ac12g14720.1 (CmABA2), encoding a cytoplasmic short chain dehydrogenase/reductase, which controls the biosynthesis of abscisic acid. This study provides valuable genetic resources for future research on melon breeding.

Pan-genome analysis sheds light on structural variation-based dissection of agronomic traits in melon crops

Sweetness and appearance of fresh fruits are key palatable and preference attributes for consumers and are often controlled by multiple genes. However, fine-mapping the key loci or genes of interest by single genome-based genetic analysis is challenging. Herein, we present the chromosome-level genome assembly of 1 landrace melon accession (Cucumis melo ssp. agrestis) with wild morphologic features and thus construct a melon pan-genome atlas via integrating sequenced melon genome datasets. Our comparative genomic analysis reveals a total of 3.4 million genetic variations, of which the presence/absence variations (PAVs) are mainly involved in regulating the function of genes for sucrose metabolism during melon domestication and improvement. We further resolved several loci that are accountable for sucrose contents, flesh color, rind stripe, and suture using a structural variation (SV)-based genome-wide association study. Furthermore, via bulked segregation analysis (BSA)-seq and map-based cloning, we uncovered that a single gene, (CmPIRL6), determines the edible or inedible characteristics of melon fruit exocarp. These findings provide important melon pan-genome information and provide a powerful toolkit for future pan-genome-informed cultivar breeding of melon.

Identification of novel sex determination loci in Japanese weedy melon

KEY MESSAGE

Japanese weedy melon exhibits unique sex expression with interactions between previously reported sex determination genes and two novel loci. Sex expression contributes to fruit quality and yield in the Cucurbitaceae. In melon, orchestrated regulation by sex determination genes explains the mechanism of sex expression, resulting in a great variety of sexual morphologies. In this study, we examined the Japanese weedy melon UT1, which does not follow the reported model of sex expression. We conducted QTL analysis using F₂ plants for flower sex on the main stem and the lateral branch and mapped "occurrence of pistil-bearing flower on the main stem" locus on Chr. 3 (Opbf3.1) and "type of pistil-bearing flower" (female or bisexual) loci on Chr. 2 (tpbf2.1) and Chr. 8 (tpbf8.1). The Opbf3.1 included the known sex determination gene CmACS11. Sequence comparison of CmACS11 between parental lines revealed three nonsynonymous SNPs. A CAPS marker developed from one of the SNPs was closely linked to the occurrence of pistil-bearing flowers on the main stem in two F₂ populations with different genetic backgrounds. The UT1 allele on Opbf3.1 was dominant in F₁ lines from crosses between UT1 and diverse cultivars and breeding lines. This study suggests that Opbf3.1 and tpbf8.1 may promote the development of pistil and stamen primordia by inhibiting CmWIP1 and CmACS-7 functions, respectively, making the UT1 plants hermaphrodite. The results of this study provide new insights into the molecular mechanisms of sex determination in melons and considerations for the application of femaleness in melon breeding.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo)

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F₂ populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

CuGenDBv2: an updated database for cucurbit genomics

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel 'Genotype' module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated 'Expression' module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

Cantaloupe melon genome reveals 3D chromatin features and structural relationship with the ancestral cucurbitaceae karyotype

Cucumis melo displays a large diversity of horticultural groups with cantaloupe melon the most cultivated type. Using a combination of single-molecule sequencing, 10X Genomics link-reads, high-density optical and genetic maps, and chromosome conformation capture (Hi-C), we assembled a chromosome scale C. melo var. cantalupensis Charentais mono genome. Integration of RNA-seq, MeDip-seq, ChIP-seq, and Hi-C data revealed a widespread compartmentalization of the melon genome, segregating constitutive heterochromatin and euchromatin. Genome-wide comparative and evolutionary analysis between melon botanical groups identified Charentais mono genome increasingly more divergent from Harukei-3 (reticulatus), Payzawat (inodorus), and HS (ssp. agrestis) genomes. To assess the paleohistory of the Cucurbitaceae, we reconstructed the ancestral Cucurbitaceae karyotype and compared it to sequenced cucurbit genomes. In contrast to other species that experienced massive chromosome shuffling, melon has retained the ancestral genome structure. We provide comprehensive genomic resources and new insights in the diversity of melon horticultural groups and evolution of cucurbits.

•

We provide a chromosome scale C. melo var. cantalupensis Charentais mono genome

•

Epigenomic analysis revealed a widespread compartmentalization of the melon genome

•

We reconstructed the ancestral Cucurbitaceae karyotype

•

Melon has retained the ancestral Cucurbitaceae genome structure

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

QTL mapping and genomic analyses of earliness and fruit ripening traits in a melon Recombinant Inbred Lines population supported by de novo assembly of their parental genomes

Earliness and ripening behavior are important attributes of fruits on and off the vine, and affect quality and preference of both growers and consumers. Fruit ripening is a complex physiological process that involves metabolic shifts affecting fruit color, firmness, and aroma production. Melon is a promising model crop for the study of fruit ripening, as the full spectrum of climacteric behavior is represented across the natural variation. Using Recombinant Inbred Lines (RILs) population derived from the parental lines "Dulce" (reticulatus, climacteric) and "Tam Dew" (inodorus, non-climacteric) that vary in earliness and ripening traits, we mapped QTLs for ethylene emission, fruit firmness and days to flowering and maturity. To further annotate the main QTL intervals and identify candidate genes, we used Oxford Nanopore long-read sequencing in combination with Illumina short-read resequencing, to assemble the parental genomes de-novo. In addition to 2.5 million genome-wide SNPs and short InDels detected between the parents, we also highlight here the structural variation between these lines and the reference melon genome. Through systematic multi-layered prioritization process, we identified 18 potential polymorphisms in candidate genes within multi-trait QTLs. The associations of selected SNPs with earliness and ripening traits were further validated across a panel of 177 diverse melon accessions and across a diallel population of 190 F1 hybrids derived from a core subset of 20 diverse parents. The combination of advanced genomic tools with diverse germplasm and targeted mapping populations is demonstrated as a way to leverage forward genetics strategies to dissect complex horticulturally important traits.

Chromosome-level genome assembly of Gynostemma pentaphyllum provides insights into gypenoside biosynthesis

Gynostemma pentaphyllum (Thunb.) Makino is an economically valuable medicinal plant belonging to the Cucurbitaceae family that produces the bioactive compound gypenoside. Despite several transcriptomes having been generated for G. pentaphyllum, a reference genome is still unavailable, which has limited the understanding of the gypenoside biosynthesis and regulatory mechanism. Here, we report a high-quality G. pentaphyllum genome with a total length of 582 Mb comprising 1,232 contigs and a scaffold N50 of 50.78 Mb. The G. pentaphyllum genome comprised 59.14% repetitive sequences and 25,285 protein-coding genes. Comparative genome analysis revealed that G. pentaphyllum was related to Siraitia grosvenorii, with an estimated divergence time dating to the Paleogene (∼48 million years ago). By combining transcriptome data from seven tissues, we reconstructed the gypenoside biosynthetic pathway and potential regulatory network using tissue-specific gene co-expression network analysis. Four UDP-glucuronosyltransferases (UGTs), belonging to the UGT85 subfamily and forming a gene cluster, were involved in catalyzing glycosylation in leaf-specific gypenoside biosynthesis. Furthermore, candidate biosynthetic genes and transcription factors involved in the gypenoside regulatory network were identified. The genetic information obtained in this study provides insights into gypenoside biosynthesis and lays the foundation for further exploration of the gypenoside regulatory mechanism.

High-quality chromosome-level genomes of Cucumis metuliferus and Cucumis melo provide insight into Cucumis genome evolution

Cucumis metuliferus (African horned cucumber), a wild relative of Cucumis sativus (cucumber) and Cucumis melo (melon), displays high-level resistance to several important plant pathogens (e.g., root-knot nematodes and several viruses). Here, we report a chromosome-level genome assembly for C. metuliferus, with a 316 Mb genome sequence comprising 29 039 genes. Phylogenetic analysis of related species in family Cucurbitaceae indicated that the divergence time between C. metuliferus and melon was 17.8 million years ago. Comparisons between the C. metuliferus and melon genomes revealed large structural variations (inversions and translocations >1 Mb) in eight chromosomes of these two species. Gene family comparison showed that C. metuliferus has the largest number of resistance-related nucleotide-binding site leucine-rich repeat (NBS-LRR) genes in Cucurbitaceae. The loss of NBS-LRR loci caused by large insertions or deletions (indels) and pseudogenization caused by small indels explained the loss of NBS-LRR genes in Cucurbitaceae. Population structure analysis suggested that C. metuliferus originated in Zimbabwe, then spread to other southern African regions where it likely underwent similar domestic selection as melon. This C. metuliferus reference sequence will accelerate the understanding of the molecular evolution of resistance-related genes and enhance cucurbit crop improvement efforts.

Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits

Most molecularly characterized plant resistance genes (R genes) belong to the nucleotide-binding-site-leucine-rich-repeat (NLR) receptor family and are prone to duplication and transposition with high sequence diversity. In this family, the Vat gene in melon is one of the few R genes known for conferring resistance to insect, i.e., Aphis gossypii, but it has been misassembled and/or mispredicted in the whole genomes of Cucurbits. We examined 14 genomic regions (about 400 kb) derived from long-read assemblies spanning Vat-related genes in Cucumis melo, Cucumis sativus, Citrullus lanatus, Benincasa hispida, Cucurbita argyrosperma, and Momordica charantia. We built the phylogeny of those genes. Investigating the paleohistory of the Vat gene cluster, we revealed a step by step process beginning from a common ancestry in cucurbits older than 50 my. We highlighted Vat exclusively in the Cucumis genera, which diverged about 20 my ago. We then focused on melon, evaluating a minimum duplication rate of Vat in 80 wild and cultivated melon lines using generalist primers; our results suggested that duplication started before melon domestication. The phylogeny of 44 Vat-CDS obtained from 21 melon lines revealed gain and loss of leucine-rich-repeat domains along diversification. Altogether, we revealed the high putative recognition scale offered in melon based on a combination of SNPs, number of leucine-rich-repeat domains within each homolog and number of homologs within each cluster that might jointly confer resistance to a large pest and pathogen spectrum. Based on our findings, we propose possible avenues for breeding programs.

Combined genomic, transcriptomic, and metabolomic analyses provide insights into chayote (Sechium edule) evolution and fruit development

Chayote (Sechium edule) is an agricultural crop in the Cucurbitaceae family that is rich in bioactive components. To enhance genetic research on chayote, we used Nanopore third-generation sequencing combined with Hi-C data to assemble a draft chayote genome. A chromosome-level assembly anchored on 14 chromosomes (N50 contig and scaffold sizes of 8.40 and 46.56 Mb, respectively) estimated the genome size as 606.42 Mb, which is large for the Cucurbitaceae, with 65.94% (401.08 Mb) of the genome comprising repetitive sequences; 28,237 protein-coding genes were predicted. Comparative genome analysis indicated that chayote and snake gourd diverged from sponge gourd and that a whole-genome duplication (WGD) event occurred in chayote at 25 ± 4 Mya. Transcriptional and metabolic analysis revealed genes involved in fruit texture, pigment, flavor, flavonoids, antioxidants, and plant hormones during chayote fruit development. The analysis of the genome, transcriptome, and metabolome provides insights into chayote evolution and lays the groundwork for future research on fruit and tuber development and genetic improvements in chayote.

The genome and transcriptome analysis of snake gourd provide insights into its evolution and fruit development and ripening

Snake gourd (Trichosanthes anguina L.), which belongs to the Cucurbitaceae family, is a popular ornamental and food crop species with medicinal value and is grown in many parts of the world. Although progress has been made in its genetic improvement, the organization, composition, and evolution of the snake gourd genome remain largely unknown. Here, we report a high-quality genome assembly for snake gourd, comprising 202 contigs, with a total size of 919.8 Mb and an N50 size of 20.1 Mb. These findings indicate that snake gourd has one of the largest genomes of Cucurbitaceae species sequenced to date. The snake gourd genome assembly harbors 22,874 protein-coding genes and 80.0% of the genome consists of repetitive sequences. Phylogenetic analysis reveals that snake gourd is closely related to sponge gourd but diverged from their common ancestor ~33-47 million years ago. The genome sequence reported here serves as a valuable resource for snake gourd genetic research and comparative genomic studies in Cucurbitaceae and other plant species. In addition, fruit transcriptome analysis reveals the candidate genes related to quality traits during snake gourd fruit development and provides a basis for future research on snake gourd fruit development and ripening at the transcript level.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

The Chromosome-Scale Genome of Melon Dissects Genetic Architecture of Important Agronomic Traits

Comparative and evolutionary genomics analyses are the powerful tools to provide mechanistic insights into important agronomic traits. Here, we completed a chromosome-scale assembly of the "neglected" but vital melon subspecies Cucumis melo ssp. agrestis using single-molecule real-time sequencing, Hi-C, and an ultra-dense genetic map. Comparative genomics analyses identified two targeted genes, UDP-sugar pyrophosphorylase and α-galactosidase, that were selected during evolution for specific phloem transport of oligosaccharides in Cucurbitaceae. Association analysis of transcriptome and the DNA methylation patterns revealed the epigenetic regulation of sucrose accumulation in developing fruits. We constructed the melon recombinant inbred lines to uncover Alkaline/Neutral Invertase (CINV), Sucrose-Phosphatase 2 (SPP2), α-galactosidase, and β-galactosidase loci related to sucrose accumulation and an LRR receptor-like serine/threonine-protein kinase associated with gummy stem blight resistance. This study provides essential genomic resources enabling functional genomics studies and the genomics-informed breeding pipelines for improving the fruit quality and disease resistance traits.

Mapping a Partial Andromonoecy Locus in Citrullus lanatus Using BSA-Seq and GWAS Approaches

The sexual expression of watermelon plants is the result of the distribution and occurrence of male, female, bisexual and hermaphrodite flowers on the main and secondary stems. Plants can be monoecious (producing male and female flowers), andromonoecious (producing male and hermaphrodite flowers), or partially andromonoecious (producing male, female, bisexual, and hermaphrodite flowers) within the same plant. Sex determination of individual floral buds and the distribution of the different flower types on the plant, are both controlled by ethylene. A single missense mutation in the ethylene biosynthesis gene CitACS4, is able to promote the conversion of female into hermaphrodite flowers, and therefore of monoecy (genotype MM) into partial andromonoecy (genotype Mm) or andromonoecy (genotype mm). We phenotyped and genotyped, for the M/m locus, a panel of 207 C. lanatus accessions, including five inbreds and hybrids, and found several accessions that were repeatedly phenotyped as PA (partially andromonoecious) in several locations and different years, despite being MM. A cosegregation analysis between a SNV in CitACS4 and the PA phenotype, demonstrated that the occurrence of bisexual and hermaphrodite flowers in a PA line is not dependent on CitACS4, but conferred by an unlinked recessive gene which we called pa. Two different approaches were performed to map the pa gene in the genome of C. lanatus: bulk segregant analysis sequencing (BSA-seq) and genome wide association analysis studies (GWAS). The BSA-seq study was performed using two contrasting bulks, the monoecious M-bulk and the partially andromonoecious PA-bulk, each one generated by pooling DNA from 20 F2 plants. For GWAS, 122 accessions from USDA gene bank, already re-sequenced by genotyping by sequencing (GBS), were used. The combination of the two approaches indicates that pa maps onto a genomic region expanding across 32.24-36.44 Mb in chromosome 1 of watermelon. Fine mapping narrowed down the pa locus to a 867 Kb genomic region containing 101 genes. A number of candidate genes were selected, not only for their function in ethylene biosynthesis and signalling as well as their role in flower development and sex determination, but also by the impact of the SNPs and indels differentially detected in the two sequenced bulks.