Next Article in Journal
Functional Conservation and Divergence of Five AP1/FUL-like Genes in Marigold (Tagetes erecta L.)
Next Article in Special Issue
Oligo-FISH Can Identify Chromosomes and Distinguish Hippophaë rhamnoides L. Taxa
Previous Article in Journal
SUMO-Based Regulation of Nuclear Positioning to Spatially Regulate Homologous Recombination Activities at Replication Stress Sites
Previous Article in Special Issue
Cloning, Expression, and Tobacco Overexpression Analyses of a PISTILLATA/GLOBOSA-like (OfGLO1) Gene from Osmanthus fragrans
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 12
10.3390/genes12122007
genes-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

Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS)

by
A. S. M. Faridul Islam
1,
Dean Sanders
2,
Amit Kumar Mishra
3 and
Vijay Joshi
1,4,*
1
Texas A&M AgriLife Research and Extension Center, Uvalde, TX 78801, USA
2
Bioinformatics Resource Center, University of Wisconsin–Madison, Madison, WI 53706, USA
3
Department of Botany, School of Life Sciences, Mizoram University, Aizawl 796004, Mizoram, India
4
Department of Horticultural Sciences, Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Genes 2021, 12(12), 2007; https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122007
Submission received: 27 November 2021 / Revised: 10 December 2021 / Accepted: 14 December 2021 / Published: 17 December 2021
(This article belongs to the Special Issue Tree Genetics and Improvement)
Browse Figures
Versions Notes

Abstract

:
Olives are one of the most important fruit and woody oil trees cultivated in many parts of the world. Olive oil is a critical component of the Mediterranean diet due to its importance in heart health. Olives are believed to have been brought to the United States from the Mediterranean countries in the 18th century. Despite the increase in demand and production areas, only a few selected olive varieties are grown in most traditional or new growing regions in the US. By understanding the genetic background, new sources of genetic diversity can be incorporated into the olive breeding programs to develop regionally adapted varieties for the US market. This study aimed to explore the genetic diversity and population structure of 90 olive accessions from the USDA repository along with six popular varieties using genotyping-by-sequencing (GBS)-generated SNP markers. After quality filtering, 54,075 SNP markers were retained for the genetic diversity analysis. The average gene diversity (GD) and polymorphic information content (PIC) values of the SNPs were 0.244 and 0.206, respectively, indicating a moderate genetic diversity for the US olive germplasm evaluated in this study. The structure analysis showed that the USDA collection was distributed across seven subpopulations; 63% of the accessions were grouped into an identifiable subpopulation. The phylogenetic and principal coordinate analysis (PCoA) showed that the subpopulations did not align with the geographical origins or climatic zones. An analysis of the molecular variance revealed that the major genetic variation sources were within populations. These findings provide critical information for future olive breeding programs to select genetically distant parents and facilitate future gene identification using genome-wide association studies (GWAS) or a marker-assisted selection (MAS) to develop varieties suited to production in the US.

1. Introduction

Olives (Olea europaea L.) are one of the economically important fruit and oil trees contributing to the Mediterranean food diet. Often referred to as ‘liquid gold’ [1], olive oil is a rich source of functional compounds such as hydroxytyrosol, oleuropein, and monounsaturated fatty acids beneficial to human health [2]. Several therapeutic studies have confirmed the utility of olive oil in alleviating the impacts of cardiovascular disease, obesity, metabolic syndrome, type 2 diabetes, and hypertension [3,4,5]. The olive is believed to have been domesticated in the Mediterranean basin about 6000 years ago, subsequently spreading through the Mediterranean countries [3]. Although most commercial olive production is confined to Mediterranean countries, more than 40 countries grow olives including Argentina, the United States (USA), Australia, Chile, and China [3,6]. California is the central oil-producing state in the USA, yielding 67,000 tons of olives at a value of USD 57,909 million [7]. Even though the US produces less than 1% of the world’s olives, it represents the third largest national market for olive oil globally, the most significant market outside the European community. Olives were believed to have been first introduced into the US by Spanish Franciscan missionaries in the late 18th century [8].
Genetic diversity is essential for any crop improvement program. Genetic improvements require in-depth screening to understand the nature and extent of the diversity within the germplasm [9]. The National Clonal Germplasm Repository (NCGR) at the University of California-Davis (UC-Davis), one of the clonal gene banks under the National Plant Germplasm System of the USDA-ARS, holds a collection of olive accessions collected from all over the world. This collection represents a valuable tool for population and evolutionary genetic studies of olives and a source of material for breeding purposes. Correctly identifying olive cultivars is challenging due to the high degree of kinship, clonal variation, mixtures with international cultivars, and exchange of plant material over the centuries. The complex distribution pattern resulted in homonymies and synonymies as well as naming errors of cultivars [10,11,12]. A detailed molecular evaluation of the olive accessions in the US repository would provide insights into the amount and organization of genetic diversity and relationships within and among different accessions. Such studies would enhance the value of regionally adapted germplasm and allow for a better utilization and management of regional challenges such as abiotic stresses (e.g., freezing, drought, and nutrient deficiencies), disease resistance (for example, cotton root rot), and oil quality traits. It is imperative to examine the genetic variation and population structure of the olive germplasm to manage the gene pool effectively, help understand the effect of domestication on genetic diversity [10], and aid in the development of new cultivars.
Several morphological traits as well as biochemical and molecular markers have been used to characterize olive germplasm resources [13]. Various types of molecular markers such as random amplified polymorphism DNA (RAPD), amplified fragment length polymorphism (AFLP), sequence-related amplified polymorphism (SRAP), simple sequence repeat (SSR), inter-simple sequence repeat (ISSR), and single-nucleotide polymorphism (SNP) have been used for genetic studies and food traceability in olives has been extensively reviewed [5,12]. However, SNPs have become a marker of choice for various genetic studies due to their unique features such as their availability throughout the genome, unbiased distribution, biallelic nature, stability (e.g., a low mutation rate), and automation in next-generation sequencing (NGS) techniques [3,10,14,15]. The use of molecular markers facilitates an understanding of the parentage in the US olive collection. In addition, the use of SNP markers allows the study of the genetic correlations among the phenotypes of interest and their heritability as well as providing an estimation of the breeding values. Such studies have already been performed in long-lived perennial plants such as almonds, apricots, and apples [14,15,16,17,18].
Genotyping-By-Sequencing (GBS) is an NGS-based genome-wide SNP discovery and genotyping technique with cost- and time-effective characteristics [19]. It can sequence multiple samples simultaneously by using NGS libraries made with methylation-sensitive restriction enzymes, avoiding genome complexities whilst better sequencing lower copy genic regions in the plant genome [11]. GBS can be applied to any species with or without a reference genome, making it the most popular approach for SNP identification [3,10]. So far, this technology has been used to study genetic characterization, association mapping, QTL mapping, and genomic selection in major crops such as wheat [20], barley [21], maize [22], and rice [23,24] as well as fruit crops such as citrus [25], peach [26], apple [27] and oil palm [28]. GBS-generated SNP markers have been utilized in olives to construct high-density genetic linkage maps, a genetic diversity analysis of the germplasm of European and Mediterranean olives, and genome-wide association mapping (GWAS) [3,5,10,11]. Aside from a single study that used microsatellite markers to examine the genetic diversity within olive accessions maintained in the US NCGR collection [29], limited attempts have been made to characterize the olive accessions in the US repository.
In the present study, we used GBS technology to genotype a collection of olive accessions assembled from the USDA-ARS National Plant Germplasm System (NPGS). The objectives of this study were to: (1) generate SNP markers by GBS technology and evaluate their characteristics; (2) determine the population structure of the USDA germplasm collection; and (3) measure the genetic relationships and sources of the genetic variations. The knowledge of the genetic diversity and relationships within and among the accessions in the US olive repository would serve as a resource for effective conservation, management, and utilization of these accessions as well as developing superior cultivars that fulfil the needs of the US market.

