Next Article in Journal
Demographic Processes Allow Echinochloa crus-galli to Compensate Seed Losses by Seed Predation
Next Article in Special Issue
Genetic Diversity of Soybeans (Glycine max (L.) Merr.) with Black Seed Coats and Green Cotyledons in Korean Germplasm
Previous Article in Journal
Potassium Application Positively Modulates Physiological Responses of Cocoa Seedlings to Drought Stress
Previous Article in Special Issue
Molecular Diversity within a Mediterranean and European Panel of Tetraploid Wheat (T. turgidum subsp.) Landraces and Modern Germplasm Inferred Using a High-Density SNP Array
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 11
Issue 3
10.3390/agronomy11030564
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genetic Diversity and Structure of Tomato Landraces from the Campania Region (Southern Italy) Uncovers a Distinct Population Identity

by
Martina Caramante
,
Youssef Rouphael
and
Giandomenico Corrado
*
Department of Agricultural Sciences, University of Naples Federico II, 80055 Portici, Italy
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(3), 564; https://doi.org/10.3390/agronomy11030564
Submission received: 9 February 2021 / Revised: 12 March 2021 / Accepted: 12 March 2021 / Published: 17 March 2021
(This article belongs to the Special Issue Old Germplasm for New Needs: Managing Crop Genetic Resources)
Browse Figures
Versions Notes

Abstract

:
Italy is one of the main producers and processors of tomato and it is considered a secondary center of diversity. In some areas, such as the Campania region (Southern Italy), a range of traditional tomato landraces is still cultivated. The distinction of this heritage germplasm is often based only on folk taxonomy and a more comprehensive definition and understanding of its genetic identity is needed. In this work, we compared a set of 15 local landraces (representative of traditional fruit types) to 15 widely used contemporary varieties, using 14 fluorescent Simple Sequence Repeat (SSR) markers. Each of the accessions possessed a unique molecular profile and overall landraces had a genetic diversity comparable to that of the contemporary varieties. The genetic diversity, multivariate, and population structure analysis separated all the genotypes according to the pre-defined groups, indicating a very reduced admixture and the presence of a differentiated (regional) population of landraces. Our work provides solid evidence for implementing conservation actions and paves the way for the creation of a premium regional brand that goes beyond the individual landrace names of the Campania region known throughout the world.

1. Introduction

