Characterization of the olive endophytic community in genotypes displaying a contrasting response to Xylella fastidiosa.

BACKGROUND: Endophytes mediate the interactions between plants and other microorganisms, and the functional aspects of interactions between endophytes and their host that support plant-growth promotion and tolerance to stresses signify the ecological relevance of the endosphere microbiome. In this work, we studied the bacterial and fungal endophytic communities of olive tree (Olea europaea L.) asymptomatic or low symptomatic genotypes sampled in groves heavily compromised by Xylella fastidiosa subsp. pauca, aiming to characterize microbiota in genotypes displaying differential response to the pathogen. RESULTS: The relationships between bacterial and fungal genera were analyzed both separately and together, in order to investigate the intricate correlations between the identified Operational Taxonomic Units (OTUs). Results suggested a dominant role of the fungal endophytic community compared to the bacterial one, and highlighted specific microbial taxa only associated with asymptomatic or low symptomatic genotypes. In addition, they indicated the occurrence of well-adapted genetic resources surviving after years of pathogen pressure in association with microorganisms such as Burkholderia, Quambalaria, Phaffia and Rhodotorula. CONCLUSIONS: This is the first study to overview endophytic communities associated with several putatively resistant olive genotypes in areas under high X. fastidiosa inoculum pressure. Identifying these negatively correlated genera can offer valuable insights into the potential antagonistic microbial resources and their possible development as biocontrol agents.

High-density linkage mapping and genetic dissection of resistance to broomrape (Orobanche crenata Forsk.) in pea (Pisum sativum L.).

Pea (Pisum sativum L.) is a widely cultivated legume of major importance for global food security and agricultural sustainability. Crenate broomrape (Orobanche crenata Forsk.) (Oc) is a parasitic weed severely affecting legumes, including pea, in the Mediterranean Basin and the Middle East. Previously, the identification of the pea line "ROR12", displaying resistance to Oc, was reported. Two-year field trials on a segregant population of 148 F7 recombinant inbred lines (RILs), originating from a cross between "ROR12" and the susceptible cultivar "Sprinter", revealed high heritability (0.84) of the "ROR12" resistance source. Genotyping-by-sequencing (GBS) on the same RIL population allowed the construction of a high-density pea linkage map, which was compared with the pea reference genome and used for quantitative trait locus (QTL) mapping. Three QTLs associated with the response to Oc infection, named PsOcr-1, PsOcr-2, and PsOcr-3, were identified, with PsOcr-1 explaining 69.3% of the genotypic variance. Evaluation of the effects of different genotypic combinations indicated additivity between PsOcr-1 and PsOcr-2, and between PsOcr-1 and PsOcr-3, and epistasis between PsOcr-2 and PsOcr-3. Finally, three Kompetitive Allele Specific PCR (KASP) marker assays were designed on the single-nucleotide polymorphisms (SNPs) associated with the QTL significance peaks. Besides contributing to the development of pea genomic resources, this work lays the foundation for the obtainment of pea cultivars resistant to Oc and the identification of genes involved in resistance to parasitic Orobanchaceae.

Detection and distribution of two dominant alleles associated with the sweet kernel phenotype in almond cultivated germplasm.

Almond [Prunus dulcis Miller (D. A. Webb), syn. Prunus amygdalus L.)] is the major tree nut crop worldwide in terms of production and cultivated area. Almond domestication was enabled by the selection of individuals bearing sweet kernels, which do not accumulate high levels of the toxic cyanogenic glucoside amygdalin. Previously, we showed that the Sweet kernel (Sk) gene, controlling the kernel taste in almond, encodes a basic helix loop helix (bHLH) transcription factor regulating the amygdalin biosynthetic pathway. In addition, we characterized a dominant allele of this gene, further referred to as Sk-1, which originates from a C1036→T missense mutation and confers the sweet kernel phenotype. Here we provide evidence indicating that the allele further referred to as Sk-2, originally detected in the cultivar "Atocha" and arising from a T989→G missense mutation, is also dominantly inherited and confers the sweet kernel phenotype in almond cultivated germplasm. The use of single nucleotide polymorphism (SNP) data from genotyping by sequencing (GBS) for population structure and hierarchical clustering analyses indicated that Sk-2 occurs in a group of related genotypes, including the widespread cultivar "Texas", descending from the same ancestral population. KASP and dual label functional markers were developed for the accurate and high-throughput selection of the Sk-1 and Sk-2 alleles, and the genotyping of a panel of 134 almond cultivars. Overall, our results provide further insights on the understanding of the almond cultivation history. In addition, molecular marker assays and genotypic data presented in this study are expected to be of major interest for the conduction of almond breeding programs, which often need to select sweet kernel individuals in segregant populations.

