Next Article in Journal
Potential and Actual Water Savings through Improved Irrigation Scheduling in Small-Scale Vegetable Production
Next Article in Special Issue
Molecular Genetic Diversity and Population Structure of Ginseng Germplasm in RDA-Genebank: Implications for Breeding and Conservation
Previous Article in Journal
A Technical-Economic Comparison between Conventional Tillage and Conservative Techniques in Paddy-Rice Production Practice in Northern Italy
Previous Article in Special Issue
Variability in Nutraceutical Lipid Content of Selected Rice (Oryza sativa L. spp. indica) Germplasms
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 9
Issue 12
10.3390/agronomy9120887
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

Multi-Level Characterization of Eggplant Accessions from Greek Islands and the Mainland Contributes to the Enhancement and Conservation of this Germplasm and Reveals a Large Diversity and Signatures of Differentiation between both Origins

by
Pietro Gramazio
1,2,†,
Eleni Chatziefstratiou
3,4,†,
Constantinos Petropoulos
5,
Vasileia Chioti
1,6,
Photini Mylona
7,
George Kapotis
4,
Santiago Vilanova
1,
Jaime Prohens
1,* and
Vasileios Papasotiropoulos
3,*
1
Instituto de Conservación y Mejora de la Agrodiversidad Valenciana, Universitat Politècnica de València, Camino de Vera 14, 46022 Valencia, Spain
2
Faculty of Life and Environmental Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 3058572, Japan
3
Department of Agriculture, University of Patras, Theodoropoulou Terma, 27200 Amaliada, Greece
4
Department of Crop Science, University of Patras, Nea Ktiria, 30200 Mesolonghi, Greece
5
Department of Mathematics, University of Patras, University Campus, 26504 Patras, Greece
6
Department of Biology, University of Patras, University Campus, 26504 Patras, Greece
7
Institute of Plant Breeding and Genetic Resources, HAO-DEMETER, Thermi, 57001 Thessaloniki, Greece
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Agronomy 2019, 9(12), 887; https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120887
Submission received: 11 November 2019 / Revised: 10 December 2019 / Accepted: 11 December 2019 / Published: 13 December 2019
(This article belongs to the Special Issue Analysis of Crop Genetic and Germplasm Diversity)
Browse Figures
Versions Notes

Abstract

:
Crop landraces are found in many inhabited islands of Greece. Due to the particularity of environment and isolation from the mainland, Greek islands represent a natural laboratory for comparing the diversity of landraces from the islands with those of the Greek mainland. A collection of 36 Greek eggplant landraces and traditional cultivars from the mainland and the islands has been phenotypically and genetically characterized using 22 morphological descriptors and 5 SSR markers. The mineral composition (K, Mg, Cu, Fe, Mn, Zn) of fruits was also determined. The objectives of this study include the multi-level characterization of eggplant local landraces and the comparison of diversity among accessions from the Greek mainland and the islands. Characterization of eggplant landraces will contribute to the enhancement and prevention of genetic erosion in this local group and will provide a resource for future investigation and breeding. PCA analysis of morphological traits explained 45.4% of the total variance revealing the formation of two clusters, one with most of the island accessions, and another with most of the mainland ones. The SSR markers used exhibited high average values for the number of alleles/locus (4.6), expected heterozygosity (0.60) and PIC (0.55), while the observed heterozygosity was low (0.13). Both STRUCTURE and PCoA analyses based on SSR data revealed two genetic clusters, one made up mainly by the mainland accessions, while the other one was mainly made up by the island accessions. Although there was considerable variation among the landraces for the concentration of minerals studied, only average Mg concentration was significantly different between mainland and island accessions. Based on our data, the Greek eggplant landraces present considerable morphological and genetic diversity with some differentiation signatures between the island and the mainland accessions. Our results have implications for conservation of Greek landraces and suggest that Greece might be considered as part of a secondary center of diversity for eggplant in the Mediterranean basin.

1. Introduction

Greece lies at the southernmost part of the Balkan peninsula, at a crossroads between Europe, Asia, and Africa. The country is noted for its high species diversity, the extensive coverage of natural and semi-natural habitats, and many island complexes in the Archipelagos of Aegean Sea and the Ionian Sea [1]. Islands are natural laboratories for the study of the evolution and ecology of wild plants and endemisms [2,3,4]. Therefore, due to their isolation and special environmental features, islands are also considered as sites of special interest to prospect and conserve plant genetic resources [5,6]. The Greek archipelago includes ca. 7600 islands and islets in the Aegean Sea (including Crete) and ca. 300 islands and islets in the Ionian Sea [7]. Among them, 116 are inhabited, whereas only 79 have over 100 permanent residents [8]. The Aegean Archipelago is one of the largest archipelagos in the world, exhibiting high environmental and topographical heterogeneity, complex geological and palaeogeographical history, as well as high diversity and endemisms, rendering it an ideal stage for biodiversity and biogeographical studies [1,9]. Among the 13 phytogeographical regions proposed by Strid [10] to occur in Greece, five were identified in the Aegean [1]. The Ionian archipelago forms a distinct phytogeographical region, hosting 2027 plant taxa, 1827 of which are native and 89 are Greek endemics [5,11,12,13].
Landraces represent a significant part of agricultural biodiversity [14]. They are characterized by high genetic diversity, local adaptation, resilience to pathogens, and good organoleptic characteristics, although they often have low yield and display a lack of uniformity [15,16,17]. Over time, landraces have been gradually displaced from market-driven production, due to the consumers’ demand for local high-quality plant products, landraces are currently recalling an increased interest [18,19]. With the shift to modern agriculture, Greece suffered a dramatic loss of its traditional agricultural germplasm, which was displaced by modern varieties produced by local breeders or imported [5,20]. Nowadays, landraces have regained an increased interest [18,19]. Due to the geographical morphology of Greece (almost 70% is hilly and/or mountainous) and the presence of small-size farms, landraces are still cultivated, especially by farmers in villages in remote regions and isolated islands, who maintain their seeds and use them mainly for their own consumption, and, to a lesser extent, for market purposes [5,13,19]. Therefore, these areas act as local in situ conservation spots and islands represent one of the most important sources of landraces and hotspots of agricultural biodiversity in Greece [5,20,21].
Eggplant (Solanum melongena L.) is the fifth most produced vegetable worldwide and has experienced a dramatic increase in yield and production in the last decade [22]. It was introduced in Greece around the 13th century and although initially used as an ornamental plant, it is now playing an important role in the Greek gastronomy, as an integral part of the Mediterranean diet. Besides commercial F1 hybrids and imported varieties, there are several traditional cultivars, e.g., Tsakoniki, and Santorini, which are adapted to the local environment, having excellent texture and cooking quality [23]. Apart from commercial varieties, several local landraces have been found during expeditions in the islands [5,24] or maintained in Greek Gene Bank (Hellenic Agricultural Organization Demeter, Thessaloniki, Greece). However, the diversity of Greek accessions has been barely explored, and so far, no studies have focused on the comparison of accessions from the mainland and the islands of Greece have been reported. Significant differences in genetic diversity of several plant species have been revealed when island populations are compared with mainland ones [25,26,27,28]. Factors such as geographic isolation, island size, plant breeding system, or life form, and in the case of crop species, human factors, are shaping the difference levels [26,29,30,31]. Characterization of the genetic diversity and population structure of landraces and local germplasm is very important to implement the most appropriate strategies for their collection, management, efficient conservation, and utilization as a source of variation in breeding schemes [32]. Several studies have focused on the characterization of eggplant genetic resources, including the cultivated species, commercial cultivars, and their most recognized wild relatives, from different regions of origin in order to identify desired genotypes for use in eggplant breeding programs [33,34,35,36]. For instance, Liu et al. [36] examined 287 worldwide accessions for genetic diversity and population structure analysis, whereas Acquadro et al. [33] used high-throughput genotyping to assess the genetic relationships of brinjal, gboma, and scarlet eggplant complexes, which represent taxa belonging to the eggplant’s primary secondary and tertiary gene pools. Moreover, Cericola et al. [34] used a combined marker-based and morphological approach to assess genetic diversity and illuminate the genetic relationship between “Occidental” and “Oriental” eggplant germplasm groups, while Hurtado et al. [35] evaluated the phenotypic and DNA-based diversity present in a collection of accessions sampled from three geographically separated centers of diversity.
Eggplant is an important source of plant-derived nutrients, valued for its composition in phytochemicals and especially minerals such as P, K, Ca, and Mg [37,38]. Raigón et al. [39] and Arivalagan et al. [38,39,40] studied the mineral composition and the nutritional value of eggplant landraces and germplasm accessions from different regions in order to identify mineral-rich germplasm for breeding purposes. Considerable differences were found in the mineral composition among varieties and germplasm accessions, revealing the existence of ample variation, which can be exploited for the selection of germplasm for nutritionally improved characteristics.
In this study, we report the assessment of the morphological and genetic diversity of eggplant germplasm from the Greek mainland and islands. In addition, the mineral content has been determined to obtain a more comprehensive view regarding the nutritional value of those accessions. The overall aim of this study is to acquire information for the efficient management and conservation of this valuable genetic material, for the identification of traits present in landraces that can be exploited for breeding purposes, and to compare the levels of diversity among the Greek mainland and island eggplant germplasm.

2. Materials and Methods

2.1. Plant Material and Cultivation Conditions

