Genome-Wide Linkage-Disequilibrium Mapping to the Candidate Gene Level in Melon (Cucumis melo)

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.

Genetic regulation of volatile production in two melon introgression line collections with contrasting ripening behavior

The importance of melon aroma in determining fruit quality has been highlighted in recent years. The fruit volatile profile is influenced by the type of fruit ripening. Non-climacteric fruits contain predominantly aldehydes, while climacteric fruits mainly produce esters. Several genes have been described to participate in volatile organic compounds (VOCs) biosynthesis pathways, but knowledge in this area is still incomplete. In this work we analysed the volatile profile of two reciprocal Introgression Line (IL) collections generated from a cross between 'Piel de Sapo' (PS) and 'Védrantais' (VED) melons, differing in their aroma profile and ripening behaviour. SPME GC-MS was performed to identify genes responsible for VOCs formation. More than 1000 QTLs for many volatiles were detected taken together both populations. Introgressions on chromosomes 3, 5, 6, 7 and 8 modified ester-aldehyde balance and were correlated to ripening changes in both genetic backgrounds. Some previously identified QTLs for fruit ripening might be involved in these phenotypes, such as ETHQV8.1 on chromosome 8 and ETHQV6.3 on chromosome 6. PS alleles on chromosomes 2, 6, 10 and 11 were found to increase ester content when introgressed in VED melons. Terpenes showed to be affected by several genomic regions not related to ripening. In addition, several candidate genes have been hypothesized to be responsible for some of the QTLs detected. The analysis of volatile compounds in two reciprocal IL collections has increased our understanding of the relationship between ripening and aroma and offers valuable plant material to improve food quality in melon breeding programs.

Molecular Marker-Assisted Mapping, Candidate Gene Identification, and Breeding in Melon (Cucumis melo L.): A Review

Melon (Cucumis melo L.) is an important crop that is cultivated worldwide for its fleshy fruit. Understanding the genetic basis of a plant's qualitative and quantitative traits is essential for developing consumer-favored varieties. This review presents genetic and molecular advances related to qualitative and quantitative phenotypic traits and biochemical compounds in melons. This information guides trait incorporation and the production of novel varieties with desirable horticultural and economic characteristics and yield performance. This review summarizes the quantitative trait loci, candidate genes, and development of molecular markers related to plant architecture, branching patterns, floral attributes (sex expression and male sterility), fruit attributes (shape, rind and flesh color, yield, biochemical compounds, sugar content, and netting), and seed attributes (seed coat color and size). The findings discussed in this review will enhance demand-driven breeding to produce cultivars that benefit consumers and melon breeders.

QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon

Melon (Cucumis melo L.) is an economically important Cucurbitaceae crop grown around the globe. The sweetness of melon is a significant factor in fruit quality and consumer appeal, and the soluble solids content (SSC) is a key index of melon sweetness. In this study, 146 recombinant inbred lines (RILs) derived from two oriental melon materials with different levels of sweetness containing 1427 bin markers, and 213 melon accessions containing 1,681,775 single nucleotide polymorphism (SNP) markers were used to identify genomic regions influencing SSC. Linkage mapping detected 10 quantitative trait loci (QTLs) distributed on six chromosomes, seven of which were overlapped with the reported QTLs. A total of 211 significant SNPs were identified by genome-wide association study (GWAS), 138 of which overlapped with the reported QTLs. Two new stable, co-localized regions on chromosome 3 were identified by QTL mapping and GWAS across multiple environments, which explained large phenotypic variance. Five candidate genes related to SSC were identified by QTL mapping, GWAS, and qRT-PCR, two of which were involved in hydrolysis of raffinose and sucrose located in the new stable loci. The other three candidate genes were involved in raffinose synthesis, sugar transport, and production of substrate for sugar synthesis. The genomic regions and candidate genes will be helpful for molecular breeding programs and elucidating the mechanisms of sugar accumulation.

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.

Fruit Morphology and Ripening-Related QTLs in a Newly Developed Introgression Line Collection of the Elite Varieties 'Védrantais' and 'Piel de Sapo'

