Next Article in Journal
Using Micro- and Macro-Level Network Metrics Unveils Top Communicative Gene Modules in Psoriasis
Next Article in Special Issue
Genetic Resources of Olea europaea L. in the Garda Trentino Olive Groves Revealed by Ancient Trees Genotyping and Parentage Analysis of Drupe Embryos
Previous Article in Journal
The Evolution of Gene Therapy in the Treatment of Metabolic Liver Diseases
Previous Article in Special Issue
Chloroplast Genome Variation and Evolutionary Analysis of Olea europaea L.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 8
10.3390/genes11080916
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

EST–SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives

by
Roberto Mariotti
1,†,
Angjelina Belaj
2,†,
Raul De La Rosa
2,
Lorenzo Leòn
2,
Federico Brizioli
1,
Luciana Baldoni
1,* and
Soraya Mousavi
1
1
CNR—Institute of Biosciences and Bioresources, Via Madonna Alta 130, 06128 Perugia, Italy
2
IFAPA—Centro Alameda del Obispo, Avda Menendez Pidal, s/n, E-14004 Cordoba, Spain
*
Author to whom correspondence should be addressed.
†
These authors have equally contributed to this work.
Genes 2020, 11(8), 916; https://doi.org/10.3390/genes11080916
Submission received: 14 July 2020 / Revised: 4 August 2020 / Accepted: 7 August 2020 / Published: 10 August 2020
(This article belongs to the Special Issue Oleaceae Genetics)
Browse Figures
Versions Notes

Abstract

:
Background: The species Olea europaea includes cultivated varieties (subsp. europaea var. europaea), wild plants (subsp. europaea var. sylvestris), and five other subspecies spread over almost all continents. Single nucleotide polymorphisms in the expressed sequence tag able to underline intra-species differentiation are not yet identified, beyond a few plastidial markers. Methods: In the present work, more than 1000 transcript-specific SNP markers obtained by the genotyping of 260 individuals were studied. These genotypes included cultivated, oleasters, and samples of subspecies guanchica, and were analyzed in silico, in order to identify polymorphisms on key genes distinguishing different Olea europaea forms. Results: Phylogeny inference and principal coordinate analysis allowed to detect two distinct clusters, clearly separating wilds and guanchica samples from cultivated olives, meanwhile the structure analysis made possible to differentiate these three groups. Sequences carrying the polymorphisms that distinguished wild and cultivated olives were analyzed and annotated, allowing to identify 124 candidate genes that have a functional role in flower development, stress response, or involvement in important metabolic pathways. Signatures of selection that occurred during olive domestication, were detected and reported. Conclusion: This deep EST–SNP analysis provided important information on the genetic and genomic diversity of the olive complex, opening new opportunities to detect gene polymorphisms with potential functional and evolutionary roles, and to apply them in genomics-assisted breeding, highlighting the importance of olive germplasm conservation.

1. Introduction

Olive (Olea europaea L.) is a long-living widespread tree species. Its genetic patrimony includes thousands of cultivated varieties, wild olives, and related subspecies, such as cuspidata, laperrinei, cerasiformis, maroccana, and guanchica, which grow in natural areas of different continents [1,2,3,4,5]. The cultivated and wild (also named oleaster) forms of olive (subsp. europaea, var. europaea and sylvestris, respectively), typical of the Mediterranean area, are diploid, with 23 chromosomes and a genome size of about 1.48 Giga base pairs [6]. Subspecies guanchica is a diploid subspecies spreading only in the Canary Islands [7].
The high genetic variability of cultivated olive [8,9,10,11,12] (mainly related to its allogamy), the lack of turnover with modern cultivars, a remarkable tree longevity, and the capacity to survive out of managed cultivation, allowed it to spread over a wide range of environments [13,14,15,16]. Inside this diversified scenario, the debate on olive domestication and the contribution of oleasters and related subspecies to this process, is still open and controversial [17,18,19,20,21,22]. Some recent analyses suggest that the olive tree domestication occurred in the Levant (Syrian–Turkish) area, followed by secondary diversifications in two opposite directions, to the west side across the shores of the Mediterranean [13,23,24,25] and up to the east in the Iranian plateau [16,26,27]. Recently it was reported that the gene flow between domesticated olives and the wild relatives allowed secondary diversification and contributed to the emergence of thousands of local varieties [15,18].
During plant domestication, wild species were converted into crop plants through artificial selection, with conscious and unconscious contribution [28]. Artificial selection might influence important loci variation, affecting traits of agronomic interest. Artificially selected SNPs might be consistently linked to important agronomic traits [29]. Artificial selection of crop varieties focused on modifying factors such as plant architecture and stress resistance, as a means of increasing yield [30]. The selection process in cereal crops left observable genetic changes, including a reduction of genetic diversity and an increased frequency of favorable seed size alleles in the cultivated lines, as compared to their wild progenitors [31,32,33]. The artificial selection could also be the cause of rapid demographic expansions, as reported for the cotton cultivation [34].
Domestication of perennial plants is expected to differ from that of annual plants, as a consequence of the evidence that trees are usually propagated vegetatively and have an extended juvenile stage, thus domestication bottleneck cannot be as severe as that for annual plants [35]. Considering what was reported previously, perennial crops retain, 95% of neutral variation from their progenitors, on average, while the annuals only retain 60% [36]. There are only few works related to perennial fruit trees. In grape, differentiation between wild and cultivated plants caused morphological changes on berry size, sugar content, seed morphology, and a shift from dioecy to hermaphroditic mating system [35,37]. Recently, a transcriptomic analysis on wild and cultivated olives indicated that the domestication of olives had moderate genomic consequences and was more effective on gene expression [38]. The studies on selection signatures might provide a “bottom-up” approach [39] on the genetic basis that differentiate the present cultivated plants from wild trees [22]. In fact, cultivated olives were introduced from the eastern side of the Mediterranean basin, as seedlings, grafting, and cuttings, strongly reducing the genetic diversity carried by autochthonous olives, such as the central-western oleasters [23].
Sequences related to expressed tags such as the long core repeat of EST–SSRs [40,41,42], developed from transcript collections, might contribute to start the study of gene differences between olive subspecies, highlighting the possible functional effect of polymorphisms between olive genotypes [42]. In the last years, special attention was paid to the development and use of SNP markers in olive [43,44,45,46,47] and, thanks to a wide set of EST–SNPs recently developed [48], the study of functional polymorphisms in olive is finally having a boost. The gene flow between oleasters, subspecies guanchica and Western Mediterranean and some cultivated varieties was already reported by using a different set of markers [3,49,50]. Therefore, based on the evidence of the previous studies, this work aimed to identify potential functional markers for differentiating cultivated from wild olive trees, by addressing the following questions—(a) what genes distinguish cultivars from oleasters and guanchica olives? (b) what are the functions of differently evolved genes? (c) which of these gene polymorphisms could be introgressed into cultivated varieties through genomic-assisted breeding, to improve their performance?

2. Materials and Methods

2.1. Plant Material and EST—SNP Genotyping

Data of 156 olive cultivars, 73 wild olives and 16 guanchica samples from the World Olive Germplasm Collection of IFAPA (Córdoba, Spain), were collected from a previous published study [47], while 15 cultivars from the International Olive Germplasm Collection of Lugnano in Teverina (Terni, Italy), were de-novo genotyped using the same approach [47]. The final selection of common markers included 1040 EST—SNPs for 260 genotypes (Table S1).

2.2. Genetic Differentiation and Population Genetic Structure