Thirty-six accessions, including landraces and cultivars from the mainland (16) and the islands (20) of Greece were selected to depict the Greek eggplant germplasm diversity (Table 1) and their geographical distribution (Figure 1). All materials are maintained at the National Greek Gene Bank (GGB; Institute of Genetic Resources and Plant Breeding Thermi, Greece), except I-6, I-7, and I-8, which are conserved at the Department of Agriculture of the University of Patras (Greece). Seeds were germinated in Petri dishes, following the Ranil et al. [41] protocol, and were subsequently transferred to seedling trays in a climatic chamber under a photoperiod and temperature regime of 16 h light (28 °C): 8 h dark (20 °C). One week later, the plantlets were transferred and kept in a heated greenhouse until transplantation. Nine plants per accession were grown in an open field plot at the University of Patras campus in Messolonghi, Greece (GPS coordinates: latitude, 38°36′ N; longitude, 21°47′; 1.5 m a.s.l.) following a completely randomized block design experiment with 3 blocks and 3 plants per accession/block and spaced at 0.5 m within the row and 1 m between rows. Plants were pruned and trellised and drip-fertigated. Standard horticultural practices for eggplant production in this area for pest and weed control were performed. Accession ANP-025/07 (I-1) displayed an abnormal development at the seedling stage and was discarded for the morphological and mineral characterization, although it was kept for the genotyping.

2.2. Morphological and Mineral Composition Characterization and Data Analysis

Plants were characterized using 22 morphological descriptors for plant (4), leaf (6), flower (1), and fruit (11) traits obtained from the EGGNET descriptor list, of the EU-RESGEN PL 98-113 program [42,43], Biodiversity International (formerly International Board for Plant Genetic Resources) [44], and Kumar et al. [45] descriptors lists (Table 2). For all the traits, four measures were made per plant, except for plant traits, in which one measurement was made. For mineral composition analysis, one plant per block and accession and between one and three fruits per plant (depending on fruit size) were collected and analyzed at commercial maturity. Fruits were washed, peeled, cut into pieces, and bulked. Subsequently, the bulks were weighed, dried at 80 °C until constant weight, and powdered with a mechanical grinder (T 25 digital ULTRA-TURRAX®, IKA®-Werke GmbH & Co. Staufen, Germany). A total of 0.5 g of the dried samples were calcined in a muffle furnace (Thermconcept) at 550 °C for 5 h. The ashes were then dissolved in 10 mL of HCl 1Μ, filtered and the extract was brought to 50 mL with distilled water [46]. K was analyzed by flame photometry (Sherwood Model 410, Cambridge, UK), while the rest of the minerals (Mg, Fe, Cu, Mn and Zn) were analyzed by atomic absorption spectrophotometry (AAS) using a Thermo Elemental (SOLAAR AA Spectrometers, Cambridge, UK) spectrometer.
Averages for accession were used to calculate means, ranges, and coefficient of variation for mainland and island accessions using IBM® SPSS® Statistics 25 (IBM, Armonk, NY, USA). Significance of differences between means of mainland and islands was performed using the appropriate test at the significance level of p < 0.05. The results were analyzed using both parametric (t-test) and nonparametric (Wilcoxon-Mann-Whitney) methods, the latter in the cases of non-normality. A principal component analysis (PCA) was performed using the function prcomp of the package stats (v3.6.1) in R [47] and represented graphically using the package ggplot2 [48].

2.3. SSR Characterization and Data Analysis

Total genomic DNA was extracted for each accession from nearly 100 mg of young leaf tissue [49]. DNA quality (230/260 and 260/280 nm ratios) and concentration were measured using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and a Qubit® 2.0 Fluorometer (Applied Biosystems, Waltham, MA, USA), respectively, while DNA integrity was evaluated by agarose gel electrophoresis (0.8%). Three genomic SSRs developed by Vilanova et al. [50] and two by Nunome et al. [51] were used due to their high polymorphism for the genetic characterization of the samples, which were organized in one multiplex reaction according to the expected allele size range (Table 3). The PCR amplification was performed in a total volume of 12 μL including 7.21 μL water, 1.2 μL 1× PCR buffer, 0.6 μL MgCl2 (50 mM), 0.24 μL dNTPs (10 mM), 0.06 μL forward primer with M13 tail (10 μM), 0.24 μL fluorochrome (FAM, VIC, NED, and PET, 10 μM), 0.3 μL reverse primer (10 μM), 0.15 μL Taq DNA Polymerase (5 U/μL), and 2 μL DNA template (20 ng/μL). The PCR was performed following the program: 95 °C for 3 min for a denaturation step, 30 cycles of 30 s at 95 °C followed by 30 s at 65 °C and 30 s at 72 °C and finally 72 °C for 5 min for the last extension step. The PCR products were subsequently diluted in formamide and sequenced by capillary electrophoresis through an ABI PRISM 3100-Avant sequencer (Applied Biosystems, Waltham, MA, USA) using a 500 LIZ GeneScan size standard (Applied Biosystems, Waltham, MA, USA). The fragments were analyzed using the GeneScan software (Applied Biosystems, Waltham, MA, USA) to obtain the electropherograms and the alleles were identified with the Genotyper DNA Fragment Analysis software (Applied Biosystems, Waltham, MA, USA).
The SSRs analysis was performed with the software PowerMarker [52] and GenAlEx 6.503 [53]. For each marker, the following parameters were calculated: number of alleles (Na), major allele frequency (f), number of effective alleles (Ne), number of genotypes (Ng), polymorphic information content (PIC) that was calculated using the following formula PIC = 1 − ∑ i = 1 n p i 2 − ∑ i = 1 n − 1 ∑ j = i + 1 n 2 p i 2 p j 2 (where n is the total number of alleles detected, pi is the frequency of the ith allele, and pj is the frequency of the jth allele) [54], observed heterozygosity (Ho) calculated as the number of heterozygous alleles/number of alleles, expected heterozygosity (He) calculated as He = 1 − ∑ i = 1 n p i 2 (where p i is the frequency of the ith allele) [55] and the inbreeding coefficient (F) calculated as F = 1 − (Ho/He) [56]. Molecular analysis of variance (AMOVA) among groups (mainland vs. islands), among accessions within a group, and within individuals was performed with 999 permutations (α = 0.05) using GenAlEx 6.503 software [53]. Correlations between morphological and SSR distance matrices were investigated using a Mantel test (α = 0.05) with 999 permutations using GenAlEx 6.503 software [53].
Principal coordinates analysis (PCoA) was performed using the function dudi.pco of the package ade4 (v1.7-13) in R [47] and represented graphically using the package ggplot2 [48]. A model-based Bayesian structure implemented in the software STRUCTURE (version 2.3.4) [57] was used to estimate the population structure. Twenty runs of STRUCTURE were performed by setting the number of clusters (K) from 1 to 10 with a length of the burn-in period of 500,000 steps followed by 150,000 Monte Carlo Markov chain (MCMC) replicates, assuming an admixture model and uncorrelated allele frequencies. No prior knowledge of the population of origin was introduced. The ΔK method [58] was used to identify the most likely number of clusters (K) using STRUCTURE HARVESTER 0.6.94 software [59]. Each accession was assigned to its corresponding group based on maximum membership probability, as indicated by Remington et al. [60].

3. Results

3.1. Morphological Characterization

A wide diversity was observed for the morphological traits in the germplasm collection of Greek eggplant (Table 4, Table S1). For traits measured in a quantitative scale, accession values encompassing a broad range of the scale were observed for most of the traits, except for LBLo (3.00 to 5.00), LP (0.00 to 0.92), FS (4.00 to 5.33), and FC (1.00 to 5.00). A wide variation was also observed for the quantitative traits, especially for PH, LPL, NOFPI, FPL, FL, FB and FW. Although the ranges of variation between accessions from the mainland and islands overlapped, island accessions showed a broader coefficient of variation (CV) for all the quantitative traits, except for NOLFF and LBW. However, significant (p < 0.05) differences between the two groups were found only for LBW, LBLo, FCP, FB, FLBR, FC, and FW.
The orthogonal transformation by PCA of the morphological traits revealed that the total variation was explained by 22 principal components (PCs), although the first five PCs explained altogether 71.5% (Supplementary data S2). The bi-dimensional representation of the first two PCs revealed that the PC1, which accounted for 32.1% of the variation, was able to separate the island from the mainland accessions (Figure 2A). In fact, most of the island accessions (13 out of 19) exhibited positives PC1 values, while most of the mainland accessions showed negative PC1 values (14 out of 16). Similarly, most of the island accessions (15 out of 19) displayed positive values for the PC2, which accounted for 13.3% of the variation, while mainland accessions (11 out of 16) displayed negative PC2 values. Interestingly, although the accessions from the same geographical area fell relatively close in the PCA, like the island accessions I-6, I-7, I-8, and I-9 from Leros or I-10 and I-11 from Lesvos and the mainland accessions M-8, M-9, and M-10 from Serres, they showed a certain degree of morphological diversity. No high correlation absolute values were found between the morphological traits and the first two PCs (Table S2). However, moderate to weak PC1 positive correlations were found for FB and FW and negative correlation for FLBR and FC, while weak PC2 positive correlations were found for LBW, LPL, LBLe, and NOFPI, and no negative PC2 correlations were found (<−0.3).

3.2. Concentration of Minerals

The concentrations of the macrominerals K and Mg were much higher than those of the microminerals Cu, Fe, Mn and Zn (Table 4). As occurred for morphological traits, a wide variation was observed among the accessions for mineral composition, with ranges of variation in the collection of over 10-fold for Cu, and ranges of the coefficient of variation of up to 67.2% for Cu in the accessions from the islands. However, for all traits, an overlap of the values was observed for the six minerals analyzed among the two groups, and significant differences between averages of mainland and insular groups were observed only for Mg, with slightly larger values in the mainland group.

3.3. SSR Characterization

All five SSRs were polymorphic and amplified between three (em311f04) and seven (csm4) alleles totaling 23 alleles, resulting in an average of 4.6 alleles per SSR locus (Table 5, Table S3). The major allele frequency was, on average, slightly higher than 0.5, ranging between 0.31 (csm4) and 0.65 (eme11f04). The average number of effective alleles (Ne; 2.74) was lower than the number of alleles (Na; 4.60), ranging from 1.91 (eme11f04) to 4.41 (csm4), while the number of genotypes (Ng) was higher than Na, ranging from 4 (emi02c21) to 10 (csm4), except emi02c21, which showed the same number of alleles and genotypes, indicating that all loci were in homozygosis. The PIC value had an average value of 0.55, with a range between 0.39 (eme11f04) and 0.75 (csm4). The Ho showed low values with an average of 0.13, ranged from 0 (emi02c21) to 0.17 (csm27), much lower than those of He (0.60 on average), which ranged from 0.47 (eme11f04) to 0.77 (csm4). On the contrary, F values were high, ranging from 0.66 (eme11f04) to 1 (emi02c21). The comparison between mainland and island accessions revealed that the latter presented generally higher genetic diversity (Table 5). In this respect, mainland accessions showed a lower number of alleles compared to the island accessions (3.2 and 4.2, respectively), as well as a lower number of effective alleles (2.22 vs. 2.79), genotypes (3.8 vs. 5.6), and PIC (0.42 vs. 0.53) but a higher frequency of the major allele (0.60 vs. 0.48). Average values for Ho and He were also much lower in the mainland accessions (0.03 and 0.48) than in the islands’ accessions (0.20 and 0.61), while F values were much higher in the mainland accessions (0.92) than in those from the islands (0.68).

3.4. SSR Molecular Analysis of Variance, Genetic Structure, and Multivariate Analysis

The AMOVA revealed that the variation among accessions accounts for the greatest proportion of molecular variance (69.62%), followed by the variation within individuals due to heterozygosity (19.25%), and finally by the variation among mainland and island groups, which accounts for just 11.12% of the total variation (Table 6). The ΔK statistic in the genetic structure analysis presented a maximum peak at K = 2, suggesting that two genetic clusters exist in the Greek materials evaluated (Figure S1). With very few exceptions, the percentage of the ascription of accessions to their cluster was over 80%, revealing a low degree of genetic admixture (Figure 3). Cluster one included 14 accessions, of which 10 (71.4%) were from the mainland and 4 (28.6%) were from islands. Cluster two was made up by 22 accessions, of which 6 (27.3%) are from the mainland, and 16 (72.7%) are from the islands.
The PCoA analysis based on SSRs data explained all the variation with 17 PCs, although the first five PCs accounted for 75.8% (Table S4). Like in the PCA analysis with morphological data, the PC1, which accounted for 27.4%, allowed to separate the island accessions, with positive values of PC1 (15 out of 20), from the accessions of the mainland, with negatives values of the PC1 (11 out 16) (Figure 2B). On the contrary, PC2, which accounted for 17.6%, cannot clearly separate the two groups that spread from positive to negative values. Like the PCA on morphological traits, the accessions from the same geographical area showed a certain degree of genetic diversity. Finally, the Mantel test showed a weak, although significant, correlation between genetic and morphological data (r = 0.244, p = 0.002) (Table S5).

4. Discussion

In a climate change scenario coupled with an unprecedented increase in food demand, due to dietary changes and relentless growing population [61], local germplasm characterization and enhancement is imperative. The genetic variability and allelic diversity of local landraces could mitigate the effects of a changing environment and may guarantee food security of domestic markets [62]. In addition, heirlooms and ecotypes can harbor allelic combinations useful for the development of new more resilient and locally adapted varieties [14]. In this respect, this study aimed to assess the morphological and genetic diversity of the eggplant germplasm collection from Greece which is distributed and cultivated in a wide variety of niches throughout the distinctive country’s topography. Moreover, the characterization of Greek eggplant landraces will prevent genetic erosion of this local group and provide a valuable resource for utilization in research and breeding.

4.1. Morphological Diversity

The morphological characterization displayed considerable diversity for most of the traits studied. The most heterogeneous traits were those related to fruit appearance, such as shape and size, followed by others related to vegetative traits, such as plant and organ size and type of growth, alongside with skin and flesh color. The solanaceous fruits exhibit considerable morphological diversity, including size, shape, and color, both within and between different species [63]. Morphological variations of fruits in terms of shape, size, and color are a result of adaptive evolution; in addition, variation observed in the vegetative traits could be due to the accumulation of mutations on neutral traits, to artificial or natural selection or to both [63,64]. These same types of morphological traits were also the most variable ones when eggplant accessions collected from Greece, Spain, Turkey, and remote centers of diversity were studied [23,34,35,65,66]. Consistent with previous works, the amount of diversity is comparable to that observed in the characterization of germplasm from other geographical areas of the Mediterranean region [34,35,64,65,66], providing further evidence that several areas of the Mediterranean region can be considered as secondary centers of diversity for eggplant [35,64].
PCA analysis revealed the existence of two discrete groups, one formed almost exclusively by the island accessions and the other mostly by those from the mainland ones. On average, the accessions from the islands have higher values for fruit weight, fruit calyx prickles, leaf blade width, leaf blade lobing, fruit breadth, but lower ones for fruit curvature and fruit length to breadth ratio than those from the mainland. The assignment of the eggplant accessions studied into groups according to their agromorphological characteristics correlated with their geographic origin has also been reported in other eggplant studies [35,67,68,69].

4.2. Mineral Composition

As occurred for morphological traits, the wide variation observed for mineral composition agrees with other works in which a significant number of eggplant accessions, commercial varieties, and hybrids were evaluated [65]. In the present study, significant differences were found only for Mg among island and mainland accessions, while non-significant differences were observed for Fe, Mn, Zn, Cu, and K. Raigón et al. [39] and Arivalagan et al. [38,39,40] studied the mineral composition and the nutritional value of eggplant landraces and germplasm accessions from different regions in order to identify mineral-rich germplasm for breeding purposes. For comparison, we transformed the values detected in the aforementioned studies, to the same units as in our study. The eggplant accessions studied had lower K and Cu concentration than those of Raigón et al. [39] and Arivalagan et al. [38,40]. Mg content in our study was similar to that of Raigón et al. [39] and higher than in Arivalagan et al. [38,40], while Zn concentration was significantly higher than in these studies. The environment, cultivation methods, as well as the cultivar used can influence the mineral composition of eggplants [70,71], and explain the differences observed in minerals’ concentrations among our study and the others. To our knowledge, this is the first time that an inclusive uniform mineral content characterization of eggplant fruits from a field-grown germplasm collection was examined.

4.3. Molecular Diversity

Many studies assessed the genetic diversity and population structure of local varieties collections from different origins using SSRs [34,35,36,45,70,71]. The SSR markers used, were highly informative and effective for genetic diversity analysis. For example, the PIC value in the whole collection ranged from 0.39 to 0.75 with a mean value of 0.55. This value is higher than those observed in previous studies [36,45,72,73], while it was lower than those calculated by Hurtado et al. [35] and Cericola et al. [34]. A PIC value greater than 0.5 indicates loci of high polymorphism, PIC values between 0.25 and 0.5 indicate loci of intermediate polymorphism, and PIC values less than 0.25 indicate loci of low polymorphism [72]. In addition, the number of alleles/locus in the whole collection was high (4.6), higher than that (3.67) revealed by Augustinos et al. [23], when they studied a limited number of Greek traditional eggplant cultivars. Taking into consideration other studies on eggplants from different geographic origin or different cultivar types based on SSR data, the number of alleles/locus calculated in our eggplant collection is lower than that detected by Hurtado et al. [35] and Liu et al. [36], but almost similar to those by Demir et al. [74], Muñoz-Falcón et al. [73], Tümbilen et al. [66], and Vilanova et al. [50,75]. Two of the SSR markers used in our study (csm27, csm32), have also been used by Hurtado et al. [35] and Vilanova et al. [50]. In their studies, they observed higher average number of alleles/locus as compared to those detected by us, an observation that might be explained by the fact that their collections originate from a broader geographic area, or belong to different cultivar types, in contrast to our collection which originates from a single and more geographically restricted area. However, when Hurtado et al. [35] calculated the average number of alleles/locus only for the Spanish eggplants included in their collection, it was almost similar to ours (5.0). Expected heterozygosity (He) values were high, with an average value of 0.6, quite similar or higher to the values observed in other similar studies, providing further evidence of the utility and highly informative nature of the SSRs used. The estimated values of observed heterozygosity (Ho) (average 0.13) were expected, since eggplant is basically an autogamous plant [76] and thus most heritage and commercial varieties are expected to be homozygous [34]. Muñoz-Falcón et al. [77] detected significantly lower values for Spanish landraces (0.02) and for non-Spanish landraces (0.02), whereas Cericola et al. [34] also found that most eggplant landraces used in their study had low heterozygosity values, and only 38 out of 238 materials had Ho > 0.10. Low Ho was also detected by Augustinos et al. [23] and Liu et al. [36] (0.04 and 0.03, respectively) and Vilanova et al. [50,75]. The high level of homozygosity in eggplant landraces shows that pure lines can easily be derived by individual selection from these materials [75]. Although genetic diversity based on molecular data is dependent on the type and number of markers used and the accessions tested [66], the relatively high PIC values and average number of alleles/locus detected and also the high levels of heterozygosity observed compared to other similar studies, indicate that considerable diversity exists in our eggplant collection. Although, the number of SSR markers we used has been limited, their high PIC values indicate that they provide meaningful information on the diversity of Greek eggplants. Consequently, Greece, together with other countries such as Spain and Turkey, might be considered as part of a secondary center of diversity for the eggplant in the Mediterranean basin. Further analysis with high-throughput methods may provide a confirmation for our hypothesis.
Remarkably, the island accessions seem to be more variable than those from the mainland, as indicated by the average number of alleles/locus (4.2 vs. 3.2) and PIC values (0.53 vs. 0.42). Higher values of Ho were observed in the island accessions when compared to the mainland ones (0.20 vs. 0.03), suggesting that a higher degree of outcrossing occurs in the island accessions than in those from the mainland. Although endemic island populations tend to have less diverse variation in comparison to continental populations [25,31,78], due to founder effect events [78,79], of limited population size and genetic drift [80], the heterozygote excess in our case could be possibly due to selective advantages enjoyed by heterozygote individuals, due to the agroecological conditions of the island environment, as in the case of dill [12]. Another possible explanation is that Greek islands in the Aegean Sea are in proximity with Asia minor and the Turkish coastline and thus genetic exchange could be favored through pollen dispersal, human influence and transportation. Environmental conditions also have a critical impact on the relationship between plants and pollinators. Rain can disrupt this pollen transfer and hinder sexual reproduction in plants [81]. The dry climate and the lack of rainfall in the islands may promote cross-pollination, therefore, increasing the proportion of heterozygotes in the island accessions. Consistent with the PCA analysis of the morphological traits, the genetic STRUCTURE analysis differentiates two genetic clusters, one made up mostly by mainland and another by the island accessions. These results are also confirmed in the PCoA analysis. The fact that some island accessions fell in the mainland genetic cluster and vice versa probably has contributed to the limited differentiation observed in the AMOVA analysis between the mainland and the island groups (11.12%). Discrimination among accessions originating from the islands and the mainland has also been reported for Greek dill by Ninou et al. [12] and for okra landraces by Kyriakopoulou et al. [82].

4.4. Comparison of Morphological and Molecular Diversity

Undoubtedly, even though genetic characterization has become predominant in the last decade due to the efficiency and the plummeting cost of sequencing and genotyping, the morphological phenotyping is still essential and complementary to the genetic assessment [83]. The correlations between morphological and molecular data usually vary depending on the eggplant genotype, the morphological descriptors, and the molecular markers used [35]. In this study, the Mantel test detected a weak, although significant, correlation between genetic and morphological data (r = 0.244, p = 0.002), in agreement with previous eggplant studies where weak or no correlation between genotypic and phenotypic data have been found [23,34,35,77]. However, these results need to be confirmed by additional analyses with larger molecular markers datasets.

5. Conclusions

In conclusion, our results show a remarkable phenotypic and genetic variation existing in the Greek eggplant landraces, suggesting that Greece, among other Mediterranean countries, such as Italy and Spain, might be considered as a secondary center of diversity for eggplant. Both morphological and genetic data revealed clustering of eggplant accessions into two major groups based on the geographic origin (mainland vs. island accessions). Some differentiation signatures in morphological, molecular, and fruit composition were detected between both groups underlying the significant impact of the environment, geographic isolation, and the agroecological conditions in the differentiation and characteristics of the eggplant landraces from the Greek islands and mainland. There was a weak but significant correlation between the morphological and genetic data, indicating that both types of data provide complementary information and therefore they should be taken into consideration for the characterization of plant material. Our results show that local Greek populations of eggplant are a promising material which can be of interest not only for conservation purposes but also for the manipulation of value-added traits through breeding programs and for commercial exploitation.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4395/9/12/887/s1, Table S1: Average values of the morphological and mineral characterization of the eggplant accessions of Greek islands (n = 19) and the mainland (n = 16). Table S2: Percentage of the variance explained by each principal component (PC) of the PCA analysis, correlations between the PCs and the morphological traits evaluated in Greek eggplant accessions and the eigenvalues of the accessions for each PCs. Table S3: SSR alleles identified in the eggplant accessions of Greek islands and the mainland. Table S4: Percentage of the variance explained by each principal component (PC) of the PCoA analysis and eigenvectors in Greek eggplant accessions for each PC. Table S5: Mantel test between morphological and genetic data. Figure S1: Delta K values for 2 to 10 genetic clusters for thirty-six accessions (cultivars and landraces) in Greek eggplant accessions. Delta K was calculated according to Evanno et al. (2005).

Author Contributions

Conceptualization, V.P., J.P., P.G., and G.K.; Investigation, E.C., P.G. and V.C.; Formal analysis, C.P. and P.G.; Methodology, J.P., S.V., and V.P.; Resources, P.M.; supervision, J.P., S.V., G.K. and V.P.; writing (original draft preparation), V.P., P.G., and C.P.; writing (review and editing), J.P., P.G., C.P., S.V., V.C., P.M. and V.P.; Visualization, P.G.

Funding

