Sustainable cultivation of melon landraces: Effects of grafting on the accumulation of flavor-related compounds

Melon landraces are highly appreciated by consumers who pay price premiums to compensate for lower yields, enabling on-farm conservation. However, they are highly susceptible to soilborne diseases. This study analyses the impact of Cucurbita and Cucumis rootstocks on the accumulation of flavor-related metabolites in Spanish landraces of the Ibericus melon group, as a strategy to promote their sustainable cultivation. Scion genotype was the main factor conditioning the accumulation of sugars and acids both under standard and saline organic farming conditions. The effects of grafting on organic acid accumulation were negligible, while the effects on sugar content were significant. The latter effects were dependent on specific scion-rootstock combinations, though wild Cucumis (e.g. Fian) rootstocks represent an alternative that should be further studied. The effect on the accumulation of volatiles was limited, and again depended on specific scion-rootstock combinations. The rootstock effect even differed between populations of the same landrace.

Insights into the early transcriptomic response against watermelon mosaic virus in melon

BACKGROUND

Watermelon mosaic virus (WMV) is one of the most prevalent viruses affecting melon worldwide. Recessive resistance to WMV in melon has previously been reported in the African accession TGR-1551. Moreover, the genomic regions associated to the resistance have also been described. Nevertheless, the transcriptomic response that might infer the resistance to this potyvirus has not been explored.

RESULTS

We have performed a comparative transcriptomic analysis using mock and WMV-inoculated plants of the susceptible cultivar "Bola de oro" (BO) and a resistant RIL (Recombinant inbred line) derived from the initial cross between "TGR-1551" and BO. In total, 616 genes were identified as differentially expressed and the weighted gene co-expression network analysis (WGCNA) detected 19 gene clusters (GCs), of which 7 were differentially expressed for the genotype x treatment interaction term. SNPs with a predicted high impact on the protein function were detected within the coding regions of most of the detected DEGs. Moreover, 3 and 16 DEGs were detected within the QTL regions previously described in chromosomes 11 and 5, respectively. In addition to these two specific genomic regions, we also observde large transcriptomic changes from genes spread across the genome in the resistant plants in response to the virus infection. This early response against WMV implied genes involved in plant-pathogen interaction, plant hormone signal transduction, the MAPK signaling pathway or ubiquitin mediated proteolysis, in detriment to the photosynthetic and basal metabolites pathways. Moreover, the gene MELO3C021395, which coded a mediator of RNA polymerase II transcription subunit 33A (MED33A), has been proposed as the candidate gene located on chromosome 11 conferring resistance to WMV.

CONCLUSIONS

The comparative transcriptomic analysis presented here showed that, even though the resistance to WMV in TGR-1551 has a recessive nature, it triggers an active defense response at a transcriptomic level, which involves broad-spectrum resistance mechanisms. Thus, this study represents a step forward on our understanding of the mechanisms underlaying WMV resistance in melon. In addition, it sheds light into a broader topic on the mechanisms of recessive resistances.

Further Molecular Diagnosis Determines Lack of Evidence for Real Seed Transmission of Tomato Leaf Curl New Delhi Virus in Cucurbits

Begomoviruses (family Geminiviridae) cause serious diseases in many crop families. Since 2013, the Spanish isolate of tomato leaf curl New Delhi virus (ToLCNDV) has been a limiting factor for cucurbits production in the Mediterranean basin, forcing farmers to adapt new management and control techniques. Although it is well-known that begomoviruses are naturally transmitted by the whitefly Bemisia tabaci, the capacity of these viruses to be vertically transmitted through seeds remains controversial. Clarifying the potential ToLCNDV seed transmission is essential to understand the epidemiology of this threating-for-cucurbits virus and to design appropriate control strategies. We assessed ToLCNDV distribution in the leaves, flowers and seeds of the infected plants of susceptible Cucumis melo accessions and toleration to the infected genotypes of Cucurbita moschata by conventional and quantitative PCR. We analyzed whether the viral particle was transmitted to offspring. We also evaluated ToLCNDV presence in commercial seeds of cucurbits (zucchini (Cucurbita pepo), melon (C. melo), cucumber (Cucumis sativus) and watermelon (Citrullus lanatus)) and in their progenies. As the assayed seedlings remained symptomless, we increased the reliability and accuracy of detection in these samples by searching for replicative forms of ToLCNDV by combining Southern blot hybridization and rolling-circle amplification (RCA). However, integral genomic DNA was not identified in the plants of offspring. Although the seedborne nature of ToLCNDV was confirmed, our results do not support the transmission of this virus from contaminated seeds to progeny.

Physicochemical, Structural, and Functional Properties of Snake Melon (Cucumis melo subsp. melo Var. flexuosus) Microencapsulated with Pea Protein and Pea Fibre by Freeze-Drying

The purpose of this study was to obtain a functional, stable powder product from Cucumis melo subsp. melo Var. flexuosus (L.) to promote its consumption and reduce waste and production losses. The melons were ground and freeze-dried with or without biopolymers (pea protein (PPSM) or pea fibre (PFSM)). The physicochemical, structural, and functional properties of the powder were studied. The water content, water activity, bulk density, porosity, Hausner ratio, Carr index, hygroscopicity, water solubility, water absorption index, particle size, colour, and microstructure of the powder were determined. In addition, vitamin C, folates, chlorophyll a, total phenols and carotenoids, antioxidant capacity, and powder encapsulation efficiency were analysed. Snake melon (SM) powders contained vitamin C, folates, carotenoids, chlorophyll a, and phenols, which contributed to their antioxidant capacity. The incorporation of PP or PF in the formulation before lyophilisation generated stable encapsulates that protected the bioactive compounds. PPSM and PFSM were less hygroscopic and more free-flowing and had lower water content and water activity compared to the SM. PFSM showed higher encapsulation efficiency and smaller particles with a smooth surface and oval shape.

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.