Italy is the leading European nation for tomatoes, considering the cultivated area, the total yield, and the output of the processing industry, which is second in the world only to that of the USA (https://tinyurl.com/y2n6pg5n; accessed 10 March 2021). Moreover, the tomato is a dominant element of the Italian gastronomy and a symbol of a culinary identity that was globally spread by Italian immigrants in the XX century [1].
The tomato was domesticated in the Americas and was introduced to Italy soon after the Columbian exchange [2]. However, its culinary use started in Italy since the second half of the 17th century when, despite the inability to contribute significantly to the caloric intake of the farmers, it was recognized the gustatory value of tomato as (cooked) sauce, condiment and in combination with other foods [3]. In that period, a regional diversification in the use of tomato has also begun, influenced by the political, cultural, and linguistic division of Italy at the time. The success of tomatoes in Italy was not uniform and even today, there are areas where relatively few tomatoes are consumed and cultivated. In the 18th century, the tomato was widespread in the Italian gastronomy especially in the Kingdom of Naples and Sicily (current Southern Italy), at that time under the Spanish rule of the House of Bourbon [4]. Tomato was consumed by commoners both fresh and as a slow cooked sauce, also because of the influence of the Spanish gastronomy [5]. It is especially in this century that varieties with different fruit shapes started to have different destinations of use, a phenomenon that has led to the selection of specific tomato types (e.g., the San Marzano) later employed for industrial canning [3]. After WWII, Italian production reached its peak (in terms of commercial product and dedicated area) during the 1970s, after the widespread diffusion of contemporary varieties [6]. This expansion was accompanied by a strong decline in landrace cultivation, which until recently are predominantly grown in amateur family gardens and in small-size farms for local markets.
In tomato, as well as in other crops, landraces are being recognized as an important genetic material for their yield stability, cultural value, adaptability to low-input and organic farming, resilience to stress, and fruit quality traits [7,8,9,10]. Tomato landraces are of increasing scientific interest also as a source of adaptive traits in the face of climate change and for sustainable fruit quality [11,12]. Moreover, landraces are of commercial importance following the rediscovery of local food systems and the promotion of short food supply chains [13]. These values led to the introduction of EU quality schemes, which promoted the cultivation of some traditional landraces (i.e., San Marzano and Pomodorino Vesuviano, for the Campania region), also because of the possibility of limiting the unauthorized production or marketing of goods using such a name. Regrettably, besides some globally recognized names, most of this germplasm is neglected and at risk of extinction.
In Italy the loss of tomato landrace has been significant especially in the regions where tomato was most diffused [14]. Until the 1950s, tomato breeding was predominantly based on mass and pedigree selection using local germplasm, while hybridization was later exploited [6]. Therefore, in germplasm collections, traditional varietal names can be associated with both landraces and old selected lines [15]. Another problem is the presence of so-called contaminant genotypes, deriving from the fortuitous or voluntary introduction in cultivation of morphologically related varieties, deriving also from crosses and/or on-farm seed propagation of hybrid cultivars [16,17,18]. Moreover, the intrinsic lack of uniformity and unequivocal descriptors of traditional landraces creates the opportunity for spurious, uncontrolled seed marketing, especially from on-line and unprofessional retailers [19,20,21]. Landraces are mainly distinguished by their geographical location and folk taxonomy. Often, tomato landraces or landrace groups are recognized by farmers in terms of fruit type [22]. Nonetheless, the presence of a clear genetic distinction to support landraces’ authenticity is not often available, increasing the need for a molecular discrimination [15].
Southern Italy is considered a secondary center of diversity for tomato [23,24], and in view of the above-mentioned reasons, it is not surprising that the Campania region, and more generally, Southern Italy is a rich reservoir of crop landraces [25]. During the transition to intensive modern agriculture, the seed industry in Southern Italy was in its infancy and it is possible to presume that most of the local germplasm has not been involved in formal breeding programs, for instance, by being engaged in the selection of (pre-)breeding material. Therefore, still today the pool of tomato landraces that are locally cultivated may be a genetically different group from contemporary varieties. In this work, to test this hypothesis, we analyzed and compared a population of local landraces (e.g., collected from farms of the Campania region) with an equivalent set of contemporary cultivars, picked among the most diffuse in cultivation. Among the different DNA molecular markers available for tomato [26], we used Simple Sequence Repeats (SSRs) because these highly variable multiallelic DNA markers are considered most suitable to reveal recent demographic events [27,28]. Previous works were dedicated to specific (fruit) types, single or national collections of tomato landraces [29,30,31,32,33,34]. Our specific aims were: (i) to evaluate the level and the distribution of the molecular diversity among and within the collected landraces and a set of common contemporary varieties, going beyond the single tomato landrace or landrace group evaluation; (ii) to quantify the genetic differentiation and highlight locus-specific differences that may contribute to the maintenance of a genetic diversity. Our work provided a solid and experimentally validated basis for the uniqueness of tomato landraces of a specific region, with the prospect of strengthening market opportunities and adding value for cultivation.

2. Materials and Methods

2.1. Plant Material and DNA Isolation

This work was carried out on a selection of 15 tomato (Solanum lycopersicum L.) landraces of the Campania region, mainly collected in the area to produce the “Pomodoro San Marzano dell’Agro Sarnese-nocerino DOP” and of the “Pomodorino del Piennolo del Vesuvio DOP”, two Protected Denomination of Origin (PDO) EU label schemes for the tomatoes of the Campania regions. Seeds were multiplied before analysis at the Department of Agricultural Sciences, University of Naples Federico II. This population was compared to an equal number of cultivars, from different seed companies, selected among the most diffused in the Campania and Puglia regions in the 2009–2014 period according to the cultivated area. As possible output for cladistic classification, we also analyzed two tomato wild relatives, namely S. habrochaites and S. neorickii. The list, code, and main characteristics of the fruit are reported in the Supplementary Table S1. Seeds were germinated on moist filter paper in Petri dishes and then plantlets were transferred to polystyrene seed trays in a growing chamber (24 °C). For DNA isolation, the first two true leaves from two different plants per accession were collected and immediately frozen in liquid nitrogen. Total DNA isolation and quantification were performed as previously described, starting from approximately 150 mg of finely ground frozen tissue [35].

2.2. SSR Analysis

To genotype the plant material, we used 14 Simple Sequence Repeat (SSR) loci that represent different repeat classes [36,37]. The loci and their main characteristics are presented in the Supplementary Table S2. Polymerase Chain Reaction (PCR) amplifications were performed independently on two plants per accession. Each reaction was assembled using 20 ng genomic DNA template, 1.5 mM MgCl2, 100 μM deoxyribonucleotide-triphosphates, 0.2 μM fluorescently labeled forward primer and unlabeled reverse primer (Supplementary Table S2), and 0.5 U GoTaq DNA polymerase (Promega, Milan, Italy), for a final volume of 25 μL. PCR reactions were performed in a Verity Pro thermal cycler (Thermo Fisher, Milan, Italy). The PCR cycle was as previously described using the annealing temperature indicated in the Supplementary Table S2 [38]. For allelic discrimination, amplicons were analyzed by fluorescent capillary electrophoresis (dye-labels are reported in Supplementary Table S2), with the POP6 polymer (Thermo Fisher) in an ABI Prism 3100 Genetic Analyzer (Thermo Fisher). Signal peak height and allele sizes were calculated using the ABI Prism GeneMapper software v. 4.1 (Thermo Fisher, Milan, Italy) calibrated on the GeneScan 500Liz dye size standard (Thermo Fisher). Values were rounded to integer, and if necessary scaled, based on the SSR core size (Supplementary Table S2), minimizing for each locus the average offset of the scaled alleles within the instrumental resolution of the DNA separation (±1 bp) (Supplementary Figure S1).

2.3. Data Analysis

For each SSR locus, we calculated the number of different alleles (for the whole collection and for S. lycopersicum), the observed heterozygosity (Ho; number of heterozygotes/N), the number of effective alleles (Ne; = 1/(Sum pi2)), the Polymorphic Information Content (PIC; 1 − Sum (pi2), equivalent to the expected heterozygosity), the Shannon’s Information Index (I; −1 × Sum (pi × Ln (pi)), and the Wright’s fixation index, where N is number of individuals, pi is the frequency of the i-th allele of a locus, and Sum (pi2) represents the sum of the squared pi. These calculations were performed with Genalex 6.5 [39]. Evenness of alleles was obtained with poppr [40]. This R-library was also employed to calculate pairwise genetic distances between genotypes according to the Prevosti coefficient [41]. Hierarchical clustering based on the unweighted pair-group method with arithmetic averages (UPGMA) algorithm, cophenetic correlation and tree visualization were performed in R using the factoextra package [42]. The Agglomerative Cluster (AC) value was computed as reported [43]. The genotype accumulation curve was built by boxplotting the number of multilocus genotypes obtained for an increasing number of SSRs (from 1 to plateau, 10 for cultivars and 12 for landraces) for each predefined population. The data distribution for each number of SSRs was obtained by random sampling loci (n = 10,000) [40]. As model-free ordination technique we used the Principal Coordinate Analysis (PCoA), carried out via covariance matrix with distance data standardization. The Analysis of Molecular Variance (AMOVA) was performed using the pairwise distance matrix with 9999 permutations to test for significance. These calculations were carried out using the Genalex 6.5 software [39]. To validate the population structure, the Fst and Gst between S. lycopersicum populations were calculated per locus and globally and statistically evaluated against a null distribution obtained from 10,000 permutations on alleles, using the MSA 4.05 software [44].

3. Results

3.1. Genetic Diversity of the Germplasm under Investigation

The analysis of the genetic diversity of the entire set of genotypes (n = 32) was carried out with 14 SSRs. Differences within each accession were not detected. All loci were polymorphic yet, the locus LEcaa001 was polymorphic only between wild species and S. lycopersicum and therefore, it was excluded for comparisons between landraces and contemporary cultivars. Main genetic parameters of the tomato collection under investigation are presented in Table 1. We detected 75 different alleles for an average of 5.35 alleles per locus. Large differences were present in the number of alleles per locus, which ranged from 2 (LELE25 and LEcaa001) to 11 (LEEF1Aa). The latter was the most diverse locus also considering the number of effective alleles, a measure of diversity weighted for allele frequencies. There was not a significant difference (p > 0.05; Student’s t-test) between the allelic richness in S. lycopersicum subgroups, with an average of 3.0 (respectively, 2.8) alleles per locus in the landraces (resp., contemporary varieties). The observed heterozygosity (Ho) was low (on average, 0.17) but, as expected, highly affected by the tomato groups, because heterozygote genotypes were present at one locus for the open-pollinated landraces (Legaa003) and at 12 SSRs for the contemporary (mostly hybrid) varieties. Nonetheless, two polymorphic loci (LELE25 and LEct001) were fixed in the whole germplasm collection and lacked heterozygotes. A significant correlation between the number of alleles per locus and its observed heterozygosity was present (p = 0.005; Spearman’s Rho). Considerable differences among SSRs were also present in the PIC, which ranged from 0.03 (LEcaa001) to 0.81 (LEEF1Aa), the locus with the highest number of alleles. However, the most diverse loci were LE20592 and LEtat002 considering the evenness of distribution of alleles for each SSR.
To evaluate the relationship between the genotypes, pairwise genetic distances were used to build a dendrogram (Figure 1). In this analysis, the two wild species served as possible outgroup. Genetic distances were summarized in an heatmap to simultaneously visualize how the genetic relatedness varies between genotypes and clusters of genotypes (Supplementary Figure S2). The cophenetic correlation was high (Pearson’s product-moment correlation: 0.80) and significant (p < 0.001). The Agglomerative Cluster value was 0.65. The dendrogram illustrated that each genotype had a distinct SSR profile and duplicate accessions were not present in our landrace collection. Moreover, the dendrogram indicated that the genotypes could be clearly separated in biologically meaningful groups at k = 4. The most taxonomically distant wild species had a distinctive position at the highest hierarchical node (k = 2). Within the S. lycopersicum germplasm, the landraces and the contemporary varieties agglomerated in two separate clusters, suggesting the presence of two genetically distinct populations (Figure 1).
Within the landraces, some small groups of genotypes with a similar fruit shape (Supplementary Table S1) could be also identified, such as the MUR-MAC-MOR subcluster (with “San Marzano”- type fruit and collected in the related PDO area), the ANO-AGN (“pomodorino” type, both collected in the valley of the Sarno river), and TLI-TLS (“tondo” type, but originating from different cultivation areas).
Finally, before testing for a possible genetic structure, we verified if the number of employed SSRs is suitable to capture the diversity in our germplasm. To this aim, we built genotype accumulation curves, considering the landrace or the contemporary varieties group as defined by our a priori classification. The result indicated that, also taking into account the possible specific polymorphisms within each group, the employed number of SSR loci is above the minimum number needed to fully capture the diversity in our sub-groups (Figure 2). For every number of randomly sampled SSRs, there was not a statistical difference in the average number of identified multilocus genotypes considering landraces or cultivars (p > 0.05, Student’s t-test).

3.2. Genetic Structure Analysis of the S. lycopersium Genotypes (Landraces and Contemporary Varieties)

To infer a population subdivision between contemporary varieties and landraces, we performed a Principal Coordinate Analysis. To achieve this goal, we excluded the wild species because of their distant relatedness and limited number of samples. The first two components explained 39.7% of the total variance, which implies that the different SSRs were useful in sampling a mostly uncorrelated, locus-specific genetic variation. The scatter plot of PC1 and PC2 values shows that the samples belonging to the two pre-defined populations are well separated along both PC1 and PC2 (Figure 3). Moreover, the samples of each population were similarly spread on the two-dimensional plane, an indication of a comparable level of diversity.
Consequently, we performed an Analysis of Molecular Variance to quantify the level of genetic differentiation between the predefined groups (Table 2). A main contribution to the genetic variance in our germplasm originated from the difference between samples in each subpopulation and not predominantly from all samples, as in the case of a panmictic population. Moreover, the differentiation between landraces and contemporary cultivars was high, with a 22% percent of the molecular variation present between the two sub-populations (Table 2). For these reasons, the AMOVA analysis provided evidence to support a population structure.
Finally, to validate the AMOVA results, we calculated two widely used measures of population differentiation, Fst and Gst (according to Nei), per locus and globally [45], using a distribution obtained by allele permutations for the statistical significance (Table 3). The data indicated that amount of the genetic differentiation between populations is significant. Interestingly, there were clear differences among loci, with five loci with statistically significant high Fst values. In the presence of a genetic structure, this suggests the possible presence of locus-specific adaptive molecular variation.

4. Discussion

The tomato is characterized by limited genetic diversity, due to its evolutionary history [23,46]. The breeding sector is dominated by hybrids that have progressively replaced local accessions. Previous molecular characterizations focused on specific landraces, fruit classes or national collections of tomato landraces [29,30,31,32,33,34]. In this work, we tested the hypothesis that a set of regional tomato landraces of Southern Italy represents a genetically differentiated group compared to contemporary varieties because of their origin and possible non-engagement of the breeding sector. The cladistic classification of the genotypes indicated that every accession had a unique profile. The polymorphic SSRs confirmed to be suitable in distinguish closely related tomato landraces [18]. The most informative SSR was LEEF1Aa, probably because of a composite core motif. This locus was also previously described as highly informative in tomato [47]. As expected, wild relatives had also SSR alleles other from S. lycopersicum [33]. The total number of alleles was similar or higher than in other published works [36,37,48,49], indicating a relevant allelic diversity for both landraces and contemporary varieties. It should be added that for landraces, the molecular variation is expected to include both human-driven selection and adaptation to local conditions [50]. Moreover, the selected varieties were chosen to represent different tomato market classes, which associate to different patterns of genetic variation [51]. Also in another paper based on landraces and contemporary tomato varieties [52], the average number of SSR alleles was judged higher compared to the tomato literature, even if the loci employed are not the same of our work. While the presence of distinctive profile for cultivars is not surprising, synonymous cases and intra-varietal variability has been reported for landraces of tomato and other species [53,54,55]. The explorative hierarchical cluster indicated also biologically coherent and highly balanced clusters of cultivars and landraces at k = 4, although the optimal number of clusters was not statistically evaluated. However, the cophenetic correlation was high (≥ 0.80), a reliable signal of a population structure [56]. The relatively low AC value of the UPGMA tree can be justified considering that even for highly structured populations, the “chaining” effect, unique genotypes, and small groups of accessions (together with two or more large groups) can strongly lower (less than 0.6) this index [57]. Finally, also considering the output group, the dendrogram analysis implied a similar level of diversity between the two pre-defined populations in terms of pairwise distances, which should be further explored.
To scan for a genetic structure, we performed a PCoA mainly because it does not require assumptions that may not be hold true for a selected germplasm collection (e.g., HW or linkage equilibrium) [51], and the intrinsic issues of model-based testing for k = 2 [58]. The analysis indicated that the S. lycopersicum genotypes are well separated according to their pre-defined population [51]. A complete distinction between landraces and contemporary varieties and a confirmation of the UPGMA clustering have not been always found [52]. It is necessary to add that work analyzed a nation-wide group of landraces from geographically diverse collections. Landraces of different secondary centers of diversity (i.e., Spain and Italy) were genetically differentiated [59], suggesting that the complete separation observed in our study could be also due to a more restricted geographic area of origin of our germplasm. The level of population differentiation was high [60], and in line with previous analysis of tomato germplasm that underwent conventional breeding [51]. An important implication is that despite the wide diffusion of contemporary varieties, the admixture between our local landraces and cultivars is very low, considering that we analyzed cultivars that were among the most cultivated in the collection sites. Farmers tend to save seeds from one year to another by choosing the most representative fruits, favoring therefore a maintenance breeding for the fruit morphology (e.g., size and shape), phenological characteristics (e.g., time of flowering and maturity), and quality traits (e.g., color and flavor) they attribute to their accession (Zeven, 2000). We also detected some significant differences among loci in highlighting a genetic differentiation. High Fst values (compared to neutral loci) are typically suggestive of genomic regions that are under divergent selection. Nonetheless, at each locus, Fst values are also influenced by heterozygosity and the mutation rate, two features that are more variable among multiallelic SSRs than in bi-allelic markers. For these reasons, the data motivate further genomic scans to pinpoint the effects of the adaptive and/or breeding processes that led to a differentiation between tomato landraces and contemporary varieties [61].

5. Conclusions

Our study provides a first comparative assessment of the diversity and population structure in a geographically specific collection of tomato landraces from a secondary center of diversity. The data indicated that the landraces constitute a genetically distinct population from common commercial varieties, in addition to having a historical origin and a locally recognized gastronomic identity. Therefore, our work provides a robust justification for implementing measures for in situ and ex situ conservation actions, as well as for creating premium regional brands that can even go beyond the worldwide known, individual landrace names of the Campania region (e.g., San Marzano, piennolo/Pomodorino Vesuviano). Moreover, multivariate analyses indicated that the landraces are characterized by a good, and to a relevant extent, specific diversity. Finally, our work encourages further multidisciplinary studies to unravel the genetic factors responsible for the adaptation to specific environments and potential quality traits of our landraces.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/11/3/564/s1, Figure S1: an example of the electropherograms at the locus LE20592 showing the allele peaks. For each genotype (indicated on the left; see Table S1 for the code), the height of the peak refers to the scale on the right-hand side (RFU). The top bar indicates the size reference range provided by the genotyping software, used to calibrate the peak data points to their DNA size according to the internal size standard. The rounded dimension of each fragment (bp) is indicated in the rectangle below. Figure S2: heatmap of the pairwise genetic distance (Prevosti) between genotypes (codes are reported in Table S1). The bar on the right side shows the color scale adopted to represent distances, from deep blue (0) to white (1). Table S1: List of the genotypes under investigation and main features of their fruit shape. Table S2: SSR loci employed in this study and their main features.

Author Contributions

Conceptualization, G.C.; formal analysis, Y.R. and G.C.; investigation, M.C. and G.C.; writing—original draft preparation, G.C.; writing—review and editing, Y.R. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the project “Laboratorio di genomica per l’innovazione e la valorizzazione della filiera pomodoro” (MIUR-art 12 del DM 593/00, Cod. DM17732).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gentilcore, D. Taste and the tomato in Italy: A transatlantic history. Food Hist. 2009, 7, 125–139. [Google Scholar] [CrossRef]
  2. Bergougnoux, V. The history of tomato: From domestication to biopharming. Biotechnol. Adv. 2014, 32, 170–189. [Google Scholar] [CrossRef]
  3. Gentilcore, D. Pomodoro! A History of the Tomato in Italy; Columbia University Press: New York, NY, USA, 2010. [Google Scholar]
  4. Astarita, T. The Italian Baroque Table: Cooking and Entertaining from the Golden Age of Naples. Medieval and Renaissance Texts and Studies; ACMRS Press: Tempe, AZ, USA, 2014; Volume 459. [Google Scholar]
  5. Latini, A. Lo Scalco Alla Moderna, 1993 ed.; Biblioteca Culinaria: Lodi, Italy, 1694; Volume 2. [Google Scholar]
  6. Acciarri, N.; Di Candilo, M.; Sanguineti, C.; Soressi, G.P. Tomato. In Italian Contribution to Plant Genetics and Breeding; Scarascia Mugnozza, G.T., Pagnotta, M.A., Eds.; Università della Tuscia: Viterbo, Italy, 1998; pp. 501–509. [Google Scholar]
  7. Sumalan, R.M.; Ciulca, S.I.; Poiana, M.A.; Moigradean, D.; Radulov, I.; Negrea, M.; Crisan, M.E.; Copolovici, L.; Sumalan, R.L. The Antioxidant Profile Evaluation of Some Tomato Landraces with Soil Salinity Tolerance Correlated with High Nutraceuticaland Functional Value. Agronomy 2020, 10, 500. [Google Scholar] [CrossRef] [Green Version]
  8. Landi, S.; De Lillo, A.; Nurcato, R.; Grillo, S.; Esposito, S. In-field study on traditional Italian tomato landraces: The constitutive activation of the ROS scavenging machinery reduces effects of drought stress. Plant Physiol. Biochem. 2017, 118, 150–160. [Google Scholar] [CrossRef]
  9. Digilio, M.C.; Corrado, G.; Sasso, R.; Coppola, V.; Iodice, L.; Pasquariello, M.; Bossi, S.; Maffei, M.E.; Coppola, M.; Pennacchio, F. Molecular and chemical mechanisms involved in aphid resistance in cultivated tomato. New Phytol. 2010, 187, 1089–1101. [Google Scholar] [CrossRef]
  10. Corrado, G.; Rao, R. Towards the genomic basis of local adaptation in landraces. Diversity 2017, 9, 51. [Google Scholar] [CrossRef] [Green Version]
  11. Fullana-Pericàs, M.; Conesa, M.À.; Douthe, C.; El Aou-ouad, H.; Ribas-Carbó, M.; Galmés, J. Tomato landraces as a source to minimize yield losses and improve fruit quality under water deficit conditions. Agric. Water Manag. 2019, 223, 105722. [Google Scholar] [CrossRef]
  12. Carillo, P.; Kyriacou, M.C.; El-Nakhel, C.; Pannico, A.; dell’Aversana, E.; D’Amelia, L.; Colla, G.; Caruso, G.; De Pascale, S.; Rouphael, Y. Sensory and functional quality characterization of protected designation of origin ‘Piennolo del Vesuvio’cherry tomato landraces from Campania-Italy. Food Chem. 2019, 292, 166–175. [Google Scholar] [CrossRef] [PubMed]
  13. Casals, J.; Rull, A.; Segarra, J.; Schober, P.; Simó, J. Participatory Plant Breeding and the Evolution of Landraces: A Case Study in the Organic Farms of the Collserola Natural Park. Agronomy 2019, 9, 486. [Google Scholar] [CrossRef] [Green Version]
  14. Grandillo, S.; Mustilli, A.C.; Parisi, M.; Morelli, G.; Giordano, I.; Bowler, C. Tecniche avanzate per la valutazione qualitativa del pomodoro: Il caso Campania. Agroindustria 2004, 3, 151–159. [Google Scholar]
  15. Villa, T.C.C.; Maxted, N.; Scholten, M.; Ford-Lloyd, B. Defining and identifying crop landraces. Plant Genet. Resour. 2005, 3, 373–384. [Google Scholar] [CrossRef] [Green Version]
  16. Cleveland, D.A.; Daniela, S.; Smith, S.E. A biological framework for understanding farmers’ plant breeding. Econ. Bot. 2000, 54, 377–394. [Google Scholar] [CrossRef]
  17. Rao, R.; Corrado, G.; Bianchi, M.; Di Mauro, A. (GATA) 4 DNA fingerprinting identifies morphologically characterized ‘San Marzano’tomato plants. Plant Breed. 2006, 125, 173–176. [Google Scholar] [CrossRef]
  18. Caramante, M.; Rao, R.; Monti, L.M.; Corrado, G. Discrimination of ‘San Marzano’accessions: A comparison of minisatellite, CAPS and SSR markers in relation to morphological traits. Sci. Hortic. 2009, 120, 560–564. [Google Scholar] [CrossRef]
  19. Muñoz-Falcón, J.E.; Vilanova, S.; Plazas, M.; Prohens, J. Diversity, relationships, and genetic fingerprinting of the Listada de Gandía eggplant landrace using genomic SSRs and EST-SSRs. Sci. Hortic. 2011, 129, 238–246. [Google Scholar] [CrossRef] [Green Version]
  20. Muthoni, J.; Nyamongo, D. Seed systems in Kenya and their relationship to on-farm conservation of food crops. J. New Seeds 2008, 9, 330–342. [Google Scholar] [CrossRef]
  21. Louwaars, N.P. Plant breeding and diversity: A troubled relationship? Euphytica 2018, 214, 1–9. [Google Scholar] [CrossRef] [Green Version]
  22. Zeven, A.C. Landraces: A review of definitions and classifications. Euphytica 1998, 104, 127–139. [Google Scholar] [CrossRef]
  23. Bauchet, G.; Causse, M. Genetic diversity in tomato (Solanum lycopersicum) and its wild relatives. Genet. Divers. Plants 2012, 8, 134–162. [Google Scholar]
  24. García-Martínez, S.; Andreani, L.; Garcia-Gusano, M.; Geuna, F.; Ruiz, J.J. Evaluation of amplified fragment length polymorphism and simple sequence repeats for tomato germplasm fingerprinting: Utility for grouping closely related traditional cultivars. Genome 2006, 49, 648–656. [Google Scholar] [CrossRef]
  25. Hammer, K.; Laghetti, G.; Perrino, P. Collection of plant genetic resources in South Italy, 1988. Die Kult. 1989, 37, 401–414. [Google Scholar] [CrossRef]
  26. Causse, M.; Giovannoni, J.; Bouzayen, M.; Zouine, M. The Tomato Genome; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  27. Zhivotovsky, L.A.; Rosenberg, N.A.; Feldman, M.W. Features of evolution and expansion of modern humans, inferred from genomewide microsatellite markers. Am. J. Hum. Genet. 2003, 72, 1171–1186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Cornuet, J.-M.; Ravigné, V.; Estoup, A. Inference on population history and model checking using DNA sequence and microsatellite data with the software DIYABC (v1. 0). BMC Bioinform. 2010, 11, 1–11. [Google Scholar] [CrossRef] [Green Version]
  29. Terzopoulos, P.; Bebeli, P. DNA and morphological diversity of selected Greek tomato (Solanum lycopersicum L.) landraces. Sci. Hortic. 2008, 116, 354–361. [Google Scholar] [CrossRef]
  30. Henareh, M.; Dursun, A.; Abdollahi-Mandoulakani, B.; Haliloğlu, K. Assessment of genetic diversity in tomato landraces using ISSR markers. Genetika 2016, 48, 25–35. [Google Scholar] [CrossRef]
  31. Parisi, M.; Aversano, R.; Graziani, G.; Ruggieri, V.; Senape, V.; Sigillo, L.; Barone, A. Phenotypic and molecular diversity in a collection of ‘Pomodoro di Sorrento’Italian tomato landrace. Sci. Hortic. 2016, 203, 143–151. [Google Scholar] [CrossRef]
  32. Castellana, S.; Ranzino, L.; Beritognolo, I.; Cherubini, M.; Luneia, R.; Villani, F.; Mattioni, C. Genetic characterization and molecular fingerprint of traditional Umbrian tomato (Solanum lycopersicum L.) landraces through SSR markers and application for varietal identification. Genet. Resour. Crop Evol. 2020, 67, 1807–1820. [Google Scholar] [CrossRef]
  33. Mazzucato, A.; Papa, R.; Bitocchi, E.; Mosconi, P.; Nanni, L.; Negri, V.; Picarella, M.E.; Siligato, F.; Soressi, G.P.; Tiranti, B. Genetic diversity, structure and marker-trait associations in a collection of Italian tomato (Solanum lycopersicum L.) landraces. Theor. Appl. Genet. 2008, 116, 657–669. [Google Scholar] [CrossRef]
  34. Cattáneo, R.A.; McCarthy, A.N.; Feingold, S.E. Evidence of genetic diversity within Solanum Lycopersicum L.‘Platense’landrace and identification of various subpopulations. Genet. Resour. Crop Evol. 2020, 67, 2057–2069. [Google Scholar] [CrossRef]
  35. Scarano, D.; Rubio, F.; Ruiz, J.J.; Rao, R.; Corrado, G. Morphological and genetic diversity among and within common bean (Phaseolus vulgaris L.) landraces from the Campania region (Southern Italy). Sci. Hortic. 2014, 180, 72–78. [Google Scholar] [CrossRef]
  36. Smulders, M.; Bredemeijer, G.; Rus-Kortekaas, W.; Arens, P.; Vosman, B. Use of short microsatellites from database sequences to generate polymorphisms among Lycopersicon esculentum cultivars and accessions of other Lycopersicon species. Theor. Appl. Genet. 1997, 94, 264–272. [Google Scholar] [CrossRef]
  37. He, C.; Poysa, V.; Yu, K. Development and characterization of simple sequence repeat (SSR) markers and their use in determining relationships among Lycopersicon esculentum cultivars. Theor. Appl. Genet. 2003, 106, 363–373. [Google Scholar] [CrossRef]
  38. Verdone, M.; Rao, R.; Coppola, M.; Corrado, G. Identification of zucchini varieties in commercial food products by DNA typing. Food Control 2018, 84, 197–204. [Google Scholar] [CrossRef]
  39. 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]
  40. 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]
  41. 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]
  42. Kassambara, A. Practical Guide to Cluster Analysis in R: Unsupervised Machine Learning; CreateSpace Independent Publishing Platform: Scotts Valley, CA, USA, 2017. [Google Scholar]
  43. Kaufman, L.; Rousseeuw, P.J. Finding Groups in Data: An Introduction to Cluster Analysis; John Wiley & Sons: Hoboken, NJ, USA, 2009; Volume 344. [Google Scholar]
  44. Dieringer, D.; Schlötterer, C. Microsatellite analyser (MSA): A platform independent analysis tool for large microsatellite data sets. Mol. Ecol. Notes 2003, 3, 167–169. [Google Scholar] [CrossRef]
  45. Verity, R.; Nichols, R.A. What is genetic differentiation, and how should we measure it—GST, D, neither or both? Mol. Ecol. 2014, 23, 4216–4225. [Google Scholar] [CrossRef]
  46. Mariette, S.; Tavaud, M.; Arunyawat, U.; Capdeville, G.; Millan, M.; Salin, F. Population structure and genetic bottleneck in sweet cherry estimated with SSRs and the gametophytic self-incompatibility locus. BMC Genet. 2010, 11, 1–13. [Google Scholar] [CrossRef] [Green Version]
  47. Miskoska-Milevska, E.; Popovski, Z.T.; Nestorovski, T. Usefulness of a locus LeEF1a in the genetic differentiation of tomato varieties. J. Agric. Food Environ. Sci. 2018, 72, 56–61. [Google Scholar]
  48. Scarano, D.; Rao, R.; Masi, P.; Corrado, G. SSR fingerprint reveals mislabeling in commercial processed tomato products. Food Control 2015, 51, 397–401. [Google Scholar] [CrossRef]
  49. Ruiz, J.J.; García-Martínez, S.; Picó, B.; Gao, M.; Quiros, C.F. Genetic variability and relationship of closely related Spanish traditional cultivars of tomato as detected by SRAP and SSR markers. J. Am. Soc. Hortic. Sci. 2005, 130, 88–94. [Google Scholar] [CrossRef] [Green Version]
  50. Corrado, G.; Piffanelli, P.; Caramante, M.; Coppola, M.; Rao, R. SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato. BMC Genom. 2013, 14, 1–14. [Google Scholar] [CrossRef] [Green Version]
  51. Sim, S.-C.; Van Deynze, A.; Stoffel, K.; Douches, D.S.; Zarka, D.; Ganal, M.W.; Chetelat, R.T.; Hutton, S.F.; Scott, J.W.; Gardner, R.G. High-density SNP genotyping of tomato (Solanum lycopersicum L.) reveals patterns of genetic variation due to breeding. PLoS ONE 2012, 7, e45520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Gonias, E.D.; Ganopoulos, I.; Mellidou, I.; Bibi, A.C.; Kalivas, A.; Mylona, P.V.; Osanthanunkul, M.; Tsaftaris, A.; Madesis, P.; Doulis, A.G. Exploring genetic diversity of tomato (Solanum lycopersicum L.) germplasm of genebank collection employing SSR and SCAR markers. Genet. Resour. Crop Evol. 2019, 66, 1295–1309. [Google Scholar] [CrossRef]
  53. Boccacci, P.; Aramini, M.; Valentini, N.; Bacchetta, L.; Rovira, M.; Drogoudi, P.; Silva, A.; Solar, A.; Calizzano, F.; Erdoğan, V. Molecular and morphological diversity of on-farm hazelnut (Corylus avellana L.) landraces from southern Europe and their role in the origin and diffusion of cultivated germplasm. Tree Genet. Genomes 2013, 9, 1465–1480. [Google Scholar] [CrossRef]
  54. Mercati, F.; Longo, C.; Poma, D.; Araniti, F.; Lupini, A.; Mammano, M.M.; Fiore, M.C.; Abenavoli, M.R.; Sunseri, F. Genetic variation of an Italian long shelf-life tomato (Solanum lycopersicon L.) collection by using SSR and morphological fruit traits. Genet. Resour. Crop. Evol. 2015, 62, 721–732. [Google Scholar] [CrossRef]
  55. Foroni, I.; Baptista, C.; Monteiro, L.; Lopes, M.S.; Mendonça, D.; Melo, M.; Carvalho, C.; Monjardino, P.; Lopes, D.J.; da Câmara Machado, A. The use of microsatellites to analyze relationships and to decipher homonyms and synonyms in Azorean apples (Malus × domestica Borkh.). Plant Syst. Evol. 2012, 298, 1297–1313. [Google Scholar] [CrossRef]
  56. Odong, T.; Van Heerwaarden, J.; Jansen, J.; van Hintum, T.J.; Van Eeuwijk, F. Determination of genetic structure of germplasm collections: Are traditional hierarchical clustering methods appropriate for molecular marker data? Theor. Appl. Genet. 2011, 123, 195–205. [Google Scholar] [CrossRef] [Green Version]
  57. Johnson, R.A.; Wichern, D.W. Applied Multivariate Statistical Analysis; Prentice Hall: Upper Saddle River, NJ, USA, 2002; Volume 5. [Google Scholar]
  58. Janes, J.K.; Miller, J.M.; Dupuis, J.R.; Malenfant, R.M.; Gorrell, J.C.; Cullingham, C.I.; Andrew, R.L. The K = 2 conundrum. Mol. Ecol. 2017, 26, 3594–3602. [Google Scholar] [CrossRef] [Green Version]
  59. García-Martínez, S.; Corrado, G.; Ruiz, J.J.; Rao, R. Diversity and structure of a sample of traditional Italian and Spanish tomato accessions. Genet. Resour. Crop Evol. 2013, 60, 789–798. [Google Scholar] [CrossRef]
  60. Frankham, R.; Ballou, S.E.J.D.; Briscoe, D.A.; Ballou, J.D. Introduction to Conservation Genetics; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  61. Kim, M.; Jung, J.-K.; Shim, E.-J.; Chung, S.-M.; Park, Y.; Lee, G.P.; Sim, S.-C. Genome-wide SNP discovery and core marker sets for DNA barcoding and variety identification in commercial tomato cultivars. Sci. Hortic. 2021, 276, 109734. [Google Scholar] [CrossRef]