“PlantUP” (MIS 5002803) which is implemented under the Action “Reinforcement of the Research and Innovation Infrastructure”, funded by the Operational Programme “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014–2020) and co-financed by Greece and the European Union (European Regional Development Fund). Funding was also received from the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement No 677379 (G2P-SOL project: Linking genetic resources, genomes and phenotypes of Solanaceous crops), and from Ministerio de Ciencia, Innovación y Universidades, Agencia Estatal de Investigacion and Fondo Europeo de Desarrollo Regional (grant RTI-2018-094592-B-100 from MCIU/AEI/FEDER, UE). Pietro Gramazio is grateful to Universitat Politècnica de València and to Japan Society for the Promotion of Science for their respective postdoctoral grants (PAID-10-18 and FY2019 JSPS Postdoctoral Fellowship for Research in Japan [Standard]).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kougioumoutzis, K.; Valli, A.T.; Georgopoulou, E.; Simaiakis, S.M.; Triantis, K.A.; Trigas, P. Network biogeography of a complex island system: The Aegean Archipelago revisited. J. Biogeogr. 2017, 44, 651–660. [Google Scholar] [CrossRef]
  2. Graham, N.R.; Gruner, D.S.; Lim, J.Y.; Gillespie, R.G. Island ecology and evolution: Challenges in the Anthropocene. Environ. Conserv. 2017, 44, 323–335. [Google Scholar] [CrossRef] [Green Version]
  3. Price, J.P.; Otto, R.; Menezes de Sequeira, M.; Kueffer, C.; Schaefer, H.; Caujapé-Castells, J.; Fernández-Palacios, J.M. Colonization and diversification shape species–area relationships in three Macaronesian archipelagos. J. Biogeogr. 2018, 45, 2027–2039. [Google Scholar] [CrossRef]
  4. Warren, B.H.; Simberloff, D.; Ricklefs, R.E.; Aguilée, R.; Condamine, F.L.; Gravel, D.; Morlon, H.; Mouquet, N.; Rosindell, J.; Casquet, J.; et al. Islands as model systems in ecology and evolution: Prospects fifty years after MacArthur-Wilson. Ecol. Lett. 2015, 18, 200–217. [Google Scholar] [CrossRef]
  5. Douma, C.; Koutis, K.; Thanopoulos, R.; Tsigou, R.; Galanidis, A.; Bebeli, P.J. Diversity of agricultural plants on Lesvos Island (Northeast Aegean, Greece) with emphasis on fruit trees. Sci. Hortic. 2016, 210, 65–84. [Google Scholar] [CrossRef]
  6. Hagenblad, J.; Leino, M.W.; Hernández Afonso, G.; Afonso Morales, D. Morphological and genetic characterization of barley (Hordeum vulgare L.) landraces in the Canary Islands. Genet. Resour. Crop. Evol. 2019, 66, 465–480. [Google Scholar] [CrossRef] [Green Version]
  7. Médail, F. The specific vulnerability of plant biodiversity and vegetation on Mediterranean islands in the face of global change. Reg. Environ. Chang. 2017, 17, 1775–1790. [Google Scholar] [CrossRef] [Green Version]
  8. Hellenic Statistical Authority. 2019. Available online: http://www.statistics.gr/en/home/ (accessed on 1 August 2019).
  9. Sfenthourakis, S.; Triantis, K.A. The Aegean archipelago: A natural laboratory of evolution, ecology and civilisations. J. Biol. Res. 2017, 24, 4. [Google Scholar] [CrossRef] [Green Version]
  10. Strid, A. Phytogeographia Aegaea and the Flora Hellenica Database. Ann. Nat. Mus. Wien. 1996, 98, 279–289. [Google Scholar]
  11. Dimopoulos, P.; Raus, T.; Bergmeier, E.; Constantinidis, T.; Iatrou, G.; Kokkini, S.; Strid, A.; Tzanoudakis, D. Vascular plants of Greece: An annotated checklist. Willdenowia 2016, 46, 301–347. [Google Scholar] [CrossRef] [Green Version]
  12. Ninou, E.G.; Mylonas, I.G.; Tsivelikas, A.L.; Ralli, P.E. Phenotypic diversity of Greek dill (Anethum graveolens L.) landraces. Acta Agric. Scand. Sect. B Soil Plant. Sci. 2017, 67, 318–325. [Google Scholar]
  13. Tsanakas, G.F.; Mylona, P.V.; Koura, K.; Gleridou, A.; Polidoros, A.N. Genetic diversity analysis of the Greek lentil (Lens culinaris) landrace “Eglouvis” using morphological and molecular markers. Plant Genet. Resour. Characterisation Util. 2018, 16, 469–477. [Google Scholar] [CrossRef]
  14. Dwivedi, S.; Goldman, I.; Ortiz, R. Pursuing the potential of heirloom cultivars to improve adaptation, nutritional, and culinary features of food crops. Agronomy 2019, 9, 441. [Google Scholar] [CrossRef] [Green Version]
  15. Bota, J.; Conesa, M.À.; Ochogavia, J.M.; Medrano, H.; Francis, D.M.; Cifre, J. Characterization of a landrace collection for Tomàtiga de Ramellet (Solanum lycopersicum L.) from the Balearic Islands. Genet. Resour. Crop. Evol. 2014, 61, 1131–1146. [Google Scholar] [CrossRef]
  16. 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]
  17. Schmidt, S.B.; George, T.S.; Brown, L.K.; Booth, A.; Wishart, J.; Hedley, P.E.; Martin, P.; Russell, J.; Husted, S. Ancient barley landraces adapted to marginal soils demonstrate exceptional tolerance to manganese limitation. Ann. Bot. 2019, 123, 831–843. [Google Scholar] [CrossRef] [Green Version]
  18. Missio, J.C.; Rivera, A.; Figàs, M.R.; Casanova, C.; Camí, B.; Soler, S.; Simó, J. A comparison of landraces vs. Modern varieties of lettuce in organic farming during the winter in the Mediterranean area: An approach considering the viewpoints of breeders, consumers, and farmers. Front. Plant Sci. 2018, 9, 1491. [Google Scholar] [CrossRef]
  19. Petropoulos, S.A.; Barros, L.; Ferreira, I.C.F.R. Editorial: Rediscovering local landraces: Shaping horticulture for the future. Front. Plant Sci. 2019, 10, 126. [Google Scholar] [CrossRef] [Green Version]
  20. Karanikolas, P.; Bebeli, P.J.; Thanopoulos, R. Farm economic sustainability and agrobiodiversity: Identifying viable farming alternatives during the economic crisis in Greece. J. Environ. Econ. Policy 2018, 7, 69–84. [Google Scholar] [CrossRef]
  21. McDorman, B.; Thomas, S. The Importance of Saving Seeds. In Promoting Biodiversity in Food Systems; Hawkins, I.W., Ed.; CRC PressI Llc: Boca Ratón, FL, USA, 2018; pp. 189–210. [Google Scholar]
  22. FAO STATISTICAL DATABASES. 2017. Available online: http://www.fao.org/faostat/ (accessed on 1 August 2019).
  23. Augustinos, A.A.; Petropoulos, C.; Karasoulou, V.; Bletsos, F.; Papasotiropoulos, V. Assessing diversity among traditional Greek and foreign eggplant cultivars using molecular markers and morphometrical descriptors. Span. J. Agric. Res. 2016, 14, e0710. [Google Scholar] [CrossRef] [Green Version]
  24. Thomas, K.; Thanopoulos, R.; Knüpffer, H.; Bebeli, P.J. Plant genetic resources of Lemnos (Greece), an isolated island in the Northern Aegean Sea, with emphasis on landraces. Genet Resour. Crop. Evol. 2012, 59, 1417–1440. [Google Scholar] [CrossRef]
  25. García-Verdugo, C.; Sajeva, M.; La Mantia, T.; Harrouni, C.; Msanda, F.; Caujapé-Castells, J. Do island plant populations really have lower genetic variation than mainland populations? Effects of selection and distribution range on genetic diversity estimates. Mol. Ecol. 2015, 24, 726–741. [Google Scholar] [CrossRef] [PubMed]
  26. Hiraoka, Y.; Tamaki, I.; Watanabe, A. The origin of wild populations of Toxicodendron succedaneum on mainland Japan revealed by genetic variation in chloroplast and nuclear DNA. J. Plant Res. 2018, 131, 225–238. [Google Scholar] [CrossRef] [PubMed]
  27. Jiménez, A.; Weigelt, B.; Santos-Guerra, A.; Caujapé-Castells, J.; Fernández-Palacios, J.M.; Conti, E. Surviving in isolation: Genetic variation, bottlenecks and reproductive strategies in the Canarian endemic Limonium macrophyllum (Plumbaginaceae). Genetica 2017, 145, 91–104. [Google Scholar] [CrossRef]
  28. Wheelwright, N.T.; Begin, E.; Ellwanger, C.; Taylor, S.H.; Stone, J.L. Minimal loss of genetic diversity and no inbreeding depression in blueflag iris (Iris versicolor) on islands in the Bay of Fundy. Botany 2016, 94, 543–554. [Google Scholar] [CrossRef] [Green Version]
  29. Hedrén, M.; Olofsson, S.N.; Paun, O. Orchid colonization: Multiple parallel dispersal events and mosaic genetic structure in Dactylorhiza majalis ssp. lapponica on the Baltic island of Gotland. Ann. Bot. 2018, 122, 1019–1032. [Google Scholar]
  30. Hufford, K.M.; Mazer, S.J.; Hodges, S.A. Genetic variation among mainland and island populations of a native perennial grass used in restoration. AoB Plants 2014, 6. [Google Scholar] [CrossRef] [Green Version]
  31. Mcglaughlin, M.E.; Wallace, L.E.; Wheeler, G.L.; Bresowar, G.; Riley, L.; Britten, N.R.; Helenurm, K. Do the island biogeography predictions of MacArthur and Wilson hold when examining genetic diversity on the near mainland California Channel Islands? Examples from endemic Acmispon (Fabaceae). Bot. J. Linn. Soc. 2014, 174, 289–304. [Google Scholar] [CrossRef] [Green Version]
  32. Idrissi, O.; Piergiovanni, A.R.; Toklu, F.; Houasli, C.; Udupa, S.M.; De Keyser, E.; Van Damme, P.; De Riek, J. Molecular variance and population structure of lentil (Lens culinaris Medik.) landraces from Mediterranean countries as revealed by simple sequence repeat DNA markers: Implications for conservation and use. Plant Genet. Resour. Characterisation Util. 2018, 16, 249–259. [Google Scholar] [CrossRef]
  33. Acquadro, A.; Barchi, L.; Gramazio, P.; Portis, E.; Vilanova, S.; Comino, C.; Plazas, M.; Prohens, J.; Lanteri, S. Coding SNPs analysis highlights genetic relationships and evolution pattern in eggplant complexes. PLoS ONE 2017, 12, e0180774. [Google Scholar] [CrossRef]
  34. Cericola, F.; Portis, E.; Toppino, L.; Barchi, L.; Acciarri, N.; Ciriaci, T.; Sala, T.; Rotino, G.L.; Lanteri, S. The population structure and diversity of eggplant from Asia and the Mediterranean Basin. PLoS ONE 2013, 8, e73702. [Google Scholar] [CrossRef] [Green Version]
  35. Hurtado, M.; Vilanova, S.; Plazas, M.; Gramazio, P.; Fonseka, H.H.; Fonseka, R.; Prohens, J. Diversity and relationships of eggplants from three geographically distant secondary centers of diversity. PLoS ONE 2012, 7, e41748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Liu, J.; Yang, Y.; Zhou, X.; Bao, S.; Zhuang, Y. Genetic diversity and population structure of worldwide eggplant (Solanum melongena L.) germplasm using SSR markers. Genet. Resour. Crop. Evol. 2018, 65, 1663–1670. [Google Scholar] [CrossRef]
  37. Rodriguez-Jimenez, J.; Amaya-Guerra, C.; Baez-Gonzalez, J.; Aguilera-Gonzalez, C.; Urias-Orona, V.; Nino-Medina, G. Physicochemical, functional, and nutraceutical properties of eggplant flours obtained by different drying methods. Molecules 2018, 23, 3210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Arivalagan, M.; Bhardwaj, R.; Gangopadhyay, K.K.; Prasad, T.V.; Sarkar, S.K. Mineral composition and their genetic variability analysis in eggplant (Solanum melongena L.) germplasm. J. Appl. Bot. Food Qual. 2013, 86, 99–103. [Google Scholar]
  39. Raigón, M.D.; Prohens, J.; Muñoz-Falcón, J.E.; Nuez, F. Comparison of eggplant landraces and commercial varieties for fruit content of phenolics, minerals, dry matter and protein. J. Food Compos. Anal. 2008, 21, 370–376. [Google Scholar] [CrossRef]
  40. Arivalagan, M.; Gangopadhyay, K.K.; Kumar, G.; Bhardwaj, R.; Prasad, T.V.; Sarkar, S.K.; Roy, A. Variability in mineral composition of Indian eggplant (Solanum melongena L.) genotypes. J. Food Compos. Anal. 2012, 26, 173–176. [Google Scholar] [CrossRef]
  41. Ranil, R.H.G.; Niran, H.M.L.; Plazas, M.; Fonseka, R.M.; Fonseka, H.H.; Vilanova, S.; Andújar, I.; Gramazio, P.; Fita, A.; Prohens, J. Improving seed germination of the eggplant rootstock Solanum torvum by testing multiple factors using an orthogonal array design. Sci. Hortic. 2015, 193, 174–181. [Google Scholar] [CrossRef] [Green Version]
  42. van der Weerden, G.M.; Barendse, G.W.M. A web-based searchable database developed for the EGGNET project and applied to the radboud university Solanaceae database. Acta Hortic. 2007, 745, 503–506. [Google Scholar] [CrossRef]
  43. Kaushik, P.; Prohens, J.; Vilanova, S.; Gramazio, P.; Plazas, M. Phenotyping of eggplant wild relatives and interspecific hybrids with conventional and phenomics descriptors provides insight for their potential utilization in breeding. Front. Plant Sci. 2016, 7, 677. [Google Scholar] [CrossRef]
  44. IBPGR. Descriptors for Eggplant; International Board for Plant Genetic Resources: Rome, Italy, 1990. [Google Scholar]
  45. Kumar, G.; Meena, B.L.; Kar, R.; Tiwari, S.K.; Gangopadhyay, K.K.; Bisht, I.S.; Mahajan, R.K. Morphological diversity in brinjal (Solanum melongena L.) germplasm accessions. Plant Genet. Res. 2008, 6, 232–236. [Google Scholar] [CrossRef]
  46. Campbell, C.R.; Plank, C.O. Preparation of Plant Tissue for Laboratory Analysis. In Handbook of Reference Methods for Plant Analysis; Kalra, Y.P., Ed.; CRC Press: Boca Ratón, FL, USA, 1998; pp. 37–50. [Google Scholar]
  47. Ihaka, R.; Gentleman, R. R: A Language for Data Analysis and Graphics. J. Comput. Graph. Stat. 1996, 5, 299–314. [Google Scholar]
  48. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; ISBN 978-3-319-24277-4. [Google Scholar]
  49. Alonso, D.; Gramazio, P.; Plazas, M.; Villanueva, G.; García-Fortea, E.; Díez, M.J.; Prohens, J.; Vilanova, S. Protocolo innovador para la obtención de ADN genómico vegetal de alta calidad. In Proceedings of the Congress “I Congreso de Jóvenes Investigadores en Ciencias Agroalimentarias”, Almeria, Spain, 20 December 2018. [Google Scholar]
  50. Vilanova, S.; Manzur, J.P.; Prohens, J. Development and characterization of genomic simple sequence repeat markers in eggplant and their application to the study of diversity and relationships in a collection of different cultivar types and origins. Mol. Breed. 2012, 30, 647–660. [Google Scholar] [CrossRef]
  51. Nunome, T.; Negoro, S.; Kono, I.; Kanamori, H.; Miyatake, K.; Yamaguchi, H.; Ohyama, A.; Fukuoka, H. Development of SSR markers derived from SSR-enriched genomic library of eggplant (Solanum melongena L.). Theor. Appl. Genet. 2009, 119, 1143–1153. [Google Scholar] [CrossRef] [PubMed]
  52. Liu, K.; Muse, S.V. PowerMaker: An integrated analysis environment for genetic maker analysis. Bioinformatics 2005, 21, 2128–2129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Peakall, P.; Smouse, R. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [Green Version]
  54. Botstein, D.; White, R.L.; Skolnick, M.; Davis, R.W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 1980, 32, 314–331. [Google Scholar]
  55. Nei, M. Genetic distance between populations. Am. Nat. 1972, 106, 283–292. [Google Scholar] [CrossRef]
  56. Wright, S. The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution 1965, 19, 395–420. [Google Scholar] [CrossRef]
  57. Pritchard, J.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar]
  58. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Earl, D.A.; VonHoldt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  60. Remington, D.L.; Thornsberry, J.M.; Matsuoka, Y.; Wilson, L.M.; Whitt, S.R.; Doebley, J.; Kresovich, S.; Goodman, M.M.; Buckler, I.V.E.S. Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc. Natl. Acad. Sci. USA 2001, 98, 11479–11484. [Google Scholar] [CrossRef] [Green Version]
  61. Beerling, D.J.; Leake, J.R.; Long, S.P.; Scholes, J.D.; Ton, J.; Nelson, P.N.; Bird, M.; Kantzas, E.; Taylor, L.L.; Sarkar, B.; et al. Farming with crops and rocks to address global climate, food and soil security. Nat. Plants 2018, 4, 138–147. [Google Scholar] [CrossRef] [PubMed]
  62. Marconi, G.; Ferradini, N.; Russi, L.; Concezzi, L.; Veronesi, F.; Albertini, E. Genetic characterization of the apple germplasm collection in central Italy: The value of local varieties. Front. Plant Sci. 2018, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Wang, L.; Li, J.; Zhao, J.; He, C. Evolutionary developmental genetics of fruit morphological variation within the Solanaceae. Front. Plant Sci. 2015, 6, 248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Prohens, J.; Blanca, J.M.; Nuez, F. Morphological and molecular variation in a collection of eggplants from a secondary center of diversity: Implications for conservation and breeding. J. Am. Soc. Hortic. Sci. 2005, 130, 54–63. [Google Scholar] [CrossRef] [Green Version]
  65. Boyaci, H.F.; Topcu, V.; Tepe, A.; Yildirim, I.K.; Oten, M.; Aktas, A. Morphological and molecular characterization and relationships of Turkish local eggplant heirlooms. Not. Bot. Horti Agrobot. Cluj Napoca 2015, 43, 100–107. [Google Scholar] [CrossRef] [Green Version]
  66. Tümbilen, Y.; Frary, A.; Mutlu, S.; Doganlar, S. Genetic diversity in Turkish eggplant (Solanum melongena) varieties as determined by morphological and molecular analyses. Int. Res. J. Biotechnol. 2011, 2, 16–25. [Google Scholar]
  67. Muñoz-Falcón, J.E.; Prohens, J.; Rodríguez-Burruezo, A.; Nuez, F. Potential of local varieties and their hybrids for the improvement of eggplant production in the open field and greenhouse cultivation. J. Food Agric. Environ. 2008, 6, 83–88. [Google Scholar]
  68. Muñoz-Falcón, J.E.; Prohens, J.; Vilanova, S.; Nuez, F. Characterization, diversity, and relationships of the Spanish striped (Listada) eggplants: A model for the enhancement and protection of local heirlooms. Euphytica 2008, 164, 405–419. [Google Scholar] [CrossRef]
  69. Polignano, G.; Uggenti, P.; Bisignano, V.; Gatta, C. Della Genetic divergence analysis in eggplant (Solanum melongena L.) and allied species. Genet. Resour. Crop. Evol. 2010, 57, 171–181. [Google Scholar] [CrossRef]
  70. Russo, V.M. Cultural methods and mineral content of eggplant Solanum melongena fruit. J. Sci. Food Agric. 1996, 71, 119–123. [Google Scholar] [CrossRef]
  71. Raigón, M.D.; Rodríguez-Burruezo, A.; Prohens, J. Effects of organic and conventional cultivation methods on composition of eggplant fruits. J. Agric. Food Chem. 2010, 58, 6833–6840. [Google Scholar] [CrossRef] [PubMed]
  72. Ge, H.; Liu, Y.; Jiang, M.; Zhang, J.; Han, H.; Chen, H. Analysis of genetic diversity and structure of eggplant populations (Solanum melongena L.) in China using simple sequence repeat markers. Sci. Hortic. 2013, 162, 71–75. [Google Scholar] [CrossRef]
  73. 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]
  74. Demir, K.; Bakir, M.; Sarikamiş, G.; Acunalp, S. Genetic diversity of eggplant (Solanum melongena) germplasm from Turkey assessed by SSR and RAPD markers. Genet. Mol. Res. 2010, 9, 1568–1576. [Google Scholar] [CrossRef]
  75. Vilanova, S.; Hurtado, M.; Cardona, A.; Plazas, M.; Gramazio, P.; Herraiz, F.J.; Andújar, I.; Prohens, J. Genetic diversity and relationships in local varieties of eggplant from different cultivar groups as assessed by genomic SSR markers. Not. Bot. Horti Agrobot. Cluj Napoca 2014, 42, 59–65. [Google Scholar] [CrossRef] [Green Version]
  76. Pessarakli, M.; Dris, R. Pollination and breeding of eggplants. J. Food Agric. Environ. 2004, 2, 218–219. [Google Scholar]
  77. Muñoz-Falcón, J.E.; Prohens, J.; Vilanova, S.; Nuez, F. Diversity in commercial varieties and landraces of black eggplants and implications for broadening the breeders’ gene pool. Ann. Appl. Biol. 2009, 154, 453–465. [Google Scholar] [CrossRef]
  78. Frankham, R. Do island populations have less genetic variation than mainland populations? Heredity 1997, 78, 311–327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Whittaker, R.; Fernández-Palacios, J. Island Biogeography: Ecology, Evolution, and Conservation; Oxford University Press: Oxford, UK, 2007; ISBN 978-0198566120. [Google Scholar]
  80. Hargreaves, S.; Maxted, N.; Hirano, R.; Abberton, M.; Skøt, L.; Ford-Lloyd, B.V. Islands as refugia of Trifolium repens genetic diversity. Conserv. Genet. 2010, 11, 1317–1326. [Google Scholar] [CrossRef]
  81. Lawson, D.A.; Rands, S.A. The effects of rainfall on plant–pollinator interactions. Arthropod Plant Interact. 2019, 13, 561–569. [Google Scholar] [CrossRef] [Green Version]
  82. Kyriakopoulou, O.G.; Arens, P.; Pelgrom, K.T.B.; Karapanos, I.; Bebeli, P.; Passam, H.C. Genetic and morphological diversity of okra (Abelmoschus esculentus [L.] Moench.) genotypes and their possible relationships, with particular reference to Greek landraces. Sci. Hortic. 2014, 171, 58–70. [Google Scholar] [CrossRef]
  83. Caguiat, X.G.I.; Hautea, D.M. Genetic diversity analysis of eggplant (Solanum melongena L.) and related wild species in the Philippines using morphological and SSR markers. Sabrao J. Breed. Genet. 2014, 46, 183–201. [Google Scholar]