Melon is an economically important crop with widely diverse fruit morphology and ripening characteristics. Its diploid sequenced genome and multiple genomic tools make this species suitable to study the genetic architecture of fruit traits. With the development of this introgression line population of the elite varieties 'Piel de Sapo' and 'Védrantais', we present a powerful tool to study fruit morphology and ripening traits that can also facilitate characterization or pyramidation of QTLs in inodorous melon types. The population consists of 36 lines covering almost 98% of the melon genome, with an average of three introgressions per chromosome and segregating for multiple fruit traits: morphology, ripening and quality. High variability in fruit morphology was found within the population, with 24 QTLs affecting six different traits, confirming previously reported QTLs and two newly detected QTLs, FLQW5.1 and FWQW7.1. We detected 20 QTLs affecting fruit ripening traits, six of them reported for the first time, two affecting the timing of yellowing of the rind (EYELLQW1.1 and EYELLQW8.1) and four at the end of chromosome 8 affecting aroma, abscission and harvest date (EAROQW8.3, EALFQW8.3, ABSQW8.3 and HARQW8.3). We also confirmed the location of several QTLs, such as fruit-quality-related QTLs affecting rind and flesh appearance and flesh firmness.

Spanish Melon Landraces: Revealing Useful Diversity by Genomic, Morphological, and Metabolomic Analysis

Spain is a secondary centre of the diversification of the melon (Cucumis melo L.), with high diversity represented in highly appreciated landraces belonging to the Flexuosus and Ibericus groups. A collection of 47 accessions of Flexuosus, Chate, Piel de Sapo, Tendral, Amarillo, Blanco, and Rochet was analysed using a genotyping-by-sequencing (GBS) approach. A total of 66,971 quality SNPs were identified. Genetic analysis differentiated Ibericus accessions and exotic materials (Ameri, Momordica, Kachri, and Agrestis), while Flexuous accessions shared ancestry between them. Within the Ibericus group, no clear genomic distinction could be identified for the different landraces evaluated, with accessions of different landraces showing high genetic similarity. The morphological characterization confirmed that the external colour and fruit shape had been used as recognition patterns for Spanish melon landraces, but variability within a landrace exists. Differences were found in the sugars and acid and volatile profiles of the materials. Flexuosus and Chate melons at the immature commercial stage accumulated malic acid and low levels of hexoses, while Ibericus melons accumulated high contents of sucrose and citric acid. Specific trends could be identified in the Ibericus landraces. Tendral accumulated low levels of sugars and citric acid and high of malic acid, maintaining higher firmness, Rochet reached higher levels of sugars, and Amarillo tended to lower malic acid contents. Interestingly, high variability was found within landraces for the acidic profile, offering possibilities to alter taste tinges. The main volatile organic compounds (VOCs) in Flexuosus and Chate were aldehydes and alcohols, with clear differences between both groups. In the Ibericus landraces, general trends for VOC accumulation could be identified, but, again, a high level of variation exists. This situation highlights the necessity to develop depuration programs to promote on-farm in situ conservation and, at the same time, offers opportunities to establish new breeding program targets and to take advantage of these sources of variation.

Fine-mapping and identification of a candidate gene controlling seed coat color in melon (Cucumis melo L. var. chinensis Pangalo)

MELO3C019554 encoding a homeobox protein (PHD transcription factor) is a candidate gene that involved in the formation of seed coat color in melon. Seed coat color is related to flavonoid content which is closely related to seed dormancy. According to the genetic analysis of a six-generation population derived from two parents (IC2508 with a yellow seed coat and IC2518 with a brown seed coat), we discovered that the yellow seed coat trait in melon is controlled by a single dominant gene, named CmBS-1. Bulked segregant analysis sequencing (BSA-Seq) revealed that the gene is located at 11,860,000-15,890,000 bp (4.03 Mb) on Chr 6. The F₂ population was genotyped using insertion-deletions (InDels), from which cleaved amplified polymorphic sequence (dCAPS) markers were derived to construct a genetic map. The gene was then fine-mapped to a 233.98 kb region containing 12 genes. Based on gene sequence analysis with two parents, we found that the MELO3C019554 gene encoding a homeobox protein (PHD transcription factor) had a nonsynonymous single nucleotide polymorphism (SNP) mutation in the coding sequence (CDS), and the SNP mutation resulted in the conversion of an amino acid (A → T) at residue 534. In addition, MELO3C019554 exhibited lower relative expression levels in the yellow seed coat than in the brown seed coat. Furthermore, we found that MELO3C019554 is related to 12 flavonoid metabolites. Thus, we predicted that MELO3C019554 is a candidate gene controlling seed coat color in melon. The study lays a foundation for further cloning projects and functional analysis of this gene, as well as marker-assisted selection breeding.