Genotyping-by-Sequencing Defines Genetic Structure within the "Acquaviva" Red Onion Landrace.

Genetic structure and distinctive features of landraces, such as adaptability to local agro-ecosystems and specific qualitative profiles, can be substantially altered by the massive introduction of allochthonous germplasm. The landrace known as "Cipolla rossa di Acquaviva" (Acquaviva red onion, further referred to as ARO) is traditionally cultivated and propagated in a small area of the Apulia region (southern Italy). However, the recent rise of its market value and cultivation area is possibly causing genetic contamination with foreign propagating material. In this work, genotyping-by-sequencing (GBS) was used to characterize genetic variation of seven onion populations commercialized as ARO, as well as one population of the landrace "Montoro" (M), which is phenotypically similar, but originates from another cultivation area and displays different qualitative features. A panel of 5011 SNP markers was used to perform parametric and non-parametric genetic structure analyses, which supported the hypothesis of genetic contamination of germplasm commercialized as ARO with a gene pool including the M landrace. Four ARO populations formed a core genetic group, homogeneous and clearly distinct from the other ARO and M populations. Conversely, the remaining three ARO populations did not display significant differences with the M population. A set of private alleles for the ARO core genetic group was identified, indicating the possibility to trace the ARO landrace by means of a SNP-based molecular barcode. Overall, the results of this study provide a framework for further breeding activities and the traceability of the ARO landrace.

Merging genotyping-by-sequencing data from two ex situ collections provides insights on the pea evolutionary history.

Pea (Pisum sativum L. subsp. sativum) is one of the oldest domesticated species and a widely cultivated legume. In this study, we combined next generation sequencing (NGS) data referring to two genotyping-by-sequencing (GBS) libraries, each one prepared from a different Pisum germplasm collection. The selection of single nucleotide polymorphism (SNP) loci called in both germplasm collections caused some loss of information; however, this did not prevent the obtainment of one of the largest datasets ever used to explore pea biodiversity, consisting of 652 accessions and 22 127 markers. The analysis of population structure reflected genetic variation based on geographic patterns and allowed the definition of a model for the expansion of pea cultivation from the domestication centre to other regions of the world. In genetically distinct populations, the average decay of linkage disequilibrium (LD) ranged from a few bases to hundreds of kilobases, thus indicating different evolutionary histories leading to their diversification. Genome-wide scans resulted in the identification of putative selective sweeps associated with domestication and breeding, including genes known to regulate shoot branching, cotyledon colour and resistance to lodging, and the correct mapping of two Mendelian genes. In addition to providing information of major interest for fundamental and applied research on pea, our work describes the first successful example of integration of different GBS datasets generated from ex situ collections - a process of potential interest for a variety of purposes, including conservation genetics, genome-wide association studies, and breeding.

Screening of Olive Biodiversity Defines Genotypes Potentially Resistant to Xylella fastidiosa.

The recent outbreak of the Olive Quick Decline Syndrome (OQDS), caused by Xylella fastidiosa subsp. pauca (Xf), is dramatically altering ecosystem services in the peninsula of Salento (Apulia Region, southeastern Italy). Here we report the accomplishment of several exploratory missions in the Salento area, resulting in the identification of thirty paucisymptomatic or asymptomatic plants in olive orchards severely affected by the OQDS. The genetic profiles of such putatively resistant plants (PRPs), assessed by a selection of ten simple sequence repeat (SSR) markers, were compared with those of 141 Mediterranean cultivars. Most (23) PRPs formed a genetic cluster (K1) with 22 Italian cultivars, including 'Leccino' and 'FS17', previously reported as resistant to Xf. The remaining PRPs displayed relatedness with genetically differentiated germplasm, including a cluster of Tunisian cultivars. Markedly lower colonization levels were observed in PRPs of the cluster K1 with respect to control plants. Field evaluation of four cultivars related to PRPs allowed the definition of partial resistance in the genotypes 'Frantoio' and 'Nocellara Messinese'. Some of the PRPs identified in this study might be exploited in cultivation, or as parental clones of breeding programs. In addition, our results indicate the possibility to characterize resistance to Xf in cultivars genetically related to PRPs.