Agronomy 09 00887 g001 550
Figure 1. Distribution map of the 36 eggplant Greek accessions from the mainland and islands used in this study. Study codes as in Table 1.
Figure 1. Distribution map of the 36 eggplant Greek accessions from the mainland and islands used in this study. Study codes as in Table 1.
Agronomy 09 00887 g001
Agronomy 09 00887 g002 550
Figure 2. Principal component analysis (PCA) similarities based on the morphological characterization of 35 eggplant accessions used in this study (A) and principal coordinates analysis (PCoA) based on the genotyping characterization of all the accessions used in this study from the mainland (16) and island (20) Greek areas (B). The first and second principal coordinates (PC) are displayed. Study codes as in Table 1.
Figure 2. Principal component analysis (PCA) similarities based on the morphological characterization of 35 eggplant accessions used in this study (A) and principal coordinates analysis (PCoA) based on the genotyping characterization of all the accessions used in this study from the mainland (16) and island (20) Greek areas (B). The first and second principal coordinates (PC) are displayed. Study codes as in Table 1.
Agronomy 09 00887 g002
Agronomy 09 00887 g003 550
Figure 3. Genetic population structure of the 36 Greek eggplant accessions from the mainland and island areas based on SSRs data (K = 2). Each individual is represented by a vertical bar that is partitioned into colored segments that represent the estimated ratio of membership of an individual to a cluster. Cluster one is in yellow while cluster 2 is in orange. Study codes as in Table 1.
Figure 3. Genetic population structure of the 36 Greek eggplant accessions from the mainland and island areas based on SSRs data (K = 2). Each individual is represented by a vertical bar that is partitioned into colored segments that represent the estimated ratio of membership of an individual to a cluster. Cluster one is in yellow while cluster 2 is in orange. Study codes as in Table 1.
Agronomy 09 00887 g003
Table
Table 1. Names/collection number, geographic origin, collection site, geographical coordinates and status of the materials used for the study of morphological and molecular variation (SSRs) in a collection of Greek eggplants.
Table 1. Names/collection number, geographic origin, collection site, geographical coordinates and status of the materials used for the study of morphological and molecular variation (SSRs) in a collection of Greek eggplants.
Accession CodeStudy CodeTerritoryGeographic RegionCollection SiteLongitudeLatitudeAltitude (m)Status
ANP-025/07I-1IslandAegean SeaAmorgos, Aegiali36°90′ Ν25°99′ Ε298Landrace
X-034/06I-2IslandAegean SeaChios, Ag. Georgios Sikousis38°19′ N26°03′ E369Landrace
IS-031/07I-3IslandAegean SeaIkaria, Droutsoulas37°36′ N26°11′ E410Landrace
LKK-094/07I-4IslandAegean SeaKos, Antimacheia36°48′ N27°05′ E133Landrace
X-116/06I-5IslandAegean SeaLemnos, Kontopouli39°55′ N25°20′ E40Landrace
White LerosI-6IslandAegean SeaLeros Island37°09′ N26°49′ E32Cultivar
Wide PurpleI-7IslandAegean SeaLeros Island37°08′ N26°51′ E6Landrace
Long PurpleI-8IslandAegean SeaLeros Island37°08′ N26°51′ E6Landrace
LKK-008/07I-9IslandAegean SeaLeros, Kamara37°09′ N26°49′ E29Landrace
M-039/06I-10IslandAegean SeaLesvos, Paleokipos39°05′ N26°45′ E58Landrace
M-069/06I-11IslandAegean SeaLesvos, Keramia39°12′ N26°42′ E14Landrace
MFS-030/07I-12IslandAegean SeaMilos, Mitakas36°44′ N24°29′ E50landrace
ANP-180/07I-13IslandAegean SeaNaxos, Agia Anna37°04′ Ν25°21′ Ε4Landrace
ANP-215/07I-14IslandAegean SeaParos, Lefkes37°07′ Ν25°12′ Ε210Landrace
HL-027/07I-15IslandAegean SeaSantorini, Vourvoulos36°26′ N25°26′ E116Landrace
SAS-078/07I-16IslandAegean SeaSkopelos, Chora39°07′ N23°43′ E41Landrace
ATS-110/06I-17IslandAegean SeaSyros Island, Chrousa37°24′ N24°55′ E320Landrace
GRC-002/08I-18IslandIonian SeaCorfu, Skripero39°42′ N19°46′ E139Landrace
HL-237/07I-19IslandCreteIraklion, Moni Savvathianon35°37′ N25°00′ E467Landrace
IK-082/06I-20IslandIonian SeaIthaki Island, Perachori38°34′ N20°71′ E343Landrace
K-153/06M-1MainlandMacedoniaGrevena, Pontini40°04′ N21°40′ E819Landrace
EMIM-2MainlandMacedoniaIPGRB/HAO DEMETER40°32′ N22°59′ Ε19Cultivar
LagkadaM-3MainlandMacedoniaIPGRB/HAO DEMETER40°32′ N22°59′ Ε19Cultivar
KD-053/07M-4MainlandMacedoniaKavala, Platanotopos40°50′ N24°03′ E244Landrace
F-154/06M-5MainlandMacedoniaKastoria, Ampelokipoi40°46′ N21°31′ E638Landrace
SK-044/066M-6MainlandMacedoniaKilkis, Eptalofos41°00′ N23°08′ E451Landrace
K-054/06M-7MainlandMacedoniaKozani, Anarachi40°29′ N21°34′ E727Landrace
VG-011/083M-8MainlandMacedoniaSerres41°05′ N23°35′ E45Landrace
SK-056/06M-9MainlandMacedoniaSerres, Platanakia41°17′ N22°56′ E313Landrace
ScoutariM-10MainlandMacedoniaSerres/Skoutari41°01′ N23°31′ E14Cultivar
P-175/06M-11MainlandPeloponneseLakonia, Lyra36°38′ N22°57′ E400Landrace
TsakonikiM-12MainlandPeloponneseLeonidio37°10′ N22°51′ E40Cultivar
P-084/06M-13MainlandPeloponneseMessinia, Kakana37°18′ Ν21°44′ E136Landrace
T-099/06M-14MainlandThessalyKarditsa, Neo Ikonio39°27′ N22°21′ E107Landrace
T-527/06M-15MainlandThessalyTrikala, Megarxis39°36′ N21°45′ E142Landrace
GRC 1430/04M-16MainlandEpirusIoannina, Pogoni Vasiliko40°00′ N20°35′ E805Landrace
Table
Table 2. Descriptors used for phenotyping of Greek eggplant accessions and their range/scale. All descriptors based on IBPGR descriptors [44], except NOLFF [42,43] and DSFS [45].
Table 2. Descriptors used for phenotyping of Greek eggplant accessions and their range/scale. All descriptors based on IBPGR descriptors [44], except NOLFF [42,43] and DSFS [45].
Descriptor CodeDescriptor NameDescriptor Scale/Unit
Plant descriptors
PGHPlant Growth Habit3–7 (3 = upright; 7 = prostrate)
PHPlant Heightcm
NOLFFNumber of Leaves to First Flowernumber
DSFSDays Since Fruit Setnumber
Leaf descriptors
LPLLeaf Petiole Lengthcm
LBLeLeaf Blade Lengthcm
LBWLeaf Blade Widthcm
LBLoLeaf Blade Lobing1–9 (1 = very weak; 9 = very strong)
LSSLeaf Surface Shape1–9 (1 = flat; 9 = very convex or bullate)
LPLeaf Prickles0–9 (0 = none; 9 = more than 20)
Flower descriptors
NOFPINumber of Flowers Per Inflorescencenumber
Fruit descriptors
FCPFruit Calyx Prickles0–9 (0 = none; 9 = more than 30)
FPLFruit Pedicel Lengthcm
FLFruit Lengthcm
FBFruit Breadthmm
FLBRFruit Length to Breadth Ratio1–9 (1 = broader than long; 9 = several times longer than broad)
FSFruit Shape3–7 (Position of widest part of fruit: 3 = ¼ way from base to tip; 5 = ½ way from base to tip; 7 = ¾ way from base to tip)
FCFruit Curvature1–9 (1 = none, fruit straight; 9 = U shaped)
FWFruit Weightg
FPCFruit Predominant Color1–9 (1 = milk white; 9 = black)
FACFruit Additional Color1–9 (1 = milk white; 9 = black)
FFCFruit Flesh Color3–7 (3 = white; 7 = green)
Table
Table 3. Characteristics of the SSR markers used for the molecular characterization of Greek eggplant accessions.
Table 3. Characteristics of the SSR markers used for the molecular characterization of Greek eggplant accessions.
SSR LocusMotifForward Primer and Reverse Primer (5′→3′)Size Range (bp)Tº AnnealingDyeSource
csm4(GA)15GCGTACCAATTCTAACCACAAG238–25460PETVilanova et al. [50]
GTAATCCGCTTCCCATTTCTC
csm27(GA)23TGTTTGGAGGTGAGGGAAAG193–21060VICVilanova et al. [50]
TCCAACTCACCGGAAAAATC
csm32(AG)23TCGAAAGTACAGCGGAGAAAG248–25460NEDVilanova et al. [50]
GGGGGTTTGATTTTCATTTTC
emi02c21(AC)13A(TA)4TGTGAGGAGAAGAATCAGAGGATCA126–13660VICNunome et al. [51]
CGCGACTAAGTTTTGTTCCTGAAA
eme11f04(TC)16ACCCCCAAATCAAATCATTTACCC88–10060FAMNunome et al. [51]
GGCATGGTTAGGGTTTTTAGCGTT
Table
Table 4. Mean, range (between brackets), and coefficient of variation (CV %) for the morphological and mineral descriptors studied in eggplant accessions of the Greek islands (n = 19) and the mainland (n = 16) and significance of mean differences among these two groups (p < 0.05).
Table 4. Mean, range (between brackets), and coefficient of variation (CV %) for the morphological and mineral descriptors studied in eggplant accessions of the Greek islands (n = 19) and the mainland (n = 16) and significance of mean differences among these two groups (p < 0.05).
MainlandIslandsp-Value
MeanCV (%)MeanCV (%)
(Range)(Range)
Plant descriptors
PGH4.3623.705.0921.500.052
(3.00–6.50)(3.00–7.00)
PH80.5916.3071.7221.900.083
(56.83–106.00)(46.50–106.17)
NOLFF4.9713.305.4112.500.060
(4.00–6.67)(4.20–6.83)
DSFS50.518.2048.269.100.129
(42.83–58.33)(41.00–55.67)
Leaf descriptors
LPL9.8317.6010.2219.400.541
(6.79–12.29)(6.53–14.63)
LBLe31.518.4031.410.500.918
(27.50–36.97)(25.57–37.48)
LBW20.4311.5022.149.800.032 *
(16.35–25.62)(18.16–25.63)
LBLo3.3823.804.3721.900.004 *
(3.00–5.00)(3.00–5.00)
LSS4.8827.905.4219.800.203
(3.00–7.00)(3.00–7.00)
LP0.18126.400.21125.700.892
(0.00–0.92)(0.00–0.92)
Flower descriptors
NOFPI1.1511.101.2216.900.388
(1.00–1.40)(1.00–1.60)
Fruit descriptors
FCP1.5183.903.3257.900.003 *
(0.00–4.00)(0.33–7.00)
FPL6.3316.706.2432.100.840
(4.91–8.22)(3.58–9.72)
FL18.4613.0018.442.800.246
(14.50–22.32)(11.2–43.85)
FB61.9316.2085.120.300.000 *
(50.91–87.78)(52.75–110.76)
FLBR7.4911.505.244.400.002 *
(5.00–8.50)(1.33–8.33)
FS4.697.004.3612.200.106
(4.00–5.33)(3.33–5.00)
FC3.332.301.8153.700.000 *
(1.00–5.00)(1.00–4.00)
FW246.3319.70329.6420.800.000 *
(177.32–368.40)(215.66–500.67)
FPC6.5222.105.8340.500.550
(2.00–8.00)(1.00–8.17)
FAC4.3148.004.7848.200.527
(1.00–7.50)(1.00–8.17)
FFC4.7315.104.4520.600.218
(3.00–5.33)(3.00–6.00)
Minerals
Κ (mg/g DW)25.0513.2123.1411.580.073
(20.07–31.50)(19.80–29.00)
Mg (mg/g DW)2.226.752.0512.190.015 *
(1.96–2.46)(1.63–2.43)
Cu (mg/Kg DW)9.5647.917.6569.670.278
(2.91–18.14)(2.35–26.93)
Fe (mg/Kg DW)31.814.6831.4815.780.850
(25.14–43.94)(22.02–41.84)
Mn (mg/Kg DW)12.258.8211.749.710.198
(9.61–13.67)(9.79–13.84)
Zn (mg/Kg DW)49.7526.8549.3931.140.942
(30.30–77.26)(29.53–76.90)
* indicates a significant difference among islands and the mainland accessions at p-value < 0.05.
Table
Table 5. Diversity statistics of SSR markers for all the Greek eggplant accessions and for the island and the mainland groups. Diversity statistics evaluated include: number of alleles (Na), major allele frequency (f), number of effective alleles (Ne), number of genotypes (Ng), polymorphic information content (PIC), observed heterozygosity (Ho), expected heterozygosity (He), and inbreeding coefficient (F).
Table 5. Diversity statistics of SSR markers for all the Greek eggplant accessions and for the island and the mainland groups. Diversity statistics evaluated include: number of alleles (Na), major allele frequency (f), number of effective alleles (Ne), number of genotypes (Ng), polymorphic information content (PIC), observed heterozygosity (Ho), expected heterozygosity (He), and inbreeding coefficient (F).
SSR LocusNafNeNgPICHoHeF
All accessions (n = 36)
csm470.324.4110.000.760.140.770.83
csm2740.502.637.000.580.170.620.75
csm3250.532.577.000.550.170.610.73
emi02c2140.612.194.000.480.000.541.00
eme11f0430.651.925.000.400.170.480.66
Mean ± SE4.6 ± 0.70.522 ± 0.052.746 ± 0.436.6 ± 1.000.552 ± 0.050.130 ± 0.030.605 ± 0.040.796 ± 0.05
Islands accessions (n = 20)
csm460.334.658.000.750.200.790.76
csm2740.532.285.000.470.250.560.57
csm3250.582.537.000.560.250.610.60
emi02c2130.502.203.000.440.000.551.00
eme11f0430.502.295.000.470.300.560.49
Mean ± SE4.2 ± 0.60.485 ± 0.042.79 ± 0.465.6 ± 0.900.538 ± 0.050.200 ± 0.050.612 ± 0.040.686 ± 0.09
Mainland accessions (n = 16)
csm440.373.445.000.660.070.710.91
csm2730.502.574.000.540.070.610.90
csm3230.472.264.000.460.060.560.89
emi02c2140.751.714.000.390.000.411.00
eme11f0420.941.132.000.110.000.121.00
Mean ± SE3.2 ± 0.40.604 ± 0.102.22 ± 0.393.8 ± 0.500.429 ± 0.090.039 ± 0.010.482 ± 0.100.923 ± 0.02
Table
Table 6. Molecular analysis of variance in 36 Greek eggplant accessions. p-value estimates are based on 999 permutations (α = 0.05). The following sources of variation were considered: between groups (mainland vs. islands), among accessions within a group, and within individuals. Acronyms: df = degrees of freedom, SS = Sums of Squares, MS = mean squared deviations, Fst = genetic differentiation among populations within the total sample, Fis = genetic differentiation among individuals within populations, Fit = genetic differentiation among individuals within the total sample.
Table 6. Molecular analysis of variance in 36 Greek eggplant accessions. p-value estimates are based on 999 permutations (α = 0.05). The following sources of variation were considered: between groups (mainland vs. islands), among accessions within a group, and within individuals. Acronyms: df = degrees of freedom, SS = Sums of Squares, MS = mean squared deviations, Fst = genetic differentiation among populations within the total sample, Fis = genetic differentiation among individuals within populations, Fit = genetic differentiation among individuals within the total sample.
SourcedfS.S.M.S.Variance ComponentsPercentage of VariationF-Statisticsp-Value
Between groups19.199.190.1911.12Fst 0.1110.001
Among accessions3489.422.631.1669.62Fis 0.7830.001
Within accessions3611.500.320.3219.25Fit 0.8070.001
Total71110.11 1.66