A cryptic variation in a member of the Ovate Family Proteins is underlying the melon fruit shape QTL fsqs8.1

The gene underlying the melon fruit shape QTL fsqs8.1 is a member of the Ovate Family Proteins. Variation in fruit morphology is caused by changes in gene expression likely due to a cryptic structural variation in this locus. Melon cultivars have a wide range of fruit morphologies. Quantitative trait loci (QTL) have been identified underlying such diversity. This research focuses on the fruit shape QTL fsqs8.1, previously detected in a cross between the accession PI 124112 (CALC, producing elongated fruit) and the cultivar 'Piel de Sapo' (PS, producing oval fruit). The CALC fsqs8.1 allele induced round fruit shape, being responsible for the transgressive segregation for this trait observed in that population. In fact, the introgression line CALC8-1, carrying the fsqs8.1 locus from CALC into the PS genetic background, produced perfect round fruit. Following a map-based cloning approach, we found that the gene underlying fsqs8.1 is a member of the Ovate Family Proteins (OFP), CmOFP13, likely a homologue of AtOFP1 and SlOFP20 from Arabidopsis thaliana and tomato, respectively. The induction of the round shape was due to the higher expression of the CALC allele at the early ovary development stage. The fsqs8.1 locus showed an important structural variation, being CmOFP13 surrounded by two deletions in the CALC genome. The deletions are present at very low frequency in melon germplasm. Deletions and single nucleotide polymorphisms in the fsqs8.1 locus could not be not associated with variation in fruit shape among different melon accessions, what indicates that other genetic factors should be involved to induce the CALC fsqs8.1 allele effects. Therefore, fsqs8.1 is an example of a cryptic variation that alters gene expression, likely due to structural variation, resulting in phenotypic changes in melon fruit morphology.

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.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Underground heterosis for yield improvement in melon

Heterosis, the superiority of hybrids over their parents, is a major genetic force associated with plant fitness and crop yield enhancement. We investigated root-mediated yield heterosis in melons (Cucumis melo) by characterizing a common variety grafted onto 190 hybrid rootstocks, resulting from crossing 20 diverse inbreds in a diallel-mating scheme. Hybrid rootstocks improved yield by more than 40% compared with their parents, and the best hybrid yield outperformed the reference commercial variety by 65% under both optimal and minimal irrigation treatments. To characterize the genetics of underground heterosis we conducted whole genome re-sequencing of the 20 founder lines, and showed that parental genetic distance was no predictor for the level of heterosis. Through inference of the 190 hybrid genotypes from their parental genomes, followed by genome-wide association analysis, we mapped multiple quantitative trait loci for root-mediated yield. Yield enhancement of the four best-performing hybrid rootstocks was validated in multiple experiments with four different scion varieties. Our grafting approach is complementary to the common roots genetic approach that focuses mainly on variation in root system architecture, and is a step towards discovery of candidate genes involved in root function and yield enhancement.

Genetic characterization of melon accessions in the U.S. National Plant Germplasm System and construction of a melon core collection

Melon (C. melo L.) is an economically important vegetable crop cultivated worldwide. The melon collection in the U.S. National Plant Germplasm System (NPGS) is a valuable resource to conserve natural genetic diversity and provide novel traits for melon breeding. Here we use the genotyping-by-sequencing (GBS) technology to characterize 2083 melon accessions in the NPGS collected from major melon production areas as well as regions where primitive melons exist. Population structure and genetic diversity analyses suggested that C. melo ssp. melo was firstly introduced from the centers of origin, Indian and Pakistan, to Central and West Asia, and then brought to Europe and Americas. C. melo ssp. melo from East Asia was likely derived from C. melo ssp. agrestis in India and Pakistan and displayed a distinct genetic background compared to the rest of ssp. melo accessions from other geographic regions. We developed a core collection of 383 accessions capturing more than 98% of genetic variation in the germplasm, providing a publicly accessible collection for future research and genomics-assisted breeding of melon. Thirty-five morphological characters investigated in the core collection indicated high variability of these characters across accessions in the collection. Genome-wide association studies using the core collection panel identified potentially associated genome regions related to fruit quality and other horticultural traits. This study provides insights into melon origin and domestication, and the constructed core collection and identified genome loci potentially associated with important traits provide valuable resources for future melon research and breeding.