Almond diversity and homozygosity define structure, kinship, inbreeding, and linkage disequilibrium in cultivated germplasm, and reveal genomic associations with nut and seed weight.

Almond [Prunus dulcis Miller (D.A. Webb)] is the main tree nut species worldwide. Here, genotyping-by-sequencing (GBS) was applied to 149 almond cultivars from the ex situ collections of the Italian Council for Agricultural Research (CREA) and the Spanish National Research Council (CSIC), leading to the detection of 93,119 single-nucleotide polymorphisms (SNPs). The study of population structure outlined four distinct genetic groups and highlighted diversification between the Mediterranean and Californian gene pools. Data on SNP diversity and runs of homozygosity (ROHs) allowed the definition of kinship, inbreeding, and linkage disequilibrium (LD) decay in almond cultivated germplasm. Four-year phenotypic observations, gathered on 98 cultivars of the CREA collection, were used to perform a genome-wide association study (GWAS) and, for the first time in a crop species, homozygosity mapping (HM), resulting in the identification of genomic associations with nut, shell, and seed weight. Both GWAS and HM suggested that loci controlling nut and seed weight are mostly independent. Overall, this study provides insights on the almond cultivation history and delivers information of major interest for almond genetics and breeding. In a broader perspective, our results encourage the use of ROHs in crop science to estimate inbreeding, choose parental combinations minimizing the risk of inbreeding depression, and identify genomic footprints of selection for specific traits.

Data on the proximate composition, bioactive compounds, physicochemical and functional properties of a collection of faba beans (Vicia faba L.) and lentils (Lens culinaris Medik.).

This dataset is referred to a collection of 41 faba bean (Vicia faba L.) and 15 lentil (Lens culinaris Medik.) accessions from the ex situ repository of the Institute of Biosciences and Bioresources of the Italian National Research Council (CNR-IBBR). All the accessions were grown at the experimental farm "P. Martucci" of the University of Bari "Aldo Moro" (41°01'22.1'' N 16°54'21.0'' E) during the growing season 2017-2018, according to a randomized block design with two replicates, each constituted by 10 individual plants. The dataset reports raw and elaborated analytical data determined on the flour produced from individual accessions, concerning proximate composition, bioactive compounds, antioxidant activity, fatty acid composition, and physicochemical and functional properties. Elaborated data might be used to understand the compositional variability within the species and, together with raw data, to highlight peculiar accessions characterized by valuable nutritional and/or technological attitude useful in research institutions and food industries. Furthermore, the data can be used for genetic studies aimed at identifying genomic regions underlying nutritional and technological traits.

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

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

Assessment of Genetic Diversity of the "Acquaviva Red Onion" (Allium cepa L.) Apulian Landrace.

Onion (Allium cepa L.) is the second most important vegetable crop worldwide and is widely appreciated for its health benefits. Despite its significant economic importance and its value as functional food, onion has been poorly investigated with respect to its genetic diversity. Herein, we surveyed the genetic variation in the "Acquaviva red onion" (ARO), a landrace with a century-old history of cultivation in a small town in the province of Bari (Apulia, Southern of Italy). A set of 11 microsatellite markers were used to explore the genetic variation in a germplasm collection consisting of 13 ARO populations and three common commercial types. Analyses of genetic structure with parametric and non-parametric methods highlighted that the ARO represents a well-defined gene pool, clearly distinct from the Tropea and Montoro landraces with which it is often mistaken. In order to provide a description of bulbs, usually used for fresh consumption, soluble solid content and pungency were evaluated, showing higher sweetness in the ARO with respect to the two above mentioned landraces. Overall, the present study is useful for the future valorization of the ARO, which could be promoted through quality labels which could contribute to limit commercial frauds and improve the income of smallholders.

Data on the chemical composition, bioactive compounds, fatty acid composition, physico-chemical and functional properties of a global chickpea collection.

