Use of targeted SNP selection for an improved anchoring of the melon (Cucumis melo L.) scaffold genome assembly

Genomic insights into Mediterranean pepper diversity using ddRADSeq

This work investigated the genetic diversity and population structure of 99 pepper lines (Capsicum annuum L.), acclimated to Mediterranean climate conditions, using double-digest restriction site-associated DNA sequencing (ddRADSeq). The aims were to understand the genetic relationships among these lines, correlate genetic clusters with botanical classifications, and provide insights into pepper domestication in the region. Obtained were 318.76 million raw sequence reads overall, averaging 3.21 million reads per sample. A total of 8475 high-quality SNPs were identified and used to assess genetic diversity and population structure. Chromosome NC_061113.1 displayed the highest amount and Chromosome NC_061118.1 the fewest of these SNPs, which were not equally spaced around the genome. Heterozygosity measures and a negative inbreeding coefficient point to the great genetic diversity seen, therefore highlighting the genetic health of the population. Different genetic clusters found by phylogenetic study and STRUCTURE analysis can be used in breeding programs to mix desired features from many genetic backgrounds. This work showed how well ddRADSeq generates high-quality SNPs for genomic research on peppers, therefore offering useful molecular tools for genomic selection and marker-assisted selection. The analysis identified significant genetic diversity and distinct genetic clusters which are valuable for breeding programs focused on crop improvement. These findings enhance our understanding of pepper domestication and provide valuable genetic resources for breeding programs aimed at improving pepper varieties.

Identification of powdery mildew resistance quantitative trait loci in melon and development of resistant near-isogenic lines through marker-assisted backcrossing

BACKGROUND

Melon (Cucumis melo L.), an important cucurbit crop, faces production limitations due to powdery mildew (PM). Developing resistant varieties offers a sustainable, genetics-based alternative to chemical treatments. Therefore, identifying PM resistance quantitative trait loci (QTL) and creating trait-associated markers are essential for efficient melon PM resistance improvement through marker-assisted backcrossing (MABC).

RESULTS

Three F2 populations, A6, B2, and C4, were generated for QTL mapping of PM resistance. Major QTL were identified on chromosome 2 in A6, chromosome 5 in B2, and chromosomes 5 and 12 in C4. A series of TaqMan® assays targeting regions on chromosomes 2, 5, and 12 were developed and validated for foreground and recombinant selection, complemented by the double digest restriction-site associated DNA genotyping system to evaluate the recurrent parent genome recovery. Three MABC programs using resistant donor parents from A6 and C4 crossed with elite susceptible recurrent parents with green and orange fruit flesh were implemented. After two to three cycles of MABC, individual QTL was successfully introgressed into elite genetic backgrounds, giving six PM resistance lines in each green- and orange-fleshed background. PM inoculation on the twelve near-isogenic lines confirmed their resistance to PM.

CONCLUSIONS

We have identified major PM resistance QTL for melon on chromosomes 2, 5, and 12 and have introgressed individual QTL to elite genetic backgrounds using MABC in three and a half years. This study demonstrates the power of combining high-throughput genotyping with breeding efforts and showcases the efficiency of molecular breeding.

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.

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.

Lateral Root Initiation in Cucumber (Cucumis sativus): What Does the Expression Pattern of Rapid Alkalinization Factor 34 (RALF34) Tell Us?

In Arabidopsis, the small signaling peptide (peptide hormone) RALF34 is involved in the gene regulatory network of lateral root initiation. In this study, we aimed to understand the nature of the signals induced by RALF34 in the non-model plant cucumber (Cucumis sativus), where lateral root primordia are induced in the apical meristem of the parental root. The RALF family members of cucumber were identified using phylogenetic analysis. The sequence of events involved in the initiation and development of lateral root primordia in cucumber was examined in detail. To elucidate the role of the small signaling peptide CsRALF34 and its receptor CsTHESEUS1 in the initial stages of lateral root formation in the parental root meristem in cucumber, we studied the expression patterns of both genes, as well as the localization and transport of the CsRALF34 peptide. CsRALF34 is expressed in all plant organs. CsRALF34 seems to differ from AtRALF34 in that its expression is not regulated by auxin. The expression of AtRALF34, as well as CsRALF34, is regulated in part by ethylene. CsTHESEUS1 is expressed constitutively in cucumber root tissues. Our data suggest that CsRALF34 acts in a non-cell-autonomous manner and is not involved in lateral root initiation in cucumber.

Perspectives and recent progress of genome-wide association studies (GWAS) in fruits

BACKGROUND

Earlier next-generation sequencing technologies are being vastly used to explore, administer, and investigate the gene space with accurate profiling of nucleotide variations in the germplasm. OVERVIEW AND PROGRESS: Recently, novel advancements in high-throughput sequencing technologies allow a genotyping-by-sequencing approach that has opened up new horizons for extensive genotyping exploiting single-nucleotide-polymorphisms (SNPs). This method acts as a bridge to support and minimize a genotype to phenotype gap allowing genetic selection at the genome-wide level, named genomic selection that could facilitate the selection of traits also in the pomology sector. In addition to this, genome-wide genotyping is a prerequisite for genome-wide association studies that have been used successfully to discover the genes, which control polygenic traits including the genetic loci, associated with the trait of interest in fruit crops.