2. Materials and Methods

2.1. Plant Materials

The study comprised 90 olive accessions obtained from the National Clonal Germplasm Repository (NCGR) at the University of California-Davis (UC-Davis) along with samples of 6 regionally popular olive varieties. For efficient rooting, the basal end of the cuttings was dipped in 1000 ppm indole butyric acid [30] for 10 s. After the IBA treatment, the cuttings were inserted in 2.5 × 14 inch Deepot tree pots containing perlite and kept under mist (80 to 90% relative humidity) at the Texas A&M AgriLife Research and Extension Center, Uvalde, TX, USA. The intermittent mist system was operated as needed to maintain uniform moisture around the cuttings. The olive accessions originated from 18 countries. The cuttings with newly sprouted leaves were transferred to new pots for the subsequent management. The 18 countries were categorized into 5 major climatic zones based on the Köppen Climate Classification [31] as tropical, dry, temperate, continental, and polar (Table 1). According to this classification, most samples (69 out of 96 accessions) belonged to the temperate climatic zone, followed by the dry zone (16), tropical zone (3), and the continental zone (2). The six olive accessions with no origin information were considered to be of an unknown origin.

2.2. DNA Extraction and Genotyping-By-Sequencing (GBS) Procedures

Leaf samples were collected from the different accessions in 2 mL centrifuge tubes and flash-frozen in liquid nitrogen. The frozen leaf tissue was homogenized to a fine powder in a Harbil model 5G-HD paint shaker (Harbil, Wheeling, IL, USA) using 3 mm Demag stainless steel balls (Abbott Ball Company, West Hartford, CT, USA). The total DNA was extracted using a DNeasy® Plant Mini Kit (QIAGEN Sciences, Germantown, MD, USA) as per the manufacturer’s protocol and treated with RNase A. The purity of the DNA was analyzed using a NanoPhotometer spectrophotometer (IMPLEN, Westlake Village, CA, USA). An ApeKI restriction enzyme was used to construct the DNA libraries for the GBS. The library construction and sequencing by NovaSeq 6000 (Illumina, San Diego, CA, USA) were performed in the Bioinformatics Resource Center, University of Wisconsin–Madison. For the sequence analysis, low-quality reads and adapter sequences were removed from the raw fastq files using computational pipelines developed at the Bioinformatics Resource Center (BRC) at the University of Wisconsin–Madison (https://www.biotech.wisc.edu/, accessed on 23 February 2021) using a trimming software, Skewer [32]. The raw GBS sequences were processed using a standard TASSEL-GBS pipeline [33]. The details of the TASSEL-GBS Pipeline Version 2 (https://tassel.bitbucket.io/, accessed on 23 February 2021) are shown in Supplementary Figure S1.

2.3. Data Analysis

2.3.1. SNP Discovery

The raw reads were quality trimmed to remove the adapters and low-quality bases (Phred ≥ 20) using Skewer software [32]. Once the quality raw fastq files were generated, the TASSEL Version 2 GBS pipeline was implemented to conduct the GBS analysis. In brief, a unique tag database was created by a GBSSeqtoTagDBPlugin that took quality controlled raw fastq files as the input data and then converted them into fastq files by the TagExportToFastqPlugin for the next alignment step. Bowtie 2 software was then used to align the exported tags against the Olea europaea var. europaea L. reference genome [34] and generate a sequence alignment map file (SAM) [35]. The SAM files were utilized by a SamToGBSdbPlugin to input the mapped genomic coordinates of each tag into the TASSEL database. The SNPs were identified using the aligned tags that were positioned at the same genomic coordinates using the DiscoverySNPCallerPluginV2, which required a MAF > 0.01 and a minimum locus coverage in all taxa of 10% (0.1). In the end, 349,851 unfiltered SNPs were discovered in the GBS analysis.

2.3.2. SNP Marker Properties

To conduct the genetic diversity analysis, 54,075 biallelic SNP markers were finally selected based on the following filtering criteria using TASSEL 5 [36]: missingness rate < 0.5, minor allele frequency (MAF) > 0.05, and heterozygosity < 0.1.
The summary statistics including the minor allele frequency (MAF), gene diversity (GD), and polymorphic information content (PIC) for all SNP markers were calculated using the snpReady package in [37].

2.3.3. Population Molecular Characterization

A model-based (Bayesian) method implemented in STRUCTURE 2.3.4 software [38] was used to infer the most probable number of clusters or subpopulations in our germplasm. The admixture model and correlated allele frequency were used to run five independent runs for each K ranging from 1 to 10 to assign a genotype into a particular subpopulation. For each run, 10,000 and 50,000 replications were used for the burn-in time and Markov Chain Monte Carlo (MCMC), respectively. The result of the STRUCTURE software was then submitted to CLUMPAK [39] to determine the best K by using ΔK values following the method of Evanno et al. [40]. A principal coordinate analysis (PCoA) was performed based on the Euclidean distance method to determine the overall genetic difference among the accessions. Both studies were conducted in the adegenet R package [30]. In addition, an unrooted phylogeny tree using the neighbor-joining method was constructed using MEGA7 [41].
An analysis of molecular variance (AMOVA) was carried out to determine the sources of the genetic variance within and among the populations detected by the STRUCTURE software. The Poppr R package [42] was used to calculate the AMOVA using the Euclidean genetic distance method with 999 permutations to declare the significance of a particular genetic variance. Furthermore, the nucleotide diversity per site and fixation index (Fst, Weir, and Cockerham’s 1984) were calculated to measure the genetic diversity within and between populations, respectively, using VCFtools [43].

3. Results

3.1. Characterization and Distribution of the GBS-Generated SNPs in the Olive Genomes

A total of 96 olives accessions were sequenced and genotyped using GBS. After filtering out the raw reads, the total demultiplexed reads for all the genotypes were 418.78 M with the average reads per accession being 4.36 M. The lowest and highest number of reads was 0.21 M and 11.35 M, respectively (Table S2). After processing the raw reads via the TASSEL-GBS pipeline and applying VCF filtering control thresholds, we were left with a subselected set of 54,075 SNPs. The dataset of 54,075 SNPs was then used for a further genetic diversity analysis. These SNPs were mapped onto 23 chromosomes along with 466 scaffolds. The highest and lowest SNPs mapped per chromosome were 22,870 and 6644 on chromosome 10 and chromosome 23, respectively, with an average of 11,906.22 SNP/chromosome (Figure 1a). On average, 163.11 SNPs were mapped to 466 scaffolds ranging from 1 to 2083 SNPs. Among the 54,075 SNPs, the transitions were more frequent (59%, 31,885 SNPs) than the transversions (41%, 22,190 SNPs) with an overall ratio of 1.44. The C/T transition had the highest frequency (30%) and the C/G transversion had the lowest (7%). The frequency of the two transition types was similar (A/G 29%, C/T 30%). The highest frequency among the transversions was found at A/T (14%). Two transversion SNP types (A/C and G/T) had the same frequency (10%) (Table 2).

3.2. Characterization of the SNP Markers

A total of 54,075 SNP markers were selected, satisfying the filtering criteria mentioned in the Method section, to conduct a genetic diversity analysis of the 96 olive germplasm used in this study. To explain the total variability of each marker, the minor allele frequency (MAF), gene diversity (GD), and polymorphic information content (PIC) were used. Although markers with a MAF < 0.05 were removed, the average MAF value was 0.160 with a minimum of 0.05 and a maximum of 0.50. About half (~44%) of the total markers (23,610 out of 54,075 markers) had a MAF less than or equal to 0.1 (Figure 1b). The mean gene diversity value was 0.244 with a maximum of 0.10 for 4389 markers and a minimum of 0.50 for 1546 markers (Figure 1c). The polymorphism information content (PIC) values ranged from 0.09 to 0.38 with an average of 0.21 (Figure 1d). Although 2451 SNPs had the lowest PIC value, 33% SNPs were found with a PIC value of half of its maximum theoretical PIC value (0.5); i.e., ≥0.25.

3.3. Characterization of the Population and the Genetic Relationships

3.3.1. Structure Analysis and the Genetic Relationships

To understand the pattern of the genetic structure, a Bayesian clustering analysis in STRUCTURE was performed (Figure 2a). A population structure analysis was conducted using K values ranging from 2 to 10 with an admixture model and five independent runs for each K value were performed. An Evanno test was then performed to determine the log-likelihood (LnP(D)) values and ΔK between each K number. From the test, the top ΔK peak was found at K = 7, indicating that the US olive germplasm could be grouped into seven subpopulations with admixture accessions (Figure 2b). With a membership probability threshold of 0.70, a total of 60 olive accessions (63%) were grouped into one of seven subpopulations and the remaining 36 accessions were considered to be an admixture group (Figure 2a). The highest number of accessions that were grouped into a particular subpopulation (Pop7) was 22, followed by 13, 8, and 7 accessions clustered into 3 different subpopulations; Pop3, Pop4, and Pop5, respectively. From the remaining groups, two were composed of three accessions (Pop1 and Pop6), and one (Pop2) had four accessions (Table S1). In terms of the proportion of genotypes per climatic zone and distribution across the seven populations, 39, 20, and 14% of the temperate accessions were grouped into admixture, Pop7, and Pop3. In contrast, 38 and 19% of the dry climatic zone accessions were clustered into Pop7 and Pop3, respectively. The accessions from the continental zone were not grouped into any specific subpopulation outside the admixture group. Similarly, 67% of the tropical accessions did not belong to any subpopulation and the remaining 33% belonged to Pop2 (Figure 2c).
The principal coordinate analysis (PCoA) (Figure 3) agreed with the relationships revealed by the structure analysis. The PCoA based on the SNPs revealed seven clusters of the 96 accessions. Among the seven clusters, Pop1, Pop4, Pop5, and Pop6 were clustered distinctively from the remaining groups, indicating a genetically distinct relationship from the other groups. In contrast, Pop2, Pop3, and Pop7 were clustered into the same ellipsis, suggesting a close genetic relationship among those accessions (Figure 3).
The genetic diversity among the subpopulations identified by the structure analysis and PCoA was computed using the average nucleotide diversity (π) and fixation index (Fst). The highest average π was observed for the whole population (0.246), followed by the admixture group and Pop5. The lowest genetic diversity was found for Pop6 (Figure 4a). We also estimated the average π per site (0.25) across the various chromosomes, which ranged from 0.27 (chromosome 8) to 0.23 (chromosomes 14 and 22), to understand the genome-wide bottleneck effects and genetic diversity (Supplementary Figure S2). Based on the Fst values, Pop2, Pop3, and Pop7 were in the same cluster, indicating their genetic relatedness in agreement with the PCoA. Pop1 and Pop5 were clustered separately whereas Pop4, Pop6, and the admixture were clustered together. Although the nucleotide diversity showed that Pop1 and Pop2 had almost similar π values (0.091 for Pop1 and 0.086 for Pop2), they were grouped into two different clusters, indicating that they were genetically distant but had less variability. Similarly, Pop3 and Pop4 had similar π values but were genetically diverse (Figure 4).

3.3.2. Cluster Analysis

We conducted a genotype-based phylogenetic analysis using the neighbor-joining method implemented in MEGA7 [41] The genotype-based cluster analysis reflected a similar population structure, resulting in seven distinct clusters but not aligning with the climatic zones. For example, all the accessions from the dry and temperate zones were clustered across all subpopulations. Among the three tropical accessions, one grouped in Pop2 and the remaining two were in the admixture. All Continental accessions were also grouped in the admixture (Figure 5).

3.3.3. Analysis of Molecular Variance (AMOVA)

An AMOVA analysis was performed to understand the underlying sources of the genetic variation in the germplasm. When the accessions were divided based upon the population structure as of the first level of stratification, the results of the AMOVA indicated that the genetic variation mainly occurred within a population (67%) whereas 33% of the variation was attributed to the difference among populations (Table 3). However, when the climatic zone was included as a second level of stratification of the population structure to group the accessions, even though the majority of genetic variations arose from within the samples (70%, as shown in Table 3), the primary source of the genetic variation (26%) was the climatic zone within a population than the among the population variation (4%) (Table 4).

3.4. Marker Characteristics across the Populations

For the minor allele frequency (MAF), gene diversity (GD), and polymorphic information content (PIC), the admixture group had the highest MAF, GD, and PIC values overall, followed by Pop5 and Pop7, among the populations. In contrast, Pop6 had the lowest value for each of the characteristics. The remaining Pop1, Pop2, Pop3, and Pop4 had almost similar MAF, GD, and PIC values (Figure 6).

4. Discussion

4.1. GBS Analysis of the Olive Genomes

Knowing the genetic variability in the collection of the available pre-adapted olive genotypes is a prerequisite for the US olive improvement program. Despite the availability of genetic studies of numerous European or Mediterranean germplasm accessions, little is known about the population structure or genetic diversity of the existing USDA collection of olives. An accurate molecular documentation is critical for germplasm curators, breeders, and geneticists as well as plant pathologists. The olive germplasm collections at several centers have been largely influenced by natural dissemination and human migration as well as multilocal selection, breeding, and propagation [29,44,45]. The genetic structure of the USDA germplasm collection consisting of 110 olive cultivars characterized using fifteen microsatellite SSRs markers [29] showed a significant diversity but low levels of differentiation among the olive cultivars within this collection. We chose to use SNPs to distinguish the germplasm collection because of the advent of next-generation sequencing technologies and genome-wide screening capabilities [5]. GBS has become the most popular SNP discovery and genotyping technique in plant species [46,47]. So far, a few studies have used the GBS technique to understand the genetic diversity of local collections of olive cultivars, leading to regional olive improvement programs [10,11,48]. A genome-wide association study (GWAS) successfully used GBS markers to map five agronomic traits using a collection of olive accessions [3]. A recent study involving 57 olive cultivars of European and Mediterranean origins showed that GBS-SNP loci effectively corrected the relationship among different cultivars, further confirming the utility of GBS markers for genetic diversity analyses [5]. Here, we performed a molecular characterization of 96 olive genotypes from the USDA core collection, including six regionally popular varieties using GBS-generated SNP markers. The accessions represented diverse cultivars from Olea europaea L. originating from 18 countries across four climatic zones.

4.2. Features of the SNP Markers

The average number of 3.99 million sequence reads per sample obtained in this study was much higher than other similar studies of olives [10,11]. We obtained 54,075 high-quality SNPs after filtering, which was much higher than previous studies of olives [3,11,49]. On the contrary, these numbers were smaller than a study [5] involving 57 olive cultivars. The results obtained in this study were mainly due to different sources of olive collections and the platform used to resolve the amplified products. However, unlike previous studies, we used an improved reference genome of Olea europaea cv.; “Farga” (version Oe9) [3] was developed by anchoring the previously used Oe6 version [48] to a publicly available genetic map [11].
Transition SNPs were more available than transversions, consistent with the previous studies of olives and other plant species [5,19,49,50,51]. As parameters for describing the variability of the SNP marker, the minor allele frequency, gene diversity (GD), and polymorphic information content (PIC) were used in the study. Understanding the GD and PIC values were critical to finding the polymorphisms among the accessions, selection pressure on the allele, and locus mutation rate over the period [19]. In the current study, the average GD and PIC values were 0.244 and 0.206, respectively. The PIC value of the SNP marker was not equal or close to the GD because of its biallelic nature that restricts an allele frequency increase [52]. Most PIC values were below 50%, consistent with previous studies of olives [5,53]. A total of 33% of the SNPs of this study had a PIC value > 0.25, which was half of its maximum theoretical PIC value (0.5) [54]. Hence, one third of the identified SNPs in this study that showed a 50% resolution ability over its maximum capacity could help in marker-assisted selections or breeding for developing new varieties.

4.3. Features of the US Repository Olive Population

The population structure analysis revealed seven subpopulations (K = 7) with admixture groups in the US olive germplasm collection based on a 0.70 membership probability. This analysis coincided well with the PCoA and genotype-based phylogenetic analysis. The structure analysis outcomes of the study were similar to a previous study [3], which identified six subpopulations in a subset of the US collection of olives. The Bayesian and distance-based clustering of the accessions in this study did not align with the geographical origin and climatic zones, suggesting that the US accessions collected from different countries and climatic zones were genetically similar, possibly because of a limited selection and crossing during early domestication [3]. Although this is in agreement with the geographic origins was consistent with previous GBS-based studies of olives [3,5,55], it was contrary to a report that used a subset of SSR markers to show a partial clustering of the US repository accessions with the geographical origin [56,57]. The difference may be due to a significantly lower number of markers used to make the conclusions described in the SSR study. Nonetheless, a low level of genetic differentiation was consistent with little diversity among the clusters identified by the population structure and phylogenetic analysis in our study. Further, in our study, two major groups with subgroups based on Fst values had a low genetic divergence among the clusters.
The US accessions were grouped into Pop1, 3, 4, 6, and 7 as well as the admixture group. Interestingly, most accessions of these subpopulations were from temperate and dry climatic zones, signifying the restrictive selection of the accessions in these regions. Similarly, the accessions from the non-native areas (Russia, Japan, and Peru) were grouped into the admixture groups with accessions from the Mediterranean regions (e.g., Cyprus, Greece, Italy, Spain, and France). It may be due to a shared ancestry, outcrossing genotypes from diverse backgrounds, or selection pressure to adapt to the local environment.
The Fst values describe the genetic differentiation between two subpopulations [3]. This study observed significant differences among the seven subpopulations except between Pop7 vs. Pop2 and Pop7 vs. Pop3, which was consistent with the phylogenetic analysis and PCoA where the genotypes were grouped based on the population. The pairwise Fst values between most of the population groups were higher than 0.15, suggesting a high genetic differentiation in the US repository collection. These results also supported the AMOVA analysis, which showed that a significant (p ≤ 0.001) genetic variation existed among and within the subpopulations. Although most of the genetic variation (67%) came from within the population, among the subpopulations was the source of a 33% genetic variation. After including the climatic zone as a sub-source of the genetic variation, the climatic zones within the subpopulations accounted for 26 and 4% of the genetic variation among the populations. In both AMOVA analyses, the genetic variation within the subpopulations was the highest and the geographical location also explained a major source of genetic variability. It is believed that the domestication of the olive started 6000 years ago in the eastern Mediterranean basin, from where it migrated to different regions across the world during classical civilization [34,56]. The selection pressure for quality and productivity traits alongside adaptation to a particular climatic area may have contributed towards the variation by geographic locations. The low level of diversity among the subpopulations was consistent with studies of other plants where the among population contributed to the low amount of genetic diversity [19,55]. In terms of allelic patterns and genetic diversity within subpopulations or clusters, Pop5 and Pop7 were more genetically diverse—as reflected by the higher gene and nucleotide diversity values than the other remaining groups, possibly because: (1) the ancestor accessions of these groups were more genetically diverse; (2) there was a selective crossing or natural outcrossings with genotypically diverse cultivars; and (3) there was a selection by the environmental conditions to be more diverse. This information is essential for parent selection in the US olive improvement program to maintain and monitor the genetic diversity required for a successful breeding program or to broaden the genetic basis of the US olive germplasm.

5. Conclusions

This study used high-throughput GBS technology for SNP genotyping to explore genetic diversity and structures in the olive accessions within the US repository. The study identified and examined the features of the SNP markers to facilitate future efforts of US olive genetic improvements. The SNP markers performed well in terms of the polymorphism, genetic diversity, and population structure analysis. A total of 33% of the total SNPs used in the study showed half of their maximum PIC value (0.5), indicating their suitability for future marker-assisted breeding. Overall, our germplasm was genetically diverse with seven subpopulations. This genetic diversity could be helpful for future breeding programs through the selection of suitable parents for developing new olive cultivars with desirable agronomical characteristics adapted to US climatic challenges. The seven subpopulations identified in this study did not align with the geographical origin or climate zones, possibly due to regional selection and domestication. The subpopulations Pop5 and Pop7were genetically diverse whereas Pop6 was less diverse. This information could be helpful to select parents from various groups to widen genetic diversity during olive improvement and breeding programs. The overall findings of this study will help to conduct genetic mapping, association mapping, genomic selection, and marker-assisted breeding.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes12122007/s1: Table S1: Olive accessions with their corresponding USDA plant ID and population groups identified in the structure analysis. Table S2: GBS-generated sequencing reads per sample. Figure S1: Tassel GBS Pipeline Version 2. Figure S2: Nucleotide diversity per site (π) for each chromosome and scaffold.

Author Contributions

Conceptualization, V.J.; methodology, V.J.; software, A.S.M.F.I.; validation, A.S.M.F.I. and D.S.; formal analysis, A.S.M.F.I. and D.S.; investigation, V.J.; resources, A.K.M.; data curation, D.S.; writing—original draft preparation, A.S.M.F.I. and V.J; writing—review and editing, A.S.M.F.I., D.S., A.K.M. and V.J.; visualization, A.S.M.F.I. and D.S.; supervision, V.J.; project administration, V.J.; funding acquisition, V.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Texas Department of Agriculture, Specialty Crop Block Grant SC-1920-13 and by the USDA National Institute of Food and Agriculture (Hatch Project No TEX09647).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data can be made available upon reasonable request.

Acknowledgments

We are thankful to Blake Greene and Madhumita Joshi for their assistance in the DNA extractions, Dalton Thompson for technical support, and Kimberly Cochran for assistance in plant maintenance. The authors thank the Bioinformatics Resource Center, University of Wisconsin–Madison, Madison, for providing the GBS facilities and services. Portions of the data analysis of this research were conducted with the advanced computing resources and consultation provided by Texas A&M High Performance Research Computing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sebastiani, L.; Busconi, M. Recent developments in olive (Olea europaea L.) genetics and genomics: Applications in taxonomy, varietal identification, traceability and breeding. Plant Cell Rep. 2017, 36, 1345–1360. [Google Scholar] [CrossRef] [PubMed]
  2. Gouvinhas, I.; Machado, N.; Sobreira, C.; Domínguez-Perles, R.; Gomes, S.; Rosa, E.; Barros, A.I.R.N.A. Critical Review on the Significance of Olive Phytochemicals in Plant Physiology and Human Health. Molecules 2017, 22, 1986. [Google Scholar] [CrossRef] [Green Version]
  3. Kaya, H.B.; Akdemir, D.; Lozano, R.; Cetin, O.; Kaya, H.S.; Sahin, M.; Smith, J.L.; Tanyolac, B.; Jannink, J.-L. Genome wide association study of 5 agronomic traits in olive (Olea europaea L.). Sci. Rep. 2019, 9, 18764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. López-Miranda, J.; Pérez-Jiménez, F.; Ros, E.; De Caterina, R.; Badimón, L.; Covas, M.; Escrich, E.; Ordovás, J.; Soriguer, F.; Abiá, R.; et al. Olive oil and health: Summary of the II international conference on olive oil and health consensus report, Jaén and Córdoba (Spain) 2008. Nutr. Metab. Cardiovasc. Dis. 2010, 20, 284–294. [Google Scholar] [CrossRef] [PubMed]
  5. Zhu, S.; Niu, E.; Shi, A.; Mou, B. Genetic Diversity Analysis of Olive Germplasm (Olea europaea L.) With Genotyping-by-Sequencing Technology. Front. Genet. 2019, 10, 755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kaniewski, D.; Van Campo, E.; Boiy, T.; Terral, J.-F.; Khadari, B.; Besnard, G. Primary domestication and early uses of the emblematic olive tree: Palaeobotanical, historical and molecular evidence from the Middle East. Biol. Rev. 2012, 87, 885–899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. National Agricultural Statistics Service (NASS); U.S.A. Department of Agriculture. Noncitrus Fruits and Nuts Summary. Available online: https://downloads.usda.library.cornell.edu/usda-esmis/files/zs25x846c/sf269213r/6t054c23t/ncit0521.pdf (accessed on 11 November 2021).
  8. Lanner, R.M.; Taylor, J.M. The Olive in California: History of an Immigrant Tree. West. Hist. Q. 2002, 33, 494. [Google Scholar] [CrossRef]
  9. Lavanya, G.R.; Srivastava, J.; Ranade, S.A. Molecular assessment of genetic diversity in mung bean germplasm. J. Genet. 2008, 87, 65–74. [Google Scholar] [CrossRef] [PubMed]
  10. D’Agostino, N.; Taranto, F.; Camposeo, S.; Mangini, G.; Fanelli, V.; Gadaleta, S.; Miazzi, M.M.; Pavan, S.; Di Rienzo, V.; Sabetta, W.; et al. GBS-derived SNP catalogue unveiled wide genetic variability and geographical relationships of Italian olive cultivars. Sci. Rep. 2018, 8, 15877. [Google Scholar] [CrossRef] [Green Version]
  11. Ipek, A.; Yılmaz, K.; Sıkıcı, P.; Tangu, N.A.; Öz, A.T.; Bayraktar, M.; Ipek, M.; Gülen, H. SNP Discovery by GBS in Olive and the Construction of a High-Density Genetic Linkage Map. Biochem. Genet. 2016, 54, 313–325. [Google Scholar] [CrossRef] [PubMed]
  12. Sion, S.; Savoia, M.A.; Gadaleta, S.; Piarulli, L.; Mascio, I.; Fanelli, V.; Montemurro, C.; Miazzi, M.M. How to Choose a Good Marker to Analyze the Olive Germplasm (Olea europaea L.) and Derived Products. Genes 2021, 12, 1474. [Google Scholar] [CrossRef]
  13. Pontikis, C.A.; Loukas, M.; Kousounis, G. The use of Biochemical Markers to Distinguish Olive Cultivars. J. Hortic. Sci. 1980, 55, 333–343. [Google Scholar] [CrossRef]
  14. Decroocq, S.; Cornille, A.; Tricon, D.; Babayeva, S.; Chague, A.; Eyquard, J.-P.; Karychev, R.; Dolgikh, S.; Kostritsyna, T.; Liu, S.; et al. Data from: New insights into the history of domesticated and wild apricots and its contribution to Plum pox virus resistance. Mol. Ecol. 2016, 25, 4712–4729. [Google Scholar] [CrossRef] [PubMed]
  15. Delplancke, M.; Alvarez, N.; Benoit, L.; Espindola, A.; Joly, H.; Neuenschwander, S.; Arrigo, N. Evolutionary history of almond tree domestication in the M editerranean basin. Mol. Ecol. 2013, 22, 1092–1104. [Google Scholar] [CrossRef] [PubMed]
  16. Richards, C.M.; Volk, G.M.; Reilley, A.A.; Henk, A.D.; Lockwood, D.R.; Reeves, P.A.; Forsline, P.L. Genetic diversity and population structure in Malus sieversii, a wild progenitor species of domesticated apple. Tree Genet. Genomes 2009, 5, 339–347. [Google Scholar] [CrossRef]
  17. Kumar, S.; Garrick, D.J.; Bink, M.C.; Whitworth, C.; Chagné, D.; Volz, R.K. Novel genomic approaches unravel genetic architecture of complex traits in apple. BMC Genom. 2013, 14, 393. [Google Scholar] [CrossRef] [Green Version]
  18. Kouassi, A.B.; Durel, C.-E.; Costa, F.; Tartarini, S.; van de Weg, E.; Evans, K.; Fernandez-Fernandez, F.; Govan, C.; Boudichevskaja, A.; Dunemann, F.; et al. Estimation of genetic parameters and prediction of breeding values for apple fruit-quality traits using pedigreed plant material in Europe. Tree Genet. Genomes 2009, 5, 659–672. [Google Scholar] [CrossRef]
  19. Luo, Z.; Brock, J.; Dyer, J.M.; Kutchan, T.; Schachtman, D.; Augustin, M.; Ge, Y.; Fahlgren, N.; Abdel-Haleem, H. Genetic Diversity and Population Structure of a Camelina sativa Spring Panel. Front. Plant Sci. 2019, 10, 184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Allen, A.M.; Winfield, M.O.; Burridge, A.J.; Downie, R.C.; Benbow, H.R.; Barker, G.L.A.; Wilkinson, P.A.; Coghill, J.; Waterfall, C.; Davassi, A.; et al. Characterization of a Wheat Breeders’ Array suitable for high-throughput SNP genotyping of global accessions of hexaploid bread wheat (Triticum aestivum). Plant Biotechnol. J. 2017, 15, 390–401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Poland, J.A.; Brown, P.J.; Sorrells, M.E.; Jannink, J.-L. Development of High-Density Genetic Maps for Barley and Wheat Using a Novel Two-Enzyme Genotyping-by-Sequencing Approach. PLoS ONE 2012, 7, e32253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Wang, N.; Yuan, Y.; Wang, H.; Yu, D.; Liu, Y.; Zhang, A.; Gowda, M.; Nair, S.K.; Hao, Z.; Lu, Y.; et al. Applications of genotyping-by-sequencing (GBS) in maize genetics and breeding. Sci. Rep. 2020, 10, 16308. [Google Scholar] [CrossRef] [PubMed]
  23. Kitony, J.; Sunohara, H.; Tasaki, M.; Mori, J.-I.; Shimazu, A.; Reyes, V.; Yasui, H.; Yamagata, Y.; Yoshimura, A.; Yamasaki, M.; et al. Development of an Aus-Derived Nested Association Mapping (Aus-NAM) Population in Rice. Plants 2021, 10, 1255. [Google Scholar] [CrossRef]
  24. McNally, K.L.; Childs, K.L.; Bohnert, R.; Davidson, R.M.; Zhao, K.; Ulat, V.J.; Zeller, G.; Clark, R.M.; Hoen, D.R.; Bureau, T.E.; et al. Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. Proc. Natl. Acad. Sci. USA 2009, 106, 12273–12278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Fujii, H.; Shimada, T.; Nonaka, K.; Kita, M.; Kuniga, T.; Endo, T.; Ikoma, Y.; Omura, M. High-throughput genotyping in citrus accessions using an SNP genotyping array. Tree Genet. Genomes 2013, 9, 145–153. [Google Scholar] [CrossRef]
  26. Micheletti, D.; Dettori, M.T.; Micali, S.; Aramini, V.; Pacheco, I.; Linge, C.D.S.; Foschi, S.; Banchi, E.; Barreneche, T.; Quilot-Turion, B.; et al. Whole-Genome Analysis of Diversity and SNP-Major Gene Association in Peach Germplasm. PLoS ONE 2015, 10, e0136803. [Google Scholar] [CrossRef]
  27. Larsen, B.; Gardner, K.; Pedersen, C.; Ørgaard, M.; Migicovsky, Z.; Myles, S.; Toldam-Andersen, T. Population structure, relatedness and ploidy levels in an apple gene bank revealed through genotyping-by-sequencing. PLoS ONE 2018, 13, e0201889. [Google Scholar] [CrossRef] [Green Version]
  28. Xia, W.; Luo, T.; Zhang, W.; Mason, A.; Huang, D.; Huang, X.; Tang, W.; Dou, Y.; Zhang, C.; Xiao, Y. Development of High-Density SNP Markers and Their Application in Evaluating Genetic Diversity and Population Structure in Elaeis guineensis. Front. Plant Sci. 2019, 10, 130. [Google Scholar] [CrossRef]
  29. Koehmstedt, A.M.; Aradhya, M.K.; Soleri, D.; Smith, J.L.; Polito, V.S. Molecular characterization of genetic diversity, structure, and differentiation in the olive (Olea europaea L.) germplasm collection of the United States Department of Agriculture. Genet. Resour. Crop. Evol. 2010, 58, 519–531. [Google Scholar] [CrossRef] [Green Version]
  30. Jombart, T.; Kamvar, Z.N.; Collins, C.; Lustrik, R.; Beugin, M.-P.; Knaus, B.J.; Jombart, M.T. Package ‘Adegenet’. Github Repository. Available online: https://github.com/thibautjombart/adegenet (accessed on 11 November 2021).
  31. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World Map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef]
  32. Jiang, H.; Lei, R.; Ding, S.-W.; Zhu, S. Skewer: A fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinform. 2014, 15, 182. [Google Scholar] [CrossRef] [PubMed]
  33. Glaubitz, J.C.; Casstevens, T.M.; Lu, F.; Harriman, J.; Elshire, R.J.; Sun, Q.; Buckler, E.S. TASSEL-GBS: A High Capacity Genotyping by Sequencing Analysis Pipeline. PLoS ONE 2014, 9, e90346. [Google Scholar] [CrossRef]
  34. Julca, I.; Marcet-Houben, M.; Cruz, F.; Gómez-Garrido, J.; Gaut, B.S.; Díez, C.M.; Gut, I.G.; Alioto, T.S.; Vargas, P.; Gabaldón, T. Genomic evidence for recurrent genetic admixture during the domestication of Mediterranean olive trees (Olea europaea L.). BMC Biol. 2020, 18, 148. [Google Scholar] [CrossRef] [PubMed]
  35. Nardi, E.P.; Evangelista, F.S.; Tormen, L.; Saint’pierre, T.D.; Curtius, A.J.; de Souza, S.S.; Barbosa, F. The use of inductively coupled plasma mass spectrometry (ICP-MS) for the determination of toxic and essential elements in different types of food samples. Food Chem. 2009, 112, 727–732. [Google Scholar] [CrossRef]
  36. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef] [PubMed]
  37. Granato, I.S.C.; Galli, G.; Couto, E.G.O.; Souza, M.B.; Mendonça, L.F.; Fritsche-Neto, R. snpReady: A tool to assist breeders in genomic analysis. Mol. Breed. 2018, 38, 102. [Google Scholar] [CrossRef]
  38. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. In Genetics; 2000; 155, pp. 945–959. Available online: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/10835412 (accessed on 11 November 2021). [CrossRef] [PubMed]
  39. Kopelman, N.M.; Mayzel, J.; Jakobsson, M.; Rosenberg, N.A.; Mayrose, I. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 2015, 15, 1179–1191. [Google Scholar] [CrossRef] [Green Version]
  40. 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] [Green Version]
  41. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kamvar, Z.N.; Tabima, J.F.; Grünwald, N.J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2014, 2, e281. [Google Scholar] [CrossRef] [Green Version]
  43. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef]
  44. Belaj, A.; Satovic, Z.; Rallo, L.; Trujillo, I. Genetic diversity and relationships in olive (Olea europaea L.) germplasm collections as determined by randomly amplified polymorphic DNA. Theor. Appl. Genet. 2002, 105, 638–644. [Google Scholar] [CrossRef]
  45. Besnard, G.; Baradat, P.; Bervillé, A. Genetic relationships in the olive (Olea europaea L.) reflect multilocal selection of cultivars. Theor. Appl. Genet. 2001, 102, 251–258. [Google Scholar] [CrossRef]
  46. He, J.; Zhao, X.; Laroche, A.; Lu, Z.-X.; Liu, H.K.; Li, Z. Genotyping-by-sequencing (GBS), an ultimate marker-assisted selection (MAS) tool to accelerate plant breeding. Front. Plant Sci. 2014, 5, 484. [Google Scholar] [CrossRef] [Green Version]
  47. Scheben, A.; Batley, J.; Edwards, D. Genotyping-by-sequencing approaches to characterize crop genomes: Choosing the right tool for the right application. Plant Biotechnol. J. 2017, 15, 149–161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Cruz, F.; Julca, I.; Gómez-Garrido, J.; Loska, D.; Marcet-Houben, M.; Cano, E.; Galán, B.; Frias, L.; Ribeca, P.; Derdak, S.; et al. Genome sequence of the olive tree, Olea europaea. GigaScience 2016, 5, 29. [Google Scholar] [CrossRef]
  49. Mantello, C.C.; Cardoso-Silva, C.B.; Da Silva, C.C.; De Souza, L.M.; Junior, E.J.S.; Gonçalves, P.D.S.; Vicentini, R.; De Souza, A.P. De Novo Assembly and Transcriptome Analysis of the Rubber Tree (Hevea brasiliensis) and SNP Markers Development for Rubber Biosynthesis Pathways. PLoS ONE 2014, 9, e102665. [Google Scholar] [CrossRef] [Green Version]
  50. Clarke, W.E.; Parkin, I.A.; Gajardo, H.A.; Gerhardt, D.J.; Higgins, E.; Sidebottom, C.; Sharpe, A.G.; Snowdon, R.J.; Federico, M.L.; Iniguez-Luy, F.L. Genomic DNA enrichment using sequence capture microarrays: A novel approach to discover sequence nucleotide polymorphisms (SNP) in Brassica napus L. PLoS ONE 2013, 8, e81992. [Google Scholar]
  51. Morton, B.R.; Bi, I.V.; McMullen, M.D.; Gaut, B.S. Variation in Mutation Dynamics Across the Maize Genome as a Function of Regional and Flanking Base Composition. Genetics 2006, 172, 569–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Shete, S.; Tiwari, H.; Elston, R.C. On Estimating the Heterozygosity and Polymorphism Information Content Value. Theor. Popul. Biol. 2000, 57, 265–271. [Google Scholar] [CrossRef] [PubMed]
  53. Belaj, A.; De La Rosa, R.; Lorite, I.J.; Mariotti, R.; Cultrera, N.G.; Beuzón, C.R.; González-Plaza, J.J.; Muñoz-Merida, A.; Trelles, O.; Baldoni, L. Usefulness of a New Large Set of High Throughput EST-SNP Markers as a Tool for Olive Germplasm Collection Management. Front. Plant Sci. 2018, 9, 1320. [Google Scholar] [CrossRef] [Green Version]
  54. Anderson, J.A.; Churchill, G.A.; Autroque, J.E.; Tanksley, S.D.; Swells, M.E. Optimising selection for plant linkage map. Genome 1993, 36, 181–186. [Google Scholar] [CrossRef] [PubMed]
  55. Eltaher, S.; Sallam, A.; Belamkar, V.; Emara, H.A.; Nower, A.; Salem, K.; Poland, J.; Baenziger, P.S. Genetic Diversity and Population Structure of F3:6 Nebraska Winter Wheat Genotypes Using Genotyping-By-Sequencing. Front. Genet. 2018, 9, 76. [Google Scholar] [CrossRef] [PubMed]
  56. Belaj, A.; del Carmen Dominguez-García, M.; Atienza, S.G.; Urdíroz, N.M.; De la Rosa, R.; Satovic, Z.; Martín, A.; Kilian, A.; Trujillo, I.; Valpuesta, V.; et al. Developing a core collection of olive (Olea europaea L.) based on molecular markers (DArTs, SSRs, SNPs) and agronomic traits. Tree Genet. Genomes 2012, 8, 365–378. [Google Scholar] [CrossRef]
  57. El Bakkali, A.; Haouane, H.; Moukhli, A.; Costes, E.; Van Damme, P.; Khadari, B. Construction of Core Collections Suitable for Association Mapping to Optimize Use of Mediterranean Olive (Olea europaea L.) Genetic Resources. PLoS ONE 2013, 8, e61265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 12 02007 g001 550