The data article refers to the paper "Nutritional, physico-chemical and functional characterization of a global chickpea collection" [1]. The data are referred to a germplasm collection of 57 chickpea accessions from the ex situ repositories of the United States Department of Agriculture (USDA), the Department of Plant, Soil and Food Science of the University of Bari, Italy (DiSSPA), and the Institute of Biosciences and Bioresources of the Italian National Research Council (CNR-IBBR). Thirty-six accessions, belonging to desi and kabuli types, were representative of the geographic distribution of chickpea global cultivation, whereas twenty-one accessions, referable to the Apulian black type, derived from different area of the Apulian region, south of Italy. All the accessions were grown at the experimental farm "P. Martucci" of the University of Bari "Aldo Moro" (41°01'22.1″ N 16°54'21.0″ E) during the growing season 2017-2018, according to a randomized block design with two replicates, each replicate formed by 30 individual plants. This article reports the data of the proximate composition, the total bioactive compounds content, the fatty acid composition and the physico-chemical and functional properties of chickpea flour. Information provided in this article can be used by food industry to develop chickpea-based foods and by geneticists for studies of association mapping aimed at the identification of genomic regions controlling the nutritional and technological traits.

Genotyping by Sequencing of Cultivated Lentil (Lens culinaris Medik.) Highlights Population Structure in the Mediterranean Gene Pool Associated With Geographic Patterns and Phenotypic Variables.

Cultivated lentil (Lens culinaris Medik.) is one of the oldest domesticated crops and one of the most important grain legumes worldwide. The Mediterranean Basin holds large part of lentil biodiversity; however, no genetic structure was defined within the Mediterranean gene pool. In this study, we used high-throughput genotyping by sequencing to resolve the genetic structure of the Mediterranean ex situ lentil collection held at the Italian National Research Council. Sequencing of a 188-plex genotyping-by-sequencing library and bioinformatics treatment of data yielded 6,693 single nucleotide polymorphisms. Analysis of nonredundant genotypes with nonparametric and parametric methods highlighted the occurrence of five highly differentiated genetic clusters. Clustering could be related to geographic patterns and phenotypic traits, indicating that post-domestication routes introducing cultivation in Mediterranean countries and selection were major forces shaping lentil population structure. The estimation of the fixation index FST at individual single nucleotide polymorphism loci allowed the identification of distinctive alleles across clusters, suggesting the possibility to set up molecular keys for the assignment of lentil germplasm to specific genetic groups. Finally, significant associations between markers and phenotypic data were identified. Overall, the results of this study are of major importance for lentil conservation genetics and breeding and provide insights on the lentil evolutionary history.

Evolutionary conservation of MLO gene promoter signatures.

BACKGROUND: Powdery mildew (PM) is a widespread fungal disease of plants in temperate climates, causing significant economic losses in agricultural settings. Specific homologs of the MLO gene family are PM susceptibility factors, as their loss-of function results in durable PM resistance (mlo resistance) in several plant species. The role of MLO susceptibility genes in plant-pathogen interactions is still elusive, however it is known that they are strongly upregulated following PM infection. RESULTS: In this study, we investigated the structure of 414 Putative Promoter Regions (PPRs) of MLO genes and highlighted motif and regulatory element patterns related to genomic relationships among species and phylogenetic distance among homologs. A TC box-like motif and a thymine-rich motif were found to be overrepresented in MLO genes transcriptionally upregulated upon infection with PM fungi. As proof of concept, we showed that the expression of a melon (Cucumis melo L.) gene enriched for the motifs above mentioned was strongly upregulated upon infection with the PM fungus Podosphaera xanthii. CONCLUSION: While identifying a candidate MLO susceptibility gene in melon, this study provides insight on the transcriptional control of MLO genes and indicates diagnostic features useful to identify MLO susceptibility genes across species affected by the PM disease.

Carbon dioxide in office diagnostic hysteroscopy: An open question. A multicenter randomized trial on 1982 procedures.