AIMS AND PROSPECTS

This review article emphasizes the role of genome-wide approaches to unlock and explore the genetic potential along with the detection of SNPs affecting the phenotype of fruit crops and highlights the prospects of genome-wide association studies in fruits.

Identification of Bacterial Wilt (Erwinia tracheiphila) Resistances in USDA Melon Collection

Bacterial wilt (BW) caused by the Gram-negative bacterium, Erwinia tracheiphila (Et.), is an important disease in melon (Cucumis melo L.). BW-resistant commercial melon varieties are not widely available. There are also no effective pathogen-based disease management strategies as BW-infected plants ultimately die. The purpose of this study is to identify BW-resistant melon accessions in the United States Department of Agriculture (USDA) collection. We tested 118 melon accessions in two inoculation trials under controlled environments. Four-week-old seedlings of test materials were mechanically inoculated with the fluorescently (GFP) labeled or unlabeled E. tracheiphila strain, Hca1-5N. We recorded the number of days to wilting of inoculated leaf (DWIL), days to wilting of whole plant (DWWP) and days to death of the plant (DDP). We identified four melon lines with high resistance to BW inoculation based on all three parameters. Fluorescent microscopy was used to visualize the host colonization dynamics of labeled bacteria from the point of inoculation into petioles, stem and roots in resistant and susceptible melon accessions, which provides an insight into possible mechanisms of BW resistance in melon. The resistant melon lines identified from this study could be valuable resistance sources for breeding of BW resistance as well as the study of cucurbit—E. tracheiphila interactions.

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.

Trimethylamine N-oxide (TMAO): a new attractive target to decrease cardiovascular risk

Cardiovascular disease (CVD) is one of the greatest disease burdens and takes the lives of many each year. There are many risk factors both modifiable and non-modifiable which contribute to the onset and progression of the disease. Trimethylamine N-oxide (TMAO) in recent years has been found to have a correlation with CVD onset. Those with increased levels of the metabolite have a markedly increased risk of future development of cardiometabolic disorders.This literature review aimed to critique past studies undertaken to find a consensus of the significance of the interrelationship between TMAO and cardiovascular risk. A definite link between TMAO levels and a CVD outcome was found. The majority of the literature stated the relationship with evidence; however, there is still some uncertainty as to why and how the correlation occurs. Further study needs to be done to further dissect and understand the relationship between TMAO and CVD risk.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps

BACKGROUND

Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown.

RESULTS

Two high-density genetic maps were constructed by whole-genome resequencing with two F₂ segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci.

CONCLUSION

The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Transcriptomic analysis of a near-isogenic line of melon with high fruit flesh firmness during ripening

BACKGROUND

A near-isogenic line (NIL) of melon (SC10-2) with introgression in linkage group X was studied from harvest (at firm-ripe stage of maturity) until day 18 of postharvest storage at 20.5 °C together with its parental control ('Piel de Sapo', PS).

RESULTS

SC10-2 showed higher flesh firmness and whole fruit hardness but lower juiciness than its parental. SC10-2 showed a decrease in respiration rate accompanied by a decrease in ethylene production during ripening, both of which fell to a greater extent than in PS. The introgression affected 11 volatile organic compounds (VOCs), the levels of which during ripening were generally higher in SC10-2 than in PS. Transcriptomic analysis from RNA-Seq revealed differentially expressed genes (DEGs) associated with the effects studied. For example, 909 DEGs were exclusive to the introgression, and only 23 DEGs were exclusive to postharvest ripening time. Major functions of the DEGs associated with introgression or ripening time were identified by cluster analysis. About 37 genes directly and/or indirectly affected the delay in ripening of SC10-2 compared with PS in general and, more particularly, the physiological and quality traits measured and, probably, the differential non-climacteric response. Of the former genes, we studied in more detail at least five that mapped in the introgression in linkage group (LG) X, and 32 outside it.

CONCLUSION

There is an apparent control of textural changes, VOCs and fruit ripening by an expression quantitative trait locus located in LG X together with a direct control on them due to genes presented in the introgression (CmTrpD, CmNADH1, CmTCP15, CmGDSL esterase/lipase, and CmHK4-like) and CmNAC18. © 2020 Society of Chemical Industry.

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.

Major QTLs for Trunk Height and Correlated Agronomic Traits Provide Insights into Multiple Trait Integration in Oil Palm Breeding

Superior oil yield is always the top priority of the oil palm industry. Short trunk height (THT) and compactness traits have become increasingly important to improve harvesting efficiency since the industry started to suffer yield losses due to labor shortages. Breeding populations with low THT and short frond length (FL) are actually available, such as Dumpy AVROS pisifera (DAV) and Gunung Melayu dura (GM). However, multiple trait stacking still remains a challenge for oil palm breeding, which usually requires 12–20 years to complete a breeding cycle. In this study, yield and height increment in the GM × GM (GM-3341) and the GM × DAV (GM-DAV-3461) crossing programs were evaluated and palms with good yield and smaller height increment were identified. In the GM-3341 family, non-linear THT growth between THT_2008 (seven years old) and THT_2014 (13 years old) was revealed by a moderate correlation, suggesting that inter-palm competition becomes increasingly important. In total, 19 quantitative trait loci (QTLs) for THT_2008 (8), oil per palm (O/P) (7) and FL (4) were localized on the GM-3341 linkage map, with an average mapping interval of 2.01 cM. Three major QTLs for THT_2008, O/P and FL are co-located on chromosome 11 and reflect the correlation of THT_2008 with O/P and FL. Multiple trait selection for high O/P and low THT (based on the cumulative effects of positive alleles per trait) identified one palm from 100 palms, but with a large starting population of 1000–1500 seedling per cross, this low frequency could be easily compensated for during breeding selection.

Long-read bitter gourd (Momordica charantia) genome and the genomic architecture of nonclassic domestication

The genetic architecture of quantitative traits is determined by both Mendelian and polygenic factors, yet classic examples of plant domestication focused on selective sweep of newly mutated Mendelian genes. Here we report the chromosome-level genome assembly and the genomic investigation of a nonclassic domestication example, bitter gourd (Momordica charantia), an important Asian vegetable and medicinal plant of the family Cucurbitaceae. Population resequencing revealed the divergence between wild and South Asian cultivars about 6,000 y ago, followed by the separation of the Southeast Asian cultivars about 800 y ago, with the latter exhibiting more extreme trait divergence from wild progenitors and stronger signs of selection on fruit traits. Unlike some crops where the largest phenotypic changes and traces of selection happened between wild and cultivar groups, in bitter gourd large differences exist between two regional cultivar groups, likely reflecting the distinct consumer preferences in different countries. Despite breeding efforts toward increasing female flower proportion, a gynoecy locus exhibits complex patterns of balanced polymorphism among haplogroups, with potential signs of selective sweep within haplogroups likely reflecting artificial selection and introgression from cultivars back to wild accessions. Our study highlights the importance to investigate such nonclassic example of domestication showing signs of balancing selection and polygenic trait architecture in addition to classic selective sweep in Mendelian factors.

Evolution through segmental duplications and losses: a Super-Reconciliation approach

The classical gene and species tree reconciliation, used to infer the history of gene gain and loss explaining the evolution of gene families, assumes an independent evolution for each family. While this assumption is reasonable for genes that are far apart in the genome, it is not appropriate for genes grouped into syntenic blocks, which are more plausibly the result of a concerted evolution. Here, we introduce the Super-Reconciliation problem which consists in inferring a history of segmental duplication and loss events (involving a set of neighboring genes) leading to a set of present-day syntenies from a single ancestral one. In other words, we extend the traditional Duplication-Loss reconciliation problem of a single gene tree, to a set of trees, accounting for segmental duplications and losses. Existency of a Super-Reconciliation depends on individual gene tree consistency. In addition, ignoring rearrangements implies that existency also depends on gene order consistency. We first show that the problem of reconstructing a most parsimonious Super-Reconciliation, if any, is NP-hard and give an exact exponential-time algorithm to solve it. Alternatively, we show that accounting for rearrangements in the evolutionary model, but still only minimizing segmental duplication and loss events, leads to an exact polynomial-time algorithm. We finally assess time efficiency of the former exponential time algorithm for the Duplication-Loss model on simulated datasets, and give a proof of concept on the opioid receptor genes.

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

Accurate reference genomes have become indispensable tools for characterization of genetic and functional variations. Here we generated a high-quality assembly of the melon Payzawat using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. The final 12 chromosome-level scaffolds cover ∼94.13% of the estimated genome (398.57 Mb). Compared with the published DHL92 genome, our assembly exhibits a 157-fold increase in contig length and remarkable improvements in the assembly of centromeres and telomeres. Six genes within STHQF12.4 on pseudochromosome 12, identified from whole-genome comparison between Payzawat and DHL92, may explain a considerable proportion of the skin thickness. In addition, our population study showed that melon domesticated at multiple times from whole-genome perspective and melons in China are introduced from different routes. Selective sweeps underlying the genes related to desirable traits, haplotypes of alleles associated with agronomic traits, and the variants from resequencing data enable efficient breeding.

•

Provides a high-quality assembly for melon genome

•

Explains a considerable proportion of epidermis thickness

•

Melons in China are introduced from different routes

•

Haplotypes of alleles associated with agronomic traits enable efficient breeding

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.

An Improved Melon Reference Genome With Single-Molecule Sequencing Uncovers a Recent Burst of Transposable Elements With Potential Impact on Genes

