Next Article in Journal
Friends in Arms: Flavonoids and the Auxin/Cytokinin Balance in Terrestrialization
Previous Article in Journal
Structural and Physiological Traits of Compound Leaves of Ceratonia siliqua Trees Grown in Urban and Suburban Ambient Conditions
Previous Article in Special Issue
A Field Collection of Indigenous Grapevines as a Valuable Repository for Applied Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morphological and Genetic Clonal Diversity within the ‘Greco Bianco’ Grapevine (Vitis vinifera L.) Variety

1
Department of Agricultural Science, University of Naples Federico II, 80055 Portici, Italy
2
Department of Environmental, Biological and Pharmaceutical Sciences and Technologies, University of Campania “Luigi Vanvitelli”, 81100 Caserta, Italy
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(3), 515; https://doi.org/10.3390/plants12030515
Submission received: 4 January 2023 / Revised: 18 January 2023 / Accepted: 19 January 2023 / Published: 23 January 2023
(This article belongs to the Special Issue Conservation of Grapevine Genetic Resources)

Abstract

:
Grapevine (Vitis vinifera L.) has been propagated vegetatively for hundreds of years. Therefore, plants tend to accumulate somatic mutations that can result in an intra-varietal diversity capable of generating distinct clones. Although it is common that winemakers request specific clones or selections for planting new vineyards, relatively limited information is available on the extent, degree, and morphological impact of the clonal diversity in traditional, highly valued grapevine varieties within production areas protected by geographical denomination of origin. Here, we present a morphological and genetic investigation of the intra-varietal diversity in ‘Greco Bianco’, the grapevine variety used to produce the DOCG and PDO “Greco di Tufo” wine. Seventeen clones from different farms (all within the allowed production area) were phenotypically characterized using ampelographic and ampelometric traits. The clones were also genotyped with Simple Sequence Repeats (SSR) and retrotransposon-based DNA markers (REMAP). The morphological analysis indicated a uniformity in the qualitatively scored traits, and a limited variability for the quantitative traits of the bunch and of the berry composition. The molecular markers also depicted variability among clones, which was more evident with the use of REMAPs. The comparison of the discriminatory information of the three analyses indicated that they provided different estimates of the level of diversity. The evaluation described herein of the clonal variability has implications for the management and protection of clonal selections in ‘Greco Bianco’ and prompts for further multidisciplinary investigations on its possible role in winemaking.

1. Introduction

