SNP genotyping in melons: genetic variation, population structure, and linkage disequilibrium

Genetic Dissection of ToLCNDV Resistance in Resistant Sources of Cucumis melo

Tomato leaf curl New Delhi virus (ToLCNDV) is a begomovirus causing significant melon (Cucumis melo) crop losses globally. This study aims to map the ToLCNDV resistance in the PI 414723 melon accession, previously identified and characterized through phenotypic studies, thereby exploring shared genomic regions with the established resistant source WM-7. In the present study, WM-7 and PI 414723 were crossed with the susceptible accessions 'Rochet' and 'Blanco' respectively, to generate F1 hybrids. These hybrids were self-pollinated to generate the populations for mapping the ToLCNDV resistance region and designing markers for marker-assisted selection. Disease evaluation included visual symptom scoring, viral-load quantification and tissue printing. Genotyping-by-sequencing and SNP markers were used for quantitative trait loci (QTL) mapping. For genetic analysis, qPCR and bulked segregant RNA-seq (BSR-seq) were performed. Gene expression was assessed using RNA-seq, and qRT-PCR was used for confirmation. The research narrows the candidate region for resistance in WM-7 and identifies overlapping QTLs on chromosome 11 in PI 414723, found in the region of the DNA primase large subunit. BSR-seq and expression analyses highlight potential regulatory roles of chromosome 2 in conferring resistance. Differential expression was confirmed for six genes in the candidate region on chromosome 2. This study confirms the existence of common resistance genes in PI 414723 and WM-7.

Simple, fast and inexpensive hot sodium hydroxide and tris DNA extraction method for genotyping tomato and melon seeds

Seed commerce is a highly profitable global market. Most commercialized seeds are hybrid seeds originating from a controlled cross between two selected parental lines. The market value of hybrid seeds depends on their hybrid genetic purity. DNA molecular markers are a reliable and widespread tool to genotype plant materials; however, DNA extraction from seeds is challenging, often laborious and expensive. With the ultimate goal of creating a tomato and melon hybrid seeds purity test, various challenges arise. To overcome these problems and with the purpose of crude DNA extraction, a simple, fast, inexpensive and easily scalable adaptation of the hot sodium hydroxide and tris method coupled to a competitive allele-specific PCR genotyping method is proposed.

A Novel Introgression Line Library Derived from a Wild Melon Gives Insights into the Genetics of Melon Domestication, Uncovering New Genetic Variability Useful for Breeding

A collection of 30 melon introgression lines (ILs) was developed from the wild accession Ames 24297 (TRI) into 'Piel de Sapo' (PS) genetic background. Each IL carried an average of 1.4 introgressions from TRI, and the introgressions represented 91.4% of the TRI genome. Twenty-two ILs, representing 75% of the TRI genome, were evaluated in greenhouse (Algarrobo and Meliana) and field (Alcàsser) trials, mainly to study traits related to domestication syndrome such as fruit weight (FW) and flesh content (FFP), as well as other fruit quality traits as fruit shape (FS), flesh firmness (FF), soluble solid concentration (SSC), rind color and abscission layer. The IL collection showed an impressive variation in size-related traits, with FW ranging from 800 to 4100 g, reflecting the strong effect of the wild genome on these traits. Most of the ILs produced smaller fruits compared with PS; however, unexpectedly, the IL TRI05-2 produced bigger fruits, likely due to new epistatic interacions with the PS genetic background. In contrast, the genotypic effect for FS was smaller, and few QTLs with notable effects were detected. Interestingly, variability was also observed for FFP, FF and SSC, rind color and abscission layer formation. Genes in these introgressions are candidates for having been involved in melon domestication and diversification as well. These results confirm that the TRI IL collection is a very powerful tool for mapping traits of agronomic interest in melon, allowing the confirmation of previously reported QTLs and the identification of new ones to better understand the domestication process of this crop.

Melon diversity on the Silk Road by molecular phylogenetic analysis in Kazakhstan melons

To uncover population structure, phylogenetic relationship, and diversity in melons along the famous Silk Road, a seed size measurement and a phylogenetic analysis using five chloroplast genome markers, 17 RAPD markers and 11 SSR markers were conducted for 87 Kazakh melon accessions with reference accessions. Kazakh melon accessions had large seed with exception of two accessions of weedy melon, Group Agrestis, and consisted of three cytoplasm types, of which Ib-1/-2 and Ib-3 were dominant in Kazakhstan and nearby areas such as northwestern China, Central Asia and Russia. Molecular phylogeny showed that two unique genetic groups, STIa-2 with Ib-1/-2 cytoplasm and STIa-1 with Ib-3 cytoplasm, and one admixed group, STIAD combined with STIa and STIb, were prevalent across all Kazakh melon groups. STIAD melons that phylogenetically overlapped with STIa-1 and STIa-2 melons were frequent in the eastern Silk Road region, including Kazakhstan. Evidently, a small population contributed to melon development and variation in the eastern Silk Road. Conscious preservation of fruit traits specific to Kazakh melon groups is thought to play a role in the conservation of Kazakh melon genetic variation during melon production, where hybrid progenies were generated through open pollination.

Genetic diversity analysis and variety identification using SSR and SNP markers in melon

Melon is an important horticultural crop with a pleasant aromatic flavor and abundance of health-promoting substances. Numerous melon varieties have been cultivated worldwide in recent years, but the high number of varieties and the high similarity between them poses a major challenge for variety evaluation, discrimination, as well as innovation in breeding. Recently, simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs), two robust molecular markers, have been utilized as a rapid and reliable method for variety identification. To elucidate the genetic structure and diversity of melon varieties, we screened out 136 perfect SSRs and 164 perfect SNPs from the resequencing data of 149 accessions, including the most representative lines worldwide. This study established the DNA fingerprint of 259 widely-cultivated melon varieties in China using Target-seq technology. All melon varieties were classified into five subgruops, including ssp. agrestis, ssp. melo, muskmelon and two subgroups of foreign individuals. Compared with ssp. melo, the ssp. agrestis varieties might be exposed to a high risk of genetic erosion due to their extremely narrow genetic background. Increasing the gene exchange between ssp. melo and ssp. agrestis is therefore necessary in the breeding procedure. In addition, analysis of the DNA fingerprints of the 259 melon varieties showed a good linear correlation (R² = 0.9722) between the SSR genotyping and SNP genotyping methods in variety identification. The pedigree analysis based on the DNA fingerprint of 'Jingyu' and 'Jingmi' series melon varieties was consistent with their breeding history. Based on the SNP index analysis, ssp. agrestis had low gene exchange with ssp. melo in chromosome 4, 7, 10, 11and 12, two specific SNP loci were verified to distinguish ssp. agrestis and ssp. melon varieties. Finally, 23 SSRs and 40 SNPs were selected as the core sets of markers for application in variety identification, which could be efficiently applied to variety authentication, variety monitoring, as well as the protection of intellectual property rights in melon.

Thorough Characterization of ETHQB3.5, a QTL Involved in Melon Fruit Climacteric Behavior and Aroma Volatile Composition

The effect of the QTL involved in climacteric ripening ETHQB3.5 on the fruit VOC composition was studied using a set of Near-Isogenic Lines (NILs) containing overlapping introgressions from the Korean accession PI 16375 on the chromosome 3 in the climacteric 'Piel de Sapo' (PS) genetic background. ETHQB3.5 was mapped in an interval of 1.24 Mb that contained a NAC transcription factor. NIL fruits also showed differences in VOC composition belonging to acetate esters, non-acetate esters, and sulfur-derived families. Cosegregation of VOC composition (23 out of 48 total QTLs were mapped) and climacteric ripening was observed, suggesting a pleiotropic effect of ETHQB3.5. On the other hand, other VOCs (mainly alkanes, aldehydes, and ketones) showed a pattern of variation independent of ETHQB3.5 effects, indicating the presence of other genes controlling non-climacteric ripening VOCs. Network correlation analysis and hierarchical clustering found groups of highly correlated compounds and confirmed the involvement of the climacteric differences in compound classes and VOC differences. The modification of melon VOCs may be achieved with or without interfering with its physiological behavior, but it is likely that high relative concentrations of some type of ethylene-dependent esters could be achieved in climacteric cultivars.

Advanced Genetic Studies on Powdery Mildew Resistance in TGR-1551

Cucurbits powdery mildew (CPM) is one of the main limiting factors of melon cultivation worldwide. Resistance to races 1, 2, and 5 has been reported in the African accession TGR-1551, whose resistance is controlled by a dominant–recessive epistasis. The dominant and recessive quantitative trail loci (QTL) have previously been located in chromosomes 5 and 12, respectively. We used several densely genotyped BC₃ families derived from the cross between TGR-1551 and the susceptible cultivar ‘Bola de Oro’ to finely map these resistance regions. The further phenotyping and genotyping of the selected BC₅, BC₅S₁, BC₅S₂, BC₄S₁, BC₄xPS, and (BC₄xPS) S₁ offspring allowed for the narrowing of the candidate intervals to a 250 and 381 kb region in chromosomes 5 and 12, respectively. Moreover, the temperature effect over the resistance provided by the dominant gene has been confirmed. High resolution melting markers (HRM) were tightly linked to both resistance regions and will be useful in marker-assisted selection programs. Candidate R genes with variants between parents that caused a potential modifier impact on the protein function were identified within both intervals. These candidate genes provide targets for future functional analyses to better understand the resistance to powdery mildew in melons.

Multi-environment genomic prediction for soluble solids content in peach (Prunus persica)

Genotype-by-environment interaction (G × E) is a common phenomenon influencing genetic improvement in plants, and a good understanding of this phenomenon is important for breeding and cultivar deployment strategies. However, there is little information on G × E in horticultural tree crops, mostly due to evaluation costs, leading to a focus on the development and deployment of locally adapted germplasm. Using sweetness (measured as soluble solids content, SSC) in peach/nectarine assessed at four trials from three US peach-breeding programs as a case study, we evaluated the hypotheses that (i) complex data from multiple breeding programs can be connected using GBLUP models to improve the knowledge of G × E for breeding and deployment and (ii) accounting for a known large-effect quantitative trait locus (QTL) improves the prediction accuracy. Following a structured strategy using univariate and multivariate models containing additive and dominance genomic effects on SSC, a model that included a previously detected QTL and background genomic effects was a significantly better fit than a genome-wide model with completely anonymous markers. Estimates of an individual's narrow-sense and broad-sense heritability for SSC were high (0.57-0.73 and 0.66-0.80, respectively), with 19-32% of total genomic variance explained by the QTL. Genome-wide dominance effects and QTL effects were stable across environments. Significant G × E was detected for background genome effects, mostly due to the low correlation of these effects across seasons within a particular trial. The expected prediction accuracy, estimated from the linear model, was higher than the realised prediction accuracy estimated by cross-validation, suggesting that these two parameters measure different qualities of the prediction models. While prediction accuracy was improved in some cases by combining data across trials, particularly when phenotypic data for untested individuals were available from other trials, this improvement was not consistent. This study confirms that complex data can be combined into a single analysis using GBLUP methods to improve understanding of G × E and also incorporate known QTL effects. In addition, the study generated baseline information to account for population structure in genomic prediction models in horticultural crop improvement.

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.

On the genetic architecture in a public tropical maize panel of the symbiosis between corn and plant growth-promoting bacteria aiming to improve plant resilience

Exploring the symbiosis between plants and plant growth-promoting bacteria (PGPB) is a new challenge for sustainable agriculture. Even though many works have reported the beneficial effects of PGPB in increasing plant resilience for several stresses, its potential is not yet widely explored. One of the many reasons is the differential symbiosis performance depending on the host genotype. This opens doors to plant breeding programs to explore the genetic variability and develop new cultivars with higher responses to PGPB interaction and, therefore, have higher resilience to stress. Hence, we aimed to study the genetic architecture of the symbiosis between PGPB and tropical maize germplasm, using a public association panel and its impact on plant resilience. Our findings reveal that the synthetic PGPB population can modulate and impact root architecture traits and improve resilience to nitrogen stress, and 37 regions were significant for controlling the symbiosis between PGPB and tropical maize. In addition, we found two overlapping SNPs in the GWAS analysis indicating strong candidates for further investigations. Furthermore, genomic prediction analysis with genomic relationship matrix computed using only significant SNPs obtained from GWAS analysis substantially increased the predictive ability for several traits endorsing the importance of these genomic regions for the response of PGPB. Finally, the public tropical panel reveals a significant genetic variability to the symbiosis with the PGPB and can be a source of alleles to improve plant resilience.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01257-6.

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.

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.

Resistance to Cucumber Green Mottle Mosaic Virus in Cucumis melo

Cucumber green mottle mosaic virus (CGMMV) is a severe threat to melon production worldwide. At present, there are no cultivars available on the market which show an effective resistance or tolerance to CGMMV infection; only wild Cucumis species were reported as resistant. Germplasm accessions of Cucumis melo, as well as C. anguria, C. ficifolius, C. myriocarpus and C. metuliferus, were mechanically infected with isolates belonging to the European and Asian strain of CGMMV and screened for resistance by scoring symptom severity and comparing the accumulation of virus by qRT-PCR. The wild species C. anguria and C. ficifolius showed no symptoms and did not accumulate CGGMV following inoculation, while C. metuliferus was highly susceptible to the isolates of both strains of CGMMV. The virus accumulated also in C. myriocarpus and the European isolate produced symptoms, but the Asian isolate did not. Thirty C. melo accessions were susceptible to CGMMV. An isolate-dependent expression of symptoms was observed in 16 melon accessions: they showed mild and severe symptoms at 14 and 21 days after inoculation with the European and Asian isolate, respectively. Freeman's Cucumber showed few or no symptoms following inoculation with the isolate of either CGMMV strain. This particular accession also showed reduced virus accumulation, whereas most other tested germplasm accessions showed significantly higher viral loads and, therefore, may well be a candidate for breeding programs aiming to reduce the losses produced by CGMMV with resistant commercial melon cultivars.

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.

Melon Genetic Resources Characterization for Rind Volatile Profile

A melon core collection was analyzed for rind volatile compounds as, despite the fact that they are scarcely studied, these compounds play an important role in consumer preferences. Gas chromatography coupled to mass spectrometry allowed the detection of 171 volatiles. The high volatile diversity found was analyzed by Hierarchical Cluster Analysis (HCA), giving rise to two major clusters of accessions. The first cluster included climacteric and aromatic types such as Cantalupensis, Ameri, Dudaim and Momordica, rich in esters; the second one mainly included non-climacteric non-aromatic types such as Inodorus, Flexuosus, Acidulus, Conomon and wild Agrestis, with low volatiles content, specifically affecting esters. Many interesting accessions were identified, with different combinations of aroma profiles for rind and flesh, such as Spanish Inodorus landraces with low aroma flesh but rind levels of esters similar to those in climacteric Cantalupensis, exotic accessions sharing high contents of specific compounds responsible for the unique aroma of Dudaim melons or wild Agrestis with unexpected high content of some esters. Sesquiterpenes were present in rinds of some Asian Ameri and Momordica landraces, and discriminate groups of cultivars (sesquiterpene-rich/-poor) within each of the two most commercial melon horticultural groups (Cantalupensis and Inodorus), suggesting that the Asian germplasm is in the origin of specific current varieties or that this feature has been introgressed more recently from Asian sources. This rind characterization will encourage future efforts for breeding melon quality as many of the characterized landraces and wild accessions have been underexploited.

A machine learning framework for genotyping the structural variations with copy number variant

Background

Genotyping of structural variation is an important computational problem in next generation sequence data analysis. However, in cancer genomes, the copy number variant(CNV) often coexists with other types of structural variations which significantly reduces the accuracy of the existing genotype methods. The bias on sequencing coverage and variant allelic frequency can be observed on a CNV region, which leads to the genotyping approaches that misinterpret the heterozygote as a homozygote. Furthermore, other data signals such as split mapped read, abnormal read will also be misjudged because of the CNV. Therefore, genotyping the structural variations with CNV is a complicated computational problem which should consider multiple features and their interactions.

Methods

Here we proposed a computational method for genotyping indels in the CNV region, which introduced a machine learning framework to comprehensively incorporate a set of data features and their interactions. We extracted fifteen kinds of classification features as input and different from the traditional genotyping problem, here the structure of variant may fall into types of normal homozygote, homozygous variant, heterozygous variant without CNV, heterozygous variant with a CNV on the mutated haplotype, and heterozygous variant with a CNV on the wild haplotype. The Multiclass Relevance Vector Machine (M-RVM) was used as a machine learning framework combined with the distribution characteristics of the features.

Results

We applied the proposed method to both simulated and real data, and compared it with the existing popular softwares include Gindel, Facets, GATK, and also compared with other machine learning cores: Support Vector Machine, Lanrange-SVM with OVO multiple classification, Naïve Bayes and BP Neural Network. The results demonstrated that the proposed method outperforms others on accuracy, stability and efficiency.

Conclusion

This work shows that the genotyping of structural variations on the CNV region cannot be solved as a traditional genotyping problem. More features should be used to efficiently complete the five-category task. According to the result, the proposed method can be a practical algorithm to correct genotype structural variations with CNV on the next generation sequence data. The source codes have been uploaded at https://github.com/TrinaZ/Mixgenotypefor academic usage only.

Melon Genome Regions Associated with TGR-1551-Derived Resistance to Cucurbit yellow stunting disorder virus