The published melon (Cucumis melo L.) reference genome assembly (v3.6.1) has still 41.6 Mb (Megabases) of sequences unassigned to pseudo-chromosomes and about 57 Mb of gaps. Although different approaches have been undertaken to improve the melon genome assembly in recent years, the high percentage of repeats (~40%) and limitations due to read length have made it difficult to resolve gaps and scaffold's misassignments to pseudomolecules, especially in the heterochromatic regions. Taking advantage of the PacBio single- molecule real-time (SMRT) sequencing technology, an improvement of the melon genome was achieved. About 90% of the gaps were filled and the unassigned sequences were drastically reduced. A lift-over of the latest annotation v4.0 allowed to re-collocate protein-coding genes belonging to the unassigned sequences to the pseudomolecules. A direct proof of the improvement reached in the new melon assembly was highlighted looking at the improved annotation of the transposable element fraction. By screening the new assembly, we discovered many young (inserted less than 2Mya), polymorphic LTR-retrotransposons that were not captured in the previous reference genome. These elements sit mostly in the pericentromeric regions, but some of them are inserted in the upstream region of genes suggesting that they can have regulatory potential. This improved reference genome will provide an invaluable tool for identifying new gene or transposon variants associated with important phenotypes.

An improved assembly and annotation of the melon (Cucumis melo L.) reference genome

We report an improved assembly (v3.6.1) of the melon (Cucumis melo L.) genome and a new genome annotation (v4.0). The optical mapping approach allowed correcting the order and the orientation of 21 previous scaffolds and permitted to correctly define the gap-size extension along the 12 pseudomolecules. A new comprehensive annotation was also built in order to update the previous annotation v3.5.1, released more than six years ago. Using an integrative annotation pipeline, based on exhaustive RNA-Seq collections and ad-hoc transposable element annotation, we identified 29,980 protein-coding loci. Compared to the previous version, the v4.0 annotation improved gene models in terms of completeness of gene structure, UTR regions definition, intron-exon junctions and reduction of fragmented genes. More than 8,000 new genes were identified, one third of them being well supported by RNA-Seq data. To make all the new resources easily exploitable and completely available for the scientific community, a redesigned Melonomics genomic platform was released at http://melonomics.net. The resources produced in this work considerably increase the reliability of the melon genome assembly and resolution of the gene models paving the way for further studies in melon and related species.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

Construction of a genome-anchored, high-density genetic map for melon (Cucumis melo L.) and identification of Fusarium oxysporum f. sp. melonis race 1 resistance QTL

Four QTLs and an epistatic interaction were associated with disease severity in response to inoculation with Fusarium oxysporum f. sp. melonis race 1 in a recombinant inbred line population of melon. The USDA Cucumis melo inbred line, MR-1, harbors a wealth of alleles associated with resistance to several major diseases of melon, including powdery mildew, downy mildew, Alternaria leaf blight, and Fusarium wilt. MR-1 was crossed to an Israeli cultivar, Ananas Yok'neam, which is susceptible to all of these diseases, to generate a recombinant inbred line (RIL) population of 172 lines. In this study, the RIL population was genotyped to construct an ultra-dense genetic linkage map with 5663 binned SNPs anchored to the C. melo genome and exhibits the overall high quality of the assembly. The utility of the densely genotyped population was demonstrated through QTL mapping of a well-studied trait, resistance to Fusarium wilt caused by Fusarium oxysporum f. sp. melonis (Fom) race 1. A major QTL co-located with the previously validated resistance gene Fom-2. In addition, three minor QTLs and an epistatic interaction contributing to Fom race 1 resistance were identified. The MR-1 × AY RIL population provides a valuable resource for future QTL mapping studies and marker-assisted selection of disease resistance in melon.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

GLABROUS (CmGL) encodes a HD-ZIP IV transcription factor playing roles in multicellular trichome initiation in melon

Map-based cloning identified CmGL that encodes a HD-ZIP type IV transcription factor that controls multicellular trichome initiation in melon. Trichomes are small hairs covering the aerial parts of plants that originate from the epidermal cells, which can protect plants against the damage by insects and pathogens. The regulatory pathway of unicellular trichomes has been well studied in the model plant Arabidopsis. Little is known about the genetic control and regulation of trichome development in melon (Cucumis melo L.) which has multicellular trichomes. In this study, we identified a melon mutant, cmgl, which showed completely glabrous on all aerial organs. A bulked segregant analysis was conducted to identify polymorphic markers for linkage analysis in a population with 256 F₂ plants, which allowed to locate the cmgl locus in melon chromosome VIII. Next-generation sequencing-aided marker discovery and fine mapping in a large population with 1536 F₂ plants narrowed the candidate gene region to 12 kb that harbored only one candidate gene for cmgl, which encoded a class IV homeodomain-associated leucine zipper transcription factor. Four SNPs in the coding region of the CmGL gene were identified between the two parental lines; a single base substitution from C to A resulted in a premature termination codon and a truncated protein in the cmgl. The SNP was converted into a dCAPS marker, which showed co-segregation in the F₂ population and 564 melon accessions. Result of this study will be helpful for better understanding of genetic control of trichome development in melon and marker-assisted selection in developing new cultivars.