Genome-Wide SNP Markers for Genotypic and Phenotypic Differentiation of Melon (Cucumis melo L.) Varieties Using Genotyping-by-Sequencing

Melon (Cucumis melo L.) is an economically important horticultural crop with abundant morphological and genetic variability. Complex genetic variations exist even among melon varieties and remain unclear to date. Therefore, unraveling the genetic variability among the three different melon varieties, muskmelon (C. melo subsp. melo), makuwa (C. melo L. var. makuwa), and cantaloupes (C. melo subsp. melo var. cantalupensis), could provide a basis for evolutionary research. In this study, we attempted a systematic approach with genotyping-by-sequencing (GBS)-derived single nucleotide polymorphisms (SNPs) to reveal the genetic structure and diversity, haplotype differences, and marker-based varieties differentiation. A total of 6406 GBS-derived SNPs were selected for the diversity analysis, in which the muskmelon varieties showed higher heterozygote SNPs. Linkage disequilibrium (LD) decay varied significantly among the three melon varieties, in which more rapid LD decay was observed in muskmelon (r² = 0.25) varieties. The Bayesian phylogenetic tree provided the intraspecific relationships among the three melon varieties that formed, as expected, individual clusters exhibiting the greatest genetic distance based on the posterior probability. The haplotype analysis also supported the phylogeny result by generating three major networks for 48 haplotypes. Further investigation for varieties discrimination allowed us to detect a total of 52 SNP markers that discriminated muskmelon from makuwa varieties, of which two SNPs were converted into cleaved amplified polymorphic sequence markers for practical use. In addition to these markers, the genome-wide association study identified two SNPs located in the genes on chromosome 6, which were significantly associated with the phenotypic traits of melon seed. This study demonstrated that a systematic approach using GBS-derived SNPs could serve to efficiently classify and manage the melon varieties in the genebank.

Comprehensive Analysis of Transcriptome and Metabolome Reveals the Flavonoid Metabolic Pathway Is Associated with Fruit Peel Coloration of Melon

The peel color is an important external quality of melon fruit. To explore the mechanisms of melon peel color formation, we performed an integrated analysis of transcriptome and metabolome with three different fruit peel samples (grey-green ‘W’, dark-green ‘B’, and yellow ‘H’). A total of 40 differentially expressed flavonoids were identified. Integrated transcriptomic and metabolomic analyses revealed that flavonoid biosynthesis was associated with the fruit peel coloration of melon. Twelve differentially expressed genes regulated flavonoids synthesis. Among them, nine (two 4CL, F3H, three F3′H, IFS, FNS, and FLS) up-regulated genes were involved in the accumulation of flavones, flavanones, flavonols, and isoflavones, and three (2 ANS and UFGT) down-regulated genes were involved in the accumulation of anthocyanins. This study laid a foundation to understand the molecular mechanisms of melon peel coloration by exploring valuable genes and metabolites.

Pairwise linkage disequilibrium estimation for polyploids

Many tasks in statistical genetics involve pairwise estimation of linkage disequilibrium (LD). The study of LD in diploids is mature. However, in polyploids, the field lacks a comprehensive characterization of LD. Polyploids also exhibit greater levels of genotype uncertainty than diploids, yet no methods currently exist to estimate LD in polyploids in the presence of such genotype uncertainty. Furthermore, most LD estimation methods do not quantify the level of uncertainty in their LD estimates. Our study contains three major contributions. (i) We characterize haplotypic and composite measures of LD in polyploids. These composite measures of LD turn out to be functions of common statistical measures of association. (ii) We derive procedures to estimate haplotypic and composite LD in polyploids in the presence of genotype uncertainty. We do this by estimating LD directly from genotype likelihoods, which may be obtained from many genotyping platforms. (iii) We derive standard errors of all LD estimators that we discuss. We validate our methods on both real and simulated data. Our methods are implemented in the R package ldsep, available on the Comprehensive R Archive Network https://cran.r-project.org/package=ldsep.