OBJECTIVE: to compare carbon dioxide and saline solution distension in diagnostic hysteroscopies with regards to patient discomfort and procedural time and in accordance with the instrument diameter (5 mm vs 3.5 mm). The secondary outcome was to evaluate the role of the patient age and the obstetrical history on perception of pain and procedural time. STUDY DESIGN: This is a prospective multicenter randomized study including 1982 patients that underwent office diagnostic hysteroscopy in: Policlinico Abano Terme, Università Cattolica del Sacro Cuore in Rome and Ente Ecclesiastico Ospedale Generale Regionale "F. Miulli" in Acquaviva delle Fonti. They were firstly randomized according to distension medium and secondly according to instrument diameter. Pain perception after the procedure was assessed by VAS (Visual Analogue Scale) score and procedural time was registered. Mann-Whitney U test was used to compare data. RESULTS: Lower pain score and procedural time were recorded with the employment of Carbon Dioxide (p < 001). Patient discomfort and procedural time were significantly influenced by the instrument diameter independent of the distension medium used, though in the subgroup where gas was employed VAS score after 1 min (VAS1) resulted lower compared to saline solution in both the traditional and mini-hysteroscopy procedures (p < 001). CONCLUSION: Carbon dioxide and saline solution are both suitable distension media for outpatient diagnostic hysteroscopy; nonetheless, carbon dioxide confers advantages in terms of pain perception and procedural time.

Synteny-Based Development of CAPS Markers Linked to the Sweet kernel LOCUS, Controlling Amygdalin Accumulation in Almond (Prunus dulcis (Mill.) D.A.Webb).

The bitterness and toxicity of wild-type seeds of Prunoideae is due to the cyanogenic glucoside amygdalin. In cultivated almond (Prunus dulcis (Mill.) D.A. Webb), a dominant mutation at the Sk locus prevents amygdalin accumulation and thus results in edible sweet kernels. Here, we exploited sequence similarity and synteny between the genomes of almond and peach (Prunus persica (L.) Batsch) to identify cleaved amplified polymorphic sequence (CAPS) molecular markers linked to the Sk locus. A segregant F1 population was used to map these markers on the Sk genomic region, together with Sk-linked simple sequence repeat (SSR) markers previously described. Molecular fingerprinting of a cultivar collection indicated the possibility to use CAPS polymorphisms identified in this study in breeding programs arising from different parental combinations. Overall, we highlight a set of codominant markers useful for early selection of sweet kernel genotypes, an aspect of primary importance in almond breeding. In addition, by showing collinearity between the physical map of peach and the genetic map of almond with respect to the Sk genomic region, we provide valuable information for further marker development and Sk positional cloning.

The Effect of Tomatine on Gene Expression and Cell Monolayer Integrity in Caco-2.

More understanding of the risk-benefit effect of the glycoalkaloid tomatine is required to be able to estimate the role it might play in our diet. In this work, we focused on effects towards intestinal epithelial cells based on a Caco-2 model in order to analyze the influence on the cell monolayer integrity and on the expression levels of genes involved in cholesterol/sterol biosynthesis (LDLR), lipid metabolism (NR2F2), glucose and amino acid uptake (SGLT1, PAT1), cell cycle (PCNA, CDKN1A), apoptosis (CASP-3, BMF, KLF6), tight junctions (CLDN4, OCLN2) and cytokine-mediated signaling (IL-8, IL1ß, TSLP, TNF-α). Furthermore, since the bioactivity of the compound might vary in the presence of a food matrix and following digestion, the influence of both pure tomatine and in vitro digested tomatine with and without tomato fruit matrix was studied. The obtained results suggested that concentrations <20 µg/mL of tomatine, either undigested or in vitro digested, do not compromise the viability of Caco-2 cells and stimulate cytokine expression. This effect of tomatine, in vitro digested tomatine or in vitro digested tomatine with tomato matrix differs slightly, probably due to variations of bioactivity or bioavailability of the tomatine. The results lead to the hypothesis that tomatine acts as hormetic compound that can induce beneficial or risk toxic effects whether used in low or high dose.

Functional Characterization of a Syntaxin Involved in Tomato (Solanum lycopersicum) Resistance against Powdery Mildew.

Specific syntaxins, such as Arabidopsis AtPEN1 and its barley ortholog ROR2, play a major role in plant defense against powdery mildews. Indeed, the impairment of these genes results in increased fungal penetration in both host and non-host interactions. In this study, a genome-wide survey allowed the identification of 21 tomato syntaxins. Two of them, named SlPEN1a and SlPEN1b, are closely related to AtPEN1. RNAi-based silencing of SlPEN1a in a tomato line carrying a loss-of-function mutation of the susceptibility gene SlMLO1 led to compromised resistance toward the tomato powdery mildew fungus Oidium neolycopersici. Moreover, it resulted in a significant increase in the penetration rate of the non-adapted powdery mildew fungus Blumeria graminis f. sp. hordei. Codon-based evolutionary analysis and multiple alignments allowed the detection of amino acid residues that are under purifying selection and are specifically conserved in syntaxins involved in plant-powdery mildew interactions. Our findings provide both insights on the evolution of syntaxins and information about their function which is of interest for future studies on plant-pathogen interactions and tomato breeding.