QTL Analyses in Multiple Populations Employed for the Fine Mapping and Identification of Candidate Genes at a Locus Affecting Sugar Accumulation in Melon (Cucumis melo L.)

Sugar content is the major determinant of both fruit quality and consumer acceptance in melon (Cucumis melo L), and is a primary target for crop improvement. Near-isogenic lines (NILs) derived from the intraspecific cross between a "Piel de Sapo" (PS) type and the exotic cultivar "Songwhan Charmi" (SC), and several populations generated from the cross of PS × Ames 24294 ("Trigonus"), a wild melon, were used to identify QTL related to sugar and organic acid composition. Seventy-eight QTL were detected across several locations and different years, with three important clusters related to sugar content located on chromosomes 4, 5, and 7. Two PS × SC NILs (SC5-1 and SC5-2) sharing a common genomic interval of 1.7 Mb at the top of chromosome 5 contained QTL reducing soluble solids content (SSC) and sucrose content by an average of 29 and 68%, respectively. This cluster collocated with QTL affecting sugar content identified in other studies in lines developed from the PS × SC cross and supported the presence of a stable consensus locus involved in sugar accumulation that we named SUCQSC5.1. QTL reducing soluble solids and sucrose content identified in the "Trigonus" mapping populations, as well as QTL identified in previous studies from other ssp. agrestis sources, collocated with SUCQSC5.1, suggesting that they may be allelic and implying a role in domestication. In subNILs derived from the PS × SC5-1 cross, SUCQSC5.1 reduced SSC and sucrose content by an average of 18 and 34%, respectively, and was fine-mapped to a 56.1 kb interval containing four genes. Expression analysis of the candidate genes in mature fruit showed differences between the subNILs with PS alleles that were "high" sugar and SC alleles of "low" sugar phenotypes for MELO3C014519, encoding a putative BEL1-like homeodomain protein. Sequence differences in the gene predicted to affect protein function were restricted to SC and other ssp. agrestis cultivar groups. These results provide the basis for further investigation of genes affecting sugar accumulation in melon.

A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo

The melon (Cucumis melo) genome and genetic maps with hundreds to thousands of single nucleotide polymorphism markers were recently released. However, a high-resolution genetic map was lacking. Gummy stem blight (Gsb) is a destructive disease responsible for considerable economic losses during melon production. We herein describe the development of an ultra-dense genetic map consisting of 12,932 recombination bin markers covering 1,818 cM, with an average distance of 0.17 cM between adjacent tags. A comparison of the genetic maps for melon, watermelon, and cucumber revealed chromosome-level syntenic relationships and recombination events among the three Cucurbitaceae species. Our genetic map was useful for re-anchoring the genome scaffolds of melon. More than 92% assembly was anchored to 12 pseudo-chromosomes and 90% of them were oriented. Furthermore, 1,135 recombination hotspots revealed an unbalanced recombination rate across the melon genome. Genetic analyses of the Gsb-resistant and -susceptible lines indicated the resistance phenotype is mediated by a single dominant gene. We identified Gsb-resistance gene candidates in a 108-kb region on pseudo-chromosome 4. Our findings verify the utility of an ultra-dense genetic map for mapping a gene of interest, and for identifying new disease resistant genes.

A mutation in the melon Vacuolar Protein Sorting 41prevents systemic infection of Cucumber mosaic virus

In the melon exotic accession PI 161375, the gene cmv1, confers recessive resistance to Cucumber mosaic virus (CMV) strains of subgroup II. cmv1 prevents the systemic infection by restricting the virus to the bundle sheath cells and impeding viral loading to the phloem. Here we report the fine mapping and cloning of cmv1. Screening of an F2 population reduced the cmv1 region to a 132 Kb interval that includes a Vacuolar Protein Sorting 41 gene. CmVPS41 is conserved among plants, animals and yeast and is required for post-Golgi vesicle trafficking towards the vacuole. We have validated CmVPS41 as the gene responsible for the resistance, both by generating CMV susceptible transgenic melon plants, expressing the susceptible allele in the resistant cultivar and by characterizing CmVPS41 TILLING mutants with reduced susceptibility to CMV. Finally, a core collection of 52 melon accessions allowed us to identify a single amino acid substitution (L348R) as the only polymorphism associated with the resistant phenotype. CmVPS41 is the first natural recessive resistance gene found to be involved in viral transport and its cellular function suggests that CMV might use CmVPS41 for its own transport towards the phloem.

A high-resolution reference genetic map positioning 8.8 K genes for the conifer white spruce: structural genomics implications and correspondence with physical distance