Resequencing of 297 melon accessions reveals the genomic history of improvement and loci related to fruit traits in melon

Domestication and improvement are two important stages in crop evolution. Melon (Cucumis melo L.) is an important vegetable crop with wide phenotypic diversity in many horticultural traits, especially fruit size, flesh thickness and aroma, which are likely the results of long-term extensive selection during its evolution. However, selective signals in domestication and improvement stages for these remarkable variations remain unclear. We resequenced 297 wild, landrace and improved melon accessions and obtained 2 045 412 high-quality SNPs. Population structure and genetic diversity analyses revealed independent and two-step selections in two subspecies of melon: ssp. melo and ssp. agrestis during melon breeding. We detected 233 (~18.35 Mbp) and 159 (~17.71 Mbp) novel potential selective signals during the improvement stage in ssp. agrestis and spp. melo, respectively. Two alcohol acyltransferase genes (CmAATs) unique to the melon genome compared with other cucurbit crops may have undergone stronger selection in ssp. agrestis for the characteristic aroma as compared with other cucurbits. Genome-wide association analysis identified eight fruit size and seven flesh thickness signals overlapping with selective sweeps. Compared with thin-skinned ssp. agrestis, thick-skinned ssp. melo has undergone a stronger selection for thicker flesh. In most melon accessions, CmCLV3 has pleiotropic effects on carpel number and fruit shape. Findings from this study provide novel insights into melon crop evolution, and new tools to advance melon breeding.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo

Melon is an important crop that exhibits broad variation for fruit morphology traits that are the substrate for genetic mapping efforts. In the post-genomic era, the link between genetic maps and physical genome assemblies is key for leveraging QTL mapping results for gene cloning and breeding purposes. Here, using a population of 164 melon recombinant inbred lines (RILs) that were subjected to genotyping-by-sequencing, we constructed and compared high-density sequence- and linkage-based recombination maps that were aligned to the reference melon genome. These analyses reveal the genome-wide variation in recombination frequency and highlight regions of disrupted collinearity between our population and the reference genome. The population was phenotyped over 3 years for fruit size and shape as well as rind netting. Four QTLs were detected for fruit size, and they act in an additive manner, while significant epistatic interaction was found between two neutral loci for this trait. Fruit shape displayed transgressive segregation that was explained by the action of four QTLs, contributed by alleles from both parents. The complexity of rind netting was demonstrated on a collection of 177 diverse accessions. Further dissection of netting in our RILs population, which is derived from a cross of smooth and densely netted parents, confirmed the intricacy of this trait and the involvement of major locus and several other interacting QTLs. A major netting QTL on chromosome 2 co-localized with results from two additional populations, paving the way for future study toward identification of a causative gene for this trait.

Comparative Metabolomics and Molecular Phylogenetics of Melon (Cucumis melo, Cucurbitaceae) Biodiversity

The broad variability of Cucumis melo (melon, Cucurbitaceae) presents a challenge to conventional classification and organization within the species. To shed further light on the infraspecific relationships within C. melo, we compared genotypic and metabolomic similarities among 44 accessions representative of most of the cultivar-groups. Genotyping-by-sequencing (GBS) provided over 20,000 single-nucleotide polymorphisms (SNPs). Metabolomics data of the mature fruit flesh and rind provided over 80,000 metabolomic and elemental features via an orchestra of six complementary metabolomic platforms. These technologies probed polar, semi-polar, and non-polar metabolite fractions as well as a set of mineral elements and included both flavor- and taste-relevant volatile and non-volatile metabolites. Together these results enabled an estimate of "metabolomic/elemental distance" and its correlation with the genetic GBS distance of melon accessions. This study indicates that extensive and non-targeted metabolomics/elemental characterization produced classifications that strongly, but not completely, reflect the current and extensive genetic classification. Certain melon Groups, such as Inodorous, clustered in parallel with the genetic classifications while other genome to metabolome/element associations proved less clear. We suggest that the combined genomic, metabolic, and element data reflect the extensive sexual compatibility among melon accessions and the breeding history that has, for example, targeted metabolic quality traits, such as taste and flavor.