Cucurbit yellow stunting disorder virus (CYSDV) is one of the main limiting factors of melon cultivation worldwide. To date, no commercial melon cultivars resistant to CYSDV are available. The African accession TGR-1551 is resistant to CYSDV. Two major quantitative trait loci (QTLs) have been previously reported, both located near each other in chromosome 5. With the objective of further mapping the gene or genes responsible of the resistance, a recombinant inbred line (RIL) population derived from the cross between TGR-1551 and the susceptible cultivar ‘Bola de Oro’ was evaluated for resistance to CYSDV in five different assays and genotyped in a genotyping by sequencing (GBS) analysis. The major effect of one of the two QTLs located on chromosome 5 was confirmed in the multienvironment RIL assay and additionally verified through the analysis of three segregating BC₁S₁ populations derived from three resistant RILs. Furthermore, progeny test using the offspring of selected BC₃ plants allowed the narrowing of the candidate interval to a 700 kb region. The SNP markers identified in this work will be useful in marker-assisted selection in the context of introgression of CYSDV resistance in elite cultivars.

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.

Microsatellite Markers Reveal Genetic Diversity and Relationships within a Melon Collection Mainly Comprising Asian Cultivated and Wild Germplasms

Melon, Cucumis melo L., is an important horticultural crop with abundant morphological variability, but the genetic diversity and relationships within wild and cultivated melons remain unclear to date. In this study, thick-skinned (TC) (cultivated subspecies melo), thin-skinned (TN) (cultivated subspecies agrestis), and wild accessions were analyzed for genetic diversity and relationships using 36 microsatellite markers. A total of 314 alleles were detected with a mean allelic number of 8.72 and polymorphism information content of 0.67. Cluster analysis of the accessions resulted in four distinct clusters (I, II, III, and IV) broadly matching with the TC, TN, and wild groups. Cluster I contained only two Indian wild accessions. Cluster II was consisted of 49 South Asian accessions, 34 wild accessions, and 15 TN accessions. Cluster III was a typical TC group including 51 multiorigin TC accessions and one wild accession. The remaining 88 accessions, including 75 TN accessions, 6 wild accessions, and 7 TC accessions, formed the cluster IV, and all the TN and wild accessions in this cluster were from China. These findings were also confirmed by Principal component analysis and STRUCTURE analysis. The South Asian subspecies agrestis accessions, wild and cultivated, had close genetic relationships with a distinctive genetic background. Chinese wild melons showed closeness to cultivated subspecies agrestis landraces and could be a return from the indigenous cultivated melons. The AMOVA and pairwise F statistics (F_ST) presented genetic differentiation among the three groups, with the strongest differentiation (F_ST = 0.380) between TC and TN melons. These results offer overall information on genetic diversity and affiliations within a variety of melon germplasms and favor efficient organization and utilization of these resources for the current breeding purpose.

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.

CmVPS41 Is a General Gatekeeper for Resistance to Cucumber Mosaic Virus Phloem Entry in Melon

Melon production is often compromised by viral diseases, which cannot be treated with chemicals. Therefore, the use of genetic resistances is the main strategy for generating crops resistant to viruses. Resistance to Cucumber mosaic virus (CMV) in melon is scarcely described in few accessions. Until recently, the only known resistant accessions were Freeman's Cucumber and PI 161375, cultivar Songwhan Charmi (SC). Resistance to CMV in melon is recessive and generally oligogenic and quantitative. However, in SC, the resistance to CMV strains of subgroup II is monogenic, depending only on one gene, cmv1, which is able to stop CMV movement by restricting the virus to the bundle sheath cells and preventing a systemic infection. This restriction depends on the viral movement protein (MP). Chimeric viruses carrying the MP of subgroup II strains, like the strain LS (CMV-LS), are restricted in the bundle sheath cells, whereas those carrying MP from subgroup I, like the strain FNY (CMV-FNY), are able to overcome this restriction. cmv1 encodes a vacuolar protein sorting 41 (CmVPS41), a protein involved in the transport of cargo proteins from the Golgi to the vacuole through late endosomes. We have analyzed the variability of the gene CmVPS41 in a set of 52 melon accessions belonging to 15 melon groups, both from the spp melo and the spp agrestis. We have identified 16 different haplotypes, encoding 12 different CmVPS41 protein variants. Challenging members of all haplotypes with CMV-LS, we have identified nine new resistant accessions. The resistance correlates with the presence of two mutations, either L348R, previously found in the accession SC and present in other three melon genotypes, or G85E, present in Freeman's Cucumber and found also in four additional melon genotypes. Moreover, the new resistant accessions belong to three different melon horticultural groups, Conomon, Makuwa, and Dudaim. In the new resistant accessions, the virus was able to replicate and move cell to cell, but was not able to reach the phloem. Therefore, resistance to phloem entry seems to be a general strategy in melon controlled by CmVPS41. Finally, the newly reported resistant accessions broaden the possibilities for the use of genetic resistances in new melon breeding strategies.

A high density SLAF-SNP genetic map and QTL detection for fibre quality traits in Gossypium hirsutum

Background

Upland Cotton (Gossypium hirsutum) is a very important cash crop known for its high quality natural fiber. Recent advances in sequencing technologies provide powerful tools with which to explore the cotton genome for single nucleotide polymorphism marker identification and high density genetic map construction toward more reliable quantitative trait locus mapping.

Results

In the present study, a RIL population was developed by crossing a Chinese high fiber quality cultivar (Yumian 1) and an American high fiber quality line (CA3084), with distinct genetic backgrounds. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to discover SNPs, and a genetic map containing 6254 SNPs was constructed, covering 3141.72 cM with an average distance of 0.5 cM between markers. A total of 95 QTL were detected for fiber quality traits in three environments, explaining 5.5-24.6% of the phenotypic variance. Fifty-five QTL found in multiple environments were considered stable QTL. Nine of the stable QTL were found in all three environments. We identified 14 QTL clusters on 13 chromosomes, each containing one or more stable QTL.

Conclusion