Figure 1. GBS-generated SNP marker characterization. (a) Number of SNPs per chromosome; (b) minor allele frequency (MAF); (c) gene diversity (GD); (d) polymorphic information content (PIC).
Figure 1. GBS-generated SNP marker characterization. (a) Number of SNPs per chromosome; (b) minor allele frequency (MAF); (c) gene diversity (GD); (d) polymorphic information content (PIC).
Genes 12 02007 g001
Genes 12 02007 g002 550
Figure 2. Population structure for 96 olive accessions from the USDA core collection using 54,075 SNP markers. (a) Subpopulation grouping inferred by the STRUCTURE software indicated in seven different colors. The y-axis values indicate the probability of the population. (b) Evanno test for the ideal population number using LnP(D)-derived ΔK from 2 to 10. At a value of K = 7, it reaches its highest value, indicating the most probable subpopulations in the germplasm. (c) Olive accession number based on the populations identified by the STRUCTURE software.
Figure 2. Population structure for 96 olive accessions from the USDA core collection using 54,075 SNP markers. (a) Subpopulation grouping inferred by the STRUCTURE software indicated in seven different colors. The y-axis values indicate the probability of the population. (b) Evanno test for the ideal population number using LnP(D)-derived ΔK from 2 to 10. At a value of K = 7, it reaches its highest value, indicating the most probable subpopulations in the germplasm. (c) Olive accession number based on the populations identified by the STRUCTURE software.
Genes 12 02007 g002
Genes 12 02007 g003 550
Figure 3. Principal coordinate analysis (PCoA). The analysis depicts the overall genetic diversity of the germplasm for 96 olive accessions using 54,075 SNP markers.
Figure 3. Principal coordinate analysis (PCoA). The analysis depicts the overall genetic diversity of the germplasm for 96 olive accessions using 54,075 SNP markers.
Genes 12 02007 g003
Genes 12 02007 g004 550
Figure 4. The genetic diversity analysis. The genetic diversity based on (a) the nucleotide diversity per site (π), (b) the fixation index (Fst), and (c) the matrix showing the pairwise Fst values between the populations.
Figure 4. The genetic diversity analysis. The genetic diversity based on (a) the nucleotide diversity per site (π), (b) the fixation index (Fst), and (c) the matrix showing the pairwise Fst values between the populations.
Genes 12 02007 g004
Genes 12 02007 g005 550
Figure 5. Phylogenetic analyses of the 96 olive cultivars using the neighbor-joining method. Different colors depict the structure analysis generated populations. Legends indicate the climatic zones from where the accessions were originated. Colors represent different subpopulations of the germplasm; red color = subpopulation1, green = subpopulation2, blue = subpopulation3, yellow = subpopulation4, pink = subpopulation5, cyan = subpopulation6, maroon = subpopulation7, black = admixture group.
Figure 5. Phylogenetic analyses of the 96 olive cultivars using the neighbor-joining method. Different colors depict the structure analysis generated populations. Legends indicate the climatic zones from where the accessions were originated. Colors represent different subpopulations of the germplasm; red color = subpopulation1, green = subpopulation2, blue = subpopulation3, yellow = subpopulation4, pink = subpopulation5, cyan = subpopulation6, maroon = subpopulation7, black = admixture group.
Genes 12 02007 g005
Genes 12 02007 g006 550
Figure 6. Population-wide marker characteristics. The marker characteristics in terms of the gene diversity (GD), polymorphic information content (PIC), and minor allele frequency (MAF) across the seven subpopulations and the whole population.
Figure 6. Population-wide marker characteristics. The marker characteristics in terms of the gene diversity (GD), polymorphic information content (PIC), and minor allele frequency (MAF) across the seven subpopulations and the whole population.
Genes 12 02007 g006
Table
Table 1. List of the countries of origin of the 96 olive accessions with their corresponding climatic zone.
Table 1. List of the countries of origin of the 96 olive accessions with their corresponding climatic zone.
OriginNo. of Genotypes
Climatic zoneCountry
ContinentalRussia2
DryAlgeria1
Argentina1
Chile1
Egypt3
Israel1
Morocco1
Pakistan1
Palestine1
Syria1
Tunisia5
TemperateAlbania2
Cyprus3
France5
Greece5
Italy22
Japan1
Spain14
US17
TropicalColombia1
Peru2
Unknown 6
Table
Table 2. Summary result of the SNP types.
Table 2. Summary result of the SNP types.
TransitionsTransversions
Type of SNPA/GC/TA/CA/TG/TC/G
Number of sites15,87516,0105513747652883913
Frequencies29.36%29.61%10.20%13.83%9.78%7.24%
Total31,885 (59.96%)22,190 (41.04%)
Table
Table 3. Analysis of molecular variance (AMOVA) based on the population as being the source of the genetic variation.
Table 3. Analysis of molecular variance (AMOVA) based on the population as being the source of the genetic variation.
Source of VariationdfSSMSVariationp-Value
Between the population7630,346.0890,049.4432.940.001
Within the population881,277,390.8414,515.8067.060.001
Total951,907,736.9120,081.44100.00
Table
Table 4. Analysis of molecular variance (AMOVA) based on the population and climatic zone as being the sources of the genetic variation.
Table 4. Analysis of molecular variance (AMOVA) based on the population and climatic zone as being the sources of the genetic variation.
Source of VariationdfSSMSVariationp-Value
Between the population7310,271.0444,324.434.290.04
Between the climatic zones within the population13493,233.7437,941.0625.820.001
Within the population751,104,232.1314,723.1069.890.001
Total951,907,736.9120,081.44100.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Islam, A.S.M.F.; Sanders, D.; Mishra, A.K.; Joshi, V. Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS). Genes 2021, 12, 2007. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122007

AMA Style

Islam ASMF, Sanders D, Mishra AK, Joshi V. Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS). Genes. 2021; 12(12):2007. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122007

Chicago/Turabian Style

Islam, A. S. M. Faridul, Dean Sanders, Amit Kumar Mishra, and Vijay Joshi. 2021. "Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS)" Genes 12, no. 12: 2007. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122007

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