Pre- and Post-Zygotic Barriers Contribute to Reproductive Isolation and Correlate with Genetic Distance in Cucumis

Hybridization between Cucumis species, including cultivated melon (C. melo), is hampered by Interspecific Reproductive Barriers (IRBs). However, the nature of IRBs in Cucumis is largely unknown. This study explores locations, timing, and contribution to reproductive isolation (RI) of pre- and post-zygotic IRBs in Cucumis. To do this, we assessed crossability among Cucumis African wild species and C. melo at the pre-zygotic level by visualizing pollen tubes under fluorescence microscopy and, post-zygotically, by evaluating fruit/seed set and F₁ hybrid fertility. Genetic distances among Cucumis species were inferred from Genotyping-by-Sequencing, and its correlation with RI stages was analyzed. Observed pre- and post-zygotic IRBs included pollen tube arrest, fruit set failure, and hybrid male sterility. Unilateral cross-incongruity/incompatibility (UCI) was detected in some hybridizations, and dominant gene action is suggested for pistil-side UCI in interspecific F₁ hybrids. Notably, the allotetraploid C. ficifolius was very fertile as a seed parent but infertile in all reciprocal crosses. Contribution to RI was found significant for both pre- and post-zygotic IRBs. Additionally, a significant positive correlation was detected between genetic distance and pre- and post-zygotic RI stages. Interestingly, UCI offers an accessible system to dissect the genetics of IRBs in Cucumis, which may facilitate the use of wild relatives in breeding.

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

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

Carotenoid fortification of zucchini fruits using a viral RNA vector

BACKGROUND

Carotenoids are health-promoting metabolites in livestock and human diets. Some important crops have been genetically modified to increase their content. Although the usefulness of transgenic plants to alleviate nutritional deficiencies is obvious, their social acceptance has been controversial.

RESULTS

Here, we demonstrate an alternative biotechnological strategy for carotenoid fortification of edible fruits in which no transgenic DNA is involved. A viral RNA vector derived from Zucchini yellow mosaic virus (ZYMV) was modified to express a bacterial phytoene synthase (crtB), and inoculated to zucchini (Cucurbita pepo L.) leaves nurturing pollinated flowers. After the viral vector moved to the developing fruit and expressed crtB, the rind and flesh of the fruits developed yellow-orange rather than green color. Metabolite analyses showed a substantial enrichment in health-promoting carotenoids, such as α- and β-carotene (provitamin A), lutein and phytoene, in both rind and flesh.

CONCLUSION

Although this strategy is perhaps not free from controversy due to the use of genetically modified viral RNA, our work does demonstrate the possibility of metabolically fortifying edible fruits using an approach in which no transgenes are involved. This article is protected by copyright. All rights reserved.

Combined Stress Conditions in Melon Induce Non-additive Effects in the Core miRNA Regulatory Network

Climate change has been associated with a higher incidence of combined adverse environmental conditions that can promote a significant decrease in crop productivity. However, knowledge on how a combination of stresses might affect plant development is still scarce. MicroRNAs (miRNAs) have been proposed as potential targets for improving crop productivity. Here, we have combined deep-sequencing, computational characterization of responsive miRNAs and validation of their regulatory role in a comprehensive analysis of response of melon to several combinations of four stresses (cold, salinity, short day, and infection with a fungus). Twenty-two miRNA families responding to double and/or triple stresses were identified. The regulatory role of the differentially expressed miRNAs was validated by quantitative measurements of the expression of the corresponding target genes. A high proportion (ca. 60%) of these families (mainly highly conserved miRNAs targeting transcription factors) showed a non-additive response to multiple stresses in comparison with that observed under each one of the stresses individually. Among those miRNAs showing non-additive response to stress combinations, most interactions were negative, suggesting the existence of functional convergence in the miRNA-mediated response to combined stresses. Taken together, our results provide compelling pieces of evidence that the response to combined stresses cannot be easily predicted from the study individual stresses.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

Resistant Sources and Genetic Control of Resistance to ToLCNDV in Cucumber

Tomato leaf curl New Delhi virus (ToLCNDV) is a severe threat for cucurbit production worldwide. Resistance has been reported in several crops, but at present, there are no described accessions with resistance to ToLCNDV in cucumber (Cucumis sativus). C. sativus var. sativus accessions were mechanically inoculated with ToLCNDV and screened for resistance, by scoring symptom severity, tissue printing, and PCR (conventional and quantitative). Severe symptoms and high load of viral DNA were found in plants of a nuclear collection of Spanish landraces and in accessions of C. sativus from different geographical origins. Three Indian accessions (CGN23089, CGN23423, and CGN23633) were highly resistant to the mechanical inoculation, as well as all plants of their progenies obtained by selfing. To study the inheritance of the resistance to ToLCNDV, plants of the CGN23089 accession were crossed with the susceptible accession BGV011742, and F₁ hybrids were used to construct segregating populations (F₂ and backcrosses), which were mechanically inoculated and evaluated for symptom development and viral load by qPCR. The analysis of the genetic control fit with a recessive monogenic inheritance model, and after genotyping with SNPs distributed along the C. sativus genome, a QTL associated with ToLCNDV resistance was identified in chromosome 2 of cucumber.

Grafting Snake Melon [Cucumis melo L. subsp. melo Var. flexuosus (L.) Naudin] in Organic Farming: Effects on Agronomic Performance; Resistance to Pathogens; Sugar, Acid, and VOC Profiles; and Consumer Acceptance

The performance of snake melon [Cucumis melo var. flexuosus (L.)] in organic farming was studied under high biotic and salt stress conditions. Soilborne diseases (mainly caused by Macrophomina phaseolina and Neocosmospora falciformis), combined with virus incidence [Watermelon mosaic virus (WMV), Zucchini yellow mosaic virus (ZYMV), and Tomato leaf curl New Delhi virus (ToLCNDV)] and Podosphaera xanthii attacks, reduced yield by more than 50%. Snake melon susceptibility to M. phaseolina and Monosporascus cannonballus was proved in pathogenicity tests, while it showed some degree of resistance to Neocosmospora keratoplastica and N. falciformis. On the contrary, salt stress had a minor impact, although a synergic effect was detected: yield losses caused by biotic stress increased dramatically when combined with salt stress. Under biotic stress, grafting onto the melon F₁Pat81 and wild Cucumis rootstocks consistently reduced plant mortality in different agroecological conditions, with a better performance compared to classic Cucurbita commercial hybrids. Yield was even improved under saline conditions in grafted plants. A negative effect was detected, though, on consumer acceptability, especially with the use of Cucurbita rootstocks. Cucumis F₁Pat81 rootstock minimized this side effect, which was probably related to changes in the profile of sugars, acids, and volatiles. Grafting affected sugars and organic acid contents, with this effect being more accentuated with the use of Cucurbita rootstocks than with Cucumis. In fact, the latter had a higher impact on the volatile organic compound profile than on sugar and acid profile, which may have resulted in a lower effect on consumer perception. The use of Cucumis rootstocks seems to be a strategy to enable organic farming production of snake melon targeted to high-quality markets in order to promote the cultivation of this neglected crop.

RNA-Seq Transcriptome Analysis Provides Candidate Genes for Resistance to Tomato Leaf Curl New Delhi Virus in Melon

Tomato leaf curl New Delhi virus (ToLCNDV) emerged in the Mediterranean Basin in 2012 as the first DNA bipartite begomovirus (Geminiviridae family), causing severe yield and economic losses in cucurbit crops. A major resistance locus was identified in the wild melon accession WM-7 (Cucumis melo kachri group), but the mechanisms involved in the resistant response remained unknown. In this work, we used RNA-sequencing to identify disease-associated genes that are differentially expressed in the course of ToLCNDV infection and could contribute to resistance. Transcriptomes of the resistant WM-7 genotype and the susceptible cultivar Piñonet Piel de Sapo (PS) (C. melo ibericus group) in ToLCNDV and mock inoculated plants were compared at four time points during infection (0, 3, 6, and 12 days post inoculation). Different gene expression patterns were observed over time in the resistant and susceptible genotypes in comparison to their respective controls. Differentially expressed genes (DEGs) in ToLCNDV-infected plants were classified using gene ontology (GO) terms, and genes of the categories transcription, DNA replication, and helicase activity were downregulated in WM-7 but upregulated in PS, suggesting that reduced activity of these functions reduces ToLCNDV replication and intercellular spread and thereby contributes to resistance. DEGs involved in the jasmonic acid signaling pathway, photosynthesis, RNA silencing, transmembrane, and sugar transporters entail adverse consequences for systemic infection in the resistant genotype, and lead to susceptibility in PS. The expression levels of selected candidate genes were validated by qRT-PCR to corroborate their differential expression upon ToLCNDV infection in resistant and susceptible melon. Furthermore, single nucleotide polymorphism (SNPs) with an effect on structural functionality of DEGs linked to the main QTLs for ToLCNDV resistance have been identified. The obtained results pinpoint cellular functions and candidate genes that are differentially expressed in a resistant and susceptible melon line in response to ToLCNDV, an information of great relevance for breeding ToLCNDV-resistant melon cultivars.

Synthetic conversion of leaf chloroplasts into carotenoid-rich plastids reveals mechanistic basis of natural chromoplast development

Plastids, the defining organelles of plant cells, undergo physiological and morphological changes to fulfill distinct biological functions. In particular, the differentiation of chloroplasts into chromoplasts results in an enhanced storage capacity for carotenoids with industrial and nutritional value such as beta-carotene (provitamin A). Here, we show that synthetically inducing a burst in the production of phytoene, the first committed intermediate of the carotenoid pathway, elicits an artificial chloroplast-to-chromoplast differentiation in leaves. Phytoene overproduction initially interferes with photosynthesis, acting as a metabolic threshold switch mechanism that weakens chloroplast identity. In a second stage, phytoene conversion into downstream carotenoids is required for the differentiation of chromoplasts, a process that involves a concurrent reprogramming of nuclear gene expression and plastid morphology for improved carotenoid storage. We hence demonstrate that loss of photosynthetic competence and enhanced production of carotenoids are not just consequences but requirements for chloroplasts to differentiate into chromoplasts.

A Major QTL Located in Chromosome 8 of Cucurbita moschata Is Responsible for Resistance to Tomato Leaf Curl New Delhi Virus

Tomato leaf curl New Delhi virus (ToLCNDV) is a bipartite whitefly transmitted begomovirus, responsible since 2013 of severe damages in cucurbit crops in Southeastern Spain. Zucchini (Cucurbita pepo) is the most affected species, but melon (Cucumis melo) and cucumber (Cucumis sativus) are also highly damaged by the infection. The virus has spread across Mediterranean basin and European countries, and integrated control measures are not being enough to reduce economic losses. The identification of resistance genes is required to develop resistant cultivars. In this assay, we studied the inheritance of the resistance to ToLCNDV previously identified in two Cucurbita moschata accessions. We generated segregating populations crossing both resistant pumpkins, an American improved cultivar Large Cheese (PI 604506) and an Indian landrace (PI 381814), with a susceptible C. moschata genotype (PI 419083). The analysis of symptoms and viral titers of all populations established the same monogenic recessive genetic control in both resistant accessions, and the allelism tests suggest the occurrence of alleles of the same locus. By genotyping with a single nucleotide polymorphism (SNP) collection evenly distributed along the C. moschata genome, a major quantitative trait locus (QTL) was identified in chromosome 8 controlling resistance to ToLCNDV. This major QTL was also confirmed in the interspecific C. moschata × C. pepo segregating populations, although C. pepo genetic background affected the resistance level. Molecular markers here identified, linked to the ToLCNDV resistance locus, are highly valuable for zucchini breeding programs, allowing the selection of improved commercial materials. The duplication of the candidate region within the C. moschata genome was studied, and genes with paralogs or single-copy genes were identified. Its synteny with the region of chromosome 17 of the susceptible C. pepo revealed an INDEL including interesting candidate genes. The chromosome 8 candidate region of C. moschata was also syntenic to the region in chromosome 11 of melon, previously described as responsible of ToLCNDV resistance. Common genes in the candidate regions of both cucurbits, with high- or moderate-impact polymorphic SNPs between resistant and susceptible C. moschata accessions, are interesting to study the mechanisms involved in this recessive resistance.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

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.

De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus

The Cucurbita genus (squashes, pumpkins and gourds) includes important domesticated species such as C. pepo, C. maxima and C. moschata. In this study, we present a high-quality draft of the zucchini (C. pepo) genome. The assembly has a size of 263 Mb, a scaffold N50 of 1.8 Mb and 34 240 gene models. It includes 92% of the conserved BUSCO core gene set, and it is estimated to cover 93.0% of the genome. The genome is organized in 20 pseudomolecules that represent 81.4% of the assembly, and it is integrated with a genetic map of 7718 SNPs. Despite the small genome size, three independent lines of evidence support that the C. pepo genome is the result of a whole-genome duplication: the topology of the gene family phylogenies, the karyotype organization and the distribution of 4DTv distances. Additionally, 40 transcriptomes of 12 species of the genus were assembled and analysed together with all the other published genomes of the Cucurbitaceae family. The duplication was detected in all the Cucurbita species analysed, including C. maxima and C. moschata, but not in the more distant cucurbits belonging to the Cucumis and Citrullus genera, and it is likely to have occurred 30 ± 4 Mya in the ancestral species that gave rise to the genus.

De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus

The Cucurbita genus (squashes, pumpkins and gourds) includes important domesticated species such as C. pepo, C. maxima and C. moschata. In this study, we present a high-quality draft of the zucchini (C. pepo) genome. The assembly has a size of 263 Mb, a scaffold N50 of 1.8 Mb and 34 240 gene models. It includes 92% of the conserved BUSCO core gene set, and it is estimated to cover 93.0% of the genome. The genome is organized in 20 pseudomolecules that represent 81.4% of the assembly, and it is integrated with a genetic map of 7718 SNPs. Despite the small genome size, three independent lines of evidence support that the C. pepo genome is the result of a whole-genome duplication: the topology of the gene family phylogenies, the karyotype organization and the distribution of 4DTv distances. Additionally, 40 transcriptomes of 12 species of the genus were assembled and analysed together with all the other published genomes of the Cucurbitaceae family. The duplication was detected in all the Cucurbita species analysed, including C. maxima and C. moschata, but not in the more distant cucurbits belonging to the Cucumis and Citrullus genera, and it is likely to have occurred 30 ± 4 Mya in the ancestral species that gave rise to the genus.

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.

An SNP-based saturated genetic map and QTL analysis of fruit-related traits in Zucchini using Genotyping-by-sequencing

Background

Cucurbita pepo is a cucurbit with growing economic importance worldwide. Zucchini morphotype is the most important within this highly variable species. Recently, transcriptome and Simple Sequence Repeat (SSR)- and Single Nucleotide Polymorphism (SNP)-based medium density maps have been reported, however further genomic tools are needed for efficient molecular breeding in the species. Our objective is to combine currently available complete transcriptomes and the Zucchini genome sequence with high throughput genotyping methods, mapping population development and extensive phenotyping to facilitate the advance of genomic research in this species.

Results

We report the Genotyping-by-sequencing analysis of a RIL population developed from the inter subspecific cross Zucchini x Scallop (ssp. pepo x ssp. ovifera). Several thousands of SNP markers were identified and genotyped, followed by the construction of a high-density linkage map based on 7,718 SNPs (average of 386 markers/linkage group) covering 2,817.6 cM of the whole genome, which is a great improvement with respect to previous maps. A QTL analysis was performed using phenotypic data obtained from the RIL population from three environments. In total, 48 consistent QTLs for vine, flowering and fruit quality traits were detected on the basis of a multiple-environment analysis, distributed in 33 independent positions in 15 LGs, and each QTL explained 1.5–62.9% of the phenotypic variance. Eight major QTLs, which could explain greater than 20% of the phenotypic variation were detected and the underlying candidate genes identified.

Conclusions

Here we report the first SNP saturated map in the species, anchored to the physical map. Additionally, several consistent QTLs related to early flowering, fruit shape and length, and rind and flesh color are reported as well as candidate genes for them. This information will enhance molecular breeding in C. pepo and will assist the gene cloning underlying the studied QTLs, helping to reveal the genetic basis of the studied processes in squash.

Electronic supplementary material

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

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.

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.

Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.)

BACKGROUND

We have identified a kind of parthenocarpy in zucchini squash which is associated with an incomplete andromonoecy, i.e. a partial conversion of female into bisexual flowers. Given that andromonoecy in this and other cucurbit species is caused by a reduction of ethylene production in the female flower, the associated parthenocarpic development of the fruit suggested the involvement of ethylene in fruit set and early fruit development.

RESULTS

We have compared the production of ethylene as well as the expression of 13 ethylene biosynthesis and signalling genes in pollinated and unpollinated ovaries/fruits of two cultivars, one of which is parthenocarpic (Cavili), while the other is non-parthenocarpic (Tosca). In the latter, unpollinated ovaries show an induction of ethylene biosynthesis and ethylene signal transduction pathway genes three days after anthesis, which is concomitant with the initiation of fruit abortion and senescence. Fruit set and early fruit development in pollinated flowers of both cultivars and unpollinated flowers of Cavili is coupled with low ethylene biosynthesis and signalling, which would also explain the partial andromonoecy in the parthenocarpic genotype. The reduction of ethylene production in the ovary cosegregates with parthenocarpy and partial andromonoecy in the selfing progeny of Cavili. Moreover, the induction of ethylene in anthesis (by ethephon treatments) reduced the percentage of bisexual parthenocarpic flowers in Cavili, while the inhibition of ethylene biosynthesis or response (by AVG and STS treatments) induces not only andromonoecy but also the parthenocarpic development of the fruit in both cultivars.

CONCLUSIONS

Results demonstrate that a reduction of ethylene production or signalling in the zucchini flower is able to induce fruit set and early fruit development, and therefore that ethylene is actively involved in fruit set and early fruit development. Auxin and TIBA treatments, inducing fruit set and early fruit development in this species, also inhibit ethylene production and the expression of ethylene biosynthesis and response genes. A model is presented that discusses the crosstalk between ethylene and auxin in the control of fruit set and early fruit development in zucchini squash.

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

Novel sequencing technologies were recently used to generate sequences from multiple melon (Cucumis melo L.) genotypes, enabling the in silico identification of large single nucleotide polymorphism (SNP) collections. In order to optimize the use of these markers, SNP validation and large-scale genotyping are necessary. In this paper, we present the first validated design for a genotyping array with 768 SNPs that are evenly distributed throughout the melon genome. This customized Illumina GoldenGate assay was used to genotype a collection of 74 accessions, representing most of the botanical groups of the species. Of the assayed loci, 91 % were successfully genotyped. The array provided a large number of polymorphic SNPs within and across accessions. This set of SNPs detected high levels of variation in accessions from this crop's center of origin as well as from several other areas of melon diversification. Allele distribution throughout the genome revealed regions that distinguished between the two main groups of cultivated accessions (inodorus and cantalupensis). Population structure analysis showed a subdivision into five subpopulations, reflecting the history of the crop. A considerably low level of LD was detected, which decayed rapidly within a few kilobases. Our results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in melon. Since many of the genotyped accessions are currently being used as the parents of breeding populations in various programs, this set of mapped markers could be used for future mapping and breeding efforts.

Root transcriptional responses of two melon genotypes with contrasting resistance to Monosporascus cannonballus (Pollack et Uecker) infection

BACKGROUND

Monosporascus cannonballus is the main causal agent of melon vine decline disease. Several studies have been carried out mainly focused on the study of the penetration of this pathogen into melon roots, the evaluation of symptoms severity on infected roots, and screening assays for breeding programs. However, a detailed molecular view on the early interaction between M. cannonballus and melon roots in either susceptible or resistant genotypes is lacking. In the present study, we used a melon oligo-based microarray to investigate the gene expression responses of two melon genotypes, Cucumis melo 'Piel de sapo' ('PS') and C. melo 'Pat 81', with contrasting resistance to the disease. This study was carried out at 1 and 3 days after infection (DPI) by M. cannonballus.

RESULTS

Our results indicate a dissimilar behavior of the susceptible vs. the resistant genotypes from 1 to 3 DPI. 'PS' responded with a more rapid infection response than 'Pat 81' at 1 DPI. At 3 DPI the total number of differentially expressed genes identified in 'PS' declined from 451 to 359, while the total number of differentially expressed transcripts in 'Pat 81' increased from 187 to 849. Several deregulated transcripts coded for components of Ca2+ and jasmonic acid (JA) signalling pathways, as well as for other proteins related to defence mechanisms. Transcriptional differences in the activation of the JA-mediated response in 'Pat 81' compared to 'PS' suggested that JA response might be partially responsible for their observed differences in resistance.

CONCLUSIONS

As a result of this study we have identified for the first time a set of candidate genes involved in the root response to the infection of the pathogen causing melon vine decline. This information is useful for understanding the disease progression and resistance mechanisms few days after inoculation.

Application of genomic tools in plant breeding

Plant breeding has been very successful in developing improved varieties using conventional tools and methodologies. Nowadays, the availability of genomic tools and resources is leading to a new revolution of plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. Next Generation Sequencing (NGS) technologies are allowing the mass sequencing of genomes and transcriptomes, which is producing a vast array of genomic information. The analysis of NGS data by means of bioinformatics developments allows discovering new genes and regulatory sequences and their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Genomic approaches include TILLING and EcoTILLING, which make possible to screen mutant and germplasm collections for allelic variants in target genes. Re-sequencing of genomes is very useful for the genome-wide discovery of markers amenable for high-throughput genotyping platforms, like SSRs and SNPs, or the construction of high density genetic maps. All these tools and resources facilitate studying the genetic diversity, which is important for germplasm management, enhancement and use. Also, they allow the identification of markers linked to genes and QTLs, using a diversity of techniques like bulked segregant analysis (BSA), fine genetic mapping, or association mapping. These new markers are used for marker assisted selection, including marker assisted backcross selection, ‘breeding by design’, or new strategies, like genomic selection. In conclusion, advances in genomics are providing breeders with new tools and methodologies that allow a great leap forward in plant breeding, including the ‘superdomestication’ of crops and the genetic dissection and breeding for complex traits.

Transcriptome sequencing for SNP discovery across Cucumis melo

Background