Agronomy 11 00564 g001 550
Figure 1. Dendrogram of the tomato genotypes based on the UPGMA cluster using the Prevosti (1975) genetic distance (the scale is reported at the bottom). Different colors denote the subclusters at k = 4 namely, red garnet: S. neorickii; purple: S. habrochaites; blue navy: S. lycopersicum landraces; azure: S. lycopersicum contemporary varieties. The codes of the genotypes are available in Table S1.
Figure 1. Dendrogram of the tomato genotypes based on the UPGMA cluster using the Prevosti (1975) genetic distance (the scale is reported at the bottom). Different colors denote the subclusters at k = 4 namely, red garnet: S. neorickii; purple: S. habrochaites; blue navy: S. lycopersicum landraces; azure: S. lycopersicum contemporary varieties. The codes of the genotypes are available in Table S1.
Agronomy 11 00564 g001
Agronomy 11 00564 g002 550
Figure 2. Genotype accumulation analysis in the two pre-defined S. lycopersium populations. In each panel, the boxplots summarize the descriptive statistics relative to the quartiles of the number of multilocus genotypes obtained by randomly sampling loci without replacement (n = 10,000). The number of sampled SSR loci (from 1 to 12) is indicated in the top dark gray bar. Dots represent outliers (i.e., values outside 1.5 times the interquartile range above the upper and below the lower quartile). Contemporary varieties (respectively, landraces) boxplots are in deep salmon (resp., cyan) color.
Figure 2. Genotype accumulation analysis in the two pre-defined S. lycopersium populations. In each panel, the boxplots summarize the descriptive statistics relative to the quartiles of the number of multilocus genotypes obtained by randomly sampling loci without replacement (n = 10,000). The number of sampled SSR loci (from 1 to 12) is indicated in the top dark gray bar. Dots represent outliers (i.e., values outside 1.5 times the interquartile range above the upper and below the lower quartile). Contemporary varieties (respectively, landraces) boxplots are in deep salmon (resp., cyan) color.
Agronomy 11 00564 g002
Agronomy 11 00564 g003 550
Figure 3. Principal Coordinate Analysis scatter plot of the tomato landraces and contemporary varieties under investigation. Each genotype is represented by a colored circle (deep salmon: contemporary cultivars; cyan: landraces). PC1 (respectively, PC2) explained 22.7 % (resp. 17.0%) of the variance.
Figure 3. Principal Coordinate Analysis scatter plot of the tomato landraces and contemporary varieties under investigation. Each genotype is represented by a colored circle (deep salmon: contemporary cultivars; cyan: landraces). PC1 (respectively, PC2) explained 22.7 % (resp. 17.0%) of the variance.
Agronomy 11 00564 g003
Table
Table 1. Main genetic indices estimated with the SSR.
Table 1. Main genetic indices estimated with the SSR.
SSRNaNaSlNeHoPICIFEvALC
LE20592533.350.250.701.300.640.888
LELE25221.070.000.060.151.000.442
LE21085532.340.190.571.110.670.668
LEEF1Aa11115.220.160.811.990.810.6613
LEta003853.390.250.711.450.650.759
LEat002421.950.130.490.850.740.716
LEga003641.780.280.440.940.360.507
LEtat002533.350.250.701.300.640.889
LEaat002532.420.130.591.040.790.787
LEcaa001211.030.030.030.08−0.020.382
LEaat007532.420.130.591.040.790.777
LEct001432.320.000.570.981.000.785
LEctt001652.050.380.511.050.270.567
LEta015763.800.190.741.540.740.779
Legend: Na: total number of alleles; NaSl: total number of alleles in S. lycopersicum; Ne: number of effective alleles; Ho: observed heterozygosity; PIC: polymorphic information content; I = Shannon’s Information Index; F: fixation index; Ev: evenness; ALC: allelic combinations.
Table
Table 2. Analysis of the Molecular Variance (AMOVA) table statistics.
Table 2. Analysis of the Molecular Variance (AMOVA) table statistics.
Source of VariationdfSSMSVar%
Between Pops131.58331.5830.88422
Between samples within Pop28141.6675.0601.93830
Within samples3035.5001.1831.18348
Total59208.750 4.006100
Legend: df, degree of freedom; SS: sum of squares; MS: Mean squares; Var: estimated variance component; %: percentage of genetic variation.
Table
Table 3. Statistical analysis of the level of genetic differentiation between landraces and contemporary cultivars according to the S. lycopersicum polymorphic SSR loci.
Table 3. Statistical analysis of the level of genetic differentiation between landraces and contemporary cultivars according to the S. lycopersicum polymorphic SSR loci.
Fstp ValueGstp Value
All loci0.220.0000.150.000
LE205920.340.0000.230.000
LE210850.080.1020.150.002
LEaat0070.050.2370.050.239
LEat0020.040.2840.050.284
LEctt0010.200.0010.600.000
LEEF1Aa0.150.0020.110.002
LEga0030.120.0310.080.031
LELE250.001.0000.031.000
LEta0030.080.0750.070.074
LEta0030.080.0750.070.074
LEta0150.150.0080.110.006
LEtat0020.050.2370.050.239
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Caramante, M.; Rouphael, Y.; Corrado, G. The Genetic Diversity and Structure of Tomato Landraces from the Campania Region (Southern Italy) Uncovers a Distinct Population Identity. Agronomy 2021, 11, 564. https://doi.org/10.3390/agronomy11030564

AMA Style

Caramante M, Rouphael Y, Corrado G. The Genetic Diversity and Structure of Tomato Landraces from the Campania Region (Southern Italy) Uncovers a Distinct Population Identity. Agronomy. 2021; 11(3):564. https://doi.org/10.3390/agronomy11030564

Chicago/Turabian Style

Caramante, Martina, Youssef Rouphael, and Giandomenico Corrado. 2021. "The Genetic Diversity and Structure of Tomato Landraces from the Campania Region (Southern Italy) Uncovers a Distinct Population Identity" Agronomy 11, no. 3: 564. https://doi.org/10.3390/agronomy11030564

APA Style

Caramante, M., Rouphael, Y., & Corrado, G. (2021). The Genetic Diversity and Structure of Tomato Landraces from the Campania Region (Southern Italy) Uncovers a Distinct Population Identity. Agronomy, 11(3), 564. https://doi.org/10.3390/agronomy11030564

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