A high-density genetic map of Gossypium hirsutum developed by using specific locus amplified fragment sequencing technology provides detailed mapping of fiber quality QTL, and identification of ‘stable QTL’ found in multiple environments. A marker-rich genetic map provides a foundation for fine mapping, candidate gene identification and marker-assisted selection of favorable alleles at stable QTL in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5294-5) contains supplementary material, which is available to authorized users.

Characterizing Palestinian snake melon (Cucumis melo var. flexuosus) germplasm diversity and structure using SNP and DArTseq markers

BACKGROUND

Crop landraces embody a source of beneficial genes potentially providing endurance to environmental stress and other agronomic qualities including yield. Our study included 88 snake melon accessions (Cucumis melo var. flexuosus) collected from 9 districts in the Palestinian West-Bank. These accessions represent four landraces of Palestinian snake melon: Green, and White Baladi, and Green, and White Sahouri.

RESULTS

This is the first report on successful application of genotyping by sequencing in snake melon. Nine thousand seven hundred fifty single-nucleotide polymorphism (SNP) and 7400 DArTseq genetic markers were employed to evaluate genetic biodiversity and population structure of Palestinian snake melon germplasm collection. Clustering based on neighbor-joining-analysis, principle coordinate and Bayesian model implemented in Structure showed that patterns of genetic diversity of snake melon landraces depends on their geographical source and unraveled the presence of two major local landraces (Sahouri, and Baladi) with accessions from each group clustering together. A significant correlation was observed between both types of markers in Mantel correlation test. A significant association between genetic and geographic matrices (P < 0.0001) was also detected. AMOVA indicated that majority of variation (90%) was due to the difference within accessions.

CONCLUSION

The Palestinian landraces seem to have unique genes that may allow the enhancement of the global snake melon gene pool and developments of the plant production worldwide. Our subsequent objective is to detect genotypes with promising qualities and to conduct association mapping studies concentrating on Fusarium-wilt resistance, yield, and environmental stresses.

Repeated domestication of melon (Cucumis melo) in Africa and Asia and a new close relative from India

PREMISE OF THE STUDY

The domestication history of melon is still unclear. An African or Asian origin has been suggested, but its closest wild relative was recently revealed to be an Australian species. The complicated taxonomic history of melon has resulted in additional confusion, with a high number of misidentified germplasm collections currently used by breeders and in genomics research.

METHODS

Using seven DNA regions sequenced for 90% of the genus and the major cultivar groups, we sort out described names and infer evolutionary origins and domestication centers.

KEY RESULTS

We found that modern melon cultivars go back to two lineages, which diverged ca. 2 million years ago. One is restricted to Asia (Cucumis melo subsp. melo), and the second, here described as C. melo subsp. meloides, is restricted to Africa. The Asian lineage has given rise to the widely commercialized cultivar groups and their market types, while the African lineage gave rise to cultivars still grown in the Sudanian region. We show that C. trigonus, an overlooked perennial and drought-tolerant species from India is among the closest living relatives of C. melo.

CONCLUSIONS

Melon was domesticated at least twice: in Africa and Asia. The African lineage and the Indian C. trigonus are exciting new resources for breeding of melons tolerant to climate change.

Fruit flesh volatile and carotenoid profile analysis within the Cucumis melo L. species reveals unexploited variability for future genetic breeding

BACKGROUND

Aroma profile and carotenoids content of melon flesh are two important aspects influencing the quality of this fruit that have been characterized using only selected genotypes. However, the extant variability of the whole species remains unknown.

RESULTS

A complete view of the volatile/carotenoid profiles of melon flesh was obtained analyzing 71 accessions, representing the whole diversity of the species. Gas chromatography-mass spectrometry and high-performance liquid chromatography were used to analyze 200 volatile compounds and five carotenoids. Genotypes were classified into two main clusters (high/low aroma), but with a large diversity of differential profiles within each cluster, consistent with the ripening behavior, flesh color and proposed evolutionary and breeding history of the different horticultural groups.

CONCLUSION

Our results highlight the huge amount of untapped aroma diversity of melon germplasm, especially of non-commercial types. Also, landraces with high nutritional value with regard to carotenoids have been identified. All this knowledge will encourage melon breeding, facilitating the selection of the genetic resources more appropriate to develop cultivars with new aromatic profiles or to minimize the impact of breeding on melon quality. The newly characterized sources provide the basis for further investigations into specific genes/alleles contributing to melon flesh quality. © 2018 Society of Chemical Industry.

Resistance to tomato leaf curl New Delhi virus in melon is controlled by a major QTL located in chromosome 11

Identification of three genomic regions and underlying candidate genes controlling the high level of resistance to ToLCNDV derived from a wild melon. SNP markers appropriated for MAS management of resistance. Tomato leaf curl New Delhi virus (ToLCNDV) is a bipartite begomovirus that severely affects melon crop (Cucumis melo) in the main production areas of Spain since 2012. In this work, we evaluated the degree of resistance of four accessions (two belonging to the subsp. agrestis var. momordica and two to the wild agrestis group) and their corresponding hybrids with a susceptible commercial melon belonging to the subsp. melo (Piel de Sapo, PS). The analysis using quantitative PCR (qPCR) allowed us to select one wild agrestis genotype (WM-7) with a high level of resistance and use it to construct segregating populations (F ₂ and backcrosses). These populations were phenotyped for symptom severity and virus content using qPCR, and genotyped with different sets of SNP markers. Phenotyping and genotyping results in the F ₂ and BC1s populations derived from the WM-7 × PS cross were used for QTL analysis. Three genomic regions controlling resistance to ToLCNDV were found, one major locus in chromosome 11 and two additional regions in chromosomes 12 and 2. The highest level of resistance (no or mild symptoms and very low viral titer) was obtained with the homozygous WM-7WM-7 genotype at the major QTL in chromosome 11, even with PSPS genotypes at the other two loci. The resistance derived from WM-7 is useful to develop new melon cultivars and the linked SNPs selected in this paper will be highly useful in marker-assisted breeding for ToLCNDV resistance in melon.

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 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.

Quantitative trait loci analysis of melon (Cucumis melo L.) domestication-related traits

Loci on LGIV, VI, and VIII of melon genome are involved in the control of fruit domestication-related traits and they are candidate to have played a role in the domestication of the crop. The fruit of wild melons is very small (20-50 g) without edible pulp, contrasting with the large size and high pulp content of cultivated melon fruits. An analysis of quantitative trait loci (QTL) controlling fruit morphology domestication-related traits was carried out using an in vitro maintained F₂ population from the cross between the Indian wild melon "Trigonus" and the western elite cultivar 'Piel de Sapo'. Twenty-seven QTL were identified in at least two out of the three field trials. Six of them were also being detected in BC1 and BC3 populations derived from the same cross. Ten of them were related to fruit morphological traits, 12 to fruit size characters, and 5 to pulp content. The Trigonus alleles decreased the value of the characters, except for the QTL at andromonoecious gene at linkage group (LG) II, and the QTL for pulp content at LGV. QTL genotypes accounted for a considerable degree of the total phenotypic variation, reaching up to 46%. Around 66% of the QTL showed additive gene action, 19% exhibited dominance, and 25% consisted of overdominance. The regions on LGIV, VI, and VIII included the QTL with more consistent and strong effects on domestication-related traits. QTLs on those regions were validated in BC2S1, BC2S2, and BC3 families, with "Trigonus" allele decreasing the fruit morphological traits in all cases. The validated QTL could represent loci involved in melon domestication, although further experiments as genomic variation studies across wild and cultivated genotypes would be necessary to confirm this hypothesis.

ETHQV6.3 is involved in melon climacteric fruit ripening and is encoded by a NAC domain transcription factor

Fruit ripening is divided into climacteric and non-climacteric types depending on the presence or absence of a transient rise in respiration rate and the production of autocatalytic ethylene. Melon is ideal for the study of fruit ripening, as both climacteric and non-climacteric varieties exist. Two introgressions of the non-climacteric accession PI 161375, encompassed in the QTLs ETHQB3.5 and ETHQV6.3, into the non-climacteric 'Piel de Sapo' background are able to induce climacteric ripening independently. We report that the gene underlying ETHQV6.3 is MELO3C016540 (CmNAC-NOR), encoding a NAC (NAM, ATAF1,2, CUC2) transcription factor that is closely related to the tomato NOR (non-ripening) gene. CmNAC-NOR was functionally validated through the identification of two TILLING lines carrying non-synonymous mutations in the conserved NAC domain region. In an otherwise highly climacteric genetic background, both mutations provoked a significant delay in the onset of fruit ripening and in the biosynthesis of ethylene. The PI 161375 allele of ETHQV6.3 is similar to that of climacteric lines of the cantalupensis type and, when introgressed into the non-climacteric 'Piel de Sapo', partially restores its climacteric ripening capacity. CmNAC-NOR is expressed in fruit flesh of both climacteric and non-climacteric lines, suggesting that the causal mutation may not be acting at the transcriptional level. The use of a comparative genetic approach in a species with both climacteric and non-climacteric ripening is a powerful strategy to dissect the complex mechanisms regulating the onset of fruit ripening.

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.

Evaluation of LD decay and various LD-decay estimators in simulated and SNP-array data of tetraploid potato

The number of SNPs required for QTL discovery is justified by the distance at which linkage disequilibrium has decayed. Simulations and real potato SNP data showed how to estimate and interpret LD decay. The magnitude of linkage disequilibrium (LD) and its decay with genetic distance determine the resolution of association mapping, and are useful for assessing the desired numbers of SNPs on arrays. To study LD and LD decay in tetraploid potato, we simulated autotetraploid genotypes and used it to explore the dependence on: (1) the number of haplotypes in the population (the amount of genetic variation) and (2) the percentage of haplotype specific SNPs (hs-SNPs). Several estimators for short-range LD were explored, such as the average r ², median r ², and other percentiles of r ² (80, 90, and 95 %). For LD decay, we looked at LD_½,90, the distance at which the short-range LD is halved when using the 90 % percentile of r ² at short range, as estimator for LD. Simulations showed that the performance of various estimators for LD decay strongly depended on the number of haplotypes, although the real value of LD decay was not influenced very much by this number. The estimator LD_½,90 was chosen to evaluate LD decay in 537 tetraploid varieties. LD_½,90 values were 1.5 Mb for varieties released before 1945 and 0.6 Mb in varieties released after 2005. LD_½,90 values within three different subpopulations ranged from 0.7 to 0.9 Mb. LD_½,90 was 2.5 Mb for introgressed regions, indicating large haplotype blocks. In pericentromeric heterochromatin, LD decay was negligible. This study demonstrates that several related factors influencing LD decay could be disentangled, that no universal approach can be suggested, and that the estimation of LD decay has to be performed with great care and knowledge of the sampled material.

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.

Whole Genome Re-Sequencing and Characterization of Powdery Mildew Disease-Associated Allelic Variation in Melon

Powdery mildew is one of the most common fungal diseases in the world. This disease frequently affects melon (Cucumis melo L.) and other Cucurbitaceous family crops in both open field and greenhouse cultivation. One of the goals of genomics is to identify the polymorphic loci responsible for variation in phenotypic traits. In this study, powdery mildew disease assessment scores were calculated for four melon accessions, 'SCNU1154', 'Edisto47', 'MR-1', and 'PMR5'. To investigate the genetic variation of these accessions, whole genome re-sequencing using the Illumina HiSeq 2000 platform was performed. A total of 754,759,704 quality-filtered reads were generated, with an average of 82.64% coverage relative to the reference genome. Comparisons of the sequences for the melon accessions revealed around 7.4 million single nucleotide polymorphisms (SNPs), 1.9 million InDels, and 182,398 putative structural variations (SVs). Functional enrichment analysis of detected variations classified them into biological process, cellular component and molecular function categories. Further, a disease-associated QTL map was constructed for 390 SNPs and 45 InDels identified as related to defense-response genes. Among them 112 SNPs and 12 InDels were observed in powdery mildew responsive chromosomes. Accordingly, this whole genome re-sequencing study identified SNPs and InDels associated with defense genes that will serve as candidate polymorphisms in the search for sources of resistance against powdery mildew disease and could accelerate marker-assisted breeding in melon.

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.