Toward the 'golden' era: The status in uncovering the regulatory control of carotenoid accumulation in plants

Carotenoids are essential pigments to plants and important natural products to humans. Carotenoids as both primary and specialized metabolites fulfill multifaceted functions in plants. As such, carotenoid accumulation (a net process of biosynthesis, degradation and sequestration) is subjected to complicated regulation throughout plant life cycle in response to developmental and environmental signals. Investigation of transcriptional regulation of carotenoid metabolic genes remains the focus in understanding the regulatory control of carotenoid accumulation. While discovery of bona fide carotenoid metabolic regulators is still challenging, the recent progress of identification of various transcription factors and regulators helps us to construct hierarchical regulatory network of carotenoid accumulation. The elucidation of carotenoid regulatory mechanisms at protein level and in chromoplast provides some insights into post-translational regulation of carotenogenic enzymes and carotenoid sequestration in plastid sink. This review briefly describes the pathways and main flux-controlling steps for carotenoid accumulation in plants. It highlights our recent understanding of the regulatory mechanisms underlying carotenoid accumulation at both transcriptional and post-translational levels. It also discusses the opportunities to expand toolbox for further shedding light upon the intrinsic regulation of carotenoid accumulation in plants.

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.

The multi-allelic APRR2 gene is associated with fruit pigment accumulation in melon and watermelon

Color and pigment contents are important aspects of fruit quality and consumer acceptance of cucurbit crops. Here, we describe the independent mapping and cloning of a common causative APRR2 gene regulating pigment accumulation in melon and watermelon. We initially show that the APRR2 transcription factor is causative for the qualitative difference between dark and light green rind in both crops. Further analyses establish the link between sequence or expression level variations in the CmAPRR2 gene and pigment content in the rind and flesh of mature melon fruits. A genome-wide association study (GWAS) of young fruit rind color in a panel composed of 177 diverse melon accessions did not result in any significant association, leading to an earlier assumption that multiple genes are involved in shaping the overall phenotypic variation in this trait. Through resequencing of 25 representative accessions and allelism tests between light rind accessions, we show that multiple independent single nucleotide polymorphisms in the CmAPRR2 gene are causative of the light rind phenotype. The multi-haplotypic nature of this gene explains the lack of detection power obtained through genotyping by sequencing-based GWAS and confirms the pivotal role of this gene in shaping fruit color variation in melon. This study demonstrates the power of combining bi- and multi-allelic designs with deep sequencing, to resolve lack of power due to high haplotypic diversity and low allele frequencies. Due to its central role and broad effect on pigment accumulation in fruits, the APRR2 gene is an attractive target for carotenoid bio-fortification of cucurbit crops.

Re-evaluation of the role of Indian germplasm as center of melon diversification based on genotyping-by-sequencing analysis

BACKGROUND

The importance of Indian germplasm as origin and primary center of diversity of cultivated melon is widely accepted. Genetic diversity of several collections from Indian has been studied previously, although an integrated analysis of these collections in a global diversity perspective has not been possible. In this study, a sample of Indian collections together with a selection of world-wide cultivars to analyze the genetic diversity structure based on Genotype by Sequence data.

RESULTS

A set of 6158 informative Single Nucleotide Polymorphism (SNP) in 175 melon accessions was generated. Melon germplasm could be classified into six major groups, in concordance with horticultural groups. Indian group was in the center of the diversity plot, with the highest genetic diversity. No strict genetic differentiation between wild and cultivated accessions was appreciated in this group. Genomic regions likely involved in the process of diversification were also found. Interestingly, some SNPs differentiating inodorus and cantalupensis groups are linked to Quantitiative Trait Loci involved in ripening behavior (a major characteristic that differentiate those groups). Linkage disequilibrium was found to be low (17 kb), with more rapid decay in euchromatic (8 kb) than heterochromatic (30 kb) regions, demonstrating that recombination events do occur within heterochromatn, although at lower frequency than in euchromatin. Concomitantly, haplotype blocks were relatively small (59 kb). Some of those haplotype blocks were found fixed in different melon groups, being therefore candidate regions that are involved in the diversification of melon cultivars.

CONCLUSIONS