General relationships within populations and amongst populations were described through a principal coordinate analysis (PCoA) performed by GenAlEx 6.5, based on a covariance matrix created from pairwise genetic distances [51], while the three-dimensional visualization was made by the online tool Cubemaker (https://tools.altiusinstitute.org/cubemaker/). Phylogenetic analyses were performed using a parametric method based on sequence alignment data. All SNP allele’s calls for the 260 genotypes were concatenated to give a specific sequence for each accession, using the IUPAC code. Multiple sequence alignment was performed using ClustalW in the Bioedit software [52], and the files were converted into MEGA format using MEGAX [53]. A phylogenetic tree was constructed with MEGAX, with a bootstrap inferred from 1000 replicates. The optimal tree had a sum of branch length equal to 6.34. The tree was drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method and were in the units of the number of base substitutions per site. All positions with less than 95% site coverage were eliminated. Alignment gaps lower than 5%, missing data, and ambiguous bases were allowed at any position. This analysis involved 260 nucleotide sequences and 1040 positions in the final dataset. The phylogenetic tree was drawn and visualized using FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/).
The genotypic data were analyzed using the Bayesian clustering algorithm of STRUCTURE v2.3 [54,55], using the admixture model and correlated allele frequencies with a burn-in of 100,000 iterations, run length of 100,000, and K = 2 to 10. The optimal value of K was determined using the delta K procedure, using STRUCTURE Harvester [56].

2.3. Genetic Diversity

The genetic diversity was calculated by using DnaSP 6 [57]. Polymorphic sites and DNA polymorphisms were detected for the overall sample set, while DNA diversity, polymorphism and divergence, gene flow, genetic differentiation, and population size changes were done within and between the three selected populations. Nucleotide diversity based on the average pairwise difference among sequences (Pi) and the average number of nucleotide differences (k) were also calculated by DnaSP 6. In order to detect if the three selected populations evolved randomly or the mutations were under selection pressure, Tajima’s D test was also performed. The obtained diversity indices were strictly related to the EST–SNPs described previously.
Heterozygosity, F statistics, Polymorphism by Population, Allele Frequencies by Population (AFP), and Allele Frequencies by Population for the Haploid Data (HAFP), were performed by GenAlEx 6.5.

2.4. Identification and Functional Analysis of EST–SNP

The set of SNPs, derived from the expressed sequence tags, and each nucleotide change could correspond to a different level of expression or protein structural variation, both for synonymous (sSNP) and non-synonymous (nsSNP) SNPs. All EST–SNPs were compared between the three olive groups and a cut-off of 45% of frequencies between two homozygous states was selected (e.g., SNP1 frequency: AAcultivated minus AAwild ≥45% yes, <45% no). All transcripts showing two different homozygosity states for the cultivated group versus oleasters and subsp. guanchica, were classified. The sequences of each contig were blasted into the NCBI database, using the BLASTN algorithm of the NCBI-blast procedure (https://blast.ncbi.nlm.nih.gov/), in order to identify nucleotide polymorphisms on the gene sequences. For the uncharacterized genes, the predicted amino acid sequences were used as a query in the same website, using BLASTP. In order to annotate the contigs that were differently allocated among the three groups of genotypes, they were used as query to blast in the NCBI database. Nucleotide changes were investigated using the ExPASy translate tool (http://web.expasy.org/translate/) for both alleles, in order to identify potential amino acid variations. Finally, for the functional classification according to Gene Ontology (GO) terms, the identified genes were used as a query for a statistical overrepresentation test, using the PANTHER online tool [58], which gave the Arabidopsis thaliana reference list with FISHER test type (http://pantherdb.org/tools/compareToRefList.jsp). In order to categorize the biological processes, the cellular component and the molecular function of each gene, the search was also performed on TAIR (https://www.arabidopsis.org/) and UniProtKB (https://www.uniprot.org/uniprot) web tools, including all plant species.

3. Results

3.1. Genetic Differentiation and Population Structure

Based on the results of the genetic differentiation and population structure analyses, the genetic diversity analyses were performed by dividing the sample set in three populations (cultivars, oleasters, and guanchica samples). The PCoA analysis (Figure 1A), based on a genetic distance matrix between the 260 olive genotypes, the three coordinates (1, 2, and 3) corresponding to the x, y, and z axes, represented the 67.57% and the 10.32% and 7.55% of total variability, respectively. It showed a clear separation between the cultivars and the group of oleaster and subsp. guanchica, with only a few mixed samples between the wild and the cultivated groups. In the first case, a supposed oleaster more likely represented a feral form (seedling from a cultivated olive), meanwhile in the latter case, the cultivated olive showed a strong wild genetic background. In fact, an interesting overlap in the crossing of axes among the cultivars and the oleasters was observed, evidencing a clear introgression of gene polymorphisms from var. sylvestris into var. europaea. The resultant guanchica samples were intermixed with the oleasters and quite distant from the cultivated genotypes.
The Unrooted Neighbor Joining tree clearly confirmed the differentiation between cultivars and oleasters (Figure 1B). This elaboration showed some guanchica plants together with oleaster, even if, almost two separated sub-clusters could be observed. Seven olive cultivars clustered within the var. sylvestris’s group, while only one oleaster, putatively a cultivar seedling, was placed within the cultivar’s group after both analyses.
When the Bayesian approach was applied, the most probable K was detected in 5 clusters (Figure 2)—the cultivated varieties were split into three sub-populations, represented by clusters 1, 2, 3, respectively, including cultivars from the West, Center, and East Mediterranean; cluster 4 included 64 out of 73 oleasters and only three cultivated varieties; and cluster 4 corresponded to the 16 guanchica samples (Table S1), together with the nine remaining oleasters (Figure 2). Allele-frequency divergence among pops (net nucleotide distance), computed using the point estimates of P, showed the highest distances between the oleasters and the cultivated varieties from Central and Eastern Mediterranean (Table S2). The subsp. guanchica group displayed similar genetic distance values as the oleasters with the cultivated varieties groups. The lowest value was observed between the oleasters and subsp. guanchica samples, followed by the distance between Western and Eastern Mediterranean cultivars. An unexpected high level of admixture was observed in most cases, especially between the cultivated varieties and the oleasters.
Based on these results, which clearly divided the three sources of materials even with few exceptions, the genetic diversity analyses able to individuate possible signs of genes differentiation were performed, which divided the sample set into three groups (cultivar, oleaster, and guanchica).

3.2. Genetic Diversity of the Cultivated Varieties, Oleasters, and Subsp. Guanchica Samples

The 1040 EST–SNPs detected on 260 samples were used to evaluate the genetic diversity overall, inside and between each group. The variable polymorphic sites (S) individuated by DnaSP, in the entire sample set, were 637. The standard deviations of heterozygosity, based on the number of variable sites among the sequences (Theta), with and without recombination, were >0.01 and >0.05, respectively. The level of total variance, with or without free recombination, was always lower than 0.01. Nucleotide diversity based on the average pairwise difference among sequences (Pi) was equal to 0.35, while the average number of nucleotide differences (k) was equal to 226.04. When the analysis was applied on each of the three selected groups, subsp. guanchica, oleaster, and cultivar, k was equal to 207.21 and Pi was 0.33 for the olive cultivars, while k = 177.67 and 137.95 and Pi equal to 0.28 and 0.22 was detected for the oleasters and subsp. guanchica, respectively (Table S3). The Pi results for each group indicated that the cultivated olive maintained a higher level of nucleotide diversity, compared to wild olives and subsp. guanchica. When the cultivated samples were matched against the oleaster group, the average number of nucleotide substitution per site (Dxy): 0.405 and number of net nucleotide substitution per site between the cultivated and wild olives (Da): 0.102. Moreover, analyzing the cultivated samples against subsp. guanchica, Dxy: 0.400 and Da: 0.129. Finally, for oleasters versus subsp. guanchica, Dxy: 0.273 and Da: 0.025 (Table 1). The calculated parameters indicated exactly what was previously reported from the above results. In fact, the main differences were between the cultivars and the other two groups, while a higher affinity was found between the oleasters and the subsp. guanchica.
The frequency analysis highlighted higher values for the observed heterozygosity than the expected ones for all three analyzed groups, but especially for the olive cultivars (+12%). Moreover, the Fixation Index F, also called the Inbreeding Coefficient, exhibited the lowest and negative value again for the cultivated olive, indicating an excess of heterozygosity, due to negative assortative mating, or a selection for heterozygotes (Table 2).
Lastly, the Tajima’s D test, applied to determine the allele frequency change, was compared between the wild, cultivated, and subsp. guanchica olives. In the cultivated olives, the Tajima’s D values were higher than those in the oleaster and subsp. guanchica samples. The Tajima’s value (0.556, p > 0.1) indicated that subsp. guanchica had a similar variation as expected, which is usually observed in populations growing in equilibrium, without selection evidence. The highest value was detected for the olive cultivars (3.244, p < 0.01), two times higher than what was observed for the oleasters. This result might suggest a lack of rare alleles or that the cultivated olives experienced a balancing selection situation or a sudden population contraction.

3.3. Identification of Genes Carrying Polymorphisms between Cultivated and Wild Olives

The allele frequencies by population for haploid data (HAFP) analysis was used by merging the two SNP bases in a single data, in order to identify genes with an opposite homozygous state, which were strictly related to each group (Table S4; Figure S1). A total of 124 genes with more than 45% differences between frequencies in each group and an opposite homozygous state were selected. Seventy-two were found between the cultivated and the wild and subsp. guanchica groups, 24 between the cultivar and only the oleaster, and 28 between the cultivar and the guanchica group.
Several genes were previously annotated in olive, meanwhile, for the genes without any correspondence in olive and for the unknown ones, the amino acid sequences were manually annotated by a BLASTP query, in order to find protein similarities in other plant species. Forty-four genes were not placed in any linkage group of the oleaster genetic map. From the total of 80 localized genes, the highest number of genes in each chromosome were eight genes in chromosome 1 and chromosome 10, seven in chromosome 18 and six in chromosome 12 (Table S5). From all 23 olive chromosomes, none of our detected genes was located in chromosomes 5, 9, and 19. Finally, by using the Expasy translate tool, alleles with different single point mutations were translated into the amino acid sequences and the possible amino acid changes were detected for each allele at each gene. From a total of 124 EST–SNPs, ninety-five were in exons and 21 of them caused amino acid changes, whereas 74 were synonymous SNPs. The other 29 EST-SNPs were located in the UTR regions (Table S5).
The GO functional classification of genes carrying among-group polymorphisms allowed us to verify their involvement in several biological processes (Tables S5 and S6). Most of them were regarding cellular processes (GO:0009987), response to stimulus and defense (GO:0050896; GO:0006952), and metabolic (GO:0008152) and biosynthetic processes (GO:0009058). The 124 genes were located in different cellular components (Table S7), mostly including membrane (GO:0016020), cytosol (GO:0005829), and cell periphery (GO:0071944). Furthermore, several molecular functions were evidenced within the 124 selected contigs (Table S8); binding (GO:0005488; GO:0043167; GO:1901363), catalytic activity (GO:0003824), and transferase activity (GO:0016740) were found.
The specific functions and processes of the 124 genes was also classified. Several genes were related to the abiotic stress response, such as salt stress—gibberellin-regulated protein 14-like GASA14 and nuclear fusion defective 4 NFD4, cold stress: catalase isozyme 3-like CAT3 and 60 S ribosomal protein L23-like, and drought stress—homeobox-leucine zipper protein ATHB-6 and ATHB-12, and aquaporin TIP2-1-like. The resultant gene V-type proton ATPase 16 were linked to cold and drought stresses. Moreover, the gene 30 S ribosomal protein 2 was related to all three abiotic stresses. Genes related to biotic stresses, such as nematode, virus, fungi, and bacteria resistance included the pectate lyase gene family A10 and 18, elongin-C-like, ADP-ribosylation factor-like protein 8a, probable pectinesterase/pectinesterase inhibitor PME34, and fatty acid amide hydrolase-like. Heterogeneous nuclear ribonucleoprotein Q-like was related to fungal and bacterial defense. Finally, cysteine protease RD19A-like was linked to both biotic and abiotic (bacterial, salt and drought) stresses. Seven genes were related to the flowering, five to seed development, whereas the FT-interacting protein 1-like and the two-component response regulator-like PRR37 were related to the photoperiod. Finally, six genes were related to fatty acid synthesis, or involved in the phenolic pathway, such as serine carboxypeptidase-like 7 SCPL and chalcone synthase J-like CHS, or the terpenoid pathway, like geraniol synthase chloroplastic-like GES, geraniol 8-hydroxylase-like, and geraniol 10-hydroxylase G10H.

4. Discussion

Plant domestication is an evolutionary process in which humans, through means of artificial selection, take advantage of the genetic diversity of a wild species and modify it, based on their needs [59]. Olive had a peculiar history in its domestication [16,22]. In fact, several cultivated varieties were exported from their original spreading area, located in the eastern side of Mediterranean basin, as seedlings and cuttings, and were diffused by grafting techniques, replacing, and decreasing the olive germplasm diversity of the Central and Western Mediterranean. This scenario played a key role, which contributed to cover and sometimes lose functional polymorphisms private to the western olive germplasm, including oleasters.
The SNP studies performed on olives up to now, were mainly focused on genetic map construction [45,60], genetic diversity studies on restricted set of cultivars [44,46], and few SNP positions were considered on the candidate genes to establish possible phenotypic/genotypic correlations [61,62]. Recently, a large set of SNPs developed on expressed sequence tags was applied for genotyping a wide set of samples, including cultivated and natural forms of olive [47]. The EST-SNP markers, deeply studied in the present work, were gene-based and, therefore, had a higher probability to represent a function-related single nucleotide variant, and were particularly valuable for characterizing genes associated with complex traits [63].
In order to identify genes carrying SNP mutations specific to cultivated and non-cultivated genotypes, in the present study, we deeply analyzed 1040 EST-SNPs on an enriched set from what was previously published [47] of 260 samples, including cultivars, oleasters, and subspecies guanchica.
Genetic diversity assessments, based on genetic distance parameters and population structure analyses, confirmed a clear separation between cultivars and oleasters, including trees of the subspecies guanchica [12,64,65,66]. The Bayesian analysis revealed that the cultivars were divided into geographically separated populations, including Western, Central, and Eastern Mediterranean groups, as already reported by the SSR genotyping studies [15,16]. Divergence between oleasters and cultivated varieties was higher than that observed among cultivars and the subsp. guanchica samples, as observed in the Net nucleotide distance analysis. An upcoming work will elucidate all these differences, especially the differences between the three phylogeographically separated varieties.
The average number of nucleotide differences and nucleotide diversity within each population demonstrated a lower level of nucleotide variation inside the oleasters and subsp. guanchica populations, with respect to the cultivated ones. This could also be explained by a lower rate of natural hybridization and introgression because of the narrow geographic environment [63], where the analyzed oleasters and guanchica trees were distributed. Conversely, the olive varieties ae cultivated across vast geographical areas [15] and, since they are mostly self-incompatible [67], the hybridization that occurred with the Eastern Mediterranean germplasm with local oleasters and subspecies might have contributed to increasing their variability [13,64,68]. In fact, the average number of differences and the nucleotide diversity were higher in the cultivated materials, with respect to the other groups. These results can reflect from one side, events of hybridization among local genotypes and, from the other, crossing between oleasters spreading in the Western Mediterranean with other olive genotypes coming from the near east [13,15]. The human artificial selection has favored heterozygous genotypes, contributing to the present genetic scenario, as already reported for other plant species [69].
On the entire set of samples, single-point mutations affecting both alleles and distinguishing the cultivated group from the others were deeply analyzed. The EST-SNPs under evaluation were particularly valuable to characterize the genes associated with functional traits [70,71]. Single nucleotide changes might lead to different levels of gene expression or variation on the protein structure, as reported for the non-synonymous as well as synonymous mutations [72,73,74]. Single-point mutations characterizing the three groups were related to 124 genes, and their molecular and biological functions, as well as their cellular component, were reported. Twenty-one EST-SNPs characterizing the cultivated group versus the oleasters and subsp. guanchica, generated nsSNPs, leading to amino acid changes. Some EST–SNPs carried by the UTR regions could play a role at the transcription level, or be linked to other functional polymorphisms. The identification of important loci and genes by accessing unique germplasm resources, such as wild populations, might allow their introgression into the cultivated varieties, in order to improve adaptation and tolerance to environmental constrains [75]. The annotation of our candidate genes revealed their involvement in important functional traits, such as stress response, flowering, and secondary metabolite synthesis. An example of genes showing interesting nsSNPs is eRF1, which binds the ribosome A site and recognizes all three stop codons in eukaryotes [76]. When the elongating ribosome encounters a termination codon, this release factor terminates translation and releases the polypeptide from the ribosome, and ribosomal subunits are recycled to participate in a new round of translation [77]. This gene is related to several protein formation and it was recently found to be essential (together with ESP1) in protein-accumulation mechanisms in rice, providing a genetic target for improving relevant agronomic traits [78].
Several genes related to terpenoid, phenol, and fatty acid synthesis were found, such as the fatty acids amide hydrolase (FAAH), where the SNP variation between oleasters and cultivars caused an amino acid change. The FAAH was recently reported to be strictly related to defoliation, mediated by the infection of plants from Verticillium dahliae. In fact, the level of FAAH gene expression considerably increased in infected cotton plants [79]. This gene is also well-known for the role in the arachinodic fatty acid synthesis, both in plants and humans, with implications for good health. Glucose-6-Phosphate Transporter 2 (GPT2), whose differential expression in olive under cold stress already showed an nsSNP mutation [80], had a key role in the promotion of cell proliferation, integrated in the growth regulatory network, to enable plant growth and development in continuously changing environmental conditions [81].
Genes related to flowering: One of the selected EST-SNPs was present on the sequence of floricaula/Leafy (LFY), a floral meristem identity gene that plays a key role in flower development [82,83], with a dual function in flower meristem identity regulation and in flowering timing [84]. The LFY promoter contains binding sites for several regulatory factors like Agamous-like MADS-box (AGL) [85], and in Citrus, the constitutive expression of LFY was found to reduce the juvenile phase [86]. An nsSNP was found between wild and cultivated olives in the XTH9 gene, whose expression was reduced in various mutants that had a defect in elongation, including acaulis pin-formed and terminal flower [87]. Moreover, under stress conditions, specific XTH isoforms were up-regulated, while others were down-regulated, suggesting a different involvement of the XTH enzymes on the effects of simultaneous exposure to dehydration and heat conditions, naturally occurring in open fields [88]. An SNP was detected on the 3’UTR region of an FT-interacting protein 1-like. Flowering time in crops represents a critical determinant of the distribution and regional adaptability of plants. In the last decade, knowledge of the florigenic function of Flowering locus T (FT) and FT-like proteins advanced substantially, but the process through which plants control the long-distance action of the FT protein remains largely unknown [89]. The accumulation of FT-like gene transcripts was documented as preceding flower initiation [90]. Expression of OeFT1/2 in olive leaves and OeFT2 in buds was recently studied together, with the expression in transgenic olive cuttings inserting this gene from Arabidopsis, giving very promising results [91]. Two-component response regulator-like gene (PRR37) played a crucial role as a floral repressor, enabling rice to be domesticated and cultivated in high-latitude regions [92]. In olive, this gene is located on chromosome 1 and the single nucleotide change modified the amino acid from aspartate to glutamate. Specific uptake inhibitor (SUI1), which is a protein translation factor, is placed on chromosome 6 in olive, and the nucleotide variation changes the amino acid from arginine to lysine. Overexpression of SUI1 in transgenic rice resulted in an increase in the number of floral organs [93]. Another SNP was detected in the 3’UTR region of the Glycerol-3-phosphate SN-2-acyltransferase 4 gene (GPAT4). A recent examination of pollen and embryo sac, suggested that down-regulation of GPAT4 caused impaired female fertility, which is the main reason for a reduced seed set [94,95]. A synonymous polymorphism was detected within the tetraketide α-pyrone reductase 1-like (TKPR1), a gene involved in the biosynthesis of hydroxylated tetraketide compounds that serve as sporopollenin precursors, the main constituents of pollen coat exine [96]. Another synonymous polymorphism was detected in the 6-phosphogluconate dehydrogenase gene (PGD), described for its important role in the growth of the male gametophytes and pollen-tube–ovule interaction in Arabidopsis [97].
Genes related to stress response: Single nucleotide mutations between wild and cultivated olives generating amino acid changes were detected in two WNK kinase genes (WNK4 and WNK9), respectively, from lysine to arginine (WNK4) and from serine to threonine (WNK9). In plants, the WNK kinase is a large family, with mechanisms to nudge intracellular proteins to their destinations, still to be recognized as a nexus for trafficking-based signal transduction [98]. A synonymous mutation was detected on a hypersensitive-induced response protein-like (HIR). Members of the HIR family are implicated in ion-channel regulation, cell cycle control, and hypersensitive reactions like programmed cell death, in response to pathogen attacks and spontaneous hypersensitive lesions, such as TaHIR3 in wheat and OsHIR1 in rice [99,100]. An SNP change from the apolar glycine to the polar serine in a Pectin Methyl Esterase Inhibitor (PMEI), catalyzes the dimethyl esterification of the homogalacturonan domains of pectin in the plant cell wall, which functions as the primary barrier against pathogens. GhPMEI3 enhances resistance to Verticillium wilt in transgenic Arabidopsis [101]. A synonymous mutation was detected on the sequence of abscisic acid receptor (PYL8), required for the ABA-mediated responses, such as stomatal closure and germination inhibition. PYL8 positively regulates ABA signaling during germination and abiotic stress responses. Furthermore, data showed that PYL8 plays essential roles in response to dark-induced senescence [102] and in drought stress [103]. An SNP variation was detected on a Cysteine protease response to dehydration (RD19A-like), which is a responsive gene to salt stress. The cysteine protease genes RD21a and RD19a of Arabidopsis thaliana are induced by water deficit and are responsive to salt stress; these play an important role in the programmed cell death pathway during stress [104]. Another SNP was found in a Stress-Related Protein gene (SRPs), found only in higher plants. These lipid droplet (LD) surface proteins are involved in post germinative growth, diurnal regulation of LDs and various stress responses independently or synergistically. However, the cellular mechanism of SRPs in response to abiotic stresses remains unclear. It is possible that SRP-mediated LDs interact with other subcellular organelles, like mitochondria, peroxisomes, and vacuoles under stress conditions [105]. Calcineurin B-like (CBLs) are EF-hand Ca2+ protein sensors, and upon Ca2+ binding, they undergo conformational changes to associate with a group of CBL-Interacting Protein Kinases (CIPKs). Different combinations of CBLs and CIPKs complexes might generate temporal and spatial specificity in Ca2+ signaling, integrating various stimuli to determine cellular responses. CBL could regulate cold-stress-induced signal transduction [106]. The SNP detected on this sequence did not determine any change on the protein sequence. These three latter genes were already tested in a recent trial of salt stress tolerance in olive with very promising results on their different levels of expression under different salt concentration and times (Mousavi et al., under submission). A Probable serine/threonine-protein kinase gene (PBL8), which regulates plant immunity, and might be involved in plant defense signaling [107], holds an inter-population SNP. A mutation distinguished a Protein RNA-directed DNA methylation like (RdDM), an essential factor for suppressing transposons, DNA damage caused by stress, and DNA nucleases. Studies of the RdDM pathway showed that the RDM family proteins are important regulatory factors with various functions that regulate this pathway [108]. A 3’UTR polymorphism was detected on a polyadenylate-binding protein (PABP), involved in salt-tolerance of halophyte plants. This gene plays critical roles in eukaryotic translation initiation and mRNA stabilization/degradation [109]. The Aquaporin-TIP2 is a gene strictly correlated with salt stress and intracellular water transport. The detected SNP is located on chromosome 22 in the olive genome. In tomato, when this is overexpressed, significant increase in fruit yield, harvest index, and plant mass are obtained, and its expression varies in relation to salt concentration in different plant tissues [110]. SlTIP2-2 is a suitable candidate gene to improve important agronomic traits, as well as stress response.
Genes related to metabolic pathways: Two SNPs were found in two different sn-glycerol-3-phosphate acyltransferases (GPAT and GPAT4), a gene family involved in acyl-lipid biosynthesis, triggering storage lipid biosynthesis and catalyzing triacylglycerol (TAG) accumulation, involved in membrane lipid and oil biosynthesis and lipid droplet production in developing pollen grains, which play a pivotal role in plant development [111]. A serine carboxypeptidase -like 7 (SCPL) was placed on chromosome 23, and carried a mutation determining an amino acid change from the apolar glycine to the polar glutamine. The gene was related to phenolic synthesis and it was detected in olive as a gene containing an EST–SSR motif [41]. Phenolic compounds are secondary metabolites involved in several processes, including resistance to biotic and abiotic stresses. The biosynthetic pathways leading to the vast diversity of plant phenolic products often include an acylation step, with phenolic compounds being the donor or acceptor molecules [112]. The biosynthesis of flavonoids from the phenylpropanoid pathway involves the chalcone synthase enzyme. CHS catalyzes the stepwise condensation of three acetate residues from naringenin chalcone, which leads to the synthesis of a variety of flavonoid derivatives [113]. As the first committed enzyme of flavonoid biosynthesis, CHS is rigorously controlled in response to a wide range of environmental and developmental stimuli. Flavonoids also function as antimicrobial compounds (phytoalexins) or insect repellents defending against phytopathogens and herbivores, and some flavonoid metabolites have pharmacological activities [114]. A chalcone synthase gene was already sequenced in olive, revealing a total of nine putative SNPs [115]. In the present study, a single sSNP was detected. An SNP variation was in charge of a Geraniol 10-hydroxylase (G10H), thought to play an important role in iridoid monoterpenoid biosynthesis. Recently, it was shown that the CrDXS and CrG10H genes isolated from Catharanthus roseus are important genes in regulating terpenoid indole alkaloid biosynthesis, filling knowledge gaps in understanding the pathway [116,117,118,119]. A 3’UTR SNP was located on a biotin carboxyl carrier protein of acetyl-CoA carboxylase 1 (BCCP), the first committed step in fatty acid synthesis, mediated by acetyl-CoA carboxylase (ACCase), a biotin-dependent enzyme that carboxylates acetyl-CoA to produce malonyl-CoA. The heteromeric ACCase is composed of four distinct subunits—biotin carboxyl transfer protein (BCCP), biotin carboxylase (BC), and a/b-carboxyltransferases (CT). These subunits are organized in two separable subcomplexes, one comprising BCCP and BC, which mediates the carboxylation of biotin, and a second a-CT/b-CT complex that catalyzes the carboxyltransferase reaction [120,121]. Olive produces a range of secondary metabolites that strongly affect the taste and nutritional properties of olive oil and fruits. Among them, the most abundant are the secoiridoids and monoterpenoids with a 3,4-dihydropyran skeleton. Some amino acid substitutions were observed in the conserved motifs of the olive 1,4-reductase [117], but the SNP identified in this work did not produce any amino acid change.
Among the EST-SNPs distinguishing cultivated and wild olive populations, 21 determining amino acid changes of enzymes directly affecting plant performance, represent an important step towards the identification of useful markers for the functional screening of olive resources or sequences useful for genetic improvement by classical breeding or via cisgenesis.

5. Conclusions

Tracing the selection footprint in combination with the identification of causative variations, would inform the future detection of important loci correlated to plant behavior and cultivation. The identification of markers specific to cultivated or wild olives might represent a new resource for the introgression of important variation into the cultivated olives through newly developed biotechniques, such as cisgenesis. The present study is an evidence of the importance of safeguarding the genetic diversity of the Olea europaea species. In fact, the availability in situ and ex situ of wild plants, together with the knowledge of variation carried by oleasters will allow the restoration or establishment of new allelic combinations, giving new opportunities for improving adaptation, stress response, or modifying the main metabolic pathways in cultivated olive, also to overcome future scenarios of climate change.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/8/916/s1, Figure S1: Genes carrying SNP polymorphisms between cultivated, oleaster, and subsp. guanchica olives. Table S1. List of olive samples analyzed, country/place of origin/diffusion. Table S2: Allele-frequency divergence among populations (Net nucleotide distance), computed using point estimates of P, using the Structure software. EM, WM, and CM, respectively, refer to East, West, and Central Mediterranean cultivars. Table S3: Genetic diversity for each population and the overall sample set, performed by the DnaSP software. Table S4: Haploid Allele Frequencies and Sample Size by Population differences in homozygosity state. Blue loci: Cultivated olives versus oleasters and subp. Guanchica. iPurple loci: Cultivars versus oleasters. Green loci: Cultivars versus subsp. Guanchica. Table S5: List of the 124 genes differentiating olive cultivars from oleasters and subsp. guanchica samples. Name of contigs, linkage group placement, NCBI accession number, gene name, and molecular process and function for each gene are reported. Table S6: List of the biological processes detected for the 124 genes differentiating olive cultivars from oleasters and subsp. guanchica samples. The analysis was performed by the PANTHER Overrepresentation Test. Table S7: List of all cellular components detected for the 124 genes differentiating olive cultivars from the oleasters and subsp. guanchica samples. The analysis was performed by the PANTHER Overrepresentation Test. Table S8: List of all molecular functions detected for the 124 genes differentiating olive cultivars from the oleasters and subsp. guanchica samples. The analysis was performed by the PANTHER Overrepresentation Test.

Author Contributions

Conceptualization: R.M. and S.M.; the acquisition of the work: R.M., A.B.; F.B., L.B., and S.M.; the performance of all the analysis of the work: R.M., F.B., A.B., L.L., and R.D.L.R.; drafting the work: R.M., S.M., A.B., and L.B.; funding acquisition: L.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Rural Development Program of the Umbria Region, 2014–2020—Measure 16.2.1, projects INNO.V.O. EU project OLIVE 4 CLIMATE—LIFE (LIFE15 CCM/IT/000141).

Conflicts of Interest

The authors declare no competing interests. Funders played no role in the study design; data collection, analysis or interpretation; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Barazani, O.; Keren-Keiserman, A.; Westberg, E.; Hanin, N.; Dag, A.; Ben-Ari, G.; Fragman-Sapir, O.; Tugendhaft, Y.; Kerem, Z.; Kadereit, J.W. Genetic variation of naturally growing olive trees in Israel: From abandoned groves to feral and wild? BMC Plant Biol. 2016, 16, 261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Di Rienzo, V.; Sion, S.; Taranto, F.; D’Agostino, N.; Montemurro, C.; Fanelli, V.; Sabetta, W.; Boucheffa, S.; Tamendjari, A.; Pasqualone, A.; et al. Genetic flow among olive populations within the Mediterranean basin. PeerJ 2018, 6, e5260. [Google Scholar] [CrossRef] [Green Version]
  3. Kassa, A.; Konrad, H.; Geburek, T. Molecular diversity and gene flow within and among different subspecies of the wild olive (Olea europaea L.): A review. Flora Morphol. Distrib. Funct. Ecol. Plants 2019, 250, 18–26. [Google Scholar] [CrossRef]
  4. Vargas, P.; Muñoz, G.; Hess, J.; Kadereit, J. Olea europaea subsp. guanchica and subsp. maroccana (Oleaceae), two new names for olive tree relatives. In Anales Del Jardín Botánico de Madrid; Instituto Botánico, Real Jardin de Botánico CSIC: Madrid, Spain, 2000; pp. 360–361. [Google Scholar]
  5. Green, P.S. A Revision of Olea L. (Oleaceae). Kew Bull. 2002, 57, 91–140. [Google Scholar] [CrossRef]
  6. Unver, T.; Wu, Z.; Sterck, L.; Turktas, M.; Lohaus, R.; Li, Z.; Yang, M.; He, L.; Deng, T.; Escalante, F.J.; et al. Genome of wild olive and the evolution of oil biosynthesis. Proc. Natl. Acad. Sci. USA 2017, 114, E9413–E9422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Lumaret, R.; Ouazzani, N. Ancient wild olives in Mediterranean forests. Nature 2001, 413, 700. [Google Scholar] [CrossRef]
  8. Hmmam, I.; Mariotti, R.; Ruperti, B.; Cultrera, N.; Baldoni, L.; Barcaccia, G. Venetian olive (Olea europaea) germplasm: Disclosing the genetic identity of locally grown cultivars suited for typical extra virgin oil productions. Genet. Resour. Crop. Evol. 2018, 65, 1733–1750. [Google Scholar] [CrossRef]
  9. Lazović, B.; Adakalić, M.; Pucci, C.; Perović, T.; Bandelj, D.; Belaj, A.; Mariotti, R.; Baldoni, L. Characterizing ancient and local olive germplasm from Montenegro. Sci. Hortic. 2016, 209, 117–123. [Google Scholar] [CrossRef]
  10. Mousavi, S.; Mazinani, S.M.H.; Arzani, K.; Ydollahi, A.; Pandolfi, S.; Baldoni, L.; Mariotti, R. Molecular and morphological characterization of Golestan (Iran) olive ecotypes provides evidence for the presence of promising genotypes. Genet. Resour. Crop. Evol. 2014, 61, 775–785. [Google Scholar] [CrossRef]
  11. Mousavi, S.; Stanzione, V.; Mencuccini, M.; Baldoni, L.; Bufacchi, M.; Mariotti, R. Biochemical and molecular profiling of unknown olive genotypes from central Italy: Determination of major and minor components. Eur. Food Res. Technol. 2019, 245, 83–94. [Google Scholar] [CrossRef]
  12. Ninot, A.; Howad, W.; Aranzana, M.J.; Senar, R.; Romero, A.; Mariotti, R.; Baldoni, L.; Belaj, A. Survey of over 4, 500 monumental olive trees preserved on-farm in the northeast Iberian Peninsula, their genotyping and characterization. Sci. Hortic. 2018, 231, 253–264. [Google Scholar] [CrossRef] [Green Version]
  13. Besnard, G.; El Bakkali, A.; Haouane, H.; Baali-Cherif, D.; Moukhli, A.; Khadari, B. Population genetics of Mediterranean and Saharan olives: Geographic patterns of differentiation and evidence for early generations of admixture. Ann. Bot. 2013, 112, 1293–1302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Mazinani, S.M.H.; Mariotti, R.; Torkzaban, B.; Sheikh-Hassani, M.; Ataei, S.; Cultrera, N.; Pandolfi, S.; Baldoni, L. High Genetic Diversity Detected in Olives beyond the Boundaries of the Mediterranean Sea. PLoS ONE 2014, 9, e93146. [Google Scholar] [CrossRef] [Green Version]
  15. Belaj, A.; Gurbuz Veral, M.; Sikaoui, H.; Moukhli, A.; Khadari, B.; Mariotti, R.; Baldoni, L. Olive Genetic Resources. In The Olive Tree Genome, Compendium of Plant Genomes; Rugini, E., Baldoni, L., Muleo, R., Sebastiani, L., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 27–54. [Google Scholar]
  16. Mousavi, S.; Mariotti, R.; Bagnoli, F.; Costantini, L.; Cultrera, N.; Arzani, K.; Pandolfi, S.; Vendramin, G.G.; Torkzaban, B.; Mazinani, S.M.H.; et al. The eastern part of the Fertile Crescent concealed an unexpected route of olive (Olea europaea L.) differentiation. Ann. Bot. 2017, 119, 1305–1318. [Google Scholar] [CrossRef]
  17. Díez, C.M.; Trujillo, I.; Martinez-Urdiroz, N.; Barranco, D.; Rallo, L.; Marfil, P.; Gaut, B.S.; Ruiz, P.F.M. Olive domestication and diversification in the Mediterranean Basin. New Phytol. 2015, 206, 436–447. [Google Scholar] [CrossRef]
  18. Besnard, G.; De Casas, R.R. Single vs multiple independent olive domestications: The jury is (still) out. New Phytol. 2016, 209, 466–470. [Google Scholar] [CrossRef] [Green Version]
  19. Dighton, A.; Fairbairn, A.S.; Bourke, S.; Faith, J.T.; Habgood, P. Bronze Age olive domestication in the north Jordan valley: New morphological evidence for regional complexity in early arboricultural practice from Pella in Jordan. Veg. Hist. Archaeobotany 2017, 26, 403–413. [Google Scholar] [CrossRef]
  20. Julca, I.; Marcet-Houben, M.; Vargas, P.; Gabaldon, T. Phylogenomics of the olive tree (Olea europaea) reveals the relative contribution of ancient allo- and autopolyploidization events. BMC Biol. 2018, 16, 15. [Google Scholar] [CrossRef] [Green Version]
  21. Langgut, D.; Cheddadi, R.; Carrión, J.; Cavanagh, M.; Colombaroli, D.; Eastwood, W.J.; Greenberg, R.; Litt, T.; Mercuri, A.M.; Miebach, A.; et al. The origin and spread of olive cultivation in the Mediterranean Basin: The fossil pollen evidence. Holocene 2019, 29, 902–922. [Google Scholar] [CrossRef]
  22. Jiménez-Ruiz, J.; Ramírez-Tejero, J.A.; Fernández-Pozo, N.; Leyva-Pérez, M.d.l.O.; Yan, H.; Rosa, R.d.l.; Belaj, A.; Montes, E.; Rodríguez-Ariza, M.O.; Navarro, F.; et al. Transposon activation is a major driver in the genome evolution of cultivated olive trees ( Olea europaea L.). Plant Genome 2020, 13, e20010. [Google Scholar] [CrossRef] [Green Version]
  23. Besnard, G.; Terral, J.; Cornille, A. On the origins and domestication of the olive: A review and perspectives. Ann. Bot. 2018, 121, 385–403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Besnard, G.; Khadari, B.; Navascués, M.; Fernández-Mazuecos, M.; El Bakkali, A.; Arrigo, N.; Baali-Cherif, D.; De Caraffa, V.B.-B.; Santoni, S.; Vargas, P.; et al. The complex history of the olive tree: From Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant. Proc. R. Soc. B Biol. Sci. 2013, 280, 20122833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Voropaeva, I.; Stika, H.-P. Plant remains from Roman period town of Obulco (today Porcuna) in Andalusia (Spain)—Distribution and domestication of olive in the west Mediterranean region. Archaeol. Anthr. Sci. 2018, 10, 865–882. [Google Scholar] [CrossRef]
  26. Naz, S.; Khan, M.R.; Awan, A.A.; Hussain, M.; Ali, S. Diversity and divergence in cultivated and wild olive germplasm collected from northern Pakistan. Int. J. Agric. Biol. 2019, 22, 1109–1115. [Google Scholar]
  27. Ashraf, U.; Peterson, A.T.; Chaudhry, M.N.; Ashraf, I.; Saqib, Z.; Ahmad, S.R.; Ali, H. Ecological niche model comparison under different climate scenarios: A case study of Oleaspp. in Asia. Ecosphere 2017, 8, e01825. [Google Scholar] [CrossRef]
  28. Qiu, J.; Zhou, Y.; Mao, L.; Ye, C.; Wang, W.; Zhang, J.; Yu, Y.; Fu, F.; Wang, Y.; Qian, F.; et al. Genomic variation associated with local adaptation of weedy rice during de-domestication. Nat. Commun. 2017, 8, 15323. [Google Scholar] [CrossRef]
  29. Lu, X.; Fu, X.; Wang, D.; Wang, J.; Chen, X.; Hao, M.; Wang, J.; Gervers, K.A.; Guo, L.; Wang, S.; et al. Resequencing of cv CRI-12 family reveals haplotype block inheritance and recombination of agronomically important genes in artificial selection. Plant Biotechnol. J. 2018, 17, 945–955. [Google Scholar] [CrossRef]
  30. Ruban, A.V. Plant science: Crops on the fast track for light. Nature 2017, 541, 36. [Google Scholar] [CrossRef]
  31. Arora, S.; Steuernagel, B.; Gaurav, K.; Chandramohan, S.; Long, Y.; Matny, O.; Johnson, R.; Enk, J.; Periyannan, S.; Singh, N.; et al. Resistance gene cloning from a wild crop relative by sequence capture and association genetics. Nat. Biotechnol. 2019, 37, 139–143. [Google Scholar] [CrossRef]
  32. Doebley, J.; Gaut, B.S.; Smith, B.D. The Molecular Genetics of Crop Domestication. Cell 2006, 127, 1309–1321. [Google Scholar] [CrossRef] [Green Version]
  33. Fan, C.; Xing, Y.; Mao, H.; Lu, T.; Han, B.; Xu, C.; Li, X.; Zhang, Q. GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor. Appl. Genet. 2006, 112, 1164–1171. [Google Scholar] [CrossRef] [PubMed]
  34. Luikart, G.; England, P.R.; Tallmon, D.; Jordan, S.; Taberlet, P. The power and promise of population genomics: From genotyping to genome typing. Nat. Rev. Genet. 2003, 4, 981–994. [Google Scholar] [CrossRef]
  35. Zhou, Y.; Massonnet, M.; Sanjak, J.S.; Cantu, D.; Gaut, B.S. Evolutionary genomics of grape (Vitis vinifera ssp. vinifera) domestication. Proc. Natl. Acad. Sci. USA 2017, 114, 11715–11720. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Miller, A.J.; Gross, B.L. From forest to field: Perennial fruit crop domestication. Am. J. Bot. 2011, 98, 1389–1414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. This, P.; Lacombe, T.; Cadle-Davidson, M.; Owens, C.L. Wine grape (Vitis vinifera L.) color associates with allelic variation in the domestication gene VvmybA1. Theor. Appl. Genet. 2007, 114, 723–730. [Google Scholar] [CrossRef]
  38. Gros-Balthazard, M.; Besnard, G.; Sarah, G.; Holtz, Y.; Leclercq, J.; Santoni, S.; Wegmann, D.; Glémain, S.; Khadari, B. Evolutionary transcriptomics reveals the origins of olives and the genomic changes associated with their domestication. Plant J. 2019, 100, 143–157. [Google Scholar] [CrossRef] [PubMed]
  39. Tao, Y.; Mace, E.S.; Tai, S.; Cruickshank, A.; Campbell, B.C.; Zhao, X.; Van Oosterom, E.J.; Godwin, I.D.; Botella, M.A.; Jordan, D.R. Whole-Genome Analysis of Candidate genes Associated with Seed Size and Weight in Sorghum bicolor Reveals Signatures of Artificial Selection and Insights into Parallel Domestication in Cereal Crops. Front. Plant Sci. 2017, 8, 8. [Google Scholar] [CrossRef] [Green Version]
  40. Arbeiter, A.B.; Hladnik, M.; Jakše, J.; Bandelj, D. Identification and validation of novel EST-SSR markers in olives. Sci. Agricola 2017, 74, 215–225. [Google Scholar] [CrossRef]
  41. De La Rosa, R.; Belaj, A.; Muñoz-Mérida, A.; Trelles, O.; Ortíz-Martín, I.; González-Plaza, J.J.; Valpuesta, V.; Beuzón, C.R. Development of EST-derived SSR Markers with Long-core Repeat in Olive and Their Use for Paternity Testing. J. Am. Soc. Hortic. Sci. 2013, 138, 290–296. [Google Scholar] [CrossRef]
  42. Mariotti, R.; Cultrera, N.G.M.; Mousavi, S.; Baglivo, F.; Rossi, M.; Albertini, E.; Alagna, F.; Carbone, F.; Perrotta, G.; Baldoni, L. Development, evaluation, and validation of new EST-SSR markers in olive (Olea europaea L.). Tree Genet. Genomes 2016, 12, 120. [Google Scholar] [CrossRef]
  43. Kyriakopoulou, C.I.; Kalogianni, D.P. Genetic Identification of the Wild Form of Olive (Olea europaea var. sylvestris) Using Allele-Specific Real-Time PCR. Foods 2020, 9, 467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. 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, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Biton, I.; Doron-Faigenboim, A.; Jamwal, M.; Mani, Y.; Eshed, R.; Rosen, A.; Sherman, A.; Ophir, R.; Lavee, S.; Avidan, B.; et al. Development of a large set of SNP markers for assessing phylogenetic relationships between the olive cultivars composing the Israeli olive germplasm collection. Mol. Breed. 2015, 35, 107. [Google Scholar] [CrossRef]
  46. Ipek, A.; Yılmaz, K.; Sıkıcı, P.; Tangu, N.A.; Öz, A.T.; Bayraktar, M.; Ipek, M.; Gulen, 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]
  47. 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]
  48. 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. [Google Scholar] [CrossRef] [Green Version]
  49. García-Verdugo, C.; Fay, M.F.; Granado-Yela, C.; Rubio de Casas, R.; Balaguer, L.; Besnard, G.; Vargas, P. Genetic diversity and differentiation processes in the ploidy series of Olea europaea L.: A multiscale approach from subsp. to insular populations. Mol. Ecol. 2009, 18, 454–467. [Google Scholar] [CrossRef]
  50. Hess, J.; Kadereit, J.W.; Vargas, P. The colonization history of Olea europaea L. in Macaronesia based on internal transcribed spacer 1 (ITS-1) sequences, randomly amplified polymorphic DNAs (RAPD), and intersimple sequence repeats (ISSR). Mol. Ecol. 2000, 9, 857–868. [Google Scholar] [CrossRef] [Green Version]
  51. Peakall, R.; Smouse, P.E. Genalex 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 2006, 6, 288–295. [Google Scholar] [CrossRef]
  52. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In Nucleic Acids Symposium Series; Information Retrieval Ltd.: London, UK, 1999; c1979–c2000; pp. 95–98. [Google Scholar]
  53. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGAX: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  54. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [PubMed]
  55. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [PubMed]
  56. Earl, D.A.; Vonholdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2011, 4, 359–361. [Google Scholar] [CrossRef]
  57. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; E Ramos-Onsins, S.; Sánchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  58. Thomas, D.C.; Campbell, M.J.; Kejariwal, A.; Mi, H.; Karlak, B.; Daverman, R.; Diemer, K.; Muruganujan, A.; Narechania, A. PANTHER: A Library of Protein Families and Subfamilies Indexed by Function. Genome Res. 2003, 13, 2129–2141. [Google Scholar] [CrossRef] [Green Version]
  59. Mastretta-Yanes, A.; Gasman, F.A.; Burgeff, C.; Ramírez, M.C.; Piñero, D.; Sarukhán, J. An Initiative for the Study and Use of Genetic Diversity of Domesticated Plants and Their Wild Relatives. Front. Plant Sci. 2018, 9, 209. [Google Scholar] [CrossRef] [Green Version]
  60. Mariotti, R.; Fornasiero, A.; Mousavi, S.; Cultrera, N.G.; Brizioli, F.; Pandolfi, S.; Passeri, V.; Rossi, M.; Magris, G.; Scalabrin, S.; et al. Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker. Front. Plant Sci. 2020, 10, 10. [Google Scholar] [CrossRef]
  61. Hedayati, V.; Mousavi, A.; Razavi, K.; Cultrera, N.; Alagna, F.; Mariotti, R.; Mazinani, S.M.H.; Baldoni, L. Polymorphisms in the AOX2 gene are associated with the rooting ability of olive cuttings. Plant Cell Rep. 2015, 34, 1151–1164. [Google Scholar] [CrossRef]
  62. Cultrera, N.G.M.; Sarri, V.; Lucentini, L.; Ceccarelli, M.; Alagna, F.; Mariotti, R.; Mousavi, S.; Ruiz, C.G.; Baldoni, L. High Levels of Variation Within Gene Sequences of Olea europaea L. Front. Plant Sci. 2019, 9, 1932. [Google Scholar] [CrossRef] [Green Version]
  63. Yang, H.; Wei, C.-L.; Liu, H.-W.; Wu, J.-L.; Li, Z.-G.; Zhang, L.; Jian, J.-B.; Li, Y.-Y.; Tai, Y.-L.; Zhang, J.; et al. Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing. PLoS ONE 2016, 11, e0151424. [Google Scholar] [CrossRef] [Green Version]
  64. Baldoni, L.; Tosti, N.; Ricciolini, C.; Belaj, A.; Arcioni, S.; Pannelli, G.; Germanà, M.A.; Mulas, M.; Porceddu, A. Genetic Structure of Wild and Cultivated Olives in the Central Mediterranean Basin. Ann. Bot. 2006, 98, 935–942. [Google Scholar] [CrossRef] [PubMed]
  65. 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]
  66. Terral, J.-F.; Alonso, N.; Capdevila, R.B.I.; Chatti, N.; Fabre, L.; Fiorentino, G.; Marinval, P.; Jordà, G.P.; Pradat, B.; Rovira, N.; et al. Historical biogeography of olive domestication (Olea europaea L.) as revealed by geometrical morphometry applied to biological and archaeological material. J. Biogeogr. 2004, 31, 63–77. [Google Scholar] [CrossRef]
  67. Saumitou-Laprade, P.; Vernet, P.; Vekemans, X.; Billiard, S.; Gallina, S.; Essalouh, L.; Mhaïs, A.; Moukhli, A.; El Bakkali, A.; Barcaccia, G.; et al. Elucidation of the genetic architecture of self-incompatibility in olive: Evolutionary consequences and perspectives for orchard management. Evol. Appl. 2017, 10, 867–880. [Google Scholar] [CrossRef]
  68. Beghè, D.; Piotti, A.; Satovic, Z.; de la Rosa, R.; Belaj, A. Pollen-mediated gene flow and fine-scale spatial genetic structure in Olea europaea subsp. europaea var. sylvestris. Ann. Bot. 2017, 119, 671–679. [Google Scholar]
  69. Wang, J.; Chu, S.; Zhang, H.; Zhu, Y.; Cheng, H.; Yu, D. Development and application of a novel genome-wide SNP array reveals domestication history in soybean. Sci. Rep. 2016, 6, 20728. [Google Scholar] [CrossRef]
  70. Bettaieb, I.; Bouktila, D. Genome-wide analysis of NBS-encoding resistance genes in the Mediterranean olive tree (Olea europaea subsp. europaea var. europaea): Insights into their molecular diversity, evolution and function. Tree Genet. Genomes 2020, 16, 23. [Google Scholar] [CrossRef]
  71. Nelson, M.R.; Marnellos, G.; Kämmerer, S.; Hoyal, C.R.; Shi, M.M.; Cantor, C.R.; Braun, A. Large-Scale Validation of Single Nucleotide Polymorphisms in Gene Regions. Genome Res. 2004, 14, 1664–1668. [Google Scholar] [CrossRef] [Green Version]
  72. Bartoszewski, R.A.; Jablonsky, M.; Bartoszewska, S.; Stevenson, L.; Dai, Q.; Kappes, J.; Collawn, J.F.; Bebok, Z. A Synonymous Single Nucleotide Polymorphism in ΔF508CFTRAlters the Secondary Structure of the mRNA and the Expression of the Mutant Protein. J. Biol. Chem. 2010, 285, 28741–28748. [Google Scholar] [CrossRef] [Green Version]
  73. Sauna, Z.E.; Kimchi-Sarfaty, C. Understanding the contribution of synonymous mutations to human disease. Nat. Rev. Genet. 2011, 12, 683–691. [Google Scholar] [CrossRef]
  74. She, R.; Jarosz, D.F. Mapping Causal Variants with Single Nucleotide Resolution Reveals the Biochemical Drivers of Phenotypic Change. Cell 2018, 172, 478–490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Qi, X.; Li, M.-W.; Xie, M.; Liu, X.; Ni, M.; Shao, G.; Song, C.; Yim, A.K.-Y.; Tao, Y.; Wong, F.-L.; et al. Identification of a novel salt tolerance gene in wild soybean by whole-genome sequencing. Nat. Commun. 2014, 5, 4340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Brogna, S.; McLeod, T.; Petric, M. The Meaning of NMD: Translate or Perish. Trends Genet. 2016, 32, 395–407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Mateyak, M.; Sasikumar, A.N.; Dunaway, S.; Kinzy, T.G. The Unique Evolutionary Distribution of Eukaryotic Elongation Factor 3. In Evolution of the Protein Synthesis Machinery and Its Regulation; Springer Science and Business Media LLC: Berlin, Germany, 2016; pp. 313–326. [Google Scholar]
  78. Elakhdar, A.; Ushijima, T.; Fukuda, M.; Yamashiro, N.; Kawagoe, Y.; Kumamaru, T. Eukaryotic peptide chain release factor 1 participates in translation termination of specific cysteine-poor prolamines in rice endosperm. Plant Sci. 2018, 281, 223–231. [Google Scholar] [CrossRef] [PubMed]
  79. Zhang, D.-D.; Wang, J.; Wang, D.; Kong, Z.-Q.; Zhou, L.; Zhang, G.-Y.; Gui, Y.-J.; Li, J.-J.; Huang, J.-Q.; Wang, B.-L.; et al. Population genomics demystifies the defoliation phenotype in the plant pathogen Verticillium dahliae. New Phytol. 2019, 222, 1012–1029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Guerra, D.; Lamontanara, A.; Bagnaresi, P.; Orrù, L.; Rizza, F.; Zelasco, S.; Beghè, D.; Ganino, T.; Pagani, D.; Cattivelli, L.; et al. Transcriptome changes associated with cold acclimation in leaves of olive tree (Olea europaea L.). Tree Genet. Genomes 2015, 11, 1–24. [Google Scholar] [CrossRef]
  81. Van Dingenen, J.; Vermeersch, M.; De Milde, L.; Hulsmans, S.; De Winne, N.; Van Leene, J.; Gonzalez, N.; Dhondt, S.; De Jaeger, G.; Rolland, F.A.; et al. The role of HEXOKINASE1 in Arabidopsis leaf growth. Plant Mol. Biol. 2018, 99, 79–93. [Google Scholar] [CrossRef]
  82. Moyroud, E.; Kusters, E.; Monniaux, M.; Koes, R.; Parcy, F. LEAFY blossoms. Trends Plant Sci. 2010, 15, 346–352. [Google Scholar] [CrossRef]
  83. Fernández-Ocaña, A.M.; García-López, M.C.; Jiménez-Ruiz, J.; Saniger, L.; Macías, D.; Navarro, F.; Oya, R.; Belaj, A.; De La Rosa, R.; Corpas, F.J.; et al. Identification of a gene involved in the juvenile-to-adult transition (JAT) in cultivated olive trees. Tree Genet. Genomes 2010, 6, 891–903. [Google Scholar] [CrossRef]
  84. Liu, C.; Thong, Z.; Yu, H. Coming into bloom: The specification of floral meristems. Development 2009, 136, 3379–3391. [Google Scholar] [CrossRef] [Green Version]
  85. Ding, F.; Zhang, S.; Chen, H.; Peng, H.; Lu, J.; He, X.; Pan, J. Functional analysis of a homologue of the FLORICAULA/LEAFY gene in litchi (Litchi chinensis Sonn.) revealing its significance in early flowering process. Genes Genom. 2018, 40, 1259–1267. [Google Scholar] [CrossRef] [PubMed]
  86. Peña, L.; Martín-Trillo, M.; Juárez, J.; Pina, J.A.; Navarro, L.; Martinez-Zapater, J.M. Constitutive expression of Arabidopsis LEAFY or APETALA1 genes in citrus reduces their generation time. Nat. Biotechnol. 2001, 19, 263–267. [Google Scholar] [CrossRef] [PubMed]
  87. Lee, Y.K.; Rhee, J.Y.; Lee, S.H.; Chung, G.C.; Park, S.J.; Segami, S.; Maeshima, M.; Choi, G. Functionally redundant LNG3 and LNG4 genes regulate turgor-driven polar cell elongation through activation of XTH17 and XTH24. Plant Mol. Biol. 2018, 97, 23–36. [Google Scholar] [CrossRef] [PubMed]
  88. Iurlaro, A.; De Caroli, M.; Sabella, E.; De Pascali, M.; Rampino, P.; De Bellis, L.; Perrotta, C.; Dalessandro, G.; Piro, G.; Fry, S.C.; et al. Drought and Heat Differentially Affect XTH Expression and XET Activity and Action in 3-Day-Old Seedlings of Durum Wheat Cultivars with Different Stress Susceptibility. Front. Plant Sci. 2016, 7, 1686. [Google Scholar] [CrossRef]
  89. Abe, M.; Kaya, H.; Watanabe-Taneda, A.; Shibuta, M.; Yamaguchi, A.; Sakamoto, T.; Kurata, T.; Ausin, I.; Araki, T.; Alonso-Blanco, C. FE, a phloem-specific Myb-related protein, promotes flowering through transcriptional activation ofFLOWERING LOCUS TandFLOWERING LOCUS T INTERACTING PROTEIN 1. Plant J. 2015, 83, 1059–1068. [Google Scholar] [CrossRef] [Green Version]
  90. Chen, Z.; Rao, P.; Yang, X.; Su, X.; Zhao, T.; Gao, K.; Yang, X.; An, X. A Global View of Transcriptome Dynamics during Male Floral Bud Development in Populus tomentosa. Sci. Rep. 2018, 8, 722. [Google Scholar] [CrossRef] [Green Version]
  91. Haberman, A.; Bakhshian, O.; Cerezo-Medina, S.; Paltiel, J.; Adler, C.; Ben-Ari, G.; Mercado, J.A.; Pliego-Alfaro, F.; Lavee, S.; Samach, A. A possible role for flowering locus T-encoding genes in interpreting environmental and internal cues affecting olive (Olea europaea L.) flower induction. Plant Cell Environ. 2017, 40, 1263–1280. [Google Scholar] [CrossRef]
  92. Koo, B.-H.; Yoo, S.-C.; Park, J.-W.; Kwon, C.-T.; Lee, B.-D.; An, G.; Zhang, Z.; Li, J.; Li, Z.; Paek, N.-C. Natural Variation in OsPRR37 Regulates Heading Date and Contributes to Rice Cultivation at a Wide Range of Latitudes. Mol. Plant 2013, 6, 1877–1888. [Google Scholar] [CrossRef] [Green Version]
  93. Song, S.; Li, L.; Li, Y.; Wang, T.; Fu, X. Overexpression of gene OsSUI1 affects floral organ development in rice (Oryza sativa L.). Mol. Breed. 2018, 38, 4. [Google Scholar] [CrossRef]
  94. Bates, P.D. Understanding the control of acyl flux through the lipid metabolic network of plant oil biosynthesis. Biochim. Biophys. Acta (BBA) Mol. Cell Biol. Lipids 2016, 1861, 1214–1225. [Google Scholar] [CrossRef]
  95. Liu, M.; Li, X.; Liu, Y.; Shi, Y.; Ma, X. Analysis of differentially expressed genes under UV-B radiation in the desert plant Reaumuria soongorica. Gene 2015, 574, 265–272. [Google Scholar] [CrossRef] [PubMed]
  96. Grienenberger, E.; Kim, S.S.; Lallemand, B.; Geoffroy, P.; Heintz, D.; Souza, C.d.A.; Heitz, T.; Douglas, C.J.; Legrand, M. Analysis of TETRAKETIDE α-PYRONE REDUCTASE Function in Arabidopsis thaliana Reveals a Previously Unknown, but Conserved, Biochemical Pathway in Sporopollenin Monomer Biosynthesis. Plant Cell 2010, 22, 4067–4083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Hölscher, C.; Lutterbey, M.-C.; Lansing, H.; Meyer, T.; Fischer, K.; Von Schaewen, A. Defects in Peroxisomal 6-Phosphogluconate Dehydrogenase Isoform PGD2 Prevent Gametophytic Interaction in Arabidopsis thaliana. Plant Physiol. 2016, 171, 192–205. [Google Scholar] [CrossRef] [Green Version]
  98. Cao-Pham, A.H.; Urano, D.; Ross-Elliott, T.J.; Jones, A.M. Nudge-nudge, WNK-WNK (kinases), say no more? New Phytol. 2018, 220, 35–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Yu, X.-M.; Zhao, W.-Q.; Yang, W.-X.; Liu, F.; Chen, J.-P.; Goyer, C.; Liu, D.-Q. Characterization of a Hypersensitive Response-Induced Gene TaHIR3 from Wheat Leaves Infected with Leaf Rust. Plant Mol. Biol. Rep. 2012, 31, 314–322. [Google Scholar] [CrossRef]
  100. Zhou, L.; Cheung, M.-Y.; Li, M.-W.; Fu, Y.; Sun, Z.; Sun, S.S.; Lam, H.-M. Rice Hypersensitive Induced Reaction Protein 1 (OsHIR1) associates with plasma membrane and triggers hypersensitive cell death. BMC Plant Biol. 2010, 10, 290. [Google Scholar] [CrossRef] [Green Version]
  101. Liu, N.; Sun, Y.; Pei, Y.; Zhang, X.; Wang, P.; Li, X.; Li, F.; Hou, Y. A Pectin Methylesterase Inhibitor Enhances Resistance to Verticillium Wilt. Plant Physiol. 2018, 176, 2202–2220. [Google Scholar] [CrossRef] [Green Version]
  102. Lee, H.-N.; Lee, K.-H.; Kim, C.S. Abscisic acid receptor PYRABACTIN RESISTANCE-LIKE 8, PYL8, is involved in glucose response and dark-induced leaf senescence in Arabidopsis. Biochem. Biophys. Res. Commun. 2015, 463, 24–28. [Google Scholar] [CrossRef]
  103. Yang, H.; Zhao, L.; Zhao, S.; Wang, J.; Shi, H. Biochemical and transcriptomic analyses of drought stress responses of LY1306 tobacco strain. Sci. Rep. 2017, 7, 1–10. [Google Scholar] [CrossRef] [Green Version]
  104. Zheng, L.; Chen, S.; Xie, L.; Lu, Z.; Liu, M.; Han, X.; Qiao, G.; Jiang, J.; Zhuo, R.; Qiu, W.; et al. Overexpression of cysteine protease gene from Salix matsudana enhances salt tolerance in transgenic Arabidopsis. Environ. Exp. Bot. 2018, 147, 53–62. [Google Scholar] [CrossRef]
  105. Park, K.Y.; Kim, W.T.; Kim, E.Y. The proper localization of RESPONSIVE TO DESICCATION 20 in lipid droplets depends on their biogenesis induced by STRESS-RELATED PROTEINS in vegetative tissues. Biochem. Biophys. Res. Commun. 2018, 495, 1885–1889. [Google Scholar] [CrossRef] [PubMed]
  106. Wu, Q.; Bai, X.; Zhao, W.; Xiang, D.; Wan, Y.; Yan, J.; Zou, L.; Zhao, G. De Novo Assembly and Analysis of Tartary Buckwheat (Fagopyrum tataricum Garetn.) Transcriptome Discloses Key Regulators Involved in Salt-Stress Response. Genes 2017, 8, 255. [Google Scholar] [CrossRef] [PubMed]
  107. Lin, Z.-J.D.; Liebrand, T.W.; Yadeta, K.A.; Coaker, G. PBL13 Is a Serine/Threonine Protein Kinase That Negatively Regulates Arabidopsis Immune Responses. Plant Physiol. 2015, 169, 2950–2962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Movahedi, A.; Sun, W.; Zhang, J.; Wu, X.; Mousavi, M.; Mohammadi, K.; Yin, T.; Zhuge, Q. RNA-directed DNA methylation in plants. Plant Cell Rep. 2015, 34, 1857–1862. [Google Scholar] [CrossRef]
  109. Dufoo-Hurtado, M.D.; Huerta-Ocampo, J.Ã.; Barrera-Pacheco, A.; Rosa, A.P.B.d.l.; Mercado-Silva, E.M. Low temperature conditioning of garlic (Allium sativum L.) “seed” cloves induces alterations in sprouts proteome. Front. Plant. Sci. 2015, 6, 332. [Google Scholar] [CrossRef] [Green Version]
  110. Xin, S.; Yu, G.; Sun, L.; Qiang, X.; Xu, N.; Cheng, X. Expression of tomato SlTIP2;2 enhances the tolerance to salt stress in the transgenic Arabidopsis and interacts with target proteins. J. Plant. Res. 2014, 127, 695–708. [Google Scholar] [CrossRef]
  111. Chen, X.; Snyder, C.L.; Truksa, M.; Shah, S.; Weselake, R.J. sn-Glycerol-3-phosphate acyltransferases in plants. Plant. Signal. Behav. 2011, 6, 1695–1699. [Google Scholar] [CrossRef] [Green Version]
  112. Bontpart, T.; Cheynier, V.; Ageorges, A.; Terrier, N. BAHD or SCPL acyltransferase? What a dilemma for acylation in the world of plant phenolic compounds. New Phytol. 2015, 208, 695–707. [Google Scholar] [CrossRef]
  113. Lijuan, C.; Huiming, G.; Yi, L.; Hongmei, C. Chalcone synthase EaCHS1 from Eupatorium adenophorum functions in salt stress tolerance in tobacco. Plant. Cell Rep. 2015, 34, 885–894. [Google Scholar] [CrossRef]
  114. Zhang, X.; Abrahan, C.; Colquhoun, T.A.; Liu, C.-J. A Proteolytic Regulator Controlling Chalcone Synthase Stability and Flavonoid Biosynthesis in Arabidopsis. Plant. Cell 2017, 29, 1157–1174. [Google Scholar] [CrossRef] [Green Version]
  115. Reale, S.; Doveri, S.; Díaz, A.; Angiolillo, A.; Lucentini, L.; Pilla, F.; Martín, A.; Donini, P.; Lee, D. SNP-based markers for discriminating olive (Olea europaea L.) cultivars. Genome 2006, 49, 1193–1205. [Google Scholar] [CrossRef] [PubMed]
  116. Kumar, K.; Kumar, S.R.; Dwivedi, V.; Rai, A.; Shukla, A.K.; Shanker, K.; Nagegowda, D.A. Precursor feeding studies and molecular characterization of geraniol synthase establish the limiting role of geraniol in monoterpene indole alkaloid biosynthesis in Catharanthus roseus leaves. Plant. Sci. 2015, 239, 56–66. [Google Scholar] [CrossRef] [PubMed]
  117. Alagna, F.; Geu-Flores, F.; Kries, H.; Panara, F.; Baldoni, L.; O’Connor, S.E.; Osbourn, A. Identification and Characterization of the Iridoid Synthase Involved in Oleuropein Biosynthesis in Olive (Olea europaea) Fruits. J. Biol. Chem. 2016, 291, 5542–5554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Servili, M.; Sordini, B.; Esposto, S.; Taticchi, A.; Urbani, S.; Sebastiani, L. Metabolomics of Olive Fruit: A Focus on the Secondary Metabolites. In The Olive Tree Genome, Compendium of Plant Genomes; Springer Science and Business Media LLC: Cham, Switzerland, 2016; pp. 123–139. [Google Scholar]
  119. Wong, D.C.J.; Pichersky, E.; Peakall, R. The Biosynthesis of Unusual Floral Volatiles and Blends Involved in Orchid Pollination by Deception: Current Progress and Future Prospects. Front. Plant. Sci. 2017, 8, 1955. [Google Scholar] [CrossRef] [Green Version]
  120. Keereetaweep, J.; Liu, H.; Zhai, Z.; Shanklin, J. Biotin Attachment Domain-Containing Proteins Irreversibly Inhibit Acetyl CoA Carboxylase. Plant. Physiol. 2018, 177, 208–215. [Google Scholar] [CrossRef] [Green Version]
  121. Liao, B.; Hao, Y.; Lu, J.-X.; Bai, H.; Guan, L.; Zhang, T. Transcriptomic analysis of Perilla frutescens seed to insight into the biosynthesis and metabolic of unsaturated fatty acids. BMC Genom. 2018, 19, 213. [Google Scholar] [CrossRef] [Green Version]