Melon (Cucumis melo L.) is a highly diverse species that is cultivated worldwide. Recent advances in massively parallel sequencing have begun to allow the study of nucleotide diversity in this species. The Sanger method combined with medium-throughput 454 technology were used in a previous study to analyze the genetic diversity of germplasm representing 3 botanical varieties, yielding a collection of about 40,000 SNPs distributed in 14,000 unigenes. However, the usefulness of this resource is limited as the sequenced genotypes do not represent the whole diversity of the species, which is divided into two subspecies with many botanical varieties variable in plant, flowering, and fruit traits, as well as in stress response. As a first step to extensively document levels and patterns of nucleotide variability across the species, we used the high-throughput SOLiD™ system to resequence the transcriptomes of a set of 67 genotypes that had previously been selected from a core collection representing the extant variation of the entire species.

Results

The deep transcriptome resequencing of all of the genotypes, grouped into 8 pools (wild African agrestis, Asian agrestis and acidulus, exotic Far Eastern conomon, Indian momordica and Asian dudaim and flexuosus, commercial cantalupensis, subsp. melo Asian and European landraces, Spanish inodorus landraces, and Piel de Sapo breeding lines) yielded about 300 M reads. Short reads were mapped to the recently generated draft genome assembly of the DHL line Piel de Sapo (inodorus) x Songwhan Charmi (conomon) and to a new version of melon transcriptome. Regions with at least 6X coverage were used in SNV calling, generating a melon collection with 303,883 variants. These SNVs were dispersed across the entire C. melo genome, and distributed in 15,064 annotated genes. The number and variability of in silico SNVs differed considerably between pools. Our finding of higher genomic diversity in wild and exotic agrestis melons from India and Africa as compared to commercial cultivars, cultigens and landraces from Eastern Europe, Western Asia and the Mediterranean basin is consistent with the evolutionary history proposed for the species. Group-specific SNVs that will be useful in introgression programs were also detected. In a sample of 143 selected putative SNPs, we verified 93% of the polymorphisms in a panel of 78 genotypes.

Conclusions

This study provides the first comprehensive resequencing data for wild, exotic, and cultivated (landraces and commercial) melon transcriptomes, yielding the largest melon SNP collection available to date and representing a notable sample of the species diversity. This data provides a valuable resource for creating a catalog of allelic variants of melon genes and it will aid in future in-depth studies of population genetics, marker-assisted breeding, and gene identification aimed at developing improved varieties.

Transcriptome sequencing for SNP discovery across Cucumis melo

BACKGROUND

Melon (Cucumis melo L.) is a highly diverse species that is cultivated worldwide. Recent advances in massively parallel sequencing have begun to allow the study of nucleotide diversity in this species. The Sanger method combined with medium-throughput 454 technology were used in a previous study to analyze the genetic diversity of germplasm representing 3 botanical varieties, yielding a collection of about 40,000 SNPs distributed in 14,000 unigenes. However, the usefulness of this resource is limited as the sequenced genotypes do not represent the whole diversity of the species, which is divided into two subspecies with many botanical varieties variable in plant, flowering, and fruit traits, as well as in stress response. As a first step to extensively document levels and patterns of nucleotide variability across the species, we used the high-throughput SOLiD™ system to resequence the transcriptomes of a set of 67 genotypes that had previously been selected from a core collection representing the extant variation of the entire species.

RESULTS

The deep transcriptome resequencing of all of the genotypes, grouped into 8 pools (wild African agrestis, Asian agrestis and acidulus, exotic Far Eastern conomon, Indian momordica and Asian dudaim and flexuosus, commercial cantalupensis, subsp. melo Asian and European landraces, Spanish inodorus landraces, and Piel de Sapo breeding lines) yielded about 300 M reads. Short reads were mapped to the recently generated draft genome assembly of the DHL line Piel de Sapo (inodorus) x Songwhan Charmi (conomon) and to a new version of melon transcriptome. Regions with at least 6X coverage were used in SNV calling, generating a melon collection with 303,883 variants. These SNVs were dispersed across the entire C. melo genome, and distributed in 15,064 annotated genes. The number and variability of in silico SNVs differed considerably between pools. Our finding of higher genomic diversity in wild and exotic agrestis melons from India and Africa as compared to commercial cultivars, cultigens and landraces from Eastern Europe, Western Asia and the Mediterranean basin is consistent with the evolutionary history proposed for the species. Group-specific SNVs that will be useful in introgression programs were also detected. In a sample of 143 selected putative SNPs, we verified 93% of the polymorphisms in a panel of 78 genotypes.

CONCLUSIONS

This study provides the first comprehensive resequencing data for wild, exotic, and cultivated (landraces and commercial) melon transcriptomes, yielding the largest melon SNP collection available to date and representing a notable sample of the species diversity. This data provides a valuable resource for creating a catalog of allelic variants of melon genes and it will aid in future in-depth studies of population genetics, marker-assisted breeding, and gene identification aimed at developing improved varieties.

Diversity in expression of phosphorus (P) responsive genes in Cucumis melo L

BACKGROUND

Phosphorus (P) is a major limiting nutrient for plant growth in many soils. Studies in model species have identified genes involved in plant adaptations to low soil P availability. However, little information is available on the genetic bases of these adaptations in vegetable crops. In this respect, sequence data for melon now makes it possible to identify melon orthologues of candidate P responsive genes, and the expression of these genes can be used to explain the diversity in the root system adaptation to low P availability, recently observed in this species.

METHODOLOGY AND FINDINGS

Transcriptional responses to P starvation were studied in nine diverse melon accessions by comparing the expression of eight candidate genes (Cm-PAP10.1, Cm-PAP10.2, Cm-RNS1, Cm-PPCK1, Cm-transferase, Cm-SQD1, Cm-DGD1 and Cm-SPX2) under P replete and P starved conditions. Differences among melon accessions were observed in response to P starvation, including differences in plant morphology, P uptake, P use efficiency (PUE) and gene expression. All studied genes were up regulated under P starvation conditions. Differences in the expression of genes involved in P mobilization and remobilization (Cm-PAP10.1, Cm-PAP10.2 and Cm-RNS1) under P starvation conditions explained part of the differences in P uptake and PUE among melon accessions. The levels of expression of the other studied genes were diverse among melon accessions, but contributed less to the phenotypical response of the accessions.