The results support the hypothesis that India is the primary center of diversity of melon, Occidental and Far-East cultivars have been developed by divergent selection. Indian germplasm is genetically distinct from African germplasm, supporting independent domestication events. The current set of traditional Indian accessions may be considered as a population rather than a standard collection of fixed landraces with high intercrossing between cultivated and wild melons.

Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits

Cucurbita pepo contains two cultivated subspecies, each of which encompasses four fruit-shape morphotypes (cultivar groups). The Pumpkin, Vegetable Marrow, Cocozelle, and Zucchini Groups are of subsp. pepo and the Acorn, Crookneck, Scallop, and Straightneck Groups are of subsp. ovifera. Recently, a de novo assembly of the C. pepo subsp. pepo Zucchini genome was published, providing insights into its evolution. To expand our knowledge of evolutionary processes within C. pepo and to identify variants associated with particular morphotypes, we performed whole-genome resequencing of seven of these eight C. pepo morphotypes. We report for the first time whole-genome resequencing of the four subsp. pepo (Pumpkin, Vegetable Marrow, Cocozelle, green Zucchini, and yellow Zucchini) morphotypes and three of the subsp. ovifera (Acorn, Crookneck, and Scallop) morphotypes. A high-depth resequencing approach was followed, using the BGISEQ-500 platform that enables the identification of rare variants, with an average of 33.5X. Approximately 94.5% of the clean reads were mapped against the reference Zucchini genome. In total, 3,823,977 high confidence single-nucleotide polymorphisms (SNPs) were identified. Within each accession, SNPs varied from 636,918 in green Zucchini to 2,656,513 in Crookneck, and were distributed homogeneously along the chromosomes. Clear differences between subspecies pepo and ovifera in genetic variation and linkage disequilibrium are highlighted. In fact, comparison between subspecies pepo and ovifera indicated 5710 genes (22.5%) with Fst > 0.80 and 1059 genes (4.1%) with Fst = 1.00 as potential candidate genes that were fixed during the independent evolution and domestication of the two subspecies. Linkage disequilibrium was greater in subsp. ovifera than in subsp. pepo, perhaps reflective of the earlier differentiation of morphotypes within subsp. ovifera. Some morphotype-specific genes have been localized. Our results offer new clues that may provide an improved understanding of the underlying genomic regions involved in the independent evolution and domestication of the two subspecies. Comparisons among SNPs unique to particular subspecies or morphotypes may provide candidate genes responsible for traits of high economic importance.

The Role of Carotenogenic Metabolic Flux in Carotenoid Accumulation and Chromoplast Differentiation: Lessons From the Melon Fruit

Carotenoids have various roles in plant physiology. Plant carotenoids are synthesized in plastids and are highly abundant in the chromoplasts of ripening fleshy fruits. Considerable research efforts have been devoted to elucidating mechanisms that regulate carotenoid biosynthesis, yet, little is known about the mechanism that triggers storage capacity, mainly through chromoplast differentiation. The Orange gene (OR) product stabilizes phytoene synthase protein (PSY) and triggers chromoplast differentiation. OR underlies carotenoid accumulation in orange cauliflower and melon. The OR's 'golden SNP', found in melon, alters the highly evolutionary conserved Arginine¹⁰⁸ to Histidine and controls β-carotene accumulation in melon fruit, in a mechanism yet to be elucidated. We have recently shown that similar carotenogenic metabolic flux is active in non-orange and orange melon fruit. This flux probably leads to carotenoid turnover but known carotenoid turnover products are not detected in non-orange fruit. Arrest of this metabolic flux, using chemical inhibitors or mutations, induces carotenoid accumulation and biogenesis of chromoplasts, regardless of the allelic state of OR. We suggest that the 'golden SNP' induces β-carotene accumulation probably by negatively affecting the capacity to synthesize downstream compounds. The accumulation of carotenoids induces chromoplast biogenesis through a metabolite-induced mechanism. Carotenogenic turnover flux can occur in non-photosynthetic tissues, which do not accumulate carotenoids. Arrest of this flux by the 'golden SNP' or other flux-arrest mutations is a potential tool for the biofortification of agricultural products with carotenoids.

QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map

BACKGROUND