Combined use of genetic and genomics resources to understand virus resistance and fruit quality traits in melon

The availability of the genome sequence of many crop species during the past few years has opened a new era in plant biology, allowing for the performance of massive genomic studies in plant species other than the classical models Arabidopsis and rice. One of these crop species is melon (Cucumis melo), a cucurbit of high economic value that has become an interesting model for the study of biological processes such as fruit ripening, sex determination and phloem transport. The recent availability of the melon genome sequence, together with a number of genetic and genomic resources, provides powerful tools that can be used to assist in the main melon breeding targets, namely disease resistance and fruit quality. In this review, we will describe recent data obtained combining the use of a melon near isogenic line (NIL) population and genomic resources to gain insight into agronomically important traits as fruit ripening, resistance to Cucumber Mosaic virus (CMV) and the accumulation of sugars in fruits.

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.

Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

Background

SNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.

Results

We construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.

Conclusions

We report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

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

Variability of candidate genes, genetic structure and association with sugar accumulation and climacteric behavior in a broad germplasm collection of melon (Cucumis melo L.)

Background

A collection of 175 melon (Cucumis melo L.) accessions (including wild relatives, feral types, landraces, breeding lines and commercial cultivars) from 50 countries was selected to study the phenotypic variability for ripening behavior and sugar accumulation. The variability of single nucleotide polymorphisms (SNPs) at 53 selected candidate genes involved in sugar accumulation and fruit ripening processes was studied, as well as their association with phenotypic variation of related traits.

Results

The collection showed a strong genetic structure, defining seven groups plus a number of accessions that could not be associated to any of the groups (admixture), which fitted well with the botanical classification of melon varieties. The variability in candidate genes for ethylene, cell wall and sugar-related traits was high and similar to SNPs located in reference genes. Variability at ripening candidate genes had an important weight on the genetic stratification of melon germplasm, indicating that traditional farmers might have selected for ripening traits during cultivar diversification. A strong relationship was also found between the genetic structure and phenotypic diversity, which could hamper genetic association studies. Accessions belonging to the ameri group are the most appropriate for association analysis given the high phenotypic and molecular diversity within the group, and lack of genetic structure.

The most remarkable association was found between sugar content and SNPs in LG III, where a hotspot of sugar content QTLs has previously been defined. By studying the differences in allelic variation of SNPs within horticultural groups with specific phenotypic features, we also detected differential variation in sugar-related candidates located in LGIX and LGX, and in ripening-related candidates located in LGII and X, all in regions with previously mapped QTLs for the corresponding traits.

Conclusions

In the current study we have found an important variability at both the phenotypic and candidate gene levels for ripening behavior and sugar accumulation in melon fruit. By combination of differences in allelic diversity and association analysis, we have identified several candidate genes that may be involved in the melon phenotypic diversity.

Electronic supplementary material

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

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

Background

The genome of the melon (Cucumis melo L.) double-haploid line DHL92 was recently sequenced, with 87.5 and 80.8% of the scaffold assembly anchored and oriented to the 12 linkage groups, respectively. However, insufficient marker coverage and a lack of recombination left several large, gene rich scaffolds unanchored, and some anchored scaffolds unoriented. To improve the anchoring and orientation of the melon genome assembly, we used resequencing data between the parental lines of DHL92 to develop a new set of SNP markers from unanchored scaffolds.

Results

A high-resolution genetic map composed of 580 SNPs was used to anchor 354.8 Mb of sequence, contained in 141 scaffolds (average size 2.5 Mb) and corresponding to 98.2% of the scaffold assembly, to the 12 melon chromosomes. Over 325.4 Mb (90%) of the assembly was oriented. The genetic map revealed regions of segregation distortion favoring SC alleles as well as recombination suppression regions coinciding with putative centromere, 45S, and 5S rDNA sites. New chromosome-scale pseudomolecules were created by incorporating to the previous v3.5 version an additional 38.3 Mb of anchored sequence representing 1,837 predicted genes contained in 55 scaffolds. Using fluorescent in situ hybridization (FISH) with BACs that produced chromosome-specific signals, melon chromosomes that correspond to the twelve linkage groups were identified, and a standardized karyotype of melon inbred line T111 was developed.

Conclusions

By utilizing resequencing data and targeted SNP selection combined with a large F2 mapping population, we significantly improved the quantity of anchored and oriented melon scaffold genome assembly. Using genome information combined with FISH mapping provided the first cytogenetic map of an inodorus melon type. With these results it was possible to make inferences on melon chromosome structure by relating zones of recombination suppression to centromeres and 45S and 5S heterochromatic regions. This study represents the first steps towards the integration of the high-resolution genetic and cytogenetic maps with the genomic sequence in melon that will provide more information on genome organization and allow for the improvement of the melon genome draft sequence.

Electronic supplementary material

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

Analysis of wild-species introgressions in tomato inbreds uncovers ancestral origins

BACKGROUND

Decades of intensive tomato breeding using wild-species germplasm have resulted in the genomes of domesticated germplasm (Solanum lycopersicum) being intertwined with introgressions from their wild relatives. Comparative analysis of genomes among cultivated tomatoes and wild species that have contributed genetic variation can help identify desirable genes, such as those conferring disease resistance. The ability to identify introgression position, borders, and contents can reveal ancestral origins and facilitate harnessing of wild variation in crop breeding.

RESULTS