CONCLUSIONS

This is the first time that these genes have been described in the context of P starvation responses in melon. There exists significant diversity in gene expression levels and P use efficiency among melon accessions as well as significant correlations between gene expression levels and phenotypical measurements.

High-throughput SNP genotyping in Cucurbita pepo for map construction and quantitative trait loci mapping

Background

Cucurbita pepo is a member of the Cucurbitaceae family, the second- most important horticultural family in terms of economic importance after Solanaceae. The "summer squash" types, including Zucchini and Scallop, rank among the highest-valued vegetables worldwide. There are few genomic tools available for this species.

The first Cucurbita transcriptome, along with a large collection of Single Nucleotide Polymorphisms (SNP), was recently generated using massive sequencing. A set of 384 SNP was selected to generate an Illumina GoldenGate assay in order to construct the first SNP-based genetic map of Cucurbita and map quantitative trait loci (QTL).

Results

We herein present the construction of the first SNP-based genetic map of Cucurbita pepo using a population derived from the cross of two varieties with contrasting phenotypes, representing the main cultivar groups of the species' two subspecies: Zucchini (subsp. pepo) × Scallop (subsp. ovifera). The mapping population was genotyped with 384 SNP, a set of selected EST-SNP identified in silico after massive sequencing of the transcriptomes of both parents, using the Illumina GoldenGate platform. The global success rate of the assay was higher than 85%. In total, 304 SNP were mapped, along with 11 SSR from a previous map, giving a map density of 5.56 cM/marker. This map was used to infer syntenic relationships between C. pepo and cucumber and to successfully map QTL that control plant, flowering and fruit traits that are of benefit to squash breeding. The QTL effects were validated in backcross populations.

Conclusion

Our results show that massive sequencing in different genotypes is an excellent tool for SNP discovery, and that the Illumina GoldenGate platform can be successfully applied to constructing genetic maps and performing QTL analysis in Cucurbita. This is the first SNP-based genetic map in the Cucurbita genus and is an invaluable new tool for biological research, especially considering that most of these markers are located in the coding regions of genes involved in different physiological processes. The platform will also be useful for future mapping and diversity studies, and will be essential in order to accelerate the process of breeding new and better-adapted squash varieties.

Towards a TILLING platform for functional genomics in Piel de Sapo melons

BACKGROUND

The availability of genetic and genomic resources for melon has increased significantly, but functional genomics resources are still limited for this crop. TILLING is a powerful reverse genetics approach that can be utilized to generate novel mutations in candidate genes. A TILLING resource is available for cantalupensis melons, but not for inodorus melons, the other main commercial group.

RESULTS

A new ethyl methanesulfonate-mutagenized (EMS) melon population was generated for the first time in an andromonoecious non-climacteric inodorus Piel de Sapo genetic background. Diverse mutant phenotypes in seedlings, vines and fruits were observed, some of which were of possible commercial interest. The population was first screened for mutations in three target genes involved in disease resistance and fruit quality (Cm-PDS, Cm-eIF4E and Cm-eIFI(iso)4E). The same genes were also tilled in the available monoecious and climacteric cantalupensis EMS melon population. The overall mutation density in this first Piel de Sapo TILLING platform was estimated to be 1 mutation/1.5 Mb by screening four additional genes (Cm-ACO1, Cm-NOR, Cm-DET1 and Cm-DHS). Thirty-three point mutations were found for the seven gene targets, six of which were predicted to have an impact on the function of the protein. The genotype/phenotype correlation was demonstrated for a loss-of-function mutation in the Phytoene desaturase gene, which is involved in carotenoid biosynthesis.

CONCLUSIONS

The TILLING approach was successful at providing new mutations in the genetic background of Piel de Sapo in most of the analyzed genes, even in genes for which natural variation is extremely low. This new resource will facilitate reverse genetics studies in non-climacteric melons, contributing materially to future genomic and breeding studies.

Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae)

Background

Cucurbita pepo belongs to the Cucurbitaceae family. The "Zucchini" types rank among the highest-valued vegetables worldwide, and other C. pepo and related Cucurbita spp., are food staples and rich sources of fat and vitamins. A broad range of genomic tools are today available for other cucurbits that have become models for the study of different metabolic processes. However, these tools are still lacking in the Cucurbita genus, thus limiting gene discovery and the process of breeding.

Results

We report the generation of a total of 512,751 C. pepo EST sequences, using 454 GS FLX Titanium technology. ESTs were obtained from normalized cDNA libraries (root, leaves, and flower tissue) prepared using two varieties with contrasting phenotypes for plant, flowering and fruit traits, representing the two C. pepo subspecies: subsp. pepo cv. Zucchini and subsp. ovifera cv Scallop. De novo assembling was performed to generate a collection of 49,610 Cucurbita unigenes (average length of 626 bp) that represent the first transcriptome of the species. Over 60% of the unigenes were functionally annotated and assigned to one or more Gene Ontology terms. The distributions of Cucurbita unigenes followed similar tendencies than that reported for Arabidopsis or melon, suggesting that the dataset may represent the whole Cucurbita transcriptome. About 34% unigenes were detected to have known orthologs of Arabidopsis or melon, including genes potentially involved in disease resistance, flowering and fruit quality. Furthermore, a set of 1,882 unigenes with SSR motifs and 9,043 high confidence SNPs between Zucchini and Scallop were identified, of which 3,538 SNPs met criteria for use with high throughput genotyping platforms, and 144 could be detected as CAPS. A set of markers were validated, being 80% of them polymorphic in a set of variable C. pepo and C. moschata accessions.

Conclusion

