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

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

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

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

A comprehensive genome variation map of melon identifies multiple domestication events and loci influencing agronomic traits.

Melon is an economically important fruit crop that has been cultivated for thousands of years; however, the genetic basis and history of its domestication still remain largely unknown. Here we report a comprehensive map of the genomic variation in melon derived from the resequencing of 1,175 accessions, which represent the global diversity of the species. Our results suggest that three independent domestication events occurred in melon, two in India and one in Africa. We detected two independent sets of domestication sweeps, resulting in diverse characteristics of the two subspecies melo and agrestis during melon breeding. Genome-wide association studies for 16 agronomic traits identified 208 loci significantly associated with fruit mass, quality and morphological characters. This study sheds light on the domestication history of melon and provides a valuable resource for genomics-assisted breeding of this important crop.

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

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

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

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

Genome of 'Charleston Gray', the principal American watermelon cultivar, and genetic characterization of 1,365 accessions in the U.S. National Plant Germplasm System watermelon collection.

Years of selection for desirable fruit quality traits in dessert watermelon (Citrullus lanatus) has resulted in a narrow genetic base in modern cultivars. Development of novel genomic and genetic resources offers great potential to expand genetic diversity and improve important traits in watermelon. Here, we report a high-quality genome sequence of watermelon cultivar 'Charleston Gray', a principal American dessert watermelon, to complement the existing reference genome from '97103', an East Asian cultivar. Comparative analyses between genomes of 'Charleston Gray' and '97103' revealed genomic variants that may underlie phenotypic differences between the two cultivars. We then genotyped 1365 watermelon plant introduction (PI) lines maintained at the U.S. National Plant Germplasm System using genotyping-by-sequencing (GBS). These PI lines were collected throughout the world and belong to three Citrullus species, C. lanatus, C. mucosospermus and C. amarus. Approximately 25 000 high-quality single nucleotide polymorphisms (SNPs) were derived from the GBS data using the 'Charleston Gray' genome as the reference. Population genomic analyses using these SNPs discovered a close relationship between C. lanatus and C. mucosospermus and identified four major groups in these two species correlated to their geographic locations. Citrullus amarus was found to have a distinct genetic makeup compared to C. lanatus and C. mucosospermus. The SNPs also enabled identification of genomic regions associated with important fruit quality and disease resistance traits through genome-wide association studies. The high-quality 'Charleston Gray' genome and the genotyping data of this large collection of watermelon accessions provide valuable resources for facilitating watermelon research, breeding and improvement.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops.

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Differential metabolism of L-phenylalanine in the formation of aromatic volatiles in melon (Cucumis melo L.) fruit.

Studies on the active pathways and the genes involved in the biosynthesis of L-phenylalanine-derived volatiles in fleshy fruits are sparse. Melon fruit rinds converted stable-isotope labeled L-phe into more than 20 volatiles. Phenylpropanes, phenylpropenes and benzenoids are apparently produced via the well-known phenylpropanoid pathway involving phenylalanine ammonia lyase (PAL) and being (E)-cinnamic acid a key intermediate. Phenethyl derivatives seemed to be derived from L-phe via a separate biosynthetic route not involving (E)-cinnamic acid and PAL. To explore for a biosynthetic route to (E)-cinnamaldehyde in melon rinds, soluble protein cell-free extracts were assayed with (E)-cinnamic acid, CoA, ATP, NADPH and MgSO4, producing (E)-cinnamaldehyde in vitro. In this context, we characterized CmCNL, a gene encoding for (E)-cinnamic acid:coenzyme A ligase, inferred to be involved in the biosynthesis of (E)-cinnamaldehyde. Additionally we describe CmBAMT, a SABATH gene family member encoding a benzoic acid:S-adenosyl-L-methionine carboxyl methyltransferase having a role in the accumulation of methyl benzoate. Our approach leads to a more comprehensive understanding of L-phe metabolism into aromatic volatiles in melon fruit.

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

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

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