Share and Cite

MDPI and ACS Style

Gramazio, P.; Chatziefstratiou, E.; Petropoulos, C.; Chioti, V.; Mylona, P.; Kapotis, G.; Vilanova, S.; Prohens, J.; Papasotiropoulos, V. Multi-Level Characterization of Eggplant Accessions from Greek Islands and the Mainland Contributes to the Enhancement and Conservation of this Germplasm and Reveals a Large Diversity and Signatures of Differentiation between both Origins. Agronomy 2019, 9, 887. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120887

AMA Style

Gramazio P, Chatziefstratiou E, Petropoulos C, Chioti V, Mylona P, Kapotis G, Vilanova S, Prohens J, Papasotiropoulos V. Multi-Level Characterization of Eggplant Accessions from Greek Islands and the Mainland Contributes to the Enhancement and Conservation of this Germplasm and Reveals a Large Diversity and Signatures of Differentiation between both Origins. Agronomy. 2019; 9(12):887. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120887

Chicago/Turabian Style

Gramazio, Pietro, Eleni Chatziefstratiou, Constantinos Petropoulos, Vasileia Chioti, Photini Mylona, George Kapotis, Santiago Vilanova, Jaime Prohens, and Vasileios Papasotiropoulos. 2019. "Multi-Level Characterization of Eggplant Accessions from Greek Islands and the Mainland Contributes to the Enhancement and Conservation of this Germplasm and Reveals a Large Diversity and Signatures of Differentiation between both Origins" Agronomy 9, no. 12: 887. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9120887

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