We present the first broad survey of gene sequences and allelic variation in C. pepo, where limited prior genomic information existed. The transcriptome provides an invaluable new tool for biological research. The developed molecular markers are the basis for future genetic linkage and quantitative trait loci analysis, and will be essential to speed up the process of breeding new and better adapted squash varieties.

Bin mapping of genomic and EST-derived SSRs in melon (Cucumis melo L.)

We report the development of 158 primer pairs flanking SSR motifs in genomic (gSSR) and EST (EST-SSR) melon sequences, all yielding polymorphic bands in melon germplasm, except one that was polymorphic only in Cucurbita species. A similar polymorphism level was found among EST-SSRs and gSSRs, between dimeric and trimeric EST-SSRs, and between EST-SSRs placed in the open reading frame or any of the 5'- or 3'-untranslated regions. Correlation between SSR length and polymorphism was only found for dinucleotide EST-SSRs located within the untranslated regions, but not for trinucleotide EST-SSRs. Transferability of EST-SSRs to Cucurbita species was assayed and 12.7% of the primer pairs amplified at least in one species, although only 5.4% were polymorphic. A set of 14 double haploid lines from the cross between the cultivar "Piel de Sapo" and the accession PI161375 were selected for the bin mapping approach in melon. One hundred and twenty-one SSR markers were newly mapped. The position of 46 SSR loci was also verified by genotyping the complete population. A final bin-map was constructed including 80 RFLPs, 212 SSRs, 3 SNPs and the Nsv locus, distributed in 122 bins with an average bin length of 10.2 cM and a maximum bin length of 33 cM. Map density was 4.2 cM/marker or 5.9 cM/SSR.

MELOGEN: an EST database for melon functional genomics

Background

Melon (Cucumis melo L.) is one of the most important fleshy fruits for fresh consumption. Despite this, few genomic resources exist for this species. To facilitate the discovery of genes involved in essential traits, such as fruit development, fruit maturation and disease resistance, and to speed up the process of breeding new and better adapted melon varieties, we have produced a large collection of expressed sequence tags (ESTs) from eight normalized cDNA libraries from different tissues in different physiological conditions.

Results

We determined over 30,000 ESTs that were clustered into 16,637 non-redundant sequences or unigenes, comprising 6,023 tentative consensus sequences (contigs) and 10,614 unclustered sequences (singletons). Many potential molecular markers were identified in the melon dataset: 1,052 potential simple sequence repeats (SSRs) and 356 single nucleotide polymorphisms (SNPs) were found. Sixty-nine percent of the melon unigenes showed a significant similarity with proteins in databases. Functional classification of the unigenes was carried out following the Gene Ontology scheme. In total, 9,402 unigenes were mapped to one or more ontology. Remarkably, the distributions of melon and Arabidopsis unigenes followed similar tendencies, suggesting that the melon dataset is representative of the whole melon transcriptome. Bioinformatic analyses primarily focused on potential precursors of melon micro RNAs (miRNAs) in the melon dataset, but many other genes potentially controlling disease resistance and fruit quality traits were also identified. Patterns of transcript accumulation were characterised by Real-Time-qPCR for 20 of these genes.

Conclusion

The collection of ESTs characterised here represents a substantial increase on the genetic information available for melon. A database (MELOGEN) which contains all EST sequences, contig images and several tools for analysis and data mining has been created. This set of sequences constitutes also the basis for an oligo-based microarray for melon that is being used in experiments to further analyse the melon transcriptome.

Quantitative detection of Cucumber vein yellowing virus in susceptible and partially resistant plants using real-time PCR

A method for the detection of Cucumber vein yellowing virus (CVYV) that combines reverse transcription with real-time PCR (SYBR((R)) Green chemistry) was developed using specific primers designed from a nucleotide sequence of the RNA polymerase gene (NIb) conserved among all the available CVYV strains. This method provided a linear assay over five to six orders of magnitude and reproducibly detected titres as low as 10(3) molecules of the target CVYV cDNA. Real-time PCR gave reproducible results for the quantification of CVYV in young leaves of susceptible and resistant cucumber landraces after mechanical inoculation. Significant differences in the starting amount of target cDNA were found between the analyzed genotypes, indicating differences in viral accumulation that correlated to their different levels of resistance. Real-time PCR results validated our previous findings using slot-blot hybridization, the dominance of the strong resistance to CVYV displayed by C.sat 10, and provided improved reliability and sensitivity of detection. This method has great potential in resistance breeding for germplasm screening, characterization of resistance mechanisms and genetic studies.

Genetic diversity of a germplasm collection of Cucurbita pepo using SRAP and AFLP markers

Cucurbita pepo is a highly polymorphic species. The cultivars can be grouped into eight morphotypes in two subspecies, ssp. pepo and ssp. ovifera. A collection of 69 accessions representative of the morphotypes and some unclassified types was used for analysing the morphological and molecular diversity of this species. This collection includes commercial cultivars and Spanish landraces, which represent the great diversification of types that have arisen in Europe after this species arrived from America. For the molecular variability studies, two PCR-based systems were employed, AFLP and SRAP, which preferentially amplify ORFs. Principal coordinates analysis and cluster analysis using the UPGMA method clearly separate the accessions into the two subspecies through the use of both markers. However, the gene diversity and the genetic identity values among morphotypes and subspecies varied between the two marker systems. The information given by SRAP markers was more concordant to the morphological variability and to the evolutionary history of the morphotypes than that of AFLP markers. In ssp. ovifera, the accessions of the different morphotypes were basically grouped according to the fruit colour. This may indicate different times of development and also the extent of breeding in the accessions used. This study has allowed identification of new types that can be employed for the development of new cultivars. The landraces of the spp. ovifera, used as ornamental in Europe, have proved to be of great interest for preserving the diversity of C. pepo.