1. Introduction
Tomato (
Solanum lycopersicum L.) is one of the most highly appraised and versatile vegetable crops globally and is widely cultivated for fresh-fruit consumption, grown as an industrial/cash crop or even for processed (concentrate, juice, paste, powder, soup or sauce) products [
1]. Globally, tomato harvests can yield more than 180 million tons [
2] indicating its significance in terms of production across countries. Moreover, tomatoes can be grown either in a field or in greenhouses, while soilless cultivation schemes enjoy increasing popularity. Such crop production systems’ flexibility allows for worldwide productivity, even in marginal or desert areas [
3]. From an alimentary supply standpoint, tomato fruits are a significant source of vitamins, fibers, sugars, and essential minerals for the human diet [
4]. Moreover, besides its universal farming and financial standing as a crop, tomato is also a pre-eminent model system for fundamental and applied plant genetic studies focusing on the ripening process [
5], secondary metabolism pathways [
6] and immunity against biotic stress [
7,
8].
Wild ancestral species are believed to be native to the Pacific shore of South America. Consequently, tomato crop wild relatives (CWRs) are adapted to an assortment of edaphoclimatic environments spanning from the sea level up to the highlands of the Andean sierra, and thrive in arid, brackish water or rainy conditions [
1]. Despite the long period of time since tomato plants were introduced into the Old World, domestication schemes and routes are still largely controversial, oscillating among a Peruvian or a Mexican center of origin. Nonetheless, historical evidence dictates that the dissemination of tomato occurred eastwards and was attributed to the conquistadors’ explorations [
9].
For that matter, early introduced genotypes (up to the nineteenth century) were mostly open-pollinated and on-farm propagated at small local scales. As a result, novel varieties emerged via spontaneous mutations and/or recombination of pre-existing genetic diversity [
10]. The bulked selection mode and the concurrent breeding from multiple farmers resulted in a plethora of tomato types and rich germplasm. Nonetheless, since tomatoes are mainly self-fertilized, intercrosses among different plants were fairly infrequent while seedlings resembled a parental phenotype. This, in turn, resulted in upholding fixed tomato populations delineated as ‘heirlooms’ which are inimitable in their size, shape, colour, taste and aroma [
1].
Nowadays, the fruit weight across modern tomato varieties range from about 10 g in acorn and cherry type tomatoes to more than 450 g in some beefsteak fresh tomatoes [
11]. To this regard, it has been established that tomato CWRs present a diverse array of characteristic small-size fruits [
12]; which nonetheless have prominent levels of valuable nutritional elements like lycopene, and elevated soluble solid content (SSC) [
13]. Still, the deliberate selection of genotypes, domestication and intensive breeding aiming at the enhancement of fruit size and durability, resulted in alternating fruit characteristics and diversification of fruit weight, shape and color [
14,
15]. Regrettably, the demand from markets for big, homogenous fruits with long postharvest life promoted the production of varieties and hybrids with such characteristics, causing a genetic bottleneck in terms of variability. Such alterations resulted in a general downgrade of nutritional properties, which are presently somewhat inferior in modern cultivars and hybrids mostly due to the increase of tomato fruit size and intensive cultivation crop schemes [
16]. Moreover, the selection of a restricted number of genotypes has caused a genetic bottleneck in terms of variability, and currently less than 10% of the total genetic diversity is present in the
S. lycopersicum gene pool [
15,
17,
18]. As a result, there is a vivid debate among consumers on the usage of recently established commercial varieties at the expense of heirloom cultivars, since the latter are considered superior in several qualitative aspects regarding flavor, aroma and phytochemicals [
4].
Currently, a significant turn towards qualitative nutriment attributes is occurring. Moreover, food science and technology have been shifting away from merely the notion of adequate calorie intake, towards an optimal and efficient nutrition scheme. Thus, currently, the breeding aim is to promote the consumption of bio-functional foods and contribute to improved physical health, therefore preventing the risk of diseases [
19]. Tomatoes are characterized by a high nutritious value and furthermore contain an extensive assortment of natural antioxidants, such as ascorbic acid, carotenoids, flavonoids and phenolic substances [
20]. The total phenolic ratios have been generally related to reactive oxygen species (ROS) scavenging and therefore are believed to be protective against cellular oxidative damage; by extension, they can be beneficial against severe illness like cardiovascular diseases, development of cancer or even neurodegenerative disorders [
21].
Tomato landraces epitomize the most basic form of existing cultivars and are largely regarded as intermediates of the breeding process [
22]. Opposed to modern-day cultivars, landraces are very heterogeneous since they were recurrently designated for their performance in adverse and low-input agricultural environments, as well as qualitative criteria e.g., aroma [
23]. Due to a unique combination of taste, tradition and functionality, heirlooms and tomato landraces are in the breeders’ spotlight and are at the epicenter of breeding efforts to re-establish nutrient and flavorsome tomato fruits. As a result, studies aiming to properly characterize the tomato germplasm are increasingly gaining attention.
When aiming to harness the phenotypic and phytochemical variation of cultivated tomato, and in order to engage this assortment in applied and basic breeding schemes, it is vital to quantify and count such traits in a precise and factual manner. Moreover, phenotypic and biochemical analyses of genetic resources are essential in order to appraise the genetic basis that connects these characters and their overall yields [
24]. However, the evaluation of such traits can be challenging and laborious, mostly due to the number of such traits [
25], and should be accompanied by genetic markers that are not affected by growing conditions. Nonetheless, while the majority of scholars emphasize the assessment of modern, as well as traditional cultivars, studies at the regional level are largely uncommon. It is generally thought that indications regarding the disparities within tomato landrace germplasm are still rather inadequate, since regionally grown traditional genotypes should not be conceived as strictly homogenous [
26].
In view of these remarks and due to the shortage of comprehensive information regarding the Cypriot tomato genetic recourses, the objectives of the current study were: (i) to perform a morphological description using scores and descriptors according to worldwide standard norms (International Plant Genetic Resources Institute (IPGRI) descriptors); (ii) to use microsatellite genetic markers and append landraces of other origins to define the genetic structure and intra-relations of Cypriot tomato varieties; (iii) to assess the variation in phytochemicals and other traits of tomato fruits from plants grown under uniform greenhouse conditions and (iv) to investigate the relations between them. The overall goal was to increase the information regarding the local tomato germplasm, which has the potential for exploitation in modern cultivation schemes and direct usage in marketplaces, breeding schemes and for future genomic surveys.
2. Materials and Methods
2.1. Plant Material
The plant material used in the current study consisted of 19 discrete landraces. Ten were Cypriot landraces, namely AR100731, ARI00732, ARI00733, ARI00734, ARI00735, ARI00736, ARI00737, ARI00872, ARI00905 and ARI00906 (
Table 1) reserved at the genebank of the Agricultural Research Institute (ARI), Cyprus. These accessions were studied in terms of morphological traits and biochemical properties, and were genotyped with microsatellites. Furthermore, eight Greek and one French heirloom varieties were genotyped in order to detect possible inter-genetic affiliations to the Cypriot tomato germplasm.
2.2. Cultivation Scheme
All Cypriot tomato genotypes were seeded at a nursery (using peat as a substrate) and transplanted at a two-leaves-stage to a greenhouse (March 2018), according to a randomized complete block design, and within row distance of 50 cm. Three discrete blocks were employed and within each plot five plants per accession were studied (in total 150 were evaluated). For fruit morphological traits, at least ten randomly selected fruits from the second and third truss were analyzed as indicated (IPGRI descriptors for Tomato). Irrigation was provided via drippers at a frequency of two days (for 30 min). Five intervals of fertilization were conducted using an all-purpose water-soluble commercial fertilizer (20-20-20), while pest management against whiteflies and Tuta absoluta was conducted via the repeated application of insecticides (Nuprid 200 SL and Bolivar 1.8 EC, respectively) when needed, before anthesis.
2.3. Analysis of Phenotypic Data
On a single tomato plant basis, 24 discrete agro-morphological traits were scored as detailed in
Supplementary Table S1. In brief, the recorded IPGRI descriptors were categorized as vegetative (6 traits), inflorescence (4 traits), fruit descriptors (12 traits) and agronomic characteristics (2 traits) and corresponded to scale, ordinal and nominal data types. Since such traits are a mixture of numerical and categorical variables, several statistical incongruities can occur when applying standard component analyses; hence for the evaluation of morphological descriptors, an R package dedicated to multivariate analysis of mixed data (PCAMIX) was employed [
27]. Euclidean distances were calculated across landraces and a distance matrix was produced. Moreover, a principal component analysis (PCA) plot was constructed from the combined morphological descriptors to depict the affiliations across the tomato landraces, and squared loadings were determined in order to detect correlations among the trait descriptors.
2.4. DNA Extraction
For nucleic acid extraction, tips from young tomato leaves were collected and kept among moist towel paper on ice, till storage at −70 °C. DNA extraction was conducted using the Dneasy Plant Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The purity and concentration of DNA were estimated by means of nanodrop spectrophotometry. Ten plants per accession were randomly sampled.
2.5. PCR and SSR Genoty**
Fifteen simple sequence repeat (SSR) markers were initially screened for the preliminary analysis of the Cypriot tomato germplasm. Based on the efficiency, reproducibility and allele size generated, ten loci (
Supplementary Table S2) were selected for the full panel of genotypes. All forward primers used were designed to have a universal M13(-21) sequence at the 5′-end and were extended by an overlap** M13 labeled primer (FAM
TM, JOE
TM ROX
TM, or TAMRA
TM). This permitted a one-tube, single-reaction nested PCR, as previously described [
28]. Based on the fragment sizes from pilot reactions, four multiplex panels (one per fluorescent dye) were determined.
For PCR reactions, each panel master mix contained 50 ng of template DNA, 10 pmole of the labeled M13 tailed forward primer, 10 pmole of the reverse and 2.5 pmole of the forward primer, 0.2 mM dNTPs, 0.5 U KAPA Taq DNA Polymerase (Kapa Biosystems, Basel, Switzerland) and a 2.5 mM final concentration of MgCl2 in a 12 μL final reaction volume. Conditions for the PCR amplification were: 94 °C (5 min for initial denaturation), followed by 39 cycles at 94/56/72 °C (60 s), and a final extension at 72 °C for 30 min.
Amplification products were verified using a standard 2% agarose electrophoresis and diluted at a 1:40 ratio with dd H2O. Diluted PCR products across panels were all mixed in one tube (per accession). One μL of the mixture was added to 10 μL deionized formamide and 0.2 μL of DNA size standard (GeneScan 500-LIZ, Applied Biosystems, Foster City, CA, USA), before denaturing at 95 °C (5 min). Allele fragments were separated by capillary electrophoresis using an Applied Biosystems 3130® Genetic Analyzer (Applied Biosystems, Foster City, CA, USA).
2.6. Molecular Data Analysis
Allele fragments were identified (de-multiplexed) by means of fluorochrome colour and/or bin size, and tandem software was utilized to verify/correct bins [
29]. Microsatellite data curation and formatting was performed via the MS Excel add-in GENALEX v. 6.501 [
30]. All genotypes were included for calculating allelic frequencies across loci studied. In order to assess the discriminating power among unique genetic profiles, a genotype accumulation curve was constructed. Additionally, genotypic diversity was assessed with several indexes (H: Shannon–Wiener Index of Multi-locus Genotype (MLG) diversity, G: Stoddart and Taylor’s index of MLG diversity, lambda Simpson’s index, E.5: Evenness of the alleles and Hexp: Nei’s unbiased gene diversity). In order to determine allelic abundance, several indexes were evaluated: Na = No. of Different Alleles; Ne = No. of Effective Alleles = 1/(Sum pi
^2); I = Shannon’s Information Index = −1* Sum (pi * Ln (pi)); Ho = Observed Heterozygosity = No. of Hets/N; He = Expected Heterozygosity = 1 − Sum pi
^2; uHe = Unbiased Expected Heterozygosity = (2N/(2N − 1)) * He; F = Fixation Index = (He − Ho)/He = 1 − (Ho/He), where pi is the frequency of the i
th allele for the population and sum pi
^2 is the sum of the squared population allele frequencies. Moreover, the proportion of polymorphic loci (PIC) and discrimination power (Dp) of each locus was determined.
The same dataset was similarly used to test for linkage disequilibrium and Hardy–Weinberg equilibrium (HWE) in the tomato accessions. Genetic relationships between individuals (MLGs) were assessed using the ‘dissimilarity’ distance algorithm and visualized as a minimum spanning network (MSN) and a discriminant analysis of principal components (DAPC). All of the above-mentioned statistics/analyses were performed using the Poppr (V. 2.8.5) package [
31] and the RStudio suite (V 1.2.5033; R V 3.6.2). The PIC and Dp indexes were calculated using the iMEC: Online Marker Efficiency Calculator (
https://irscope.shinyapps.io/iMEC/ (accessed on 1 August 2021)).
A phylogenetic tree was also constructed using the binary template (converted from allele size) using the R package polysat [
32]. An approximate likelihood-ratio test (aLRT) for branch support was achieved by means of the SH-like parameter as previously described [
33]. The Newick-formatted tree was displayed and manipulated using the iTOL v4 server [
34].
A Bayesian statistic-employing method for estimating genetic kinship was performed using Structure 2.3.4 [
35]. The admixture model was selected and 10 independent repeats per K value (extending from 1 to 20) were run. Each run involved 100,000 iterations of the burning period and 500,000 post-burning simulations. Validation of the most probable number of K-clusters and visualization was achieved using the Clumpak server (
http://clumpak.tau.ac.il/).
In order to examine the possible correlation of mean morphological traits with the genetic distance of landraces, a mantel test (999 permutations) was computed using the distances matrix of morphological descriptors and the genetic distances matrix calculated via the GENALEX v. 6.501 software [
30].
2.7. Tomato Fruits Sample Preparation
Harvesting (approximately 30 to 40 days after anthesis) was performed at the same ripeness level (red ripe; more than 90% of the surface had red colour) across all tomato fruits from the second and third trusses. Tomatoes were washed with dH2O, dried with absorbent paper and immediately placed at −80 °C to stop possible metabolic processes. For the phytochemical analyses, fruits were cut in quarters and at least five discrete fruits were bulked in order to create one sample (biological replicate). Bulked samples were blended and juice was lyophilized. In total, three discrete biological repetitions were estimated across phytochemical analyses.
2.8. Quantification of Total Phenols (TP)
Total phenols were estimated following the Folin–Ciocalteu reagent methodology revised for a microplate reader, as previously reported [
36]. To attain tomato extracts, one mL of absolute ethanol was added to 0.05 g of lyophilized tomato-fruit powder. Samples were sonicated for 10 min and centrifuged at 13,000 rpm for five min at room temperature. Forty μL aliquots of extract (supernatant) were added to an equal volume of 0.1 M Folin–Ciocalteau reagent and incubated for five min at constant stirring. Forty μL of 0.5% Na
2CO
3 were further added, incubated at 40 °C for 30 min, agitated for one min, and absorbance was estimated at a 750 nm wavelength (Tecan Infinite 200 PRO, Männedorf, Switzerland). For the standard curve, serial dilutions of gallic acid (GA) were employed, and results were expressed as mg GA equivalents/100 g FW. All trials were performed in biological triplicates.
2.9. Vitamin C Assessment
Determination of the vitamin C content followed the reduction of the 2,6-dichloroindophenol sodium salt (DCIP) method revised for a microplate reader according to Ochoa-Velasco and co-workers [
36]. In general, 0.1 g of lyophilized tomato fruit samples were mixed with 0.1% oxalic acid, incubated for five min at room temperature and centrifuged at 13,000 rpm. Forty μL of extracts were transferred into a 96-well microplate (Eppendorf, Hamburg, Germany). Equal volumes of acetate buffer and 2,6-dichloroindophenol solution (30 mg/L) were added and mixed for one min. Absorbance was estimated at a 515 nm wavelength. For the standard curve, serial dilutions of ascorbic acid (AA) were employed and results were expressed as mg AA/100 g FW. All trials were performed in biological triplicates.
2.10. Lycopene and β-Carotene Quantification
Quantification of lycopene and
β-carotene was conducted according to a method previously described [
37]. All pigments were extracted from 0.05 g of lyophilized tomato-fruits powder using as a solvent 10 mL of a 4:6 (
v/
v) acetone-hexane mixture. Samples were vigorously vortexed and incubated in dark for 20 min at 4 °C. A five min centrifugation followed and 200 μL of extracts were transferred to a 96-well polypropylene microplate (resistant to organic solvents). Absorbance readings were performed at the following wavelengths: 663 nm, 645 nm, 505 nm and 453 nm, and pigments determination was calculated using the referred equations:
All trials were performed in biological triplicates.
2.11. Total Soluble Solids (TSS), pH and Titratable Acidity
The juice from three biological replicates was used to estimate the total soluble solids (TSS) via a portable digital refractometer (DR103L, Sun Instruments Corp., Torrance, CA, USA), and results were expressed in Brix. pH was measured with a benchtop pH-meter (Hanna, Woonsocket, RI, USA). Titratable acidity (TA) was assessed employing the potentiometric titration of 0.1 mol/L NaOH up to pH 8.1, using five mL of undiluted juice. Estimations were performed on a DL22 Mettler Toledo titrator (Mettler–Toledo, Inc., Columbus, OH, USA) and were expressed as citric acid g in 1 L of tomato juice.
2.12. Macro and Micro Nutrient Content in Tomato Fruits
Tissues were ground into a fine powder and passed through a 30-mesh screen. Each sample (0.5 g) was dry-ashed in a muffle furnace at 515 °C for 5 h. Then, the ash was digested in 3 mL of 6 N HCl and diluted with double-distilled water up to 50 mL. The concentrations of P, K, Ca, Mg, Mn, Zn and Cu were determined by ICP (Perkin Elmer-Optical Emission Spectrometer, OPTIMA 2100 DV, Waltham, MA, USA). Nitrogen was determined by the Kjeldahl method (BUCHI, digest automat K-439 and Distillation Kjelflex K-360, Switzerland).
2.13. Generalized Procrustes Analysis (GPA)
In order to combine morphological, genetic, mineral and phytochemical analyses, we employed a Generalized Procrustes Analysis (GPA) employing all quantitative traits. The FactoMineR, an R package dedicated to multivariate Exploratory Data Analysis was used under the ‘Commandeur’ algorithm (
http://factominer.free.fr/).
4. Discussion
Cyprus, the largest island in the eastern Mediterranean basin is located at the crossroad point of three continents (Europe, Asia and Africa) and historically it has been on the map of millenia-old trade routes. Many discrete civilizations have shaped the Cypriot agricultural tradition throughout the centuries, each contributing to the introduction of non-indigenous species; such as tomato. Moreover, Cyprus has a complex edaphoclimatic background that is divided into four discrete geological zones: (a) the Pentadaktylos (Keryneia) zone, (b) the Troodos Ophiolite, (c) the Mamonia zone and (d) the zone of the autochthonous sedimentary rocks [
38]. Complimentary to soil type diversity, the agricultural zones of Cyprus range from the sea level to the highlands of the Troodos Sierra (1952 m). In general, conditions across Cyprus are very harsh since during summertime (tomato growing season) temperatures can easily pass and remain above the 40 °C threshold, while soil is mainly characterized as poor, alkalic and rich in calcium. Taken together, the Cypriot germplasm holds a unique global place and the landraces acclimated in such distinct conditions are of fundamental importance and factually represent the Cypriot local identity. Moreover, such lines possibly have a breeding value (irrespectively of their nutrient value) since they could be crossed with more intensive varieties.
Nowadays, a gradual detachment of consumers from established tomato cultivars towards locally-grown traditional cultivars is ongoing. This attitude is revitalizing the interest in heirloom and rustic tomato landraces and is leading to the rediscovery of local agricultural traditions [
19]. Hence, several studies are increasingly focusing on the description of agro-morphological traits, diversity analysis and phytochemical characterization of local and antique tomato germplasm [
26,
39,
40,
41,
42,
43,
44,
45]. In the current study, ten accessions of Cypriot tomato landraces were characterized using 24 agro-morphological descriptors, microsatellite genetic markers (across ten loci) as well as several phytochemical parameters in order to characterize, for the first time in literature, the local tomato germplasm and further evaluate its breeding potential.
4.1. Morphological Parameters
All landraces presented an intermediate growth trait and had dense foliage as frequently reported across Mediterranean genotypes [
46,
47]. As a result, these genotypes are not determinate and thus require pruning, formation and support, in order to sustain fruit production and achieve the optimal nutritional equilibrium. Across landraces, a rich diversity was established at almost all morphological traits. Nonetheless, the intra-variability of Cypriot accessions was lower and genotypes were highly uniform; at least in the cases of ordinal and nominal traits. Still, the coefficient of variation (CV) indexes across quantitative parameters for the collection were elevated (
Supplementary Table S3); stressing the diversity of the panel and the possibility of selecting preferable characteristics in future breeding schemes. The heterogeneity of the Cypriot panel is in accordance with similar studies focusing on the phenotypic diversity of Greek [
47], Italian [
42] and Spanish tomato landraces [
43]; although Cypriot genotypes did not reveal distinct intra-varietal morphotypes as previously reported [
39].
The present collection included genotypes producing fruits mostly appropriate for local trading and self-consumption, since tomato seeds were initially collected directly from farmers who grew tomatoes for such purposes, rather than having long shelf life for wider distribution, thus, the predominant fruit shapes and colour hues mostly reflected local preferences. We also found a considerable level of variation for fruit size and weight, while in several cases, the fruits exceeded the 500 g threshold, having an extremely large size (
Supplementary Table S3). Parisi and coworkers also reported size variation and several morphotypes in the ‘Sorrento’ tomato Italian landrace [
44] while Terzopoulos et al. reported that approximately 70% of the Greek tomato landraces studied had fruits with low weight [
47]. Such outstanding levels of measured traits have been frequently attributed to both the unique genotypic potential, as well as the environmental adaptability and the capacity of system production [
46]. Among morphological traits, the fruit shape and size, ribbing and colour are essential criteria for the definition of a tomato type [
48]. Four major types were distinguished across Cypriot landraces (
Figure 1). Large fruiting tomatoes having a flattened (ARI00731, ARI00732, ARI00737 and ARI00872) or heart-shaped outline (ARI00733, ARI00905 and ARI00906) were predominate, while two accessions had a more typical rounded figure (ARI00735 and ARI00736). Only one highly rounded type was asserted (ARI00734). In a recent study [
44], the analysis of morphological traits in Italian landraces divided the germplasm in two distinctive major groups (flattened, obcordate or oblate fruit shaped and heart-shaped or circular fruits showing angular or circular shape in cross-sections.
The majority of accessions had also distinguished greenback (green shoulders). The shape of the pistil scar was also a significant breeding trait that showed variation across genotypes. In general, flattened/heart-shaped tomato fruits presented an irregular shape of pistil scar (
Supplementary Table S3) with the exception of ARI00732, which generally developed a small and round-shaped scar. As previously reported [
49], a large and irregular size of the pistil scar can reduce the commercial value of the tomato fruit and the post-harvest shelf-life due to increased water loss and pathogen susceptibility. Hence it can be recognized that promising traits can be found within the Cypriot collection.
All morphological traits were converted to eigenvalues and a PCA was constructed (
Figure 2). The PCA analysis of mixed type data revealed that the Cypriot panel was mainly classified into three clusters while one accession remained unrelated (ARI00734). The first two axes explained less than 50% of the total variation. The rather low contribution for each eigenvalue when explaining variation in tomato landraces has been attributed to either a high intra-population variation across genotypes or due to low inter-genetic variation [
46]. Additionally, the close geographical proximity (Cyprus has an area of less than 9500 km
2) may cause a geographical bias that further restricts landrace diversity and narrows the genetic basis.
4.2. Genetic Variation
Several molecular marker systems have been employed to survey the genetic diversity of tomato collections such as amplified fragment length polymorphisms (AFLPs), randomly amplified polymorphic DNA (RAPD), sequence-characterized amplified regions (SCARs) and single nucleotide polymorphisms (SNPs) [
50]. Nonetheless, microsatellites are still considered as among the most reliable and reproducible interlaboratory techniques across scholars and are also co-dominant markers. Still, despite extreme morphological variations across tomato types, it has been established that the level of genetic polymorphisms detected by genetic markers is rather low [
39]. Recently, several studies focusing on simple sequence region germplasm characterization [
51,
52,
53,
54] reported a low mean number of alleles across the loci studied ranging from 1.8 [
54] to 9.5 [
53]. In the current study, 190 genotypes (100 of Cypriot origin) were analyzed and a mean number of 3.2 alleles per locus was established (
Table 3). The rather low genetic diversity in the European tomato germplasm seems not to be correlated to the marker system used for fingerprinting. Parisi et al. [
44] scored on average 120 SNPs in each tomato chromosome, and less than 20% resulted in polymorphisms. In the current study, the PIC index revealed that on average 0.347 loci are polymorphic across accessions. Recently in Greek tomato landraces [
53], a PIC value of approximately 0.7 was reported, while 50% of the SSR loci employed for the characterization of the Bulgarian germplasm were found to be polymorphic [
54]. Hence, it can be ascertained that SSRs are a useful marker system for tomato breeding.
Both expected and observed heterozygosity indexes (
Table 2) varied but had a generally low range (0.100–0.400 for Ho; 0.098–0.200 for He) indicative of self-pollinated species [
55]. Additionally, the majority of Cypriot landraces were found genetically uniform since only one genotype (MLG) was detected, in contrast to Greek heirlooms that presented a more heterogenous composition. This is in accordance with the findings of Terzopoulos and colleagues [
39,
47] that reported a substantial amount of intra-population heterogeneity and the concurrent presence of several morphotypes within Greek tomato landraces.
All genetic data analyses (hierarchical clustering/dendrogram, MSN, and Bayesian inference) indicated that there is neither a geographical structure nor does clustering correspond with fruit type. This is accordance with several studies that describe a genetic grou** uncorrelated to several morphological traits and fruit characteristics [
52,
54]. Hence, it can be established that the primer set used in the current study does not correspond to any functional marker linked to phenotypic traits. Moreover, a Mantel test among morphological and genetic distances revealed a significant but rather low-to-moderate correlation. Nonetheless, several significant outcomes could be established that correlate to the agro-morphological analysis. Landraces ARI00735 and ARI00736 were highly affiliated and differed at one locus out of ten hence probably represent different clones of the same landrace. The structure analysis indicated that two major genetic lineages could be distinguished (
Figure 5). Most Greek genotypes (from Aegean islands) where clustered in the second genetic group alongside Cypriot landraces ARI00735 and ARI00736. Interestingly, both genotypes from Santorini island (‘Katsari Santorinis’ and ‘Leia Santorinis’) despite having a different fruit type were highly genetically affiliated; hence the possibility of intercrossing among these genotypes cannot be uncritically ruled out. In fact, it has been reported that tomato landraces from Santorini island have a very heterogenous lineage. Tomato farmers in Santorini reserve three discrete types of fruit shape; thus, fulfilling diverse needs: rounded fruits are used for preserves and tomato juice, while slightly flattened or flattened fruits are employed for sun-dried tomatoes [
47].
Overall, the current germplasm seems to be genetically distinct and thus could serve as a valuable addition for breeding schemes targeting the genetic enhancement of modern cultivars. Moreover, the employment of heirloom varieties and tomato landraces could reverse the loss of genetic diversity, the farming of allochthonous varieties and ameliorate the ongoing genetic erosion. Moreover, since these landraces have evolved under low input agricultural systems in the semi-arid Cypriot environment, these genotypes might serve as valuable germplasm in sustainable farming systems.
4.3. Phytochemical Characterization
Taking into consideration a typical proteinic, lipidic and sugar content in order to describe the nutritional value (in terms of calories), it seems that tomato fruits do not offer a high nutritional value [
1]. Still, tomatoes offer an important basis of minerals and nutrients that are defining components for human health promotion, such as antioxidants, ascorbic acid (vitamin C), lycopene and vitamin A (
β-carotene);
Table 5. Thus, tomatoes are the foremost source of lycopene, which presents antioxidant capacity and is assumed to be a protective agent against cardiovascular diseases and malignant tumors [
56]. Tomatoes are further considered a vital and outstanding source of ascorbic acid. Comparing with modern cultivars, wild tomato genotypes and landraces are richer in ascorbic acid and may present up to five times more ascorbic acid than their cultivated counterparts [
57]. Indeed, Cypriot genotypes can be characterized as highly nutrient tomatoes since they were found to possess significant levels of minerals (21.5 ± 1.38 g/kg DW (N), 5.96 ± 1.35 g/kg DW (P), 2.32 ± 1.00 g/kg DW (Ca), 30.17 ± 5.92 mg/kg DW (Zn), 15.53 ± 3.20 mg/kg DW (Mn) etc.) and vitamins (33.28 ± 8.78 mg/100 g FW vitamin C, 3.55 ± 1.67 mg/100 g FW lycopene and 0.65 ± 0.16 mg/100 g FW
β-carotene). In that respect, rustic tomato Tuscan varieties were reported to have a three-fold increase in antioxidant phytochemicals, compared to the commercial counterparts [
19]. All Cypriot tomato landraces were found to produce fruits of high nutritional value and thus have the potential for commercial exploitation in breeding schemes and trait selection. A clear correlation of fruit shape and phytochemicals was not profound; but in some instances, round tomatoes (landraces ARI00735 and ARI00736) were found to be richer in vitamin C and total phenolics compared to flattened and heart-shaped genotypes. Nonetheless, lycopene levels were somewhat lower than the average. A correlation of fruit type and metabolite levels is not unprecedented, and has been described previously. Carli and coworkers [
58,
59] and Figas et al. [
60] reported that round and elongated fruit tomato types have higher levels of sugars and dry weight compared to other fruit types studied. Recently, a close affiliation of alkaloids to specific fruit types was also reported [
61]. Moreover, the authors suggested that high levels of these metabolites could serve as a dual-purpose element for pathogen resilience, as well as flavor enhancement and thus may present a desirable future breeding objective.
Cypriot tomato fruits also contained significant amounts of phenolic compounds (averaging 7 mg GAE/100 g FW). It has been reported that phenolics are a highly desirable feature since these moieties encompass the major contribution of antioxidant activities in tomato fruits [
62]. Flattened and heart-shaped tomato fruits were found to be rather homogeneous for total phenolics but landrace ARI00734 (the only one having a high-rounded type) presented substantially lower levels (
Table 5). A superiority of flattened tomato fruits in phenolic concentration has also been reported [
61]. Interestingly, the levels of total phenolics naturally present in the Cypriot tomato collection are comparable to the ones reported in greenhouse-cultivated tomato fruits (4.3 to 8.5 GAE/100 g FW) after the addition of biofertilizer (
Bacillus licheniformis) on a commercial cultivar [
36]. Nonetheless, high levels of phenolics are not a strictly genotypic regulated trait, but are rather the outcome of a genotype X environment interaction. Stressor factors that can induce ROS in tissues have been implicated in the increase of phenolics, in order to ameliorate cellular damage [
63].
A critical element of consumers’ demand for tomato fruits is the TSS and acidity characteristics. The pH parameter ranged from 4.55 (ARI00732) to 4.68 (ARI00735); hence Cypriot tomatoes do not present an elevated acidity. Scarano and colleagues analyzed the phytochemical composition of an Italian landrace tomato germplasm collection under an elevated temperature and established a very narrow pH range (4.3 to 4.6) that corresponds to the present study [
64]. Nonetheless, the TSS content of the Cypriot landraces has an overall lower value (4.11 °Brix on average and a maximum of 5 °Brix) compared to other local germplasms where values larger than 5 °Brix are reported [
42,
45,
64]. The sugar content is heavily influenced by both temperature and light spectra; hence the different environmental conditions and cultivation practices across diverse studies can hamper direct comparisons among landraces grown at distant locations and timelines. Moreover, postharvest treatments can severely distort TSS levels. Kasim and Kasim reported that the fructose, glucose, and the TSS content of tomatoes treated with UVB light at the red ripe stage were found to be positively affected [
65]. Moreover, earlier studies reported that TSS content for large beefsteak tomatoes fluctuated from 3 to 5 °Brix, for medium-sized tomato fruits values ranged from 5 to 7%, and for small cherry tomatoes the TSS content was estimated to be above 9% [
66,
67]. As a result, it can be established that the fruit size is negatively correlated to sugar content. In the current study, Cypriot landraces were found to produce rather large fruits; hence the relatively low sugar levels may be attributed to augmented size.