Over the last decade, extensive genetic and genomic resources have been developed for the conifer white spruce (Picea glauca, Pinaceae), which has one of the largest plant genomes (20 Gbp). Draft genome sequences of white spruce and other conifers have recently been produced, but dense genetic maps are needed to comprehend genome macrostructure, delineate regions involved in quantitative traits, complement functional genomic investigations, and assist the assembly of fragmented genomic sequences. A greatly expanded P. glauca composite linkage map was generated from a set of 1976 full-sib progeny, with the positioning of 8793 expressed genes. Regions with significant low or high gene density were identified. Gene family members tended to be mapped on the same chromosomes, with tandemly arrayed genes significantly biased towards specific functional classes. The map was integrated with transcriptome data surveyed across eight tissues. In total, 69 clusters of co-expressed and co-localising genes were identified. A high level of synteny was found with pine genetic maps, which should facilitate the transfer of structural information in the Pinaceae. Although the current white spruce genome sequence remains highly fragmented, dozens of scaffolds encompassing more than one mapped gene were identified. From these, the relationship between genetic and physical distances was examined and the genome-wide recombination rate was found to be much smaller than most estimates reported for angiosperm genomes. This gene linkage map shall assist the large-scale assembly of the next-generation white spruce genome sequence and provide a reference resource for the conifer genomics community.

In Silico identification and annotation of non-coding RNAs by RNA-seq and De Novo assembly of the transcriptome of Tomato Fruits

The complexity of the tomato (Solanum lycopersicum) transcriptome has not yet been fully elucidated. To gain insights into the diversity and features of coding and non-coding RNA molecules of tomato fruits, we generated strand-specific libraries from berries of two tomato cultivars grown in two open-field conditions with different soil type. Following high-throughput Illumina RNA-sequencing (RNA-seq), more than 90% of the reads (over one billion, derived from twelve dataset) were aligned to the tomato reference genome. We report a comprehensive analysis of the transcriptome, improved with 39,095 transcripts, which reveals previously unannotated novel transcripts, natural antisense transcripts, long non-coding RNAs and alternative splicing variants. In addition, we investigated the sequence variants between the cultivars under investigation to highlight their genetic difference. Our strand-specific analysis allowed us to expand the current tomato transcriptome annotation and it is the first to reveal the complexity of the poly-adenylated RNA world in tomato. Moreover, our work demonstrates the usefulness of strand specific RNA-seq approach for the transcriptome-based genome annotation and provides a resource valuable for further functional studies.

Genome-Wide Single Nucleotide Polymorphism Discovery and the Construction of a High-Density Genetic Map for Melon (Cucumis melo L.) Using Genotyping-by-Sequencing

Although genotyping-by-sequencing (GBS) enables the efficient and low-cost generation of large numbers of markers, the utility of resultant genotypes are limited, because they are enormously error-prone and contain high proportions of missing data. In this study, we generated single nucleotide polymorphism (SNP) markers for 109 recombinant inbred lines of melon (Cucumis melo L.) using the GBS approach and ordered them according to their physical position on the draft double haploid line DHL92 genome. Next, by investigating associations between these SNPs, we discovered that some segments on the physical map conflict with linkage relationships. Therefore, to filter out error-prone loci, 4,110 SNPs in which we have a high degree of confidence were selected as anchors to test independence with respect to unselected markers, and the resultant dataset was then analyzed using the Full-Sib Family Haplotype (FSFHap) algorithm in the software TASSEL 5.2. On the basis of this analysis, 22,933 loci that have an average rate of missing data of 0.281% were used to construct a genetic map, which spans 1,088.3 cM across 12 chromosomes and has a maximum spacing of 6.0 cM. Use of this high-quality linkage map enabled the identification of several quantitative trait loci (QTL) known to control traits in fruit and validated our approach. This study highlights the utility of GBS markers for the identification of trait-associated QTLs in melon and facilitates further investigation of genome structure.

Genotyping-by-sequencing of a melon (Cucumis melo L.) germplasm collection from a secondary center of diversity highlights patterns of genetic variation and genomic features of different gene pools

BACKGROUND

Melon (Cucumis melo L.) is one of the most important horticultural species, which includes several taxonomic groups. With the advent of next-generation sequencing, single nucleotide polymorphism (SNP) markers are widely used in the study of genetic diversity and genomics.

RESULTS

We report the first successful application of genotyping-by-sequencing (GBS) technology in melon. We detected 25,422 SNPs by the analysis of 72 accessions collected in Apulia, a secondary centre of diversity in Southern Italy. Analyses of genetic structure, principal components, and hierarchical clustering support the identification of three distinct subpopulations. One of them includes accessions known with the folk name of 'carosello', referable to the chate taxonomic group. This is one of the oldest domesticated forms of C. melo, once widespread in Europe and now exposed to the risk of genetic erosion. The second subpopulation contains landraces of 'barattiere', a regional vegetable production that was never characterized at the DNA level and we show was erroneously considered another form of chate melon. The third subpopulation includes genotypes of winter melon (C. melo var. inodorus). Genetic analysis within each subpopulation revealed patterns of diversity associated with fruit phenotype and geographical origin. We used SNP data to describe, for each subpopulation, the average linkage disequilibrium (LD) decay, and to highlight genomic regions possibly resulting from directional selection and associated with phenotypic variation.

CONCLUSIONS

We used GBS to characterize patterns of genetic diversity and genomic features within C. melo. We provide useful information to preserve endangered gene pools and to guide the use of germplasm in breeding. Finally, our findings lay a foundation for molecular breeding approaches and the identification of genes underlying phenotypic traits.

A new genomic library of melon introgression lines in a cantaloupe genetic background for dissecting desirable agronomical traits

BACKGROUND

Genomic libraries of introgression lines (ILs) consist of collections of homozygous lines with a single chromosomal introgression from a donor genotype in a common, usually elite, genetic background, representing the whole donor genome in the full collection. Currently, the only available melon IL collection was generated using Piel de sapo (var. inodorus) as the recurrent background. ILs are not available in genetic backgrounds representing other important market class cultivars, such as the cantalupensis. The recent availability of genomic tools in melon, such as SNP collections and genetic maps, facilitates the development of such mapping populations.

RESULTS

We have developed a new genomic library of introgression lines from the Japanese cv. Ginsen Makuwa (var. makuwa) into the French Charentais-type cv. Vedrantais (var. cantalupensis) genetic background. In order to speed up the breeding program, we applied medium-throughput SNP genotyping with Sequenom MassARRAY technology in early backcross generations and High Resolution Melting in the final steps. The phenotyping of the backcross generations and of the final set of 27 ILs (averaging 1.3 introgressions/plant and covering nearly 100 % of the donor genome), in three environments, allowed the detection of stable QTLs for flowering and fruit quality traits, including some that affect fruit size in chromosomes 6 and 11, others that change fruit shape in chromosomes 7 and 11, others that change flesh color in chromosomes 2, 8 and 9, and still others that increase sucrose content and delay climacteric behavior in chromosomes 5 and 10.

CONCLUSIONS

A new melon IL collection in the Charentais genetic background has been developed. Genomic regions that consistently affect flowering and fruit quality traits have been identified, which demonstrates the suitability of this collection for dissecting complex traits in melon. Additionally, pre-breeding lines with new, commercially interesting phenotypes have been observed, including delayed climacteric ripening associated to higher sucrose levels, which is of great interest for Charentais cultivar breeding.

Updated sesame genome assembly and fine mapping of plant height and seed coat color QTLs using a new high-density genetic map

Background

Sesame is an important high-quality oil seed crop. The sesame genome was de novo sequenced and assembled in 2014 (version 1.0); however, the number of anchored pseudomolecules was higher than the chromosome number (2n = 2x = 26) due to the lack of a high-density genetic map with 13 linkage groups.

Results

We resequenced a permanent population consisting of 430 recombinant inbred lines and constructed a genetic map to improve the sesame genome assembly. We successfully anchored 327 scaffolds onto 13 pseudomolecules. The new genome assembly (version 2.0) included 97.5 % of the scaffolds greater than 150 kb in size present in assembly version 1.0 and increased the total pseudomolecule length from 233.7 to 258.4 Mb with 94.3 % of the genome assembled and 97.2 % of the predicted gene models anchored. Based on the new genome assembly, a bin map including 1,522 bins spanning 1090.99 cM was generated and used to identified 41 quantitative trait loci (QTLs) for sesame plant height and 9 for seed coat color. The plant height-related QTLs explained 3–24 % the phenotypic variation (mean value, 8 %), and 29 of them were detected in at least two field trials. Two major loci (qPH-8.2 and qPH-3.3) that contributed 23 and 18 % of the plant height were located in 350 and 928-kb spaces on Chr8 and Chr3, respectively. qPH-3.3, is predicted to be responsible for the semi-dwarf sesame plant phenotype and contains 102 candidate genes. This is the first report of a sesame semi-dwarf locus and provides an interesting opportunity for a plant architecture study of the sesame. For the sesame seed coat color, the QTLs of the color spaces L*, a*, and b* were detected with contribution rates of 3–46 %. qSCb-4.1 contributed approximately 39 % of the b* value and was located on Chr4 in a 199.9-kb space. A list of 32 candidate genes for the locus, including a predicted black seed coat-related gene, was determined by screening the newly anchored genome.

Conclusions

This study offers a high-density genetic map and an improved assembly of the sesame genome. The number of linkage groups and pseudomolecules in this assembly equals the number of sesame chromosomes for the first time. The map and updated genome assembly are expected to serve as a platform for future comparative genomics and genetic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2316-4) contains supplementary material, which is available to authorized users.

Genome-Wide Differentiation of Various Melon Horticultural Groups for Use in GWAS for Fruit Firmness and Construction of a High Resolution Genetic Map

Melon (Cucumis melo L.) is a phenotypically diverse eudicot diploid (2n = 2x = 24) has climacteric and non-climacteric morphotypes and show wide variation for fruit firmness, an important trait for transportation and shelf life. We generated 13,789 SNP markers using genotyping-by-sequencing (GBS) and anchored them to chromosomes to understand genome-wide fixation indices (Fst) between various melon morphotypes and genomewide linkage disequilibrium (LD) decay. The F_ST between accessions of cantalupensis and inodorus was 0.23. The F_ST between cantalupensis and various agrestis accessions was in a range of 0.19-0.53 and between inodorus and agrestis accessions was in a range of 0.21-0.59 indicating sporadic to wide ranging introgression. The EM (Expectation Maximization) algorithm was used for estimation of 1436 haplotypes. Average genome-wide LD decay for the melon genome was noted to be 9.27 Kb. In the current research, we focused on the genome-wide divergence underlying diverse melon horticultural groups. A high-resolution genetic map with 7153 loci was constructed. Genome-wide segregation distortion and recombination rate across various chromosomes were characterized. Melon has climacteric and non-climacteric morphotypes and wide variation for fruit firmness, a very important trait for transportation and shelf life. Various levels of QTLs were identified with high to moderate stringency and linked to fruit firmness using both genome-wide association study (GWAS) and biparental mapping. Gene annotation revealed some of the SNPs are located in β-D-xylosidase, glyoxysomal malate synthase, chloroplastic anthranilate phosphoribosyltransferase, and histidine kinase, the genes that were previously characterized for fruit ripening and softening in other crops.

A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.)

Background

Pumpkin (Cucurbita maxima Duch.) is an economically important crop belonging to the Cucurbitaceae family. However, very few genomic and genetic resources are available for this species. As part of our ongoing efforts to sequence the pumpkin genome, high-density genetic map is essential for anchoring and orienting the assembled scaffolds. In addition, a saturated genetic map can facilitate quantitative trait locus (QTL) mapping.

Results

A set of 186 F₂ plants derived from the cross of pumpkin inbred lines Rimu and SQ026 were genotyped using the genotyping-by-sequencing approach. Using the SNPs we identified, a high-density genetic map containing 458 bin-markers was constructed, spanning a total genetic distance of 2,566.8 cM across the 20 linkage groups of C. maxima with a mean marker density of 5.60 cM. Using this map we were able to anchor 58 assembled scaffolds that covered about 194.5 Mb (71.7 %) of the 271.4 Mb assembled pumpkin genome, of which 44 (183.0 Mb; 67.4 %) were oriented. Furthermore, the high-density genetic map was used to identify genomic regions highly associated with an important agronomic trait, dwarf vine. Three QTLs on linkage groups (LGs) 1, 3 and 4, respectively, were recovered. One QTL, qCmB2, which was located in an interval of 0.42 Mb on LG 3, explained 21.4 % phenotypic variations. Within qCmB2, one gene, Cma_004516, encoding the gibberellin (GA) 20-oxidase in the GA biosynthesis pathway, had a 1249-bp deletion in its promoter in bush type lines, and its expression level was significantly increased during the vine growth and higher in vine type lines than bush type lines, supporting Cma_004516 as a possible candidate gene controlling vine growth in pumpkin.

Conclusions

A high-density pumpkin genetic map was constructed, which was used to successfully anchor and orient the assembled genome scaffolds, and to identify QTLs highly associated with pumpkin vine length. The map provided a valuable resource for gene cloning and marker assisted breeding in pumpkin and other related species. The identified vine length QTLs would help to dissect the underlying molecular basis regulating pumpkin vine growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2312-8) contains supplementary material, which is available to authorized users.

Transposon Insertions, Structural Variations, and SNPs Contribute to the Evolution of the Melon Genome

The availability of extensive databases of crop genome sequences should allow analysis of crop variability at an unprecedented scale, which should have an important impact in plant breeding. However, up to now the analysis of genetic variability at the whole-genome scale has been mainly restricted to single nucleotide polymorphisms (SNPs). This is a strong limitation as structural variation (SV) and transposon insertion polymorphisms are frequent in plant species and have had an important mutational role in crop domestication and breeding. Here, we present the first comprehensive analysis of melon genetic diversity, which includes a detailed analysis of SNPs, SV, and transposon insertion polymorphisms. The variability found among seven melon varieties representing the species diversity and including wild accessions and highly breed lines, is relatively high due in part to the marked divergence of some lineages. The diversity is distributed nonuniformly across the genome, being lower at the extremes of the chromosomes and higher in the pericentromeric regions, which is compatible with the effect of purifying selection and recombination forces over functional regions. Additionally, this variability is greatly reduced among elite varieties, probably due to selection during breeding. We have found some chromosomal regions showing a high differentiation of the elite varieties versus the rest, which could be considered as strongly selected candidate regions. Our data also suggest that transposons and SV may be at the origin of an important fraction of the variability in melon, which highlights the importance of analyzing all types of genetic variability to understand crop genome evolution.