Genes 11 00916 g001 550
Figure 1. The Principal Coordinate (PCoA), 2- and 3-dimensional graphs, analysis based on a genetic distance matrix between the 260 olive genotypes (A); unrooted neighbor joining tree constructed with MEGA X of 260 olive genotypes (B). Colors representation, light blue—cultivated varieties; purple—oleaster; green—subsp. guanchica.
Figure 1. The Principal Coordinate (PCoA), 2- and 3-dimensional graphs, analysis based on a genetic distance matrix between the 260 olive genotypes (A); unrooted neighbor joining tree constructed with MEGA X of 260 olive genotypes (B). Colors representation, light blue—cultivated varieties; purple—oleaster; green—subsp. guanchica.
Genes 11 00916 g001
Genes 11 00916 g002 550
Figure 2. Estimate of genetic diversity of 260 Olea europaea accessions analyzed by 1040 SNP markers. Bar-plot describes the population structure estimated by the Bayesian clustering. According to the most informative K values, five clusters were distinguished. Each individual was represented by a horizontal colored line, with a length proportional to the estimated membership coefficient. Clusters 1, 2, and 3 represent the three sub-populations of the cultivated varieties: Cluster 1 (light blue)—almost all cultivated varieties from West Mediterranean, together with nine oleasters; cluster 2 (greyish-blue)—includes cultivated varieties mostly from Center Mediterranean; cluster 3 (dark-blue)—almost all cultivated varieties from East Mediterranean; cluster 4, purple—includes 64 out of 73 oleasters (and three cultivated varieties); cluster 5, green—the group corresponding to the 16 samples of subsp. guanchica.
Figure 2. Estimate of genetic diversity of 260 Olea europaea accessions analyzed by 1040 SNP markers. Bar-plot describes the population structure estimated by the Bayesian clustering. According to the most informative K values, five clusters were distinguished. Each individual was represented by a horizontal colored line, with a length proportional to the estimated membership coefficient. Clusters 1, 2, and 3 represent the three sub-populations of the cultivated varieties: Cluster 1 (light blue)—almost all cultivated varieties from West Mediterranean, together with nine oleasters; cluster 2 (greyish-blue)—includes cultivated varieties mostly from Center Mediterranean; cluster 3 (dark-blue)—almost all cultivated varieties from East Mediterranean; cluster 4, purple—includes 64 out of 73 oleasters (and three cultivated varieties); cluster 5, green—the group corresponding to the 16 samples of subsp. guanchica.
Genes 11 00916 g002
Table
Table 1. Genetic differentiation parameters between cultivars, oleasters, and subsp. guanchica samples detected by DnaSP 6.
Table 1. Genetic differentiation parameters between cultivars, oleasters, and subsp. guanchica samples detected by DnaSP 6.
Group 1Group 2KsKxyGstDeltaStGammaStNstFstDxyDa
CultivarsWild olives19.837 ***25.770 ***0.0000.0430.1230.3300.2530.4050.102
Cultivarssubsp. guanchica20.128 ***25.492 ***0.0060.0200.0620.4020.3230.4000.129
Wild olivessubsp. guanchica17.052 ***17.386 ***0.0040.0080.0320.1070.0920.2730.025
Ks and Kxy—genetic differentiation parameters by the PM test; probability obtained by the permutation test with 1000 replicates, *** p < 0.001. Gst—Gene flow parameters, based on haplotype data information; DeltaSt, GammaSt, Nst, and Fst, based on sequence data information. Dxy—average number of nucleotide substitutions per site; Da—number of net nucleotide substitutions per site.
Table
Table 2. Heterozygosity, F-statistics, and polymorphism by population.
Table 2. Heterozygosity, F-statistics, and polymorphism by population.
Group NNaNeHoHeUHeF
CultivarMean170.432.0131.5360.3690.3240.325−0.120
SE0.0310.0040.0090.0050.0040.0040.005
var. sylvestrisMean72.921.9761.4480.2950.2750.276−0.055
SE0.0110.0050.0100.0060.0050.0050.006
subsp. guanchicaMean15.921.7381.3450.2140.2060.212−0.027
SE0.0120.0140.0110.0070.0060.0060.009
TotalMean86.431.9091.4430.2930.2680.271−0.071
SE1.1420.0050.0060.0040.0030.0030.004
SE—Standard Error; N—Sample Size; Na—number of alleles; Ne—effective alleles; Ho—observed heterozygosity; He—expected heterozygosity, UHe—unbiased expected heterozygosity; and F—Fixation Index.

Share and Cite

MDPI and ACS Style

Mariotti, R.; Belaj, A.; De La Rosa, R.; Leòn, L.; Brizioli, F.; Baldoni, L.; Mousavi, S. EST–SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives. Genes 2020, 11, 916. https://doi.org/10.3390/genes11080916

AMA Style

Mariotti R, Belaj A, De La Rosa R, Leòn L, Brizioli F, Baldoni L, Mousavi S. EST–SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives. Genes. 2020; 11(8):916. https://doi.org/10.3390/genes11080916

Chicago/Turabian Style

Mariotti, Roberto, Angjelina Belaj, Raul De La Rosa, Lorenzo Leòn, Federico Brizioli, Luciana Baldoni, and Soraya Mousavi. 2020. "EST–SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives" Genes 11, no. 8: 916. https://doi.org/10.3390/genes11080916

APA Style

Mariotti, R., Belaj, A., De La Rosa, R., Leòn, L., Brizioli, F., Baldoni, L., & Mousavi, S. (2020). EST–SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives. Genes, 11(8), 916. https://doi.org/10.3390/genes11080916

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