- 1Department of Agricultural Sciences, University of Naples Federico II, Portici, Italy
- 2Department of Soil, Plant and Food Sciences, University of Bari Aldo Moro, Bari, Italy
- 3Institute of Biosciences and Bioresources (CNR-IBBR), Bari, Italy
- 4Support Unit Bari, Institute for Sustainable Plant Protection, National Research Council of Italy (CNR), Bari, Italy
- 5SINAGRI S.r.l., Spin Off of the University of Bari Aldo Moro, Bari, Italy
Sexual reproduction has contributed to a significant degree of variability in cultivated grapevine populations. However, the additional influence of spontaneous somatic mutations has played a pivotal role in shaping the diverse landscape of grapevine agrobiodiversity. These naturally occurring selections, termed 'clones,' represent a vast reservoir of potentially valuable traits and alleles that hold promise for enhancing grape quality and bolstering plant resilience against environmental and biotic challenges. Despite their potential, many of these clones remain largely untapped.In light of this context, this study aims to delve into the population structure, genetic diversity, and distinctive genetic loci within a collection of 138 clones derived from six Campanian and Apulian grapevine varieties, known for their desirable attributes in viticulture and winemaking. Employing two reduced representation sequencing methods, we extracted Single-Nucleotide Polymorphism (SNP) markers. Population structure analysis and fixation index (FST) calculations were conducted both between populations and at individual loci. Notably, varieties originating from the same geographical region exhibited pronounced genetic similarity.The resulting SNP dataset facilitated the identification of approximately two hundred loci featuring divergent markers (FST ≥ 0.80) within annotated exons. Several of these loci exhibited associations with essential traits like phenotypic adaptability and environmental responsiveness, offering compelling opportunities for grapevine breeding initiatives. By shedding light on the genetic variability inherent in these treasured traditional grapevines, our study contributes to the broader understanding of their potential. Importantly, it underscores the urgency of preserving and characterizing these valuable genetic resources to safeguard their intra-varietal diversity and foster future advancements in grapevine cultivation.
Introduction
According to the data from the Food and Agriculture Organization (FAOSTAT, 2023), Italy holds a significant position among the leading grape-producing nations. This prominence can be attributed to several unique factors, all encapsulated in the concept of “terroir”. The Italian wine industry is variegated due to the distinct characteristics of the regions, including their physical, climatic, and cultural aspects. However, the biggest contributing factor is the availability and use of a wide range of grapevine varieties cultivated in different part of the Country (Sardaro et al., 2017; Pomarici et al., 2021).
To date, based on the Italian Vitis database (talian Vitis Database, 2023; D’Onofrio and Scalabrelli, 2010) and the National Register of Grapevine Varieties (http://catalogoviti.politicheagricole.it) more than 1,100 varieties have been collected and described; however, the actual number of grape varieties cultivated and kept in germplasm collections in Italy is much higher and cannot be estimated with certainty. The invaluable germplasm of Italian grapes includes many autochthonous varieties. Although sexual reproduction has led to a great variability among varieties (i.e., inter-varietal diversity), spontaneous somatic mutations have further contributed to the kaleidoscope of grapevine diversity with the so-called clones (Salmaso et al., 2005). As a result, autochthonous varieties are made up of groups of different clones with genotypic, morphological, and physiological characteristics slightly different from those of the original mother plant (i.e., intra-varietal diversity). An effective approach to explore the inter- and intra-varietal diversity relies on the use of molecular markers. Microsatellites are the most used markers in grapevine, and a standard set of nine loci is currently used for international cataloguing (This et al., 2004; Laucou et al., 2011; Vitis International Variety Catalogue VIVC, 2023). Single nucleotide polymorphisms (SNPs) are another class of markers widely used for grape genotyping, which promise a higher map resolution, higher throughput, lower cost, and lower error rate than microsatellites (Vezzulli et al., 2008; Kumar et al., 2012; Vishwakarma et al., 2022). Analyses of SNPs variations can only be conducted in a small fraction of the genome with Reduced Representation Sequencing (RRS) (Pavan et al., 2020). The most popular techniques use restriction enzymes to prepare DNA for sequencing (Restriction-site-Associated DNA: RAD). Many methods based on the RAD approach have been developed in recent years, differing in the number of enzymes used or in additional library preparation steps (Andrews et al., 2016). In plant, the most used are GBS (genotyping-by-sequencing), which generally relies on a single restriction enzyme (Elshire et al., 2011) and ddRADseq (double-digest RAD sequencing), in which DNA is digested with two restriction enzymes (Peterson et al., 2012; O'Leary et al., 2018). Both GBS and ddRADseq techniques have been successfully used in Vitis genotyping to study intra-varietal diversity (Miazzi et al., 2020; Calderón et al., 2021; D’Onofrio et al., 2021).
Southern Italy is recognised as the oldest wine-growing area in Italy where many traditional and autochthonous varieties are still cultivated (De Lorenzis et al., 2019). They represent a vast reservoir of traits/alleles that could be useful for improving the quality of grape as well as plant tolerance to environmental and biotic stresses. However, most of them are still underexploited. In this framework, the present study is aimed to analyse population structure, genetic diversity, and divergent loci in a panel of 138 clones belonging to six Campanian and Apulian grapevine varieties characterised by attractive traits from a viticultural and oenological point of view. To genotype this plant material, we used ddRADseq and GBS techniques on Campanian and Apulian clones, respectively. The analysis revealed the complex genetic structure of the varieties under investigation at the clonal level, their relationships, and the presence of divergent SNP loci within genes involved in grape phenology and adaptation to the environment.
Materials and methods
Plant material
A total of 138 grapevine clones belonging to 6 autochthonous varieties of Southern Italy were sampled in Campania and in Apulia (Supplementary Table S1). They are traditionally grown in environments with different climate conditions: plains with hot, dry summers and mild winters for the Apulian varieties; hills with warm, humid summers and cold winters for the Campanian varieties. The Campanian clones include two red-berry varieties, camaiola (formerly barbera del sannio) (n. 20) and aglianico lasco (n.64), and the white-berry greco b. (n. 24). aglianico lasco has loose clusters able to reduce the development of mold during ripening. camaiola has intense ruby red colour berries suitable for violet hues wines, and greco b. fruits give rise to high quality wines appreciated for the complex aroma profiles (Antonacci et al., 2017). The Apulian clones belong to two red-berry varieties, nero di troia (synonym uva di troia; n. 16) and malvasia nera di brindisi (hereafter referred as malvasia nera; n. 8), and to the white-berry minutolo b. (n. 6). nero di troia is a high-quality, flavourful variety cultivated mainly in northern Apulia to make an ever-increasing number of monovarietal wines or blends. malvasia nera is a variety of Greek origin, widely cultivated in the southern part of Apulia, which produces fruity wines often together with the local variety negramaro. minutolo b. is an aromatic white variety with relatively small and loose clusters, which in the past was widely grown in the specific area of Apulia region called Itria Valley. The map of the sampling sites is shown in Figure 1. No major morphological and physiological differences were found among the clones of the six varieties analysed. In addition, four native (aglianico del vulture, aglianico taurasi, aglianico del taburno and sangiovese) and two international (chardonnay and merlot) varieties were used for microsatellite (SSR) analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1201287/fpls-14-1201287-HTML/image_m/fpls-14-1201287-g001.jpg)
Figure 1 Map of southern Italy with the sampling sites for the Campanian and Apulian clones of the six varieties under investigation.
DNA isolation
As for Campanian clones, the total DNA was extracted from 1 gr of young leaf tissue following the method by Japelaghi et al. (2011) with a few modifications, such as a ratio of 0.2 g tissue/1 mL extraction buffer, an increase in the concentration of soluble PVP from 2% (w/v) to 4% (w/v) in the extraction buffer, and an additional washing step with 70% (v/v) ethanol (EtOH) before elution.
As for Apulian samples, the DNA was extracted following the procedure described in Spadoni et al. (2019). For all samples, the yield, quantity, and quality of extracted DNA were estimated using a NanoDrop 1000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA), a 0.8% agarose gel, and a fluorimeter (Qubit 2.0, Thermo Fisher Scientific) with the Qubit dsDNA BR Assay kit (Invitrogen, CA, USA).
SNP calling
Campanian genetic material was genotyped using ddRADseq (Peterson et al., 2012). The ddRADseq tags were aligned to the reference genome GCA_000003745.2 (Pinot Noir, PN40024) using BWA-MEM (Li and Durbin, 2009). The gstacks and populations (with options -R = 0.75 and –max-obs-het = 0.8) utilities included in Stacks v2.53 (Catchen et al., 2013) were used for SNP calling and for applying filtering options. Genotyping-by-sequencing and SNP calling on Apulian samples were performed by Elshire group Ltd. (https://www.elshiregroup.co.nz/) as described in Miazzi et al. (2020), using the same reference genome as above.
SNP filtering and statistics
ddRADseq and GBS produced two Variant Call Format (VCF) files respectively for Campanian and Apulian varieties, which underwent SNP filtering using VCFtools v. 0.1.16 (Danecek et al., 2011). For these two datasets, SNP markers with minor allele frequency (MAF) < 5%, minimum site count < 15%, and a minimum depth of coverage of 5 were filtered out. Instead, for the dataset of each variety, SNP markers with a call rate of 100%, minor allele frequency (MAF) < 5%, and a minimum depth of coverage of 5 were filtered out. The choice of not having missing data was necessary for the intra-varietal analysis, as the missing data affected the Identity-by-State (IBS) values. VCFtools were also used to generate various statistics on the datasets under investigation and to add gene annotations to VCF files.
Identity-by-State and linkage disequilibrium analysis
PLINK v.1.90 (Purcell et al., 2007) was used to obtain the IBS distance matrix for each variety. Duplicate individuals were identified by setting IBS value ≥ 0.99 as the threshold. Only one individual was retained among duplicates and used for downstream analyses.
For each dataset, linkage disequilibrium (LD) analysis was performed using the Golden Helix SNP and Variation Suite (SVS) v.8.8.3 (Golden Helix Inc.). The “Nonlinear Regression” function was used to plot the r2 values against the physical distance of the markers, and the LD decay was estimated at the critical level of r2 = 0.20.
Merging SNP data points
After removing the duplicates, bcftools +fixref (Li, 2011) was used to fix the inconsistencies in reference (REF) – alternative (ALT) alleles between the two VCF files. Then the vcf-merge utility of VCFtools (Danecek et al., 2011) was used. VCFtools (Danecek et al., 2011) and PLINK (Purcell et al., 2007) were used to filter out individuals with SNP markers with a missingness per individual > 85%, MAF < 5% and missingness per marker >20%. The kinship coefficient was calculated using the –relatedness2 option in VCFtools. The dataset was LD pruned (r2 = 0.50) using SVS.
Genetic diversity and population structure
Allele frequency and ancestry estimation were performed on filtered and pruned dataset using ADMIXTURE v. 1.3.0 (Alexander and Lange, 2011), with 10-fold cross validation (CV) for sub-populations (K) ranging from 1 to 10, and 1,000 bootstrap replicates. CV scores were used to estimate the optimal K value. A membership coefficient (qi) >0.55 was used to separate the individuals into sub-populations. Discriminant Analysis of Principal Components (DAPC) was used as an exploratory data analysis to investigate population structure (Jombart et al., 2010; Jombart and Ahmed, 2011). The optimal number of principal components (PCs) to maintain was determined using a value ≥ 1:200. The optimal number of k-means was determined using the Bayesian information criterion (BIC) as a statistical measure of goodness-of-fit. Principal component analysis (PCA) was performed using SVS and Neighbor-joining trees (bootstrap = 1,000) were built with MEGA v.11. FigTree v.1.4 was used to visualise trees (Rambaut, 2014; http://tree.bio.ed.ac.uk/software/figtree; Tamura et al., 2021). Haplotype diversity and analysis of molecular variance (AMOVA) was conducted using GenALEx version 6.5 (Peakall and Smouse, 2006) to calculate variance components and their levels of statistical significance for variation between and within populations.
Microsatellite analysis
Total genomic DNA was extracted from young leaves of Aglianico Lasco clones using the Qiagen Plant DNeasy Maxi Kit (Qiagen, Valencia, CA, USA), following the manufacturer’s procedure. Microsatellite analysis was carried out with seven nuclear markers (VVMD27, VVlb01, Vvln16, VVIp60, VVIq52, VrZAG79, and VVS2) from Laucou et al. (2011). PCR amplification, size calibration, alleles detection was performed as reported by Villano et al. (2014). Validation of results was performed with three biological and technical replicates. Allele sizes were normalised using SSR data reported in the Vitis International Variety Catalog (http://www.vivc.de/). The phylogenetic clustering tree was constructed with the neighbour joining method using MEGA X (Tamura et al., 2021). The robustness of the clusters was tested by bootstrap resampling (n = 1,000) with the Darwin software (Perrier et al., 2003).
Signature of divergence
Per-site FSTcalculation was performed on the filtered dataset using SVS and comparing pairs of populations each time as described in Taranto et al. (2020). Only divergent SNP markers (FST ≥ 0.80) within annotated exons were considered for the search for candidate genes. Grapevine gene annotations were retrieved at https://urgi.versailles.inra.fr/Species/Vitis/Annotations.
Results
Intra-varietal diversity: genome-wide SNPs discovery and the extent of linkage disequilibrium
A total of 114,465 and 180,840 SNP markers were scored in the 108 Campanian and in the 30 Apulian clones, respectively. After filtering, 59,149 and 17,296 SNPs were retained in the Campanian and Apulian datasets, respectively (Table 1; Supplementary Figures 1A, B). Transitions (Ti) were more abundant (81%) in the Campanian dataset than in Apulian (63.33%), with a Ti/Tv ratio of 4.22 and 1.72, respectively. Taking advantage of the genomic coordinates of the reference gene models, 15,611 (26.4%) and 8,877 (51.3%) SNPs fell within annotated exons for Campanian and Apulian datasets, respectively. The number of SNPs before and after filtering and the statistics of high-quality SNPs are shown in Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1201287/fpls-14-1201287-HTML/image_m/fpls-14-1201287-t001.jpg)
Table 1 Number of SNPs before and after filtering procedure and statistics of high-quality SNPs calculated for each dataset.
The clones of each variety were separated into six VCF files and each dataset was subjected to the filtering procedure. After filtering, 9,963 and 9,461 (on average, avg) SNPs were retained in the Campanian and Apulian datasets, respectively.
Of these SNPs, an average of 3,086 fell within annotated genes in the Campanian datasets and 4,544 in the Apulian ones. Malvasia showed the greatest number of SNPs within genes (N = 7,644), while Greco B. showed the least (N = 2,966). In all the varieties there were more transitions (Ts) (7,079 avg) than transversions (Tv) (2,572 avg) and the average Ts/Tv ratio was higher in the Campanian varieties (4.60) than in the Apulian ones (1.80) (Table 1).
The high-quality SNP datasets were then used to calculate the pairwise IBS distance between clones of each variety. Clone pairs with an allele sharing rate > 99% were considered duplicates and only one of them was retained for downstream analysis. No duplicate clones were found in Camaiola and Apulian germplasm, whereas 47 and 13 duplicates were removed in Aglianico Lasco and Greco B., respectively, resulting in a final number of 78 clones (Supplementary Table S1). The LD decay was estimated for each dataset. At r2 = 0.20, Aglianico Lasco showed the slower LD decay of 254.3 kb. In contrast, Camaiola, Greco B., Minutolo, Malvasia and Nero di Troia were characterised by very rapid LD decay (< 0.093kb) (Supplementary Figure S2).
Inter-varietal diversity: population structure and FST
The 78 clones retained after IBS analysis were merged into a single dataset of 2,235 SNP markers. After filtering, 1,091 high-quality SNPs and 72 individuals were retained (Supplementary Table S1). The distribution of SNPs along the 19 chromosomes is shown in Supplementary Figure S3. The dataset was further pruned for LD (r2 = 0.50) resulting in 500 SNPs, out of which 280 fall within genes. Population diversity was analysed using different approaches. The population structure indicated K = 6 and K = 7 as the best numbers of sub-populations, based on cross-validation error (Supplementary Figure S4). At K7 (Figure 2A), all the varieties grouped their own clones into separate sub-populations, with the exception of Aglianico Lasco which was separated into two sub-populations (henceforth AL1 and AL2), comprising nine and four clones, respectively. Individuals belonging to the sub-populations Minutolo, Camaiola, AL2, Greco B. and Nero di Troia shared a co-ancestry coefficient (qi) > of 0.99, whereas the clones of AL1 and Malvasia had greater variability with qi ranging from 0.99 to 0.60 and from 0.99 to 0.57, respectively. The clones AL1-62, AL1-1, Malvasia-100 and Malvasia-101 resulted admixed. In addition, clone AL1-18 was grouped with clones belonging to Camaiola instead of Aglianico Lasco. DAPC analysis was then performed which confirmed the existence of seven clusters (Supplementary Figure S5) albeit with a slightly different clone distribution (Figure 2B). The Nero Di Troia and Camaiola clones were alone in the first and second quadrant, respectively. The Greco B. clones were grouped in the centre of the plot, while those of Aglianico Lasco (AL1 and AL2) were scattered in the fourth quadrant together with Malvasia Nera and Minutolo. Inter-varietal structuring within clones was supported by the neighbour-joining tree (Figure 2C). The phylogenetic analysis highlighted the geographical differentiation of the clones. The Apulian and Campanian clones formed well distinct clusters within which the clones belonging to each variety formed well-supported subclades. Also, the clustering of Aglianico Lasco confirmed the results of ADMIXTURE. These clones were divided into two subclusters, except for AL1-1, AL1-18, and AL1-62 which fell respectively in AL2, Camaiola and Greco B.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1201287/fpls-14-1201287-HTML/image_m/fpls-14-1201287-g002.jpg)
Figure 2 (A) Bar-plot describing the population structure estimated by ADMIXTURE. Each bar is divided into K coloured segments each representing the proportion of ancestry (qi) in each individual. (B) PCA plot of 72 grapevine clones obtained using 1,091 SNPs. (C) Neighbour-Joining tree. (D) Heatmap of fixation indices (FST) between subpopulations.
The analysis of haplotype diversity corroborated these results, in fact Camaiola and Greco B. showed lower values (h =0.041 and 0.058, respectively) than AL1 (h = 0.427), AL2 (h = 0.456), Minutolo (h = 0.173), Nero Di Troia (h = 0.236) and Malvasia Nera (h = 0.350). This was supported by the within-population sum of squares (SSWP) from the AMOVA analysis, where SSWP was 349.636 and 423.789 for Greco B. and Camaiola, respectively. The highest SSWP was found in AL1 (2097.001). Additionally, the AMOVA revealed high genetic variability within (41%) and between (59%) populations. Further investigation of genetic distances was performed using the genetic fixation index (FST) between sub-populations as emerged from clustering of ADMIXTURE and DAPC (Figure 2D). The heatmap in Figure 2D shows that all varieties are genetically distinct. The greatest distance was found between Camaiola and all the other populations (FST mean value > of 0.52), with the highest value detected between Camaiola and Minutolo/Nero di troia (FST =0.58), followed by Nero di Troia with FST mean value > 0.48. Clones belonging to AL1 showed the lowest mean values (FST = 0.27) with lower FST (0.16) detected with AL2.
Genetic distance of AL1 and AL2 sub-populations based on IBS and microsatellite analyses
An additional intra-variety IBS and Kinship analysis was performed for AL1 and AL2 individuals. The pairwise genetic distance range was different within the two populations, 0.61-0.70 for AL1 and 0.78-0.99 for AL2. IBS values calculated between AL1 and AL2 individuals showed that IBS ranged from 0.53 to 0.73.
The separation of Aglianico Lasco into two populations (namely AL1 and AL2) was further investigated using a set of seven discriminating microsatellite markers and the varieties Aglianico del Vulture, Aglianico Taurasi, Aglianico del Taburno, Sangiovese, Chardonnay and Merlot. The neighbour joining tree differentiated two main clusters (Figure 3). The first included all international varieties and Sangiovese, the other included AL1, AL2 and all Aglianico biotypes. In the latter, the AL1 clones were grouped together with Aglianico del Vulture, Taurasi and Aglianico del Taburno, while the AL2 clones were clearly grouped in a separate sub-cluster where no other varieties were included.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1201287/fpls-14-1201287-HTML/image_m/fpls-14-1201287-g003.jpg)
Figure 3 Dendrogram generated by neighbour-joining clustering of SSR markers for clones AL1, AL2, and Aglianico, as well as Sangiovese, Chardonnay, and Merlot. Bootstrap values greater than 60% are indicated.
Identification of divergent loci
Genetic differentiation among the seven populations described above was further investigated by Wright Fixation Index (FST) analysis at individual loci. All comparisons revealed numerous loci with an FST value > 0.25 (high differentiation) (Figure 4). To look for strongly divergent loci, FST ≥ 0.80 was used as a threshold, resulting in 1,016 divergent SNPs, out of which 285 SNPs were non redundant. Among the latter, 200 were in gene regions, and their distribution along the chromosomes was reported in Figure S6 and Supplementary Table S3. The largest number of divergent SNPs was found in Camaiola (N. = 453) followed by Nero di troia (N. = 404) and Minutolo (N. = 317), while AL1 showed the lowest number (N. = 117). An upset plot showing the co-occurrence of pairwise divergent loci is shown in Figure 5. As indicated by the “set size” in the Figure 5, Camaiola had the largest number of divergent genes between the other varieties, 100 against AL2, 91 against Minutolo, 85 against Nero di Troia, 81 against Malvasia, 57 against Greco B. and 39 against AL1 (Supplementary Table S3). In total, 81% of divergent genes were found in Camaiola comparisons. On the counterpart, AL1 and AL2 showed only two divergent genes (Supplementary Table S3), which represented the lowest number among all comparisons. To further our understanding of the 200 gene-associated divergent SNPs, we annotated the corresponding loci and categorised them into 14 major groups (e.g., development, hormone signalling, primary and secondary metabolism, stress response, etc.) (Supplementary Table S4). The most represented group was the “primary metabolism” with 62 genes, followed by “cellular component organisation and biogenesis” (28 genes), and “regulation of gene expression” (24 genes). The largest number of divergent SNPs (N=5) were found in VIT_214s0060g02480 (mainly identified in Malvasia comparisons), VIT_215s0021g01230 and VIT_212s0059g02300 (only in Nero di Troia vs MINUTOLO) and VIT_205s0029g01285 (mostly found in Nero di Troia comparisons), which are involved in “cellular component organisation and biogenesis”, “primary metabolism”, “regulation of gene expression” and “stress response”, respectively. Additional information associated with the 200 gene-associated divergent SNPs was retrieved from published data and allowed the identification of 63 SNPs associated with genes involved in grape phenology and adaptation to the environment (Supplementary Table S4). Among the latter genes, VIT_204s0008g01910 (mainly in Malvasia and Nero di Troia comparisons) and VIT_215s0048g02790 (mostly in Camaiola comparisons) were the most enriched in divergent SNPs. They are annotated as ferredoxin-related gene and PMT26 methyltransferase, respectively, and known to be linked to the initiation of berry ripening (Supplementary Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1201287/fpls-14-1201287-HTML/image_m/fpls-14-1201287-g004.jpg)
Figure 4 Number of SNPs with FST > 0.25 (yellow), FST > 0.5 (orange), FST > 0.8 (red) and FST =1 (darkred) detected in each comparison by computing pairwise FST index between all loci.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1201287/fpls-14-1201287-HTML/image_m/fpls-14-1201287-g005.jpg)
Figure 5 UpSet plot showing the number of the co-occurence of divergent SNPs (FST ≥ 0.80) shared between the seven populations. The dot-and-line chart in the bottom combination matrix indicates intersections between the seven populations. The upper bar chart (black bars) indicated the numbers of divergent loci in each set; The lower left horizontal bar chart (coloured bars) indicated the number of divergent loci for each intersection. Different colours represent comparisons of Apulia vs Apulia varieties (red), Campania vs Campania varieties (yellow) and Apulia vs Campania varieties (blue).
Discussion
The combination of ddRADSeq and GBS leads to a single high-quality SNP dataset
Six varieties were genotyped with two RRS strategies; GBS was used for the Apulian clones, while ddRADSeq for the Campanian ones. These methods differ in the number and type of restriction enzymes (REs) used to digest and access the genome and affected the total number of SNPs spanning on the reference genome. As expected, the ddRAD-seq dataset has more SNPs compared to fewer SNPs from GBS data, due to more restriction sites that are more polymorphic than GBS (Ruperao et al., 2023). The RE affected also the number of SNPs detected in genic and intergenic regions, indeed, the RE used for the Apulian dataset was ApeK1, a methylation-sensitive RE showing increased cleavage frequency in single-copy genomic regions that are enriched with genes (Elshire et al., 2011). Otherwise, a double digestion with a pair of rare and frequent cutter REs (SphI and MobI) performed for samples from Campania, according to the method proposed by Peterson et al. (2012) and modified by Scaglione et al. (2015), showed higher number of SNPs in intergenic regions. To address this challenging task, we searched for and found common SNP markers between the Campania and Apulia datasets and corrected for discrepancies using an in silico approach. Remarkably, after removing duplicate clones from each variety, we successfully merged all clones into a single dataset consisting of 2,235 SNP markers that was then used for cross-varietal analysis. To date, only Guillardín-Calvo et al. (2019) used both genotyping methods in Mediterranean evergreen oaks for technical comparison purposes. To the authors’ knowledge, this is the first time that ddRADSeq and GBS datasets have been merged and used to analyse the structure and genetic diversity of grapevine populations. Following this approach, we provide a consistent and repeatable framework that not only streamlines computational analyses, but also facilitates more general comparisons between different datasets, adhering to ‘FAIR’ (Findability, Accessibility, Interoperability, and Reusability) principles (Wilkinson et al., 2016).
RRS-based methods describe a complex genetic structure of grape varieties at the clonal level
As acknowledged by the International Organization of Vine and Wine in resolutions OIV-VITI 424/2010 and OIV-VITI 564B-2019, a crucial conservation strategy to protect the world heritage of grapevine is based on the assessment of intra-varietal diversity and polyclonal selection (Resolution OIV-VITI 424/2010, 2010; Resolution OIV-VITI 564B-2019, 2019). In this study, the use of RRS-based analysis allowed the identification of a reliable set of clonal genetic variants to be employed in our genotyping experiment and the elimination of identical individuals within each variety studied. The presence of duplicated samples within fields, nurseries and germplasm repositories is quite common in clonally propagated species (Ipek et al., 2008; Ciftci et al., 2017; Orek et al., 2022; Selga et al., 2022). An exceptionally large number of identical genotypes was observed within the Aglianico Lasco population, where more than 70% of collected samples was duplicated. This redundancy is probably linked to the propagation history of clones, likely deriving from few homogenous plots of mother plants which were multiplied by local growers in a restricted cultivation area, under similar environmental pressures and selection criteria. After elimination of duplicated individuals, we were left with a set of 78 individuals that we employed in downstream investigations. We found that Camaiola, Greco B., Minutolo, Malvasia, and Nero di Troia split into separated populations. This finding indicates that these varieties are represented by clones with a genetically uniform profile. An exception was found in the population of Aglianico Lasco which separated into two groups called AL1 and AL2. To shed lights on the genotypic identity of each sub-population of Aglianico Lasco, we performed a comparative analysis based on microsatellites using, among others, the main red variety grown in Campania (Aglianico) with which Aglianico Lasco is generally matched and confused. The AL1 population included individuals genetically similar to all Aglianico biotypes (Taurasi, Taburno and Vulture), whereas AL2 grouped individuals with microsatellite profiles divergent from those of AL1. The high IBS values and the Kinship coefficient within AL2 (Table S2) suggest that these individuals have a more preserved genetic diversity, and therefore can be considered the true Aglianico Lasco, whereas AL1 individuals showed a greater variability, indicating that they may be homonyms (different cultivars named alike) of the most renowned cultivar Aglianico, whichhas at least three different biotypes inside.
Misattribution of names is frequently reported in grapevine (Crespan and Milani, 2001; De Lorenzis et al., 2014; Sunseri et al., 2018; Fanelli et al., 2021) and is considered the main cause of varietal confusion (Villano et al., 2022).
To explain the genetic diversity of the varieties under study, it is necessary to consider their origin and ancestry. Despite PCA plot did not reflect the geographical grouping, neighbour-joining tree clearly separated the Campanian from the Apulian varieties, indicating that the geographical origin could be one of the main drivers in defining varietal relationships in grapevine species, as observed in other studies (Fanelli et al., 2021).
Nero di Troia takes its name from the city of Troia in the Daunia region, where it is mainly cultivated. It derives from the spontaneous crossing between the French Bouteillan and the Bombino bianco, giving origin also to other Apulian varieties such as Bombino nero and Impigno (Miazzi et al., 2020; D’Onofrio et al., 2021). The reduced genetic variability between clones of Nero di Troia and their high genetic distance from the other Apulian varieties could further support the hypothesis that the province of Foggia represents a biodiversity hotspot for different crops (Conversa et al., 2018; Pisani et al., 2021); native varieties such as Nero di Troia could have been marginalised in those areas and therefore have underexplored gene pools that deserve to be studied to search for new and beneficial alleles. By contrast, the populations of Malvasia nera, Minutolo, AL1 and AL2 partially overlapped. This could be explained considering that these four varieties share a common ancestor, namely Visparola. Indeed, it has been reported that Malvasia Nera derives from a cross between Negroamaro and Malvasia bianca lunga, where Negroamaro is a descendent of Maiolica, an offspring of Visparola. In addition, the aromatic variety Minutolo derives directly from Visparola, also related to Aglianico (putative AL1 in our study) and its offspring; this confirms the centrality of Visparola in the origin of many grapevine varieties of Southern Italy (D’Onofrio et al., 2021). As the AL2 and AL1 populations partially overlap in the PCA plot, we hypothesise a parent-offspring relationship between Aglianico and Aglianico Lasco. Regarding Camaiola, PCA and FST results showed that it differed more than anyone else from all the other varieties analysed. On the origin of this genetic distance, some hypotheses can be made based on the scarce information available. For example, the geographical distribution of Camaiola is restricted to the surroundings of Castelvenere city, where it is mostly cultivated, and has probably maintained a specific and uncontaminated genetic makeup. It should also be noted that DNA profiling provided evidence that Camaiola does not share close genetic links with any other Campania cultivar (Costantini et al., 2005; Villano et al., 2014). This suggests that Camaiola has been introduced relatively recently in Campania and supports the historical research by Selvaggio and Carlo (2015).
A subset of divergent SNP loci is related to phenology and plasticity
In this work, the search for genetic differentiation at single loci yielded nearly 200 divergent SNPs and detected putative genes under selection. Camaiola displayed the greatest number of divergent loci, which supports its separation from all the other varieties in the FST analysis. About 30% of the 200 divergent loci were in genes/QTLs involved in important phenological processes. For example, 12 divergent loci have been described as “switch genes” actively involved in the shift of berry developmental from immature to mature growth (Palumbo et al., 2014). This is a crucial physiological event that marks the ripening onset (called veraison) where numerous molecular, biochemical, and physiological changes occur that strongly impact the quality of wine (Conde et al., 2007; Fasoli et al., 2012; Castellarin et al., 2016; Fasoli et al., 2018). Among the divergent “switch genes” identified, VIT_203s0088g01250 and VIT_218s0001g14270 fell within phenology-related QTLs (Delfino et al., 2019) on chromosomes 3 and 18 responsible for phenotypes related to budding, flowering, soluble solid concentration, brix, and ripening. We also found signature of divergence in some genes involved in the regulation of multiple stress response and, therefore, putative players of plant adaptation to adverse biotic and abiotic conditions (Lee et al., 2010). For example, the arachidonic acid-induced DEA1 protein (VIT_202s0154g00280) in grapevine is involved in the phospholipid signalling processes and in the regulation of programmed cell death (PCD) in the scion/rootstock joining (Cookson and Ollat, 2013; Assunção et al., 2019). Notably, regulation of PCD is of paramount importance in plant-microbe interactions and is observed in many host resistance responses (Dickman and Fluhr, 2013). Greco B. and Nero di Troia were different at this locus. The drought tolerance reported for Nero di Troia (Antonacci, 2004) and its high susceptibility to powdery mildew compared to the low sensitivity to downy and powdery mildew of Greco B. (Rauscedo, 2011) are in agreement with these data.
We also found 14 divergent rootstock-responsive loci according to Corso et al. (2016). Among them, the protein CYCLIC NUCLEOTIDE-GATED CHANNEL 14 (CNGC14, VIT_204s0069g00790) and VIT_203s0038g02140. The former is known to be essential for the first step of auxin-induced Ca2+ signalling and growth inhibition in Arabidopsis root (Shih et al., 2015; Dindas et al., 2018); the latter is an auxin influx carrier protein, which is differentially expressed in graft interface tissues and possibly involved in directing the reconnection of vascular tissues (Yin et al., 2012; Cookson and Ollat, 2013). Finally, we identified an FKBP12-rapamycin complex-associated protein (VIT_203s0088g00450), with many divergent SNPs in Nero di Troia, which probably underwent selective pressure acting on this locus. The expression of VIT_203s0088g00450 was previously associated with bud dormancy (Shangguan et al., 2020), an essential adaptation process that allows temperate woody perennials to survive adverse environmental conditions during winter (Arora et al., 2003). Taken together, our data showed that the selection process shaped the genetic diversity of viticulturally attractive loci involved in fruit quality, growth and reproductive processes, the fine-tuning of which determines the fitness of fruit trees in changing climate conditions. We hypothesise that the molecular signatures we found could be the result of the selection operated by growers. It is possible that ripening time, sugar accumulation and stress response could have played a key role in determining whether a given cultivar was adapted to their specific local climate conditions. Interestingly, these functions appear strikingly similar to those in domesticated apple, peach, apricot and pear trees in which selective sweeps pointed to genes associated with growth cycle (Groppi et al., 2021) and fruit quality (Khan et al., 2014; Wu et al., 2018; Guan et al., 2021).
Our results provided evidence that a number of divergent loci overlapped with previously identified differentially methylated regions (DMR) in grapevine samples grown under changing environmental conditions (Xie et al., 2017). Among these we found the aforementioned VIT_203s0088g00450, and VIT_208s0007g07380, an HSP40 protein (also known as the DnaJ protein), which is known to contribute to cellular protein homeostasis under various biotic and abiotic stresses (Chen et al., 2021). Mounting evidence suggests that DNA methylation can affect local adaptation or plasticity by shaping the phenotypes that allow organisms to respond to their local environments (Richards et al., 2010). However, an active debate surrounds the question of whether (and to what extent) such epigenetic variations may either directly or indirectly affect local adaptation. Further experimental work is needed to corroborate these hypotheses.
Conclusions
In this study, we used genome-wide SNP datasets generated by the GBS and ddRADSeq methods to assess the clonal diversity of six traditional grapevine varieties. We have shown that merging different SNP datasets is possible and valuable to study the inter- and intra-specific genetic diversity of grape populations. This provided that the same reference genome is used. Through such an approach, we have provided a repeatable framework to streamline future computational studies based on the retrieval of information from partial analyses performed at different times and with different techniques. Our results also demonstrated the value of advanced genomic methods in the study of population structure and synonymy/homonymy detection (as exemplified by Aglianico Lasco) as well in identifying possible recent introduced outgroup variety such as Camaiola for Campania. This is a crucial task in any germplasm management and conservation strategy. Finally, we identified several divergent SNP loci within genes involved in grapevine phenology and environmental adaptation. This evidence emphasises that some traits, such as those related to budding, flowering, and fruit quality could have played a key role in assessing whether a given cultivar was adapted to specific local climatic conditions.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author contributions
FT, CM, CV, and RA designed the experiment. CM, CV, PN, and RA established the grapevine collection. ES, SP, MM, FT, and ND’A carried out part of the bioinformatic and the genetic diversity analyses and detected loci under selection. VF and CV performed genetic diversity indices and SSR analysis. All authors were involved in the data interpretation. CV, FT, MM, and RA wrote the draft manuscript. CV, MM, and FT critically revised the manuscript for important intellectual content. DC, PN, CM, RA, ND’A, and FT finalised the latest version of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the INDIGENA (CUP G12C21000190003), RINNOVALA (CUP J92C20001030009), and DICOVALE (CUP B24I19000450009) projects, all funded by the Campania Region, under the Rural Development Program 2014/2020, Measure 16.1, Action 2. The research was also supported by the Apulia Region within the PSR Puglia 2014-2020 REGEVIP project, (CUP B17H22003530009); Misura 10-Pagamenti agro-climatico-ambientali; Sottomisura 10.2-Sostegno per la conservazione, l’uso e lo sviluppo sostenibili delle risorse genetiche in agricoltura; Operazione 10.2.1-Progetti per la conservazione e valorizzazione. This study was also carried out within the Agritech National Research Center and received funding from the European Union Next- GenerationEU (PIANO NAZIONALE DI RIPRESA E RESILIENZA 6 (PNRR) – MISSIONE 4 COMPONENTE 2, INVESTIMENTO 1.4 – D.D. 1032 17/06/2022, CN00000022).
Acknowledgments
We thank Dr. R. Garramone and Dr. A. Gallucci for technical assistance.
Conflict of interest
Author CM was employed by the company SINAGRI S.r.l.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author disclaimer
This paper reflects only the authors’ views and opinions, and neither the European Union or the European Commission can be considered responsible for them.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1201287/full#supplementary-material
Supplementary Figure 1 | SNP density plot showing the number of filtered SNPs in 1 Mb size windows for the 19 Vitis vinifera chromosomes in (A) Campanian and (B) Apulian datasets.
Supplementary Figure 2 | Scatter plot showing the decay of linkage disequilibrium (r2 = 0.20) calculated for each of the six varieties.
Supplementary Figure 3 | SNP density plot showing the number of SNPs in 1 Mb size windows for the 19 chromosomes of Vitis vinifera.
Supplementary Figure 4 | Cross-validation error estimates for each value of K (i.e., number of sub-populations) tested.
Supplementary Figure 5 | (A) Percentage cumulative variance for retained PCA eigen vectors. (B) Bayesian information criterion (BIC). (C) Scatterplot of the DAPC and DAPC F statistic.
Supplementary Figure 6 | Distribution of divergent SNPs along the 19 chromosomes of Vitis vinifera.
Supplementary Table 1 | List of clones under study. For each clone, the region of origin, the name of the variety, the clone identifier and the sampling site are reported. The clones selected after IBS analysis are in bold. The asterisks indicate the individual removed after filtering on merged dataset.
Supplementary Table 2 | Identity-by-state (IBS) values and Kinship coefficients (Relatedness-PHI) calculated for the AL1 and AL2 populations.
Supplementary Table 3 | List of divergent SNPs (FST ≥ 0.80) falling into gene regions.
Supplementary Table 4 | List of divergent SNPs derived from pairwise comparisons between clusters identified by ADMIXTURE and DAPC. Gene identifier, number of SNPs identified following pairwise comparison, gene function and functional group were reported.
References
Alexander, D. H., Lange, K. (2011). Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinf. 12 (1), 1–6. doi: 10.1186/1471-2105-12-246
Andrews, K., Good, J., Miller, M., Luikart, G., Hohenlohe, P. A. (2016). Harnessing the power of RADseq for ecological and evolutionary genomics. Nat. Rev. Genet. 17, 81–92. doi: 10.1038/nrg.2015.28
Antonacci, D., Velenosi, M., Rocco, P., Crupi, P., Ventura, M., Bergamini, C. (2017). Utilizzo di vitigni di antica coltivazione del Meridione d'Italia per la produzione di vini passiti di pregio In BIO Web of Conferences (Vol. 9, p. 04001). EDP Sciences.
Arora, R., Rowland, L. J., Tanino, K. (2003). Induction and release of bud dormancy in woody perennials: a science comes of age. HortScience 38 (5), 911–921. doi: 10.21273/HORTSCI.38.5.911
Assunção, M., Pinheiro, J., Cruz, S., Brazão, J., Queiroz, J., Dias, J. E. E., et al. (2019). Gallic acid, sinapic acid and catechin as potential chemical markers of Vitis graft success. Scientia Hortic. 246, 129–135. doi: 10.1016/j.scienta.2018.10.056
Calderón, L., Mauri., N., Muñoz, C., Carbonell-Bejerano, P., Bree, L., Bergamin, D., et al. (2021). Whole genome resequencing and custom genotyping unveil clonal lineages in ‘Malbec’grapevines (Vitis vinifera L.). Sci. Rep. 11 (1), 1–11. doi: 10.1038/s41598-021-87445-y
Castellarin, S. D., Gambetta, G. A., Wada, H., Krasnow, M. N., Cramer, G. R., Peterlunger, E., et al. (2016). Characterization of major ripening events during softening in grape: turgor, sugar accumulation, abscisic acid metabolism, colour development, and their relationship with growth. J. Exp. Bot. 67 (3), 709–722. doi: 10.1093/jxb/erv483
Catalogo nazionale delle varietà di vite. The National Register of Grapevine Varieties. Available at: http://catalogoviti.politicheagricole.it (Accessed February 6, 2023).
Catchen, J., Hohenlohe, P. A., Bassham, S., Amores, A., Cresko, W. A. (2013). Stacks: an analysis tool set for population genomics. Mol. Ecol. 22 (11), 3124–3140. doi: 10.1111/mec.12354
Chen, T., Xu, T., Zhang, T., Liu, T., Shen, L., Chen, Z., et al. (2021). Genome-wide identification and characterization of dnaJ gene family in grape (Vitis vinifera L.). Horticulturae 7 (12), 589. doi: 10.3390/horticulturae7120589
Ciftci, A., Karatay, H., Kücükosmanoğlu, F., Karahan, A., Kaya, Z. (2017). Genetic differentiation between clone collections and natural populations of European black poplar (Populus nigra L.) in Turkey. Genet. Genomes 13, 69. doi: 10.1007/s11295-017-1154-8
Conde, C., Silva, P., Fontes, N., Dias, A. C. P., Tavares, R. M., Sousa, M. J., et al. (2007). Biochemical changes throughout grape berry development and fruit and wine quality. Food 1, 1–22.
Conversa, G., Bonasia, A., Elia, A. (2018). Collecting landrace vegetable crops in Puglia region: case study of Foggia Province. Acta Hortic. 1215, 453–458. doi: 10.17660/ActaHortic.2018.1215.82
Cookson, S. J., Ollat, N. (2013). Grafting with rootstocks induces extensive transcriptional re-programming in the shoot apical meristem of grapevine. BMC Plant Biol. 13 (1), 1–14. doi: 10.1186/1471-2229-13-147
Corso, M., Vannozzi, A., Ziliotto, F., Zouine, M., Maza, E., Nicolato, T., et al. (2016). Grapevine rootstocks differentially affect the rate of ripening and modulate auxin-related genes in cabernet sauvignon berries. Front. Plant Sci. 7, 69. doi: 10.3389/fpls.2016.00069
Costantini, L., Monaco, A., Vouillamoz, J., Forlani, M., Grando, M. S. (2005). Genetic relationship among local Vitis vinifera cultivars from Campania (Italy). Vitis 44 (1), 25–34.
Crespan, M., Milani, N. (2001). The Muscats: A molecular analysis of synonyms, homonyms and genetic relationships within a large family of grapevine cultivars. VITIS-GEILWEILERHOF- 40 (1), 23–30.
Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., et al. (2011). The variant call format and VCFtools. Bioinformatics 27 (15), 2156–2158. doi: 10.1093/bioinformatics/btr330
Delfino, P., Zenoni, S., Imanifard, Z., Tornielli, G. B., Bellin, D. (2019). Selection of candidate genes controlling veraison time in grapevine through integration of meta-QTL and transcriptomic data. BMC Genomics 20, 1–19. doi: 10.1186/s12864-019-6124-0
De Lorenzis, G., Las Casas, G., Brancadoro, L., Scienza, A. (2014). Genotyping of Sicilian grapevine germplasm resources (V. vinifera L.) and their relationships with Sangiovese. Scientia Hortic. 169, 189–198. doi: 10.1016/j.scienta.2014.02.028
De Lorenzis, G., Mercati, F., Bergamini, C., Cardone, M. F., Lupini, A., Mauceri, A., et al. (2019). SNP genotyping elucidates the genetic diversity of Magna Graecia grapevine germplasm and its historical origin and dissemination. BMC Plant Biol. 19, 7. doi: 10.1186/s12870-018-1576-y
Dickman, M. B., Fluhr, R. (2013). Centrality of host cell death in plant-microbe interactions. Annu. Rev. Phytopathol. 51, 543–570. doi: 10.1146/annurev-phyto-081211-173027
Dindas, J., Scherzer, S., Roelfsema, M. R. G., von Meyer, K., Müller, H. M., Al-Rasheid, K. A. S., et al. (2018). AUX1-mediated root hair auxin influx governs SCFTIR1/AFB-type Ca2+ signaling. Nat. Commun. 9 (1), 1174. doi: 10.1038/s41467-018-03582-5
D’Onofrio, C., Scalabrelli, G. (2010). Un database viticolo universale. 2° Convegno Nazionale di Viticoltura (Marsala. Italus Hortus) Marsala 14-19 Luglio 2008 (pp. 102–102). Convegno Nazionale di Viticoltura.
D’Onofrio, C., Tumino, G., Gardiman, M., Crespan, M., Bignami, C., De Palma, L., et al. (2021). Parentage atlas of Italian grapevine varieties as inferred from SNP genotyping. Front. Plant Sci. 11, 605934. doi: 10.3389/fpls.2020.605934
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PloS One 6 (5), e19379. doi: 10.1371/journal.pone.0019379
Fanelli, V., Roseti, V., Savoia, M. A., Miazzi, M. M., Venerito, P., Savino, V. N., et al. (2021). New insight into the identity of italian grapevine varieties: the case study of calabrian germplasm. Agronomy. 11 (8), 1538. doi: 10.3390/agronomy11081538
FAOSTAT (2023) (Accessed February 12, 2023). Available at: https://www.fao.org/faostat.
Fasoli, M., Dal Santo, S., Zenoni., S., Tornielli, G. B., Farina, L., Zamboni, A., et al. (2012). The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 24 (9), 3489–3505. doi: 10.1105/tpc.112.100230
Fasoli, M., Richter, C. L., Zenoni, S., Bertini, E., Vitulo, N., Dal Santo, S., et al. (2018). Timing and order of the molecular events marking the onset of berry ripening in grapevine. Plant Physiol. 178 (3), 1187–1206. doi: 10.1104/pp.18.00559
Groppi, A., Liu, S., Cornille, A., Decroocq, S., Bui, Q. T., Tricon, D., et al. (2021). Population genomics of apricots unravels domestication history and adaptive events. Nat. Commun. 12 (1), 3956. doi: 10.1038/s41467-021-24283-6
Guan, J., Xu, Y., Yu, Y., Fu, J., Ren, F., Guo, J., et al. (2021). Genome structure variation analyses of peach reveal population dynamics and a 1.67 Mb causal inversion for fruit shape. Genome Biol. 22, 1–25. doi: 10.1186/s13059-020-02239-1
Guillardín-Calvo, L., Mora-Márquez, F., Soto, Á., López de Heredia, U. (2019). RADdesigner: a workflow to select the optimal sequencing methodology in genotyping experiments on woody plant species. Tree Genet. Genomes 15 (4), 1–14. doi: 10.1007/s11295-019-1372-3
Ipek, M., İpek, A., Simon, P. W. (2008). Rapid characterization of garlic clones with locus-specific DNA markers. Turkish J. Agric. Forestry 32 (5), 357–362.
Italian Vitis Database. (2023). (Accessed March 2, 2023). Available at: http://www.vitisdb.it/.
Japelaghi, R. H., Haddad, R., Garoosi, G. A. (2011). Rapid and efficient isolation of high-quality nucleic acids from plant tissues rich in polyphenols and polysaccharides. Mol. Biotechnol. 49 (2), 129–137. doi: 10.1007/s12033-011-9384-8
Jombart, T., Ahmed, I. (2011). Adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics 27 (21), 3070–3071. doi: 10.1093/bioinformatics/btr521
Jombart, T., Devillard, S., Balloux, F. (2010). Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11 (1), 1–15. doi: 10.1186/1471-2156-11-94
Khan, M. A., Olsen, K. M., Sovero, V., Kushad, M. M., Korban, S. S. (2014). Fruit quality traits have played critical roles in domestication of the apple. Plant Genome 7 (3), 2014–2004. doi: 10.3835/plantgenome2014.04.0018
Kumar, S., Banks, T. W., Cloutier, S. (2012). SNP discovery through next-generation sequencing and its applications. Int. J. Plant Genomics 2012, 15. doi: 10.1155/2012/831460
Laucou, V., Lacombe, T., Dechesne, F., Siret, R., Bruno, J. P., Dessup, M., et al. (2011). High throughput analysis of grape genetic diversity as a tool for germplasm collection management. Theor. Appl. Genet. 122 (6), 1233–1245. doi: 10.1007/s00122-010-1527-y
Lee, J. M., Kubota, C., Tsao, S. J., Bie, Z., Echevarria, P. H., Morra, L., et al. (2010). Current status of vegetable grafting: Diffusion, grafting techniques, automation. Scientia Hortic. 127 (2), 93–105. doi: 10.1016/j.scienta.2010.08.003
Li, H. (2011). A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27 (21), 2987–2993. doi: 10.1093/bioinformatics/btr509
Li, H., Durbin, R. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25 (14), 1754–1760. doi: 10.1093/bioinformatics/btp324
Miazzi, M. M., D’Agostino, N., di Rienzo, V., Venerito, P., Savino, V. N., Fucilli, V., et al. (2020). Marginal grapevine germplasm from Apulia (Southern Italy) represents an unexplored source of genetic diversity. Agronomy 10 (4), 563. doi: 10.3390/agronomy10040563
O'Leary, S. J., Puritz, J. B., Willis, S. C., Hollenbeck, C. M., Portnoy, D. S. (2018). These aren’t the loci you’re looking for: Principles of effective SNP filtering for molecular ecologists. Mol. Ecol. 27, 1–14. doi: 10.1111/mec.14792
Orek, C., Kyallo, M., Oluwaseyi, S., Yao, N. (2022). Genotyping by sequencing reveals genetic relatedness and duplicates amongst local cassava (Manihot esculenta Crantz) landraces and improved genotypes in Kenya. Res. Square. doi: 10.21203/rs.3.rs-1295398/v1
Palumbo, M. C., Zenoni, S., Fasoli, M., Massonnet, M., Farina, L., Castiglione, F., et al. (2014). Integrated network analysis identifies fight-club nodes as a class of hubs encompassing key putative switch genes that induce major transcriptome reprogramming during grapevine development. Plant Cell 26 (12), 4617–4635. doi: 10.1105/tpc.114.133710
Pavan, S., Delvento, C., Ricciardi, L., Lotti, C., Ciani, E., D’Agostino, N. (2020). Recommendations for choosing the genotyping method and best practices for quality control in crop genome-wide association studies. Front. Genet. 11, 447. doi: 10.3389/fgene.2020.00447
Peakall, R. O. D., Smouse, P. E. (2006). GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 6 (1), 288–295. doi: 10.1111/j.1471-8286.2005.01155.x
Perrier, X., Flori, A., Bonnot, F. (2003). “Data analysis methods,” in Genetic diversity of cultivated tropical plants. Eds. Hamon, P., Seguin, M., Perrier, X., Glaszmann, J.C. (Montpellier: Enfield, Science Publishers), 43–76.
Peterson, B. K., Weber, J. N., Kay, E. H., Fisher, H. S., Hoekstra, H. E. (2012). Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PloS One 7 (5), e37135. doi: 10.1371/journal.pone.0037135
Pisani, D., Pazienza, P., Perrino, E. V., Caporale, D., De Lucia, C. (2021). The economic valuation of ecosystem services of biodiversity components in protected areas: a review for a framework of analysis for the Gargano national park. Sustainability 13 (21), 11726. doi: 10.3390/su132111726
Pomarici, E., Corsi, A., Mazzarino, S., Sardone, R. (2021). The Italian wine sector: evolution, structure, competitiveness and future challenges of an enduring leader. Ital. Economic J. 7 (2), 259–295. doi: 10.1007/s40797-021-00144-5
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81 (3), 559–575. doi: 10.1086/519795
Rambaut, A. (2014) FigTree v1. 4.2, a graphical viewer of phylogenetic trees. Available at: http://tree.bio.ed.ac.uk/software/figtree/right.
Rauscedo (2011). Available at: https://vivairauscedo.com/contributi/it/schede/greco-di-tufo.pdf (Accessed February, 2023).
Resolution OIV-VITI 424/2010. (2010). Conservation of Vine’s Genetic Resources. International Organisation of Vine and Wine: Paris, France.
Resolution OIV-VITI 564B-2019. (2019). Process for the Recovery and Conservation of the Intravarietal Diversity and the Polyclonal Selection of the Vine in Grape Varieties with Wide Genetic Variability. International Organisation of Vine and Wine: Paris, France.
Richards, C. L., Bossdorf, O., Pigliucci, M. (2010). What role does heritable epigenetic variation play in phenotypic evolution? BioScience 60 (3), 232–237. doi: 10.1525/bio.2010.60.3.9
Ruperao, P., Bajaj, P., Subramani, R., Yadav, R., Reddy Lachagari, V. B., Lekkala, S. P., et al. (2023). A pilot-scale comparison between single and double-digest RAD markers generated using GBS strategy in sesame (Sesamum indicum L.). PloS One 18 (6), e0286599. doi: 10.1371/journal.pone.0286599
Salmaso, M., Faes, G., Segala, C., Stefanini, M., Salakhutdinov, I., Zyprian, E., et al. (2005). Genome diversity and gene haplotypes in the grapevine (Vitis vinifera L.), as revealed 12 by single nucleotide polymorphisms. Mol. Breed. 14 (4), 385–395. doi: 10.1007/s11032-005-0261-7
Sardaro, R., Bozzo, F., Petrillo, F., Fucilli, V. (2017). Measuring the financial sustainability of vine landraces for better conservation programmes of Mediterranean agro-biodiversity. Land Use Policy 68, 160–167. doi: 10.1016/j.landusepol.2017.07.045
Scaglione, D., Fornasiero, A., Cattonaro, F., Mesa, K., Pinto, C., Shiny, P., et al. (2015). A RAD-based linkage map of kiwifruit (Actinidia chinensis Pl.) and map position of genetic determinants of traits of interest to breeders. Acta Hortic. 1096, 191–198. doi: 10.17660/ActaHortic.2015.1096.20
Selga, C., Chrominski, P., Carlson-Nilsson, U., Andersson, M., Chawade, A., Ortiz, R. (2022). Diversity and population structure of Nordic potato cultivars and breeding clones. BMC Plant Biol. 22, 350. doi: 10.1186/s12870-022-03726-2
Selvaggio, M., Carlo, P. (2015). Catelvenere valdese: storia e personaggi nell'epocale trasformazione da borgo a paese (Benevento, Italy: Realtà Sannita).
Shangguan, L., Wang, Z., Fang, X., Chen, C., Liu, T., Chen, M., et al. (2020). Regulation mechanism of bud dormancy by DAM gene in perennial fruit trees. J. Nanjing Agric. Univ. 43 (5), 790–799. doi: 10.7685/jnau.201912069
Shih, H. W., DePew, C. L., Miller, N. D., Monshausen, G. B. (2015). The cyclic nucleotide-gated channel CNGC14 regulates root gravitropism in Arabidopsis thaliana. Curr. Biol. 25 (23), 3119–3125. doi: 10.1016/j.cub.2015.10.025
Spadoni, A., Sion, S., Gadaleta, S., Savoia, M. A., Piarulli, L., Fanelli, V., et al. (2019). A simple and rapid method for genomic DNA extraction and microsatellite analysis in tree plants. J. Agric. Sci. Technol. 21 (5), 1215–1226.
Sunseri, F., Lupini, A., Mauceri, A., De Lorenzis, G., Araniti, F., Brancadoro, L., et al. (2018). Single nucleotide polymorphism profiles reveal an admixture genetic structure of grapevine germplasm from Calabria, Italy, uncovering its key role for the diversification of cultivars in the Mediterranean Basin. Aust. J. Grape Wine 40 Res. 24 (3), 345–359. doi: 10.1111/ajgw.12339
Tamura, K., Stecher, G., Kumar, S. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38 (7), 3022–3027. doi: 10.1093/molbev/msab120
Taranto, F., D’ Agostino, N., Rodriguez, M., Pavan, S., Minervini, A. P., Pecchioni, N., et al. (2020). Whole genome scan reveals molecular signatures of divergence and selection related to important traits in durum wheat germplasm. Front. Genet. 217. doi: 10.3389/fgene.2020.00217
This, P., Jung, A., Boccacci, P., Borrego, J., Botta, R., Costantini, L., et al. (2004). Development of a standard set of microsatellite reference alleles for identification of grape 15 cultivars. Theor. Appl. Genet. 109 (7), 1448–1458. doi: 10.1007/s00122-004-1760-3
Vezzulli, S., Troggio, M., Coppola, G., Jermakow, A., Cartwright, D., Zharkikh, A., et al. (2008). A reference integrated map for cultivated grapevine (Vitis vinifera L.) from three crosses, based on 283 SSR and 501 SNP-based markers. Theor. Appl. Genet. 117, 499–511. doi: 10.1007/s00122-008-0794-3
Villano, C., Aiese Cigliano, R., Esposito, S., D’Amelia, V., Iovene, M., Carputo, D., et al. (2022). DNA-based technologies for grapevine biodiversity exploitation: state of the art and future perspectives. Agronomy 12 (2), 491. doi: 10.3390/agronomy12020491
Villano, C., Carputo, D., Frusciante, L., Santoro, X., Aversano, R. (2014). Use of SSR and retrotransposon-based markers to interpret the population structure of native grapevines from southern Italy. Mol. Biotechnol. 56 (11), 1011–1020. doi: 10.1007/s12033-014-9780-y
Vishwakarma, M. K., Yadav, P. S., Rai, V. P., Kumar, U., Joshi, A. K. (2022). Relationship between microbes and the environment for sustainable ecosystem services Volume 1: Microbial Products for Sustainable Ecosystem Services. (Elsevier).
Vitis International Variety Catalogue VIVC (2023). Available at: https://www.vivc.de/ (Accessed March 7, 2023).
Wilkinson, M., Dumontier, M., Aalbersberg, I., Appleton, G., Axton, M., Baak, A., et al. (2016). The FAIR Guiding Principles for scientific data management and stewardship. Sci. Data 3, 160018. doi: 10.1038/sdata.2016.18
Wu, J., Wang, Y., Xu, J., Korban, S. S., Fei, Z., Tao, S., et al. (2018). Diversification and independent domestication of Asian and European pears. Genome Biol. 19 (1), 1–16. doi: 10.1186/s13059-018-1452-y
Xie, H., Konate, M., Sai, N., Tesfamicael, K. G., Cavagnaro, T., Gilliham, M., et al. (2017). Global DNA methylation patterns can play a role in defining terroir in grapevine (Vitis vinifera cv. Shiraz). Front. Plant Sci. 8, 1860. doi: 10.3389/fpls.2017.01860
Keywords: Vitis vinifera, intra-varietal diversity, genotyping by sequencing, double digest restriction associated DNA, molecular markers, divergent loci
Citation: Villano C, Procino S, Blaiotta G, Carputo D, D’Agostino N, Di Serio E, Fanelli V, La Notte P, Miazzi MM, Montemurro C, Taranto F and Aversano R (2023) Genetic diversity and signature of divergence in the genome of grapevine clones of Southern Italy varieties. Front. Plant Sci. 14:1201287. doi: 10.3389/fpls.2023.1201287
Received: 06 April 2023; Accepted: 21 August 2023;
Published: 13 September 2023.
Edited by:
Lucian Copolovici, Aurel Vlaicu University of Arad, RomaniaReviewed by:
Ezio Portis, University of Turin, ItalyRadu Liviu Sumalan, Banat University of Agricultural Sciences and Veterinary Medicine, Romania
Annarita Marrano, Phoenix Bioinformatics, United States
Copyright © 2023 Villano, Procino, Blaiotta, Carputo, D’Agostino, Di Serio, Fanelli, La Notte, Miazzi, Montemurro, Taranto and Aversano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cinzia Montemurro, cinzia.montemurro@uniba.it; Francesca Taranto, francesca.taranto@cnr.it
†These authors have contributed equally to this work