Cucumis melo is highly diverse for fruit traits providing wide breeding and genetic research opportunities, including genome-wide association (GWA) analysis. We used a collection of 177 accessions representing the two C. melo subspecies and 11 horticultural groups for detailed characterization of fruit traits variation and evaluation of the potential of GWA for trait mapping in melon. Through genotyping-by-sequencing, 23,931 informative SNPs were selected for genome-wide analyses. We found that linkage-disequilibrium decays at ~100 Kb in this collection and that population structure effect on association results varies between traits. We mapped several monogenic traits to narrow intervals overlapping with known causative genes, demonstrating the potential of diverse collections and GWA for mapping Mendelian traits to a candidate-gene level in melon. We further report on mapping of fruit shape quantitative trait loci (QTLs) and comparison with multiple previous QTL studies. Expansion of sample size and a more balanced representation of taxonomic groups might improve efficiency for simple traits dissection. But, as in other plant species, integrated linkage-association multi-allelic approaches are likely to produce better combination of statistical power, diversity capture and mapping resolution in melon. Our data can be utilized for selection of the most appropriate accessions for such approaches.

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.

Distinct Mechanisms of the ORANGE Protein in Controlling Carotenoid Flux.

ß-Carotene adds nutritious value and determines the color of many fruits, including melon (Cucumis melo). In melon mesocarp, ß-carotene accumulation is governed by the Orange gene (CmOr) golden single-nucleotide polymorphism (SNP) through a yet to be discovered mechanism. In Arabidopsis (Arabidopsis thaliana), OR increases carotenoid levels by posttranscriptionally regulating phytoene synthase (PSY). Here, we identified a CmOr nonsense mutation (Cmor-lowß) that lowered fruit ß-carotene levels with impaired chromoplast biogenesis. Cmor-lowß exerted a minimal effect on PSY transcripts but dramatically decreased PSY protein levels and enzymatic activity, leading to reduced carotenoid metabolic flux and accumulation. However, the golden SNP was discovered to not affect PSY protein levels and carotenoid metabolic flux in melon fruit, as shown by carotenoid and immunoblot analyses of selected melon genotypes and by using chemical pathway inhibitors. The high ß-carotene accumulation in golden SNP melons was found to be due to a reduced further metabolism of ß-carotene. This was revealed by genetic studies with double mutants including carotenoid isomerase (yofi), a carotenoid-isomerase nonsense mutant, which arrests the turnover of prolycopene. The yofi F2 segregants accumulated prolycopene independently of the golden SNP Moreover, Cmor-lowß was found to inhibit chromoplast formation and chloroplast disintegration in fruits from 30 d after anthesis until ripening, suggesting that CmOr regulates the chloroplast-to-chromoplast transition. Taken together, our results demonstrate that CmOr is required to achieve PSY protein levels to maintain carotenoid biosynthesis metabolic flux but that the mechanism of the CmOr golden SNP involves an inhibited metabolism downstream of ß-carotene to dramatically affect both carotenoid content and plastid fate.

The biosynthetic pathway of the nonsugar, high-intensity sweetener mogroside V from Siraitia grosvenorii.

The consumption of sweeteners, natural as well as synthetic sugars, is implicated in an array of modern-day health problems. Therefore, natural nonsugar sweeteners are of increasing interest. We identify here the biosynthetic pathway of the sweet triterpenoid glycoside mogroside V, which has a sweetening strength of 250 times that of sucrose and is derived from mature fruit of luo-han-guo (Siraitia grosvenorii, monk fruit). A whole-genome sequencing of Siraitia, leading to a preliminary draft of the genome, was combined with an extensive transcriptomic analysis of developing fruit. A functional expression survey of nearly 200 candidate genes identified the members of the five enzyme families responsible for the synthesis of mogroside V: squalene epoxidases, triterpenoid synthases, epoxide hydrolases, cytochrome P450s, and UDP-glucosyltransferases. Protein modeling and docking studies corroborated the experimentally proven functional enzyme activities and indicated the order of the metabolic steps in the pathway. A comparison of the genomic organization and expression patterns of these Siraitia genes with the orthologs of other Cucurbitaceae implicates a strikingly coordinated expression of the pathway in the evolution of this species-specific and valuable metabolic pathway. The genomic organization of the pathway genes, syntenously preserved among the Cucurbitaceae, indicates, on the other hand, that gene clustering cannot account for this novel secondary metabolic pathway.

A bulk segregant transcriptome analysis reveals metabolic and cellular processes associated with Orange allelic variation and fruit ß-carotene accumulation in melon fruit.