Melon shows a broad diversity in fruit morphology and quality, which is still underexploited in breeding programs. The knowledge of the genetic basis of fruit quality traits is important for identifying new alleles that may be introduced in elite material by highly efficient molecular breeding tools.

RESULTS

In order to identify QTLs controlling fruit quality, a recombinant inbred line population was developed using two commercial cultivars as parental lines: "Védrantais", from the cantalupensis group, and "Piel de Sapo", from the inodorus group. Both have desirable quality traits for the market, but their fruits differ in traits such as rind and flesh color, sugar content, ripening behavior, size and shape. We used a genotyping-by-sequencing strategy to construct a dense genetic map, which included around five thousand variants distributed in 824 bins. The RIL population was phenotyped for quality and morphology traits, and we mapped 33 stable QTLs involved in sugar and carotenoid content, fruit and seed morphology and major loci controlling external color of immature fruit and mottled rind. The median confidence interval of the QTLs was 942 kb, suggesting that the high density of the genetic map helped in increasing the mapping resolution. Some of these intervals contained less than a hundred annotated genes, and an integrative strategy combining gene expression and resequencing data enabled identification of candidate genes for some of these traits.

CONCLUSION

Several QTLs controlling fruit quality traits in melon were identified and delimited to narrow genomic intervals, using a RIL population and a GBS-based genetic map.

Dissecting the Genetic Architecture of Melon Chilling Tolerance at the Seedling Stage by Association Mapping and Identification of the Elite Alleles

Low temperature is an important abiotic stress that negatively affects morphological growth and fruit development in melon (Cucumis melo L.). Chilling stress at the seedling stage causes seedling injury and poor stand establishment, prolonging vegetative growth and delaying fruit harvest. In this study, association mapping was performed for chilling tolerance at the seedling stage on an expanded melon core collection containing 212 diverse accessions by 272 SSRs and 27 CAPSs. Chilling tolerance of the melon seedlings was evaluated by calculating the chilling injury index (CII) in 2016 and 2017. Genetic diversity analysis of the whole accession panel presented two main groups, which corresponded to the two subspecies of C. melo, melo, and agrestis. Both the subspecies were sensitive to chilling but with agrestis being more tolerant. Genome-wide association study (GWAS) was conducted, respectively, on the whole panel and the two subspecies, totally detecting 51 loci that contributed to 74 marker-trait associations. Of these associations, 35 were detected in the whole panel, 21 in melo, and 18 in agrestis. About half of the associations identified in the two subspecies were also observed in the whole panel, and seven associations were shared by both the subspecies. CMCT505_Chr.1 was repeatedly detected in different populations with high phenotypic contribution and could be a key locus controlling chilling tolerance in C. melo. Nine loci were selected for evaluation of the phenotypic effects related to their alleles, which identified 11 elite alleles contributing to seedling chilling tolerance. Four such alleles existed in both the subspecies and six in either of the two subspecies. Analysis of 20 parental combinations for their allelic status and phenotypic values showed that the elite alleles collectively contributed to enhancement of the chilling tolerance. Tagging the loci responsible for chilling tolerance may simultaneously favor dissecting the complex adaptability traits and elevate the efficiency to improve chilling tolerance using marker-assisted selection in melon.

Linkage Disequilibrium Estimation in Low Coverage High-Throughput Sequencing Data

High-throughput sequencing methods provide a cost-effective approach for genotyping and are commonly used in population genetics studies. A drawback of these methods, however, is that sequencing and genotyping errors can arise...

High-throughput sequencing methods that multiplex a large number of individuals have provided a cost-effective approach for discovering genome-wide genetic variation in large populations. These sequencing methods are increasingly being utilized in population genetic studies across a diverse range of species. Two side-effects of these methods, however, are (1) sequencing errors and (2) heterozygous genotypes called as homozygous due to only one allele at a particular locus being sequenced, which occurs when the sequencing depth is insufficient. Both of these errors have a profound effect on the estimation of linkage disequilibrium (LD) and, if not taken into account, lead to inaccurate estimates. We developed a new likelihood method, GUS-LD, to estimate pairwise linkage disequilibrium using low coverage sequencing data that accounts for undercalled heterozygous genotypes and sequencing errors. Our findings show that accurate estimates were obtained using GUS-LD, whereas underestimation of LD results if no adjustment is made for the errors.