A Distinct Genetic Cluster in Cultivated Chickpea as Revealed by Genome-wide Marker Discovery and Genotyping.

The accurate description of plant biodiversity is of utmost importance to efficiently address efforts in conservation genetics and breeding. Herein, we report the successful application of a genotyping-by-sequencing (GBS) approach in chickpea ( L.), resulting in the characterization of a cultivated germplasm collection with 3187 high-quality single nucleotide polymorphism (SNP) markers. Genetic structure inference, principal component analysis, and hierarchical clustering all indicated the identification of a genetic cluster corresponding to black-seeded genotypes traditionally cultivated in Southern Italy. Remarkably, this cluster was clearly distinct at both genetic and phenotypic levels from germplasm groups reflecting commercial chickpea classification into and seed types. Fixation index estimates for individual polymorphisms pointed out loci and genomic regions that might be of significance for the diversification of agronomic and commercial traits. Overall, our findings provide information on genetic relationships within cultivated chickpea and highlight a gene pool of great interest for the scientific community and chickpea breeding, which is limited by the low genetic diversity available in the primary gene pool.

Functional characterization of the powdery mildew susceptibility gene SmMLO1 in eggplant (Solanum melongena L.).

Eggplant (Solanum melongena L.) is one of the most important vegetables among the Solanaceae and can be a host to fungal species causing powdery mildew (PM) disease. Specific homologs of the plant Mildew Locus O (MLO) gene family are PM susceptibility factors, as their loss of function results in a recessive form of resistance known as mlo resistance. In a previous work, we isolated the eggplant MLO homolog SmMLO1. SmMLO1 is closely related to MLO susceptibility genes characterized in other plant species. However, it displays a peculiar non-synonymous substitution that leads to a T â†’ M amino acid change at protein position 422, in correspondence of the MLO calmodulin-binding domain. In this study, we performed the functional characterization of SmMLO1. Transgenic overexpression of SmMLO1 in a tomato mlo mutant compromised resistance to the tomato PM pathogen Oidium neolycopersici, thus indicating that SmMLO1 is a PM susceptibility factor in eggplant. PM susceptibility was also restored by the transgenic expression of a synthetic gene, named s-SmMLO1, encoding a protein identical to SmMLO1, except for the presence of T at position 422. This indicates that the T â†’ M polymorphism does not affect the protein role as PM susceptibility factor. Overall, the results of this work are of interest for the functional characterization of MLO proteins and the introduction of PM resistance in eggplant using reverse genetics.

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

BACKGROUND: Melon (Cucumis melo L.) is one of the most important horticultural species, which includes several taxonomic groups. With the advent of next-generation sequencing, single nucleotide polymorphism (SNP) markers are widely used in the study of genetic diversity and genomics. RESULTS: We report the first successful application of genotyping-by-sequencing (GBS) technology in melon. We detected 25,422 SNPs by the analysis of 72 accessions collected in Apulia, a secondary centre of diversity in Southern Italy. Analyses of genetic structure, principal components, and hierarchical clustering support the identification of three distinct subpopulations. One of them includes accessions known with the folk name of 'carosello', referable to the chate taxonomic group. This is one of the oldest domesticated forms of C. melo, once widespread in Europe and now exposed to the risk of genetic erosion. The second subpopulation contains landraces of 'barattiere', a regional vegetable production that was never characterized at the DNA level and we show was erroneously considered another form of chate melon. The third subpopulation includes genotypes of winter melon (C. melo var. inodorus). Genetic analysis within each subpopulation revealed patterns of diversity associated with fruit phenotype and geographical origin. We used SNP data to describe, for each subpopulation, the average linkage disequilibrium (LD) decay, and to highlight genomic regions possibly resulting from directional selection and associated with phenotypic variation. CONCLUSIONS: We used GBS to characterize patterns of genetic diversity and genomic features within C. melo. We provide useful information to preserve endangered gene pools and to guide the use of germplasm in breeding. Finally, our findings lay a foundation for molecular breeding approaches and the identification of genes underlying phenotypic traits.