Here we present the whole-genome sequences of two tomato inbreds, Gh13 and BTI-87, both carrying the begomovirus resistance locus Ty-3 introgressed from wild tomato species. Introgressions of different sizes on chromosome 6 of Gh13 and BTI-87, both corresponding to the Ty-3 region, were identified as from a source close to the wild species S. chilense. Other introgressions were identified throughout the genomes of the inbreds and showed major differences in the breeding pedigrees of the two lines. Interestingly, additional large introgressions from the close tomato relative S. pimpinellifolium were identified in both lines. Some of the polymorphic regions were attributed to introgressions in the reference Heinz 1706 genome, indicating wild genome sequences in the reference tomato genome.

CONCLUSIONS

The methods developed in this work can be used to delineate genome introgressions, and subsequently contribute to development of molecular markers to aid phenotypic selection, fine mapping and discovery of candidate genes for important phenotypes, and for identification of novel variation for tomato improvement. These universal methods can easily be applied to other crop plants.

Mapping and introgression of QTL involved in fruit shape transgressive segregation into ‘piel de sapo’ melon (cucumis melo l.) [corrected]

A mapping F₂ population from the cross ‘Piel de Sapo’ × PI124112 was selectively genotyped to study the genetic control of morphological fruit traits by QTL (Quantitative Trait Loci) analysis. Ten QTL were identified, five for FL (Fruit Length), two for FD (Fruit Diameter) and three for FS (Fruit Shape). At least one robust QTL per character was found, flqs8.1 (LOD = 16.85, R² = 34%), fdqs12.1 (LOD = 3.47, R² = 11%) and fsqs8.1 (LOD = 14.85, R² = 41%). flqs2.1 and fsqs2.1 cosegregate with gene a (andromonoecious), responsible for flower sex determination and with pleiotropic effects on FS. They display a positive additive effect (a) value, so the PI124112 allele causes an increase in FL and FS, producing more elongated fruits. Conversely, the negative a value for flqs8.1 and fsqs8.1 indicates a decrease in FL and FS, what results in rounder fruits, even if PI124112 produces very elongated melons. This is explained by a significant epistatic interaction between fsqs2.1 and fsqs8.1, where the effects of the alleles at locus a are attenuated by the additive PI124112 allele at fsqs8.1. Roundest fruits are produced by homozygous for PI124112 at fsqs8.1 that do not carry any dominant A allele at locus a (PiPiaa). A significant interaction between fsqs8.1 and fsqs12.1 was also detected, with the alleles at fsqs12.1 producing more elongated fruits. fsqs8.1 seems to be allelic to QTL discovered in other populations where the exotic alleles produce elongated fruits. This model has been validated in assays with backcross lines along 3 years and ultimately obtaining a fsqs8.1-NIL (Near Isogenic Line) in ‘Piel de Sapo’ background which yields round melons.

The genetic basis of fruit morphology in horticultural crops: lessons from tomato and melon

Fruits represent an important part of the human diet and show extensive variation in size and shape between and within cultivated species. The genetic basis of such variation has been studied most extensively in tomato, where currently six quantitative trait loci (QTLs) involving these traits have been fine-mapped and the genes underlying the QTLs identified. The genes responsible for the cloned QTLs belong to families with a few to many members. FASCIATED is encoded by a member of the YABBY family, CNR/FW2.2 by a member of the Cell Number Regulator family, SlKLUH/FW3.2 by a cytochrome P450 of the 78A class (CYP78A), LOCULE NUMBER by a member of the WOX family including WUSCHEL, OVATE by a member of the Ovate Family Proteins (OFP), and SUN by a member of the IQ domain family. A high portion of the history and current diversity in fruit morphology among tomato cultivars can be explained by modifications at four of these cloned QTLs. In melon, a number of QTLs involved in fruit morphology have been mapped, but the molecular basis for these QTLs is unknown. In the present review, we examine the current knowledge on the molecular basis of fruit morphology in tomato and transfer that information in order to define candidate genes of melon fruit shape and size QTLs. We hypothesize that different members of the gene families identified in tomato may have a role in the regulation of fruit morphology in other species. We anchored the published melon QTL map on the genome sequence and identified the melon family members of the six cloned tomato QTLs in the genome. We investigated the co-localization of melon fruit morphology QTLs and the candidate genes. We found that QTLs for fruit weight co-localized frequently with members of the CNR/FW2.2 and KLUH/FW3.2 families, as well as co-localizations between OFP family members and fruit-shape QTLs, making this family the most suitable to explain fruit shape variation among melon accessions.

Comparison of SSR and SNP markers in estimation of genetic diversity and population structure of Indian rice varieties

Simple sequence repeat (SSR) and Single Nucleotide Polymorphic (SNP), the two most robust markers for identifying rice varieties were compared for assessment of genetic diversity and population structure. Total 375 varieties of rice from various regions of India archived at the Indian National GeneBank, NBPGR, New Delhi, were analyzed using thirty six genetic markers, each of hypervariable SSR (HvSSR) and SNP which were distributed across 12 rice chromosomes. A total of 80 alleles were amplified with the SSR markers with an average of 2.22 alleles per locus whereas, 72 alleles were amplified with SNP markers. Polymorphic information content (PIC) values for HvSSR ranged from 0.04 to 0.5 with an average of 0.25. In the case of SNP markers, PIC values ranged from 0.03 to 0.37 with an average of 0.23. Genetic relatedness among the varieties was studied; utilizing an unrooted tree all the genotypes were grouped into three major clusters with both SSR and SNP markers. Analysis of molecular variance (AMOVA) indicated that maximum diversity was partitioned between and within individual level but not between populations. Principal coordinate analysis (PCoA) with SSR markers showed that genotypes were uniformly distributed across the two axes with 13.33% of cumulative variation whereas, in case of SNP markers varieties were grouped into three broad groups across two axes with 45.20% of cumulative variation. Population structure were tested using K values from 1 to 20, but there was no clear population structure, therefore Ln(PD) derived Δk was plotted against the K to determine the number of populations. In case of SSR maximum Δk was at K=5 whereas, in case of SNP maximum Δk was found at K=15, suggesting that resolution of population was higher with SNP markers, but SSR were more efficient for diversity analysis.