BACKGROUND: Melon fruit flesh color is primarily controlled by the "golden" single nucleotide polymorhism of the "Orange" gene, CmOr, which dominantly triggers the accumulation of the pro-vitamin A molecule, ß-carotene, in the fruit mesocarp. The mechanism by which CmOr operates is not fully understood. To identify cellular and metabolic processes associated with CmOr allelic variation, we compared the transcriptome of bulks of developing fruit of homozygous orange and green fruited F3 families derived from a cross between orange and green fruited parental lines. RESULTS: Pooling together F3 families that share same fruit flesh color and thus the same CmOr allelic variation, normalized traits unrelated to CmOr allelic variation. RNA sequencing analysis of these bulks enabled the identification of differentially expressed genes. These genes were clustered into functional groups. The relatively enriched functional groups were those involved in photosynthesis, RNA and protein regulation, and response to stress. CONCLUSIONS: The differentially expressed genes and the enriched processes identified here by bulk segregant RNA sequencing analysis are likely part of the regulatory network of CmOr. Our study demonstrates the resolution power of bulk segregant RNA sequencing in identifying genes related to commercially important traits and provides a useful tool for better understanding the mode of action of CmOr gene in the mediation of carotenoid accumulation.

A Kelch Domain-Containing F-Box Coding Gene Negatively Regulates Flavonoid Accumulation in Muskmelon.

The flavonoids are phenylpropanoid-derived metabolites that are ubiquitous in plants, playing many roles in growth and development. Recently, we observed that fruit rinds of yellow casaba muskmelons (Cucumis melo 'Inodorous Group') accumulate naringenin chalcone, a yellow flavonoid pigment. With RNA-sequencing analysis of bulked segregants representing the tails of a population segregating for naringenin chalcone accumulation followed by fine mapping and genetic transformation, we identified a Kelch domain-containing F-box protein coding (CmKFB) gene that, when expressed, negatively regulates naringenin chalcone accumulation. Additional metabolite analysis indicated that downstream flavonoids are accumulated together with naringenin chalcone, whereas CmKFB expression diverts the biochemical flux toward coumarins and general phenylpropanoids. These results show that CmKFB functions as a posttranscriptional regulator that diverts flavonoid metabolic flux.

Characterization of endophytic bacteria from cucurbit fruits with potential benefits to agriculture in melons (Cucumis melo L.).

Endophytes are microorganisms that mainly colonize vegetative parts, but are also found in reproductive and disseminating organs, and may have beneficial characteristics. To identify microorganisms associated with the agriculturally important family, Cucurbitaceae, endophytes were initially determined in fruits of Cucumis melo Reticulatus Group 'Dulce' by a cultivation-independent approach based on fluorescence in situ hybridization using double labeling of oligonucleotide probes. Alpha-, Beta-, Gammaproteobacteria, Firmicutes and Actinobacteria were localized inside the fruits. Culturable bacteria were further isolated and identified from fruit tissues of 'Dulce', from fruits of other cultivated and wild-field-grown Cucurbitaceae, and from wild fruits growing under natural conditions. Low densities of culturable bacteria were detected in the investigated fruits, especially in four out of the five wild species, regardless of their growing environment. Substantial differences were observed between the wild and cultivated cucurbit taxa in regard to the number of colonized fruits as well as the type of endophytes. Bacillus was the most dominant genus of endophytes colonizing fruits of Cucurbitaceae. The antagonistic effects of isolated endophytes were assessed against cucurbit disease agents in dual-culture assays. Several bacterial isolates exhibited antagonistic properties against the tested plant pathogens. The identified bacteria may be useful for protecting plants not only in the field, but also for post-harvest.

A Single Amino Acid Substitution in an ORANGE Protein Promotes Carotenoid Overaccumulation in Arabidopsis.

Carotenoids are crucial for plant growth and human health. The finding of ORANGE (OR) protein as a pivotal regulator of carotenogenesis offers a unique opportunity to comprehensively understand the regulatory mechanisms of carotenoid accumulation and develop crops with enhanced nutritional quality. Here, we demonstrated that alteration of a single amino acid in a wild-type OR greatly enhanced its ability to promote carotenoid accumulation. Whereas overexpression of OR from Arabidopsis (Arabidopsis thaliana; AtOR) or from the agronomically important crop sorghum (Sorghum bicolor; SbOR) increased carotenoid levels up to 2-fold, expression of AtOR(His) (R90H) or SbOR(His) (R104H) variants dramatically enhanced carotenoid accumulation by up to 7-fold in the Arabidopsis calli. Moreover, we found that AtOR(Ala) (R90A) functioned similarly to AtOR(His) to promote carotenoid overproduction. Neither AtOR nor AtOR(His) greatly affected carotenogenic gene expression. AtOR(His) exhibited similar interactions with phytoene synthase (PSY) as AtOR in posttranscriptionally regulating PSY protein abundance. AtOR(His) triggered biogenesis of membranous chromoplasts in the Arabidopsis calli, which shared structures similar to chromoplasts found in the curd of the orange cauliflower (Brassica oleracea) mutant. By contrast, AtOR did not cause plastid-type changes in comparison with the controls, but produced plastids containing larger and electron-dense plastoglobuli. The unique ability of AtOR(His) in mediating chromoplast biogenesis is responsible for its induced carotenoid overproduction. Our study demonstrates OR(His/Ala) as powerful tools for carotenoid enrichment in plants, and provides insights into the mechanisms underlying OR(His)-regulated carotenoid accumulation.

Systems approach for exploring the intricate associations between sweetness, color and aroma in melon fruits.

BACKGROUND: Melon (Cucumis melo) fruits exhibit phenotypic diversity in several key quality determinants such as taste, color and aroma. Sucrose, carotenoids and volatiles are recognized as the key compounds shaping the above corresponding traits yet the full network of biochemical events underlying their synthesis have not been comprehensively described. To delineate the cellular processes shaping fruit quality phenotypes, a population of recombinant inbred lines (RIL) was used as a source of phenotypic and genotypic variations. In parallel, ripe fruits were analyzed for both the quantified level of 77 metabolic traits directly associated with fruit quality and for RNA-seq based expression profiles generated for 27,000 unigenes. First, we explored inter-metabolite association patterns; then, we described metabolites versus gene association patterns; finally, we used the correlation-based associations for predicting uncharacterized synthesis pathways. RESULTS: Based on metabolite versus metabolite and metabolite versus gene association patterns, we divided metabolites into two key groups: a group including ethylene and aroma determining volatiles whose accumulation patterns are correlated with the expression of genes involved in the glycolysis and TCA cycle pathways; and a group including sucrose and color determining carotenoids whose accumulation levels are correlated with the expression of genes associated with plastid formation. CONCLUSIONS: The study integrates multiple processes into a genome scale perspective of cellular activity. This lays a foundation for deciphering the role of gene markers associated with the determination of fruit quality traits.

Ultrahigh-density linkage map for cultivated cucumber (Cucumis sativus L.) using a single-nucleotide polymorphism genotyping array.

Genotyping arrays are tools for high-throughput genotyping, which is beneficial in constructing saturated genetic maps and therefore high-resolution mapping of complex traits. Since the report of the first cucumber genome draft, genetic maps have been constructed mainly based on simple-sequence repeats (SSRs) or on combinations of SSRs and sequence-related amplified polymorphism (SRAP). In this study, we developed the first cucumber genotyping array consisting of 32,864 single-nucleotide polymorphisms (SNPs). These markers cover the cucumber genome with a median interval of ~2 Kb and have expected genotype calls in parents/F1 hybridizations as a training set. The training set was validated with Fluidigm technology and showed 96% concordance with the genotype calls in the parents/F1 hybridizations. Application of the genotyping array was illustrated by constructing a 598.7 cM genetic map based on a '9930' × 'Gy14' recombinant inbred line (RIL) population comprised of 11,156 SNPs. Marker collinearity between the genetic map and reference genomes of the two parents was estimated at R2 = 0.97. We also used the array-derived genetic map to investigate chromosomal rearrangements, regional recombination rate, and specific regions with segregation distortions. Finally, 82% of the linkage-map bins were polymorphic in other cucumber variants, suggesting that the array can be applied for genotyping in other lines. The genotyping array presented here, together with the genotype calls of the parents/F1 hybridizations as a training set, should be a powerful tool in future studies with high-throughput cucumber genotyping. An ultrahigh-density linkage map constructed by this genotyping array on RIL population may be invaluable for assembly improvement, and for mapping important cucumber QTLs.