Cultivated grapevine, Vitis vinifera L., is one of the most widely grown fruit crops in the world [1]. Italy ranks first in the EU as a wine producer, with 44.500 million of liters of wine produced annually. This is a consequence of the gastronomy, traditions, and pedo-climatic conditions, but also of a rich and heterogeneous cultivated germplasm. Currently, 2072 cultivars are registered in the Italian Catalogue of Grapevine Varieties (http://catalogoviti.politicheagricole.it; accessed on the 1st of December 2022). This number is probably an underestimation because the real consistency of the varieties is difficult to establish and could be higher. In the cultivated grapevine germplasm, as in other clonally propagated fruit trees, somatic mutations add genetic variability within a variety (the so-called inter-varietal diversity) [2]. For instance, bud sports are usually agamically propagated by growers and are often considered and labeled as clones. Therefore, traditional varieties should be viewed as a group of clones with genotypic, morphological, physiological, and stress-resistance traits that may not necessarily be identical to those of the original mother plant [3,4,5].
The assessment of grape diversity is an essential component of germplasm characterization and conservation, which in turn are necessary to sustain and improve crop production. Morphological traits of the fruits have often been used for the characterization of fruit trees because of their obvious economic importance and discriminating power [6]. Specifically for grape, morphological traits combined with multivariate statistical methods have been long used to assess the genetic variation and relationships among genotypes and cultivars [7,8,9,10,11]. In the last decades, grapevine discrimination methods have taken advantage of the increasing economic and technical affordability of DNA-based techniques. These can also overcome the limitations of phenotypic-based diversity analysis [12]. Among them, microsatellites (also known as Simple Sequence Repeats—SSRs) have emerged as highly versatile and informative markers for plants because they are codominant, multiallelic, and amenable to multiplexing and automatic sizing [13]. Briefly, SSR typing is based on the amplification of short, tandemly repeated sequences by using locus-specific flanking primers. For example, SSRs have been applied in studies of Italian grapes belonging to various regions [14,15,16,17,18,19,20,21,22]. These studies estimated genetic diversity in each collection and demonstrated the usefulness of these markers in discriminating varieties. Retrotransposon-based markers represent another class of anonymous and highly informative DNA markers used in plants. Briefly, retrotransposons (RTs) are dynamic and highly dispersed DNA elements, hundreds to thousands of nucleotides long, and able to translocate and change their genomic location. The various retrotransposon-based marker systems that have been developed exploit size polymorphisms (RTs cause large insertions by their transpositional activity) and the presence of conserved domains, which allow the design of specific PCR primers. For instance, in the REtrotransposon–Microsatellite Amplified Polymorphism (REMAP) system, the amplification is carried out with a primer matching the conserved end-sequence of a long terminal repeat (LTR) RT, and a primer binding to the proximity of the tandemly repeated and highly variable core sequence of an SSR locus [23]. REMAPs are cheap, easy to perform and highly polymorphic [24]. In grapevine, three retro-elements have been identified, namely Vine-1, Tvv1, and Gret1 (the only one fully sequenced) [5,25,26]. The latter is associated with mutations causing most white-fruited V. vinifera genotypes due to its insertion into the promoter of VvMybA1, which codes for a transcription factor controlling the final step in anthocyanin biosynthesis during ripening [27].
Among the rich Italian grapevine germplasm, ‘Greco Bianco’ (hereafter ‘Greco’) is a white grape variety from the Campania region used to produce the “Greco di Tufo”, an Italian wine known for its pleasant, intense, elegant, and distinctive aroma. The “Greco di Tufo” benefits from the Italian quality label “Denomination of Controlled and Guaranteed Origin” (DOCG) (DM 18.07.2003 G.U. 180, 05-08-2003), and it is also protected in the EU and the UK by the Protected Designation of Origin (PDO) geographical indication EU label (PDO-IT-A0236). It is believed that this wine has an ancient origin. One of the earliest pieces of evidence of a ‘Greco’ variety in the scientific ampelographic literature is probably the “Ampélographie (Traité général de viticulture)”, edited by Viala and Vermorel in the first decade of the last century [28]. The authors stated the hypothesis that the “Greco Bianco”, which is cultivated in the Avellino Province, is probably the same variety of the Aminaea Gemella described by ancient authors. Pliny the Elder (AD 23–79) in his Naturalis Historia reported the cultivation, which occurred on Mount Vesuvius and on the hills around Sorrento of Aminaea. The grapevines he described were characterized by twin bunches (gemella), which is a typical morphological trait of the “Greco Bianco” cultivated in modern times. From those areas, the cultivar later spread to other areas of the Campania region. Currently, the DOCG production is limited to a few towns in the Irpinia area (Avellino province), with a total area of 900 ha dedicated to viticulture (https://tinyurl.com/2p93eyyf; accessed on 1 December 2022). The wine produced in this area took the current name “Greco di Tufo”, which is probably from the tuff rock widely present in the subsoil. Interestingly, the tuff is also typical of the areas around Mount Vesuvius, as it is a common effusive magmatic rock.
The long history of cultivation of the ‘Greco’ variety suggests complex and diverse evolutionary traces following agamic propagation by farmers. Their extent and relevance at the phenotypic level is largely unknown. In this work, we examined a representative selection of varieties used by farmers to produce the DOGC “Greco di Tufo” in a long-term effort to understand the possible adaptive potential of clonal diversity in “Greco di Tufo” production and vinification. Specifically, this study aimed to analyze and quantify the intra-varietal diversity of the ‘Greco’ using DNA analysis based on two marker systems (microsatellite and REMAP markers) and a morphological characterization mainly focusing on pomological traits. Moreover, in a first attempt to understand the relevance of the observed diversity, we correlated genetic and morphological distances. The gathered information helps design germplasm maintenance strategies and, ultimately, would guide selection and breeding.

2. Results

2.1. Diversity in Morphological and Berry Composition Traits

The 17 clones did not differ in terms of the ten qualitative morphological traits analyzed in this study. In summary, they all were characterized by (a) shoots with two consecutive tendrils; (b) mature leaves having three lobes, a cordate/pentagonal blade shape, teeth straight on both sides, a U-shaped base of both the petiole and the upper lateral sinuses; (c) flowers with fully developed stamen and gynoecium; (d) conical bunch shape; (e) one wing in the primary bunch; and (f) a broad-ellipsoid berry shape.
Bud potential fertility index (PFI) did not differ among the 17 clones (with an average PFI of 1.2 bunches per shoot), whereas, at harvest, bunch morphology varied (Table 1). Large variability was found among the clones in the components of bunch compactness. The number of berries per bunch, the total rachis length, and the bunch compactness index ranged, respectively, between 77 (G16) and 157 berries per bunch (G14), between 10.6 (G12) and 21.1 cm (G1), and between 4.77 (G5) and 9.90 berries/cm (G11), respectively. Clones G11 and G14 had the most compact bunches mainly because of the high number of berries per bunch (133 and 157 berries per bunch, respectively), which is also associated, in the case of G11, with a short rachis length (14 cm). However, G1, G5, and G16 had the loosest bunches (an average IC of 4.92 berries per bunch), mainly because of a long rachis (G1 and G16) or a low number of berries per bunch (G5). The other clones had bunches with intermediate compactness (Table 1).
Berry composition at harvest differed significantly between clones (Table 2). G1 and G17 had the highest juice TSS (around 24.3 °Brix) and intermediate TA values (10.4 and 9.6 g/L tartaric acid, respectively), whereas G9 berries were characterized by relatively low TSS (20.2 °Brix) and TA (7.95 g/L tartaric acid). Juice TA at harvest was lowest in G12 (7.2 g/L tartaric acid). The other clones had intermediate TSS (between 21.0 and 23.6 °Brix) and TA (between 8.55 and 12 g/L tartaric acid). Juice pH at harvest also significantly varied between clones ranging between 2.96 (G2) and 3.21 (G6 and G12). G12 berries had the highest TSS/TA ratio (3.16), followed by G7 with TSS/TA values of 2.62. The lowest values of the TSS/TA ratio were measured in G2 and G3 (Table 2). The other clones had TSS/TA values between 1.87 (G14) and 2.55 (G15).
For multivariate exploratory analysis, we removed from the initial 18 parameters those that were not statistically different among clones, namely the ten qualitative traits and PFI. The PCA extracted a total of six principal components, with the first three having eigenvalues higher than 1 and explaining together 82.6% of the total variance (Table 3). The first principal component (Dim. 1) was positively correlated with rachis length and TA, and negatively with pH (Table 3 and Figure 1). The second principal component (Dim. 2) correlated positively with the number of berries per bunch and bunch compactness index. Dim. 3 correlated positively with rachis length (Table S1 and Figure 1).
The bi-dimensional spaces delineated combining the first two principal components (Figure 1) separated the 17 clones well. Specifically, they were uniformly scattered, as indicated, for instance, by the almost equal number of clones in each quadrant, and many were distant from the plot origin. Moreover, clusters or similar clones were not evident even when a priori knowledge and possible groups related to the geographical locations were taken into account. Finally, the PCA bi-plot indicated that the patterns of morphological variance depicted by the PCs may have different, un-orthogonally correlated causes given the uniform distribution and the correlation of the variable loadings.

2.2. Genetic Diversity

2.2.1. Microsatellite Analysis

The DNA of the 17 clones under investigation was analyzed with ten highly polymorphic SSR loci described in the literature (Table S1). All loci were polymorphic in our population (presenting more than one allele), and co-dominance was confirmed by the presence of a maximum of two alleles per clone. Main genetic indexes are reported in Table 4.
We obtained 29 alleles in total, with an average of 2.9 ± 0.31 (mean ± s.e.) alleles per locus, and 14 different multilocus genotypes. The size of the alleles ranged from 112 bp (VRZAG29) to a maximum of 320 bp (VVIp60) and was consistent with the literature. As expected for clonally propagated crops, the observed heterozygosity was very high (0.87 ± 0.10), and only two loci had homozygous genotypes; although the Ho was high (i.e., above 50%) for the VVMD7 locus, the other locus (VVS5) was almost fixed. The Expected heterozygosity (He), Fixation index (F), and Shannon’s Index were, on average, 0.50 ± 0.03, −0.68 ± 0.18, and 0.78 ± 0.05, respectively. The He was greater than or equal to 0.5 in all markers except VVS5. Values were always below 0.6, suggesting that the high level of the Ho is not due to a highly uniform distribution of the alleles. This was confirmed by the estimation of the Evenness, a measure related to the ratio between the more abundant and the rarer genotypes, whose high values were overall consistent with clonal diversity. Consequently, the F value was negative for all markers except VVS5 (0.848). The Shannon’s index was close to 1.0 in VVMD7 and lower than 0.8 in the other markers.
It was not possible to find a reference ‘Greco’ genotype considering all the markers mainly because of the variability at VVMD7. Specifically, for the biallelic VRZAG21, VVS4, VRZAG29, and VVLN16 loci, we detected only one heterozygote profile in the 17 clones. A highly predominant allelic combination (i.e., 94%) was found in VVS2 and VVIC05, having four and two alleles, respectively. For VVLH54, VVS5, and VVLP60, a common heterozygote profile was present in at least 75% of the genotypes. Six genotypes were present at the tetra allelic VVMD7 locus yet, among them, two were predominant, with a 29% frequency each. Overall, the genotypic analysis is consistent with a highly heterozygote population of clonal origin. The analysis of the VVMD7 locus suggests the presence of two main groups which may have derived from an early differentiation. Finally, the SSR analysis indicated intra-clonal variability at six of the ten employed SSRs.
To illustrate the genetic resemblance among clones, we built a dendrogram using Provesti’s genetic distance and the UPGMA algorithm (Figure 2). The correlation between the distance and the cophenetic matrices was high (r = 0.85) and significant (p ≤ 0.001; Mantel test), indicating high goodness of fit. The average distance was very low (0.10 ± 0.01), and three groups of clones were indistinguishable. Specifically, the clones G16 and G13, G8, G9, and G11, and G4, G5, G12, and G14 had an identical SSR profile.

2.2.2. REMAP Analysis

To complement and extend the analysis of the genetic diversity among clones, we also used the highly polymorphic REMAP marker system based on the Gret1 LTR elements. We obtained 26 easily scorable polymorphic fragments in our population, with an average of 11.77 bands per clone. All ‘Greco’ clones had a unique molecular profile. The relationship among clones is illustrated by the hierarchical cluster obtained using the Jaccard distance and the UPGMA algorithm (Figure 3). The linear correlation between the distance and the cophenetic distance matrices was very high (0.93) and significant (p ≤ 0.001; Mantel test), indicating very high goodness of fit of the dendrogram with the dissimilarity matrix on which it was based.

2.3. Relationships between the Classification Obtained with the Different Marker Systems

To test for agreement between the hierarchical classification obtained with morphological, SSR, and REMAP data, we compared matrices using the matrix correlation coefficient as their index of concordance (by analyzing corresponding values in two half-matrices, diagonal excluded) [29]. Moreover, a comparison of cophenetic matrices was made possible using matrices obtained with the same clustering algorithm. To this aim, we first calculated a distance matrix for the quantitative data using the Gower’s index on scaled (Z-score) data and the related cophenetic matrix based on the UPGMA algorithm. The goodness of fit between the matrices relative to the quantitative morphological traits was not high 0.71 (p ≤ 0.001; Mantel test). Afterwards, we performed the concordance test on distance and cophenetic matrices (Table 5).
The pairwise correlation between the two DNA marker systems and the quantitative traits was always not significant, both in terms of described differences among clones and at the hierarchical classification level, indicating that the different systems depicted different sources of variability.

3. Discussion

Knowledge of the differences between plants belonging to the same grape variety is useful to manage germplasm, ensure homogeneity of new vineyards, and select superior clones even when considering the limitation of genetic exchange-based breeding approaches in grapevine. In this work, we evaluated the on-farm intra-varietal diversity of the traditional ‘Greco’ grapevine in the PDO area. Previous studies investigated the clonal diversity present in plants collected in ample areas [30,31,32]. Although the long history of cultivation and farmer-to-farmer exchange (as an alternative to the use of certified material) may allow for a more diverse population [33], it is also expected that the propagation of a relatively reduced number of mother plants, along with the more recent pressure to comply with production and processing rules, will restrain the clonal diversity in the PDO areas.
For the estimation of phenotypic variability, leaf and fruit descriptors have been widely used for both grapevine cultivars and clones [34,35,36,37]. Clones typically need multivariate biometric analysis to be distinguished [38]. Our work indicated heterogeneity in some morphological traits, such as rachis length and bunch compactness, which are under strong genetic control [39]. The multivariate data analysis did not reveal highly correlated influential variables, whereas a rather uniform distribution of the plants was observed, all of which suggest that separate, independent sources are responsible for the pattern of morphological variability among the clones. Although the qualitative traits were not polymorphic, differences were detected not only for the bunch but also for the berry composition. This is consistent with previous studies that indicate the high discriminating power of berry composition in distinguishing grape genotypes [7,9]. Although the variation of this trait is often limited (considering the CV value), it remains to be tested if the clonal diversity may be associated with different specific compounds and, ultimately, wine flavors [40].
We also investigated the genetic diversity of two DNA-based markers, SSR and REMAP. The diverse nature of these markers revealed a different level of variability within clones. For instance, SSRs were unable to distinguish all the genotypes, consistent with previous knowledge indicating that this marker does not often show allelic variation among grape clones [41,42]. The analysis of the co-dominant SSRs revealed a high genetic similarity, fixed heterozygosity at several loci, and the frequent presence of predominant locus-specific profiles. Taken together, these can be considered evidence of the common origin of the material. Therefore, it is possible to propose for further evaluation that the ‘Greco’ of the PDO area is of monoclonal origin [32,38]. This is also supported by the fact that divergent genotypes usually differed for very few loci. Moreover, the analysis of the most variable SSR locus (VVMD7) implied the possibility of two main groups of clones. On the other hand, REMAP markers separated all the genotypes under investigation. Considering the lack of genetic recombination and the abundance of transposable elements in the grape genome, insertion polymorphisms due to the activity of mobile elements are considered the largest source of clonal diversity in cultivated grapevine [43,44]. For instance, retrotransposon-based markers can identify polymorphisms within a single grapevine clone even when SSRs cannot resolve the existing genetic diversity [25,45,46]. Although the REMAP marker system proved to be a useful tool for fingerprinting individual plants, its traditional disadvantages are the lack of co-dominance and association with phenotypic traits. Overall, our work confirmed the usefulness of using different DNA typing methods for the analysis of the diversity among clones [30] and that their combination is probably necessary for a better understanding of their evolutionary history.
From a molecular perspective, the polymorphism of SSRs and REMAPs originate from independent mechanisms, but the agamic propagation of the germplasm should not allow different evolutionary trajectories, making possible a potential correlation. The comparison of the classification provided by the two DNA marker systems indicated that they provided different information, also considering the analysis of the cophenetic values (i.e., the similarity levels at which two entities group together during clustering). Other works based on the diversity among grapevine and varieties indicated that the cophenetic matrix of REMAP significantly correlated with that of the Inter-Single Sequence Repeats (ISSR) markers, but not with that of SSR, AFLP, SSAP, or IRAP [25]. Likewise, also the relationship between the morphological and each of the two molecular classifications was very low and not significant. The literature on agamically propagated fruit trees reports usually low or non-significant correlations between DNA-based and morphological classifications [22,47,48,49]. This is typically explained by the fact that the employed marker systems sampled a different kind of mostly uncorrelated diversity. It is very likely that the two neutral DNA marker systems principally tested a non-adaptive diversity. Nonetheless, it is also necessary to add that interesting agronomic features are fixed in grapevine through clonal selection and agamic propagation, thus ensuring the concurrent conservation of both adaptive and neutral genetic diversity. An implication that should be tested analyzing an ample panel of grapevine varieties is that it will be unlikely to infer the presence of phenotypic differences among grape clones according to the of resemblance provided by anonymous and/or neutral DNA markers. Therefore, the classification, distinction, and descent of the clones as well as the identification of a reference genotype for the ‘Greco’ variety should be based on a multidisciplinary integration of different tools.
In conclusion, this work assessed the clonal diversity of a traditional grapevine. The data indicated a limited phenotypic variability and supported a monoclonal origin of the ‘Greco’ population present the PDO area. Although our aim did not include the identification of clone-specific variants and demographic assessments (e.g., parentage reconstruction), the SSR and REMAP analysis suggested a possible division in two groups and the presence of region-specific genomic differences given the highly uncorrelated classifications. Further validations, which should also include material outside the PDO area, will be able to identify a reference molecular profile for clone authenticity and for the future identification of adaptive genetic variation.

4. Materials and Methods

4.1. Plant Material and Experimental Site

The grapevine germplasm under investigation included 17 clones of Vitis vinifera cv. ‘Greco Bianco’ sampled in different areas of Tufo village (Avellino), Italy (Table S1). During the experiment, air temperature was measured hourly at a weather station located within 2 km from the experimental sites (Figure S1).

4.2. Morphological and Berry Composition Analysis

Ten qualitative and three quantitative morphological traits were assessed on the vines of the 17 selected clones in 2019. The qualitative morphological descriptors were categorically scored according to the indications of the OIV (available at https://www.oiv.int/public/medias/2274/code-2e-edition-finale.pdf, accessed on 1 December 2022): number of consecutive tendrils on the shoots (descriptor identification code: OIV-016), blade shape in mature leaves (OIV-067), number of lobes in mature leaves (OIV-068), teeth shape in mature leaves (OIV-76), shape of base of petiole sinus in mature leaves (OIV-080), shape of the base of upper lateral sinuses in mature leaves (OIV-083-1), morphological development of sexual organs in flowers (OIV-151), bunch shape (OIV-208), number of wings in the primary bunch (OIV-209), and berry shape (OIV-223). The quantitative traits measured in this study were the number of berries per bunch, the total rachis length (main axis + wing axes), and the compactness index (calculated as the number of berries per bunch to total rachis length ratio). These parameters were measured at harvest on five bunches per clone.
Right after fruit set, the number of bunches located on eight shoots per vine was counted separately and was used to calculate the bud fertility potential index (mean number of bunches per shoot). Grape harvest was carried out when total soluble solids (TSS) in berry juice reached around 21 °Brix. TSS was measured on the extracted juice with digital refractometer (HI96811, Hanna Instruments, TX, USA). For clones G1, G2, G3, G4, G5, G6, G8, G10, G11, G13, G14, and G16, grapes were harvested on October 2 whole for G7, G9, G12, G15, and G17 6 days later (8 October). At harvest, three 20-berry samples were collected per vine and transported to the laboratory for measuring total soluble solids (TSS), pH, and titratable acidity (TA) of berry juice. TSS was measured as previously described, whereas juice TA was measured by titration with a 0.1 M NaOH solution until reaching a pH value of 8.2 and expressed as g/L tartaric acid. The pH was monitored continuously with a pH-meter (GLP 21; Crison, Alella, Barcelona, Spain). The TSS/A ratio was calculated for each sample.
The significance of the differences among clones in the quantitative morphological traits was assessed with one-way ANOVA using the Duncan as a post-hoc test for mean separation (p ≤ 0.05). Moreover, differences between clones in these traits were studied with a Principal Component Analysis (PCA) using, as input, seven original variables: no berries per bunch, rachis length, bunch compactness index, berry juice TTS, pH, TA, and the TSS/TA ratio. Statistical analyses and biplot representation were achieved using R 4.2.

4.3. Microsatellite Analysis

For DNA analysis, six young healthy leaves per plant (in three sampled replicates per clone) were harvested and immediately frozen with liquid nitrogen. Total genomic DNA was extracted using the Plant DNeasy Maxi Kit (Qiagen, Milan, Italy) as per the manufacturer’s procedure. Microsatellite analyses were carried out with ten nuclear markers (Table S2) selected from different sources [42,50,51,52]. The amplification profile was as reported [21]. PCR reactions were performed in a 20 μL volume containing 1 × reaction buffer with 1.6 mM MgCl2, 0.2 mM of each dNTP, 30 pM FAM-labelled forward SSR primer, 30 pM reverse SSR primer (Table S2), 1 unit of Taq polymerase (Promega, Milan, Italy), and 30 ng of genomic DNA. After agarose gel electrophoresis to verify the correct amplification and estimate product size, amplicons were separated with the ABI PRISM 3130 DNA Analyzer system (Applied Biosystems, Milan, Italy). Size calibration was performed with the DNA marker size ladder GenScan 500 ROX dye Size Standard (Applied Biosystems). Binning and sizing of the SSR alleles were performed using the Peak Scanner software v 2.0 (Applied Biosystems). Experiments were carried in triplicate per clone.

4.4. REMAP Analysis

For REMAP markers amplification, three primers were chosen: one designed on the Gret1 LTR region 3′-LTR (5′- GCATTTAGAAGGATTTAGCTT-3′) and two anchored primers designed on microsatellite repeats, Microsat-GA ((GA)9C) and Microsat-CT ((CT)9G). We set up two REMAP markers that combine the LTR primer with the microsatellite primers as described (Pereira et al., 2012). PCR reactions were performed in a 20 μL reaction mixture containing 20 μg of genomic DNA template, as previously described [24]. PCR products were separated via 2% agarose gel electrophoresis at 50 V for 30 min, 70V for 30 min, and then 100 V for 3 h in 1 × TAE buffer. The bands were detected using the GelRed Nucleic Acid Gel Stain (Biotium, Hayward, CA, USA). For the image acquisition and the identification of band size, the Quantity One 1-D Analysis Software (Bio-Rad Laboratories, Hercules, CA, USA) was used. Experiments were carried in triplicate per clone.

4.5. Molecular Data Analysis

The GenAlex 6.5 software [53] was used to obtain the main indices of genetic diversity based on the SSR profiles, namely the number of alleles (Na); the allelic size range in bp; the observed heterozygosity (Ho); the expected heterozygosity (He): 1 − Σpi2; the Shannon’s index of diversity (I): −1 × Σ(pi  ×  ln(pi)); the Wright’s fixation index (F): (He − Ho)/He, where for each locus, pi is the frequency of the i-th allele, Σpi2 is the sum of the squared population allele frequencies. The Poppr library [54] was used to calculate the Evenness (E) as [(1/λ) − 1]/[(eH) − 1], where 1/λ is Stoddart and Taylor’s index, and H is the Shannon diversity.
Genetic distances were calculated using the Prevosti distance [55], which ranges from 0 to 1. REMAP fragments, scored as binary markers and were used to calculate the Jaccard distance [56], which ranges from 0 to 1. Cluster analysis was performed using the unweighted pair group method with arithmetic mean (UPGMA) algorithm [56].

4.6. Concordance Test between Marker Systems

For each marker type (i.e., morphological data, SSR and REMAP), ordered parallel triangular matrices were obtained using the matsbyname R package. From independent data, pairwise comparison of distance matrices and the cophenetic matrices obtained with the same algorithm (UPGMA) was performed in R 4.2 with the Mantel test (9999 permutations) using the ADE4 package [57].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12030515/s1, Figure S1: Minimum, mean, and maximum air daily temperature (2 m); each dot represents a daily measurement; lines are smoothed conditional means surrounded by standard error of the means (grey areas). Table S1. Plant material under investigation and their location; Table S2. Microsatellite markers used for grape genotyping and their main features (linkage group, primer sequences (5′–3′), core motif, and annealing temperatures (Ta) in °C).

Author Contributions

Conceptualization, C.V., B.B., and R.A.; methodology, C.V., B.B., G.C., and R.A.; resources, R.A.; investigation: C.V., E.D.S., A.M., and E.F.; formal analysis: C.V., B.B., and G.C., writing—original draft preparation, C.V., G.C., B.B., and R.A.; writing—review and editing, C.V., G.C., B.B., and R.A.; visualization, G.C. and B.B.; supervision, C.V, B.B., and R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data that are not already contained within the article or supplementary material will be available from the corresponding author (B.B.) upon reasonable request.

Acknowledgments

We thank Raffaele Garramone and Alessia Gallucci for technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reisch, B.I.; Owens, C.L.; Cousins, P.S. Grape. In Fruit Breeding; Springer: Berlin/Heidelberg, Germany, 2012; pp. 225–262. [Google Scholar]
  2. McKey, D.; Elias, M.; Pujol, B.; Duputié, A. The evolutionary ecology of clonally propagated domesticated plants. New Phytol. 2010, 186, 318–332. [Google Scholar] [CrossRef]
  3. Riaz, S.; Garrison, K.E.; Dangl, G.S.; Boursiquot, J.-M.; Meredith, C.P. Genetic divergence and chimerism within ancient asexually propagated winegrape cultivars. J. Am. Soc. Hortic. Sci. 2002, 127, 508–514. [Google Scholar] [CrossRef] [Green Version]
  4. Vondras, A.M.; Minio, A.; Blanco-Ulate, B.; Figueroa-Balderas, R.; Penn, M.A.; Zhou, Y.; Seymour, D.; Ye, Z.; Liang, D.; Espinoza, L.K. The genomic diversification of grapevine clones. BMC Genom. 2019, 20, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Pelsy, F. Molecular and cellular mechanisms of diversity within grapevine varieties. Heredity 2010, 104, 331–340. [Google Scholar] [CrossRef] [PubMed]
  6. Chessa, I.; Nieddu, G. Analysis of diversity in the fruit tree genetic resources from a Mediterranean island. Genet. Resour. Crop Evol. 2005, 52, 267–276. [Google Scholar] [CrossRef]
  7. Leão, P.C.D.S.; Cruz, C.D.; Motoike, S.Y. Genetic diversity of table grape based on morphoagronomic traits. Sci. Agric. 2011, 68, 42–49. [Google Scholar] [CrossRef]
  8. Vafaee, Y.; Ghaderi, N.; Khadivi, A. Morphological variation and marker-fruit trait associations in a collection of grape (Vitis vinifera L.). Sci. Hortic. 2017, 225, 771–782. [Google Scholar] [CrossRef]
  9. Abiri, K.; Rezaei, M.; Tahanian, H.; Heidari, P.; Khadivi, A. Morphological and pomological variability of a grape (Vitis vinifera L.) germplasm collection. Sci. Hortic. 2020, 266, 109285. [Google Scholar] [CrossRef]
  10. Kupe, M.; Sayinci, B.; Demir, B.; Ercisli, S.; Aslan, K.A.; Gundesli, M.A.; Baron, M.; Sochor, J. Multivariate Analysis Approaches for Dimension and Shape Discrimination of Vitis vinifera Varieties. Plants 2021, 10, 1528. [Google Scholar] [CrossRef]
  11. Muganu, M.; Dangl, G.; Aradhya, M.; Frediani, M.; Scossa, A.; Stover, E. Ampelographic and DNA characterization of local grapevine accessions of the Tuscia area (Latium, Italy). Am. J. Enol. Vitic. 2009, 60, 110–115. [Google Scholar] [CrossRef]
  12. Villano, C.; Aiese Cigliano, R.; Esposito, S.; D’Amelia, V.; Iovene, M.; Carputo, D.; Aversano, R. DNA-Based Technologies for Grapevine Biodiversity Exploitation: State of the Art and Future Perspectives. Agronomy 2022, 12, 491. [Google Scholar] [CrossRef]
  13. Mondini, L.; Noorani, A.; Pagnotta, M.A. Assessing plant genetic diversity by molecular tools. Diversity 2009, 1, 19–35. [Google Scholar] [CrossRef] [Green Version]
  14. Grassi, F.; Labra, M.; Imazio, S.; Spada, A.; Sgorbati, S.; Scienza, A.; Sala, F. Evidence of a secondary grapevine domestication centre detected by SSR analysis. Theor. Appl. Genet. 2003, 107, 1315–1320. [Google Scholar] [CrossRef] [PubMed]
  15. Labra, M.; Failla, O.; Fossati, T.; Castiglione, S.; Scienza, A.; Sala, F. Phylogenetic analysis of grapevine cv. Ansonica growing on the island of Giglio, Italy, by AFLP and SSR markers. Vitis 1999, 38, 161–166. [Google Scholar]
  16. Pastore, C.; Fontana, M.; Raimondi, S.; Ruffa, P.; Filippetti, I.; Schneider, A. Genetic characterization of grapevine varieties from Emilia-Romagna (northern Italy) discloses unexplored genetic resources. Am. J. Enol. Vitic. 2020, 71, 334–343. [Google Scholar] [CrossRef]
  17. De Lorenzis, G.; Imazio, S.; Biagini, B.; Failla, O.; Scienza, A. Pedigree Reconstruction of the Italian Grapevine Aglianico (Vitis vinifera L.) from Campania. Mol. Biotechnol. 2012, 54, 634–642. [Google Scholar] [CrossRef]
  18. Crespan, M.; Calò, A.; Giannetto, S.; Sparacio, A.; Storchi, P.; Costacurta, A. Sangiovese’ and ‘Garganega’ are two key varieties of the Italian grapevine assortment evolution. Vitis 2008, 47, 97–104. [Google Scholar]
  19. Raimondi, S.; Tumino, G.; Ruffa, P.; Boccacci, P.; Gambino, G.; Schneider, A. DNA-based genealogy reconstruction of Nebbiolo, Barbera and other ancient grapevine cultivars from northwestern Italy. Sci. Rep. 2020, 10, 15782. [Google Scholar] [CrossRef]
  20. Imazio, S.; Labra, M.; Grassi, F.; Scienza, A.; Failla, O. Chloroplast microsatellites to investigate the origin of grapevine. Genet. Resour. Crop Evol. 2006, 53, 1003–1011. [Google Scholar] [CrossRef]
  21. Villano, C.; Carputo, D.; Frusciante, L.; Santoro, X.; Aversano, R. Use of SSR and retrotransposon-based markers to interpret the population structure of native grapevines from southern Italy. Mol. Biotechnol. 2014, 56, 1011–1020. [Google Scholar] [CrossRef]
  22. Muccillo, L.; Gambuti, A.; Frusciante, L.; Iorizzo, M.; Moio, L.; Raieta, K.; Rinaldi, A.; Colantuoni, V.; Aversano, R. Biochemical features of native red wines and genetic diversity of the corresponding grape varieties from Campania region. Food Chem. 2014, 143, 506–513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kalendar, R.; Schulman, A.H. IRAP and REMAP for retrotransposon-based genotyping and fingerprinting. Nat. Protoc. 2006, 1, 2478–2484. [Google Scholar] [CrossRef] [PubMed]
  24. Kalendar, R.; Grob, T.; Regina, M.; Suoniemi, A.; Schulman, A. IRAP and REMAP: Two new retrotransposon-based DNA fingerprinting techniques. Theor. Appl. Genet. 1999, 98, 704–711. [Google Scholar] [CrossRef]
  25. D’Onofrio, C.; De Lorenzis, G.; Giordani, T.; Natali, L.; Cavallini, A.; Scalabrelli, G. Retrotransposon-based molecular markers for grapevine species and cultivars identification. Tree Genet. Genomes 2010, 6, 451–466. [Google Scholar] [CrossRef]
  26. Pereira, H.S.; Barão, A.; Delgado, M.; Morais-Cecílio, L.; Viegas, W. Genomic analysis of grapevine retrotransposon 1 (Gret1) in Vitis vinifera. Theor. Appl. Genet. 2005, 111, 871–878. [Google Scholar] [CrossRef]
  27. Kobayashi, S.; Goto-Yamamoto, N.; Hirochika, H. Retrotransposon-induced mutations in grape skin color. Science 2004, 304, 982. [Google Scholar] [CrossRef]
  28. Viala, P.; Vermorel, V. Traité Général de Viticulture: Ampélographie; Masson et Cie: Paris, France, 1901–1910; Volume 7. [Google Scholar]
  29. Rohlf, F.J. Consensus indices for comparing classifications. Math. Biosci. 1982, 59, 131–144. [Google Scholar] [CrossRef]
  30. Zombardo, A.; Meneghetti, S.; Morreale, G.; Calò, A.; Costacurta, A.; Storchi, P. Study of Inter-and Intra-varietal Genetic Variability in Grapevine Cultivars. Plants 2022, 11, 397. [Google Scholar] [CrossRef]
  31. Gonçalves, E.; Martins, A. Efficient Assessment and Large-Scale Conservation of Intra-Varietal Diversity of Ancient Grapevine Varieties: Case Study Portugal. Plants 2022, 11, 1917. [Google Scholar] [CrossRef]
  32. Vignani, R.; Bowers, J.; Meredith, C. Microsatellite DNA polymorphism analysis of clones of Vitis vinifera ‘Sangiovese’. Sci. Hortic. 1996, 65, 163–169. [Google Scholar] [CrossRef]
  33. Banilas, G.; Korkas, E.; Kaldis, P.; Hatzopoulos, P. Olive and grapevine biodiversity in Greece and Cyprus–a review. Clim. Chang. Intercropping Pest Control Benef. Microorg. 2009, 2, 401–428. [Google Scholar]
  34. Santiago, J.L.; Boso, S.; Gago, P.; Alonso-Villaverde, V.; Martínez, M.C. Molecular and ampelographic characterisation of Vitis vinifera L. “Albariño”, “Savagnin Blanc” and “Caíño Blanco” shows that they are different cultivars. Span. J. Agric. Res. 2007, 5, 333. [Google Scholar] [CrossRef]
  35. Sabir, A.; Tangolar, S.; Buyukalaca, S.; Kafkas, S. Ampelographic and molecular diversity among grapevine (Vitis spp.) cultivars. Czech J. Genet. Plant Breed. 2009, 45, 160–168. [Google Scholar] [CrossRef]
  36. Gago, P.; Santiago, J.-L.; Boso, S.; Alonso-Villaverde, V.; Grando, M.S.; Martínez, M.C. Biodiversity and Characterization of Twenty-two Vitis vinifera L. Cultivars in the Northwestern Iberian Peninsula. Am. J. Enol. Vitic. 2009, 60, 293–301. [Google Scholar] [CrossRef]
  37. Meneghetti, S.; Poljuha, D.; Frare, E.; Costacurta, A.; Morreale, G.; Bavaresco, L.; Calò, A. Inter-and intra-varietal genetic variability in Malvasia cultivars. Mol. Biotechnol. 2012, 50, 189–199. [Google Scholar] [CrossRef]
  38. Razi, M.; Amiri, M.E.; Darvishzadeh, R.; Doulati Baneh, H.; Alipour, H.; Martínez-Gómez, P. Assessment of genetic diversity of cultivated and wild Iranian grape germplasm using retrotransposon-microsatellite amplified polymorphism (REMAP) markers and pomological traits. Mol. Biol. Rep. 2020, 47, 7593–7606. [Google Scholar] [CrossRef]
  39. Tello, J.; Ibáñez, J. What do we know about grapevine bunch compactness? A state-of-the-art review. Aust. J. Grape Wine Res. 2018, 24, 6–23. [Google Scholar] [CrossRef]
  40. Šuklje, K.; Antalick, G.; Buica, A.; Langlois, J.; Coetzee, Z.A.; Gouot, J.; Schmidtke, L.M.; Deloire, A. Clonal differences and impact of defoliation on Sauvignon blanc (Vitis vinifera L.) wines: A chemical and sensory investigation. J. Sci. Food Agric. 2016, 96, 915–926. [Google Scholar] [CrossRef]
  41. Cipriani, G.; Marrazzo, M.T.; Di Gaspero, G.; Pfeiffer, A.; Morgante, M.; Testolin, R. A set of microsatellite markers with long core repeat optimized for grape (Vitis spp.) genotyping. BMC Plant Biol. 2008, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  42. Laucou, V.; Lacombe, T.; Dechesne, F.; Siret, R.; Bruno, J.-P.; Dessup, M.; Dessup, T.; Ortigosa, P.; Parra, P.; Roux, C. High throughput analysis of grape genetic diversity as a tool for germplasm collection management. Theor. Appl. Genet. 2011, 122, 1233–1245. [Google Scholar] [CrossRef]
  43. Carrier, G.; Le Cunff, L.; Dereeper, A.; Legrand, D.; Sabot, F.; Bouchez, O.; Audeguin, L.; Boursiquot, J.-M.; This, P. Transposable Elements Are a Major Cause of Somatic Polymorphism in Vitis vinifera L. PLoS ONE 2012, 7, e32973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. The French–Italian Public Consortium for Grapevine Genome Characterization. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 2007, 449, 463–467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Hocquigny, S.; Pelsy, F.; Dumas, V.; Kindt, S.; Heloir, M.C.; Merdinoglu, D. Diversification within grapevine cultivars goes through chimeric states. Genome 2004, 47, 579–589. [Google Scholar] [CrossRef]
  46. Moncada, X.; Pelsy, F.; Merdinoglu, D.; Hinrichsen, P. Genetic diversity and geographical dispersal in grapevine clones revealed by microsatellite markers. Genome 2006, 49, 1459–1472. [Google Scholar] [CrossRef]
  47. Corrado, G.; La Mura, M.; Ambrosino, O.; Pugliano, G.; Varricchio, P.; Rao, R. Relationships of Campanian olive cultivars: Comparative analysis of molecular and phenotypic data. Genome 2009, 52, 692–700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Basile, B.; Mataffo, A.; Forlani, M.; Corrado, G. Diversity in Morphometric, Pomological, and Fruit-Quality Traits of Apricot (Prunus armeniaca) Traditional Varieties: Implications for Landrace Differentiation at Regional Scale. Diversity 2022, 14, 608. [Google Scholar] [CrossRef]
  49. Manco, R.; Basile, B.; Capuozzo, C.; Scognamiglio, P.; Forlani, M.; Rao, R.; Corrado, G. Molecular and phenotypic diversity of traditional European plum (Prunus domestica L.) germplasm of Southern Italy. Sustainability 2019, 11, 4112. [Google Scholar] [CrossRef] [Green Version]
  50. Sefc, K.M.; Regner, F.; Turetschek, E.; Glössl, J.; Steinkellner, H. Identification of microsatellite sequences in Vitis riparia and their applicability for genotyping of different Vitis species. Genome 1999, 42, 367–373. [Google Scholar] [CrossRef]
  51. Thomas, M.; Scott, N. Microsatellite repeats in grapevine reveal DNA polymorphisms when analysed as sequence-tagged sites (STSs). Theor. Appl. Genet. 1993, 86, 985–990. [Google Scholar] [CrossRef]
  52. Merdinoglu, D.; Butterlin, G.; Bevilacqua, L.; Chiquet, V.; Adam-Blondon, A.-F.; Decroocq, S. Development and characterization of a large set of microsatellite markers in grapevine (Vitis vinifera L.) suitable for multiplex PCR. Mol. Breed. 2005, 15, 349–366. [Google Scholar] [CrossRef]
  53. Peakall, R.; Smouse, P.E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 2006, 6, 288–295. [Google Scholar] [CrossRef]
  54. Kamvar, Z.N.; Tabima, J.F.; Grünwald, N.J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2014, 2, e281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Prevosti, A.; Ocana, J.; Alonso, G. Distances between populations of Drosophila subobscura, based on chromosome arrangement frequencies. Theor. Appl. Genet. 1975, 45, 231–241. [Google Scholar] [CrossRef]
  56. Romesburg, C. Cluster Analysis for Researchers; Lulu. com: Morrisville, NC, USA, 2004. [Google Scholar]
  57. Thioulouse, J.; Chessel, D.; Olivier, J.-M. ADE-4: A multivariate analysis and graphical display software. Stat. Comput. 1997, 7, 75–83. [Google Scholar] [CrossRef]
Figure 1. Bi-plot of the 17 ‘Greco’ clones (from G1 to G17) for the first principal components extracted by the PCA. The vectors indicate the original variable loadings. BN—number of berries per bunch; RL—Rachis length; CI—bunch compactness index; TSS—berry juice total soluble solids; pH—berry juice pH; TA—berry juice titratable acidity.
Figure 1. Bi-plot of the 17 ‘Greco’ clones (from G1 to G17) for the first principal components extracted by the PCA. The vectors indicate the original variable loadings. BN—number of berries per bunch; RL—Rachis length; CI—bunch compactness index; TSS—berry juice total soluble solids; pH—berry juice pH; TA—berry juice titratable acidity.
Plants 12 00515 g001
Figure 2. Dendrogram of the ‘Greco’ clones (G) based on their SSR profile. The dendrogram was built with the UPGMA algorithm using Provesti’s genetic distance (reported in the scale bar on the bottom).
Figure 2. Dendrogram of the ‘Greco’ clones (G) based on their SSR profile. The dendrogram was built with the UPGMA algorithm using Provesti’s genetic distance (reported in the scale bar on the bottom).
Plants 12 00515 g002
Figure 3. Dendrogram of the ‘Greco’ clones (G) based on the REMAP profile. The dendrogram was built with the UPGMA algorithm using the Jaccard genetic distance (reported in the scale bar on the bottom).
Figure 3. Dendrogram of the ‘Greco’ clones (G) based on the REMAP profile. The dendrogram was built with the UPGMA algorithm using the Jaccard genetic distance (reported in the scale bar on the bottom).
Plants 12 00515 g003
Table 1. Bud potential fertility and bunch morphological traits. In each column, different statistical letters indicate statistically significant groups according to a one-way ANOVA followed by Duncan post-hoc test (p ≤ 0.05). The Coefficient of Variation (C.V.) among clones was reported for each measured parameter.
Table 1. Bud potential fertility and bunch morphological traits. In each column, different statistical letters indicate statistically significant groups according to a one-way ANOVA followed by Duncan post-hoc test (p ≤ 0.05). The Coefficient of Variation (C.V.) among clones was reported for each measured parameter.
CloneBud Potential Fertility IndexComponents of Bunch Compactness
No Berries
Per Bunch
Rachis Length
(cm)
CI
(No Berries/
cm Rachis)
G11.10 ± 0.14 a108 ± 23 ab21.1 ± 2.6 a4.99 ± 0.62 e
G21.00 ± 0.14 a114 ± 19 ab16.3 ± 0.8 abcd6.86 ± 0.85 bcde
G31.00 ± 0.15 a86 ± 13 b14.4 ± 2.0 abcd6.33 ± 1.42 bcde
G41.18 ± 0.18 a78 ± 13 b12.4 ± 2.0 cd6.53 ± 0.97 bcde
G51.27 ± 0.18 a77 ± 14 b16.0 ± 1.3 abcd4.77 ± 0.79 e
G61.19 ± 0.14 a110 ± 33 ab14.8 ± 2.7 abcd7.07 ± 1.14 bcde
G71.33 ± 0.19 a137 ± 27 ab17.8 ± 3.5 abc7.71 ± 0.03 abcd
G81.29 ± 0.167 a129 ± 15 ab19.9 ± 0.8 ab6.43 ± 0.55 bcde
G91.29 ± 0.17 a91 ± 20 b17.3 ± 2.4 abcd5.15 ± 0.47 de
G101.29 ± 0.19 a134 ± 1 ab19.3 ± 1.3 ab7.07 ± 0.45 bcde
G111.17 ± 0.14 a133 ± 21 ab14.0 ± 2.9 bcd9.90 ± 0.77 a
G121.20 ± 0.17 a87 ± 12 b10.6 ± 1.8 d8.28 ± 0.72 abc
G131.12 ± 0.17 a118 ± 15 ab15.9 ± 0.8 abcd7.37 ± 0.69 bcde
G141.13 ± 0.17 a157 ± 16 a18.0 ± 0.9 abc8.64 ± 0.50 ab
G151.29 ± 0.167 a91 ± 14 b15.9 ± 2.1 abcd5.91 ± 0.91 cde
G161.38 ± 0.18 a96 ± 5 b20.0 ± 1.9 ab4.98 ± 0.65 e
G171.19 ± 0.16 a91 ± 13 b15.9 ± 1.9 abcd5.76 ± 0.51 cde
C.V. (%)9.121.916.921.1
Table 2. Berry composition of the different clones at harvest. The harvest date for clones G7, G9, G12, G15, and G17 (8 October 2019) was six days later than the other clones (2 October 2019). In each column, different statistical letters indicate statistically significant groups according to a one-way ANOVA followed by Duncan post-hoc test (p ≤ 0.05). The Coefficient of Variation (C.V.) among clones was reported for each measured parameter.
Table 2. Berry composition of the different clones at harvest. The harvest date for clones G7, G9, G12, G15, and G17 (8 October 2019) was six days later than the other clones (2 October 2019). In each column, different statistical letters indicate statistically significant groups according to a one-way ANOVA followed by Duncan post-hoc test (p ≤ 0.05). The Coefficient of Variation (C.V.) among clones was reported for each measured parameter.
CloneTSS
(°Brix)
pHTA
(g/L Tartaric Acid)
TSS/TA
G124.3 ± 0.3 a3.00 ± 0.06 cd10.4 ± 0.1 abcd2.35 ± 0.03 de
G222.6 ± 0.2 bcd2.96 ± 0.04 d14.1 ± 0.1 a1.60 ± 0.01 i
G321.6 ± 0.2 de3.01 ± 0.02 cd13.4 ± 2.2 ab1.62 ± 0.01 i
G422.3 ± 0.3 cd3.15 ± 0.01 abc8.9 ± 0.1 cd2.52 ± 0.04 bc
G521.1 ± 0.4 ef3.07 ± 0.01 abcd10.1 ± 0.4 abcd2.10 ± 0.04 f
G622.8 ± 0.5 bc3.21 ± 0.09 a9.3 ± 2.1 bcd2.45 ± 0.06 cd
G723.6 ± 0.2 ab3.07 ± 0.01 abcd9.0 ± 0.3 cd2.62 ± 0.02 b
G822.2 ± 0.3 cd3.03 ± 0.04 bcd11.0 ± 0.4 abcd2.03 ± 0.03 fg
G920.2 ± 0.4 f3.13 ± 0.04 abc8.0 ± 1.3 cd2.55 ± 0.05 bc
G1021.2 ± 0.2 ef3.16 ± 0.01 ab8.6 ± 0.1 cd2.47 ± 0.026 c
G1123.2 ± 0.2 bc3.15 ± 0.03 abc9.9 ± 0.9 abcd2.34 ± 0.02 de
G1222.8 ± 0.3 bcd3.21 ± 0.04 a7.2 ± 0.1 d3.16 ± 0.046 a
G1321.1 ± 0.6 ef3.16 ± 0.07 ab10.8 ± 2.7 abcd1.95 ± 0.06 gh
G1422.5 ± 0.1 cd3.13 ± 0.07 abc12 ± 2.4 abc1.87 ± 0.01 h
G1522.5 ± 0.7 bcd3.17 ± 0.02 ab8.9 ± 0.1 cd2.55 ± 0.08 bc
G1621.1 ± 0.4 ef3.10 ± 0.01 abc9.3 ± 0.3 bcd2.27 ± 0.04 e
G1724.4 ± 0.3 a3.13 ± 0.04 abc9.6 ± 0.9 bcd2.54 ± 0.03 bc
C.V. (%)5.22.318.117.2
Table 3. Eigenvalue, percent of explained variance, and correlation with the eight original variables of the six principal components (Dim.) extracted. BN—number of berries per bunch; RL—Rachis length, CI—bunch compactness index; TSS—berry juice total soluble solids; pH—berry juice pH; TA—berry juice titratable acidity; TSS/TA—TSS to TA ratio.
Table 3. Eigenvalue, percent of explained variance, and correlation with the eight original variables of the six principal components (Dim.) extracted. BN—number of berries per bunch; RL—Rachis length, CI—bunch compactness index; TSS—berry juice total soluble solids; pH—berry juice pH; TA—berry juice titratable acidity; TSS/TA—TSS to TA ratio.
Principal
Component
EigenvalueVariance Explained
(%)
Correlation
BNRLCITSSpHTATSS/TA
Dim. 12.1035.00.28 NS0.72 **−0.34 NS0.05 NS−0.91 ***0.75 ***0.28 NS
Dim. 21.7929.90.85 ***−0.04 NS0.89***0.46 NS0.11 NS0.23 NS0.85 NS
Dim. 31.0617.70.44 NS0.64 **−0.10 NS−0.40 NS0.30 NS−0.44 NS0.44 NS
Dim. 40.8814.7−0.06 NS0.25 NS−0.24 NS0.79 **−0.01 NS−0.37 NS−0.06 NS
Dim. 50.152.50.01 NS0.04 NS−0.12 NS0.06 NS0.28 NS0.23 NS0.01 NS
Dim. 60.010.2−0.07 NS0.06 NS0.07 NS0.01 NS0.01 NS0.01 NS−0.07 NS
**, ***, and NS indicate significant differences at p ≤ 0.01, p ≤ 0.001, and not significant (p > 0.05) correlation.
Table 4. Main genetic indexes of diversity in the analyzed population of ‘Greco’ clones. Na—number of alleles; ASR—allelic size range; Ho—observed heterozygosity: He—expected heterozygosity; I—Shannon’s index of diversity; F—Wright’s fixation index; and E—evenness.
Table 4. Main genetic indexes of diversity in the analyzed population of ‘Greco’ clones. Na—number of alleles; ASR—allelic size range; Ho—observed heterozygosity: He—expected heterozygosity; I—Shannon’s index of diversity; F—Wright’s fixation index; and E—evenness.
LocusNaASRHoHeIFE
VVS24128–1261.000.560.92−0.800.83
VVMD74244–2500.650.581.05−0.120.75
VRZAG212188–2001.000.500.69−1.001.00
VVS42166–1721.000.500.69−1.001.00
VRZAG292112–1181.000.500.69−1.001.00
VVLN162149–1571.000.500.69−1.001.00
VVLH543163–1671.000.530.82−0.880.90
VVS54147–1670.060.280.600.780.47
VVIC052162–1661.000.500.69−1.001.00
VVLP604312–3201.000.560.94−0.780.83
Table 5. Pairwise comparison of molecular and morphological data. Numbers in parentheses indicate p-value (Mantel test). The results below (respectively, above) the diagonal indicated the correlation coefficient between the distance matrices (resp. the cophenetic value matrices).
Table 5. Pairwise comparison of molecular and morphological data. Numbers in parentheses indicate p-value (Mantel test). The results below (respectively, above) the diagonal indicated the correlation coefficient between the distance matrices (resp. the cophenetic value matrices).
QuantitativeSSRREMAP
Quantitative-−0.034 (0.54)−0.074 (0.64)
SSR−0.120 (0.76)-−0.084 (0.64)
REMAP−0.150 (0.82)0.028 (0.43)-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Villano, C.; Corrado, G.; Basile, B.; Di Serio, E.; Mataffo, A.; Ferrara, E.; Aversano, R. Morphological and Genetic Clonal Diversity within the ‘Greco Bianco’ Grapevine (Vitis vinifera L.) Variety. Plants 2023, 12, 515. https://doi.org/10.3390/plants12030515

AMA Style

Villano C, Corrado G, Basile B, Di Serio E, Mataffo A, Ferrara E, Aversano R. Morphological and Genetic Clonal Diversity within the ‘Greco Bianco’ Grapevine (Vitis vinifera L.) Variety. Plants. 2023; 12(3):515. https://doi.org/10.3390/plants12030515

Chicago/Turabian Style

Villano, Clizia, Giandomenico Corrado, Boris Basile, Ermanno Di Serio, Alessandro Mataffo, Elvira Ferrara, and Riccardo Aversano. 2023. "Morphological and Genetic Clonal Diversity within the ‘Greco Bianco’ Grapevine (Vitis vinifera L.) Variety" Plants 12, no. 3: 515. https://doi.org/10.3390/plants12030515

APA Style

Villano, C., Corrado, G., Basile, B., Di Serio, E., Mataffo, A., Ferrara, E., & Aversano, R. (2023). Morphological and Genetic Clonal Diversity within the ‘Greco Bianco’ Grapevine (Vitis vinifera L.) Variety. Plants, 12(3), 515. https://doi.org/10.3390/plants12030515

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop