Transposon activation is a major driver in the genome evolution of cultivated olive trees (Olea europaea L.)

Genetic markers of olive fruit weight selected to be used in breeding experiments

Olive fruit weight is a crucial trait to consider in olive breeding programs due to its impact on final yield and its relevance for mechanical harvesting and fruit processing. Although environmental conditions influence this trait, fruit weight is primarily determined by genetic factors and exhibits a high degree of heritability in breeding progenies. Despite several studies identifying potential markers associated with fruit weight, these markers have not been validated. In this study, we analyzed 40 genetic markers linked to fruit weight using a dataset comprising 73 cultivars (including 33 newly sequenced varieties) and 10 wild olives with a wide range of phenotypic characteristics, spanning from very light (0.41 g) to very heavy fruits (8.57 g). By examining the phenotype distribution for each genotype of the newly sequenced varieties, we successfully validated 16 genetic markers. Additionally, machine learning tools demonstrated that 9 out of the 16 validated markers have a high predictive ability for fruit weight. As a result, our work provides, for the first time, a set of 9 well-validated genetic markers suitable for use in marker-assisted selection during the early stages of olive breeding programs.

Multi-Omic Advances in Olive Tree (Olea europaea subsp. europaea L.) Under Salinity: Stepping Towards 'Smart Oliviculture'

Soil salinisation is threatening crop sustainability worldwide, mainly due to anthropogenic climate change. Molecular mechanisms developed to counteract salinity have been intensely studied in model plants. Nevertheless, the economically relevant olive tree (Olea europaea subsp. europaea L.), being highly exposed to soil salinisation, deserves a specific review to extract the recent genomic advances that support the known morphological and biochemical mechanisms that make it a relative salt-tolerant crop. A comprehensive list of 98 olive cultivars classified by salt tolerance is provided, together with the list of available olive tree genomes and genes known to be involved in salt response. Na+ and Cl- exclusion in leaves and retention in roots seem to be the most prominent adaptations, but cell wall thickening and antioxidant changes are also required for a tolerant response. Several post-translational modifications of proteins are emerging as key factors, together with microbiota amendments, making treatments with biostimulants and chemical compounds a promising approach to enable cultivation in already salinised soils. Low and high-throughput transcriptomics and metagenomics results obtained from salt-sensitive and -tolerant cultivars, and the future advantages of engineering specific metacaspases involved in programmed cell death and autophagy pathways to rapidly raise salt-tolerant cultivars or rootstocks are also discussed. The overview of bioinformatic tools focused on olive tree, combined with machine learning approaches for studying plant stress from a multi-omics perspective, indicates that the development of salt-tolerant cultivars or rootstocks adapted to soil salinisation is progressing. This could pave the way for 'smart oliviculture', promoting more productive and sustainable practices under salt stress.

First note of QTL mapping of low vigor traits using the updated F2 'Koroneiki' linkage map of olive

The olive tree (Olea europaea L.), which characterizes the agriculture of the Mediterranean basin, faces challenges adapting to high-density orchards and mechanized cultivation. This study addresses a key issue: controlling tree size to enhance efficiency and manageability in olive cultivation. Utilizing genetic mapping methods, we have identified significant Quantitative Trait Loci (QTL) and candidate genes associated with low-vigor traits in olive trees. Our research on the 'Koroneiki' F2 progeny, which exhibits low vigor traits but remains underutilized in breeding programs, has pinpointed a QTL linked to trunk basal diameter-a trait correlated with plant height based on morphological measurements. Results underscore a strong genetic control of these traits, with a consistent correlation observed over time. We identified two candidate genes - Acid Phosphatase 1, Shikimate O-hydroxycinnamoyltransferase, and a SNP Marker likely associated with Calcium Responsive Proteins - each potentially interacting with plant hormones to influence growth. Controlling olive tree size presents several challenges, including the genetic complexity of polygenic traits like size and vigor, and limited rootstock options. By integrating reference genomes with our genetic analysis, we offer a conceptual advancement that could substantially accelerate breeding timelines compared to traditional approaches. Although genome editing is still a future possibility due to the complexity of olive genetics and the species' recalcitrance to transformation, our study lays a foundational understanding to guide future breeding programs. By targeting the identified candidate genes, this research represents a pivotal step toward selecting new low-vigor genotypes and rootstocks, contributing to innovations in olive cultivation.

Circadian- and Light-Driven Rhythmicity of Interconnected Gene Networks in Olive Tree

A circadian clock (CC) has evolved in plants that synchronizes their growth and development with daily and seasonal cycles. A properly functioning circadian clock contributes to increasing plant growth, reproduction, and competitiveness. In plants, continuous light treatment has been a successful approach for obtaining novel knowledge about the circadian clock. The olive tree (Olea europaea L.) is one of the most important crops in the Mediterranean area, and, so far, limited information is available on its CC gene network. Here, we studied the behavior of circadian rhythm genes under LD (light/darkness) and LL (light/light) conditions, the relationships in this network, and the ability of the treatments to modulate gene expression in the photoprotective pigment and lipid biosynthesis pathways. One month of LL conditions increased olive growth performance, but LL exposure also caused reductions in vegetative growth and chlorophyll accumulation. A panel was designed for a study of the transcription expression levels of the genes involved in light perception, the CC, and secondary metabolite and fatty acid biosynthesis. Our results revealed that the levels of 78% of the transcripts exhibited intraday differences under LD conditions, and most of them retained this rhythmicity after exposure to one and two months of LL conditions. Furthermore, co-regulation within a complex network among genes of photoreceptors, anthocyanidins, and fatty acids biosynthesis was orchestrated by the transcription factor HY5. This research enriches our knowledge on olive trees grown under prolonged irradiation, which may be attractive for the scientific community involved in breeding programs for the improvement of this species.

Genome-wide association analysis of flowering date in a collection of cultivated olive tree

Flowering date in perennial fruit trees is an important trait for fruit production. Depending on the winter and spring temperatures, flowering of olive may be advanced, delayed, or even suppressed. Deciphering the genetic control of flowering date is thus key to help selecting cultivars better adapted to the current climate context. Here, we investigated the genetic determinism of full flowering date stage in cultivated olive based on capture sequencing data of 318 genotypes from the worldwide olive germplasm bank of Marrakech, Morocco. The genetic structure of this collection was organized in three clusters that were broadly attributed to eastern, central, and western Mediterranean regions, based on the presumed origin of genotypes. Flowering dates, collected over 7 years, were used to estimate the genotypic best linear unbiased predictors, which were then analyzed in a genome-wide association study. Loci with small effects were significantly associated with the studied trait, by either a single- or a multi-locus approach. The three most robust loci were located on chromosomes 01 and 04, and on a scaffold, and explained 7.1%, 6.2%, and 6.5% of the trait variance, respectively. A significantly higher accuracy in the best linear unbiased predictors of flowering date prediction was reported with Ridge- compared to LASSO-based genomic prediction model. Along with genomic association results, this suggests a complex polygenic determinism of flowering date, as seen in many other fruit perennials. These results and the screening of associated regions for candidate genes open perspectives for further studies and breeding programs targeting flowering date.

Transcriptomic Analysis During Olive Fruit Development and Expression Profiling of Fatty Acid Desaturase Genes

The olive fruit is a drupe whose development and ripening takes several months from flowering to full maturation. During this period, several biochemical and physiological changes occur that affect the skin color, texture, composition, and size of the mesocarp. The final result is a fruit rich in fatty acids, phenolic compounds, tocopherols, pigments, sterols, terpenoids, and other compounds of nutritional interest. In this work, a transcriptomic analysis was performed using flowers (T0) and mesocarp tissue at seven different stages during olive fruit development and ripening (T1-T7) of the 'Picual' cultivar. A total of 1755 genes overexpressed at any time with respect to the flowering stage were further analyzed. These genes were grouped into eight clusters based on their expression profile. The gene enrichment analysis revealed the most relevant biological process of every cluster. Highlighting the important role of hormones at very early stages of fruit development (T1, Cluster 1), whereas genes involved in fatty acid biosynthesis were relevant throughout the fruit developmental process. Hence, genes coding for different fatty acid desaturase (SAD, FAD2, FAD3, FAD4, FAD5, FAD6, and FAD7) enzymes received special attention. In particular, 26 genes coding for different fatty acid desaturase enzymes were identified in the 'Picual' genome, contributing to the improvement of the genome annotation. The expression pattern of these genes during fruit development corroborated their role in determining fatty acid composition.

Integrating genetic analysis of germplasm wealth for enhanced selection and improvement in olive (Olea europaea L.): insights from leaves

KEY MESSAGE

High-throughput next-generation sequencing of 161 olive germplas. 33 samples were selected as core olive germplasm and Fingerprints were constructed. After GWAS analysis of olive leaf shape, 14 candidate genes were localized. Olive (Olea europaea L.) has been introduced to China since the 1960s. After a prolonged period of variation and domestication, there is a lack of comprehensive research on its genetics. The olive oil directly extracted from Olea europaea L. is recognized as 'liquid gold', nevertheless, people constantly overlook the valuable wealth of olive leaves. High-throughput next-generation sequencing was performed on 161 olive germplasm to analyze the kinship, genetic structure and diversity of olives, and the core germplasm of olives were selected and fingerprints were constructed. Meanwhile, Genome-wide association analysis (GWAS) was performed to locate the gene for regulating olive leaf shape. Herein, the results parsed that most of the Chinese olive germplasm was more closely related to the Italian germplasm. A wealth of hybridized germplasm possessed high genetic diversity and had the potential to be used as superior parental material for olive germplasm. A total of 33 samples were selected and characterized as core germplasm of olive and Fingerprints were also constructed. A total of 14 candidate genes were localized after GWAS analysis of four olive leaf shape phenotypes, including leaf shape, leaf curvature shape, leaf tip and leaf base shape. Collectively, this study revealed the genetic basis of olives in China and also succeeded in constructing the core germplasm that stands for the genetic diversity of olives, which can contribute to the scientific and effective collection and preservation of olive germplasm resources, and provide a scientific basis for the in-depth excavation and utilization of genes regulating olive leaf shape.

Evaluation of genetic diversity among olive trees (Olea europaea L.) from Jordan

This study aimed to identify and evaluate the genetic diversity of olive trees in Jordan, a country located in the eastern Mediterranean, where olive domestication originated. For this purpose, a total of 386 olive trees were analyzed, including 338 collected from two surveys (JOCC-1 and JOCC-2) across seven regions, and 48 selected accessions from the Olive Germplasm Bank of Jordan (JGBOC). These trees underwent comprehensive phenotypic and molecular characterization using different tools. Significant differences in morphological traits were detected among tested regions using the Chi-square test. Principal components analysis revealed that fruit color change and growth habit as the most discriminating traits, segregating the trees into two groups, with the first group including the Kanabisi cultivar and the second group including the Kfari Baladi cultivar. Utilizing Kompetitive Allele Specific PCR assay, two sets of informative SNPs were used for the genetic diversity analysis. Cladograms were constructed using the maximum likelihood method, revealing a consistent pattern where two clades containing identical genotypes were observed to cluster with the Kfari Baladi or Kanabisi. In addition, the SNP data was used to perform a comparative analysis with the Worldwide Olive Germplasm Bank of Córdoba, which revealed 73 unreported olive genotypes from Jordan. Genetic structure analyses using Discriminant Analysis of Principal Components (DAPC) identified four clusters with distinctive patterns of relatedness among 149 unique accessions, including 52 olive accessions from various Mediterranean countries (IOCC-3). ADMIXTURE analysis revealed four genetic clusters, consistent with the clustering observed in DAPC and cladogram analysis, indicating a high level of genetic admixture among Jordanian olive germplasm. In conclusion, the results show that olive trees in Jordan are highly diverse, providing valuable information for future conservation and management plans.

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

New insights in the Spanish gene pool of olive (Olea europaea L.) preserved ex situ and in situ based on high-throughput molecular markers

In Spain, several local studies have highlighted the likely presence of unknown olive cultivars distinct from the approximately 260 ones previously described in the literature. Furthermore, recent advancements in identification techniques have significantly enhanced in terms of efficacy and precision. This scenario motivated a new nationwide prospecting effort aimed at recovering and characterizing new cultivated germplasm using high-throughput molecular markers. In the present study, the use of 96 EST-SNP markers allowed the identification of a considerable amount of new material (173 new genotypes) coming from areas with low intensification of production in different regions of Spain. As a result, the number of distinct national genotypes documented in the World Olive Germplasm Bank of IFAPA, Córdoba (WOGBC-ESP046) increased to 427. Likewise, 65 and 24 new synonymy and homonymy cases were identified, respectively. This rise in the number of different national cultivars allowed to deepen the knowledge about the underlying genetic structure. The great genetic variability of Spanish germplasm was confirmed, and a new hot spot of diversity was identified in the northern regions of La Rioja and Aragon. Analysis of the genetic structure showed a clear separation between the germplasm of southern and northern-northeastern Spain and indicated a significantly higher level of admixture in the latter. Given the expansion of modern olive cultivation with only a few cultivars, this cryptic germplasm is in great danger of disappearing. This underlines the fact that maintaining as many cultivars as possible will increase the genetic variability of the olive gene pool to meet the future challenges of olive cultivation.

RNA Polymerases IV and V Are Involved in Olive Fruit Development

Transcription is carried out in most eukaryotes by three multimeric complexes (RNA polymerases I, II and III). However, plants contain two additional RNA polymerases (IV and V), which have evolved from RNA polymerase II. RNA polymerases II, IV and V contain both common and specific subunits that may specialise some of their functions. In this study, we conducted a search for the genes that putatively code for the specific subunits of RNA polymerases IV and V, as well as those corresponding to RNA polymerase II in olive trees. Based on the homology with the genes of Arabidopsis thaliana, we identified 13 genes that putatively code for the specific subunits of polymerases IV and V, and 16 genes that code for the corresponding specific subunits of polymerase II in olives. The transcriptomic analysis by RNA-Seq revealed that the expression of the RNA polymerases IV and V genes was induced during the initial stages of fruit development. Given that RNA polymerases IV and V are involved in the transcription of long non-coding RNAs, we investigated their expression and observed relevant changes in the expression of this type of RNAs. Particularly, the expression of the intergenic and intronic long non-coding RNAs tended to increase in the early steps of fruit development, suggesting their potential role in this process. The positive correlation between the expression of RNA polymerases IV and V subunits and the expression of non-coding RNAs supports the hypothesis that RNA polymerases IV and V may play a role in fruit development through the synthesis of this type of RNAs.

Wild and cultivated olive tree genetic diversity in Greece: a diverse resource in danger of erosion

The genetic relationships between Greek wild olive tree populations and cultivars were investigated. A total of 219 wild genotypes and 67 cultivar genotypes were analyzed by employing 10 SSR markers. Data evidenced that the wild populations exhibited high levels of genetic diversity and exclusively host 40% of the total number of alleles detected. Inbreeding was observed within populations, probably as a consequence of their fragmented spatial distribution. The genetic differentiation between cultivars and wild individuals, as well as within wild populations, was low. Nevertheless, three gene pools of wild trees were detected, corresponding to the geographical areas of Northeastern Greece, Peloponnese-Crete and Epirus. Most cultivars clustered in a separate group, while the rest of them formed a heterogenous group with membership coefficients akin to the three wild olive clusters. Regarding the history of olive cultivation in Greece, bidirectional gene flow was detected between populations of Peloponnese-Crete and the gene pool that composes some of Greece's most important cultivars, such as "Koroneiki" and "Mastoidis", which is inferred as an indication of a minor domestication event in the area. A strategy for the protection of Greek-oriented olive genetic resources is proposed, along with suggestions for the utilization of the genetically diverse wild resources with regard to the introgression of traits of agronomical interest to cultivars.

One hundred single-copy nuclear sequence markers for olive variety identification: a case of fingerprinting database construction in China

Olive is an ancient oil-producing tree, widely cultivated in Mediterranean countries, and now spread to other areas of the world, including China. Recently, several molecular databases were constructed in different countries and platforms for olive identification using simple sequence repeats (SSRs) or single-nucleotide polymorphisms (SNPs). However, comparing their results across laboratories was difficult. Herein, hundreds of polymorphic single-copy nuclear sequence markers were developed from the olive genome. Using the advantage of multiplex PCR amplification and high-throughput sequencing, a fingerprint database was constructed for the majority of olives cultivated in China. We used 100 high-quality sequence loci and estimated the genetic diversity and structure among all these varieties. We found that compared with that based on SSRs, the constructed fingerprint database based on these 100 sequences or a few of them, could provide a reliable olive variety identification platform in China, with high discrimination among different varieties using the principle of BLAST algorithm. An example of such identification platform based on this study was displayed on the web for the olive database in China (http://olivedb.cn/jianding). After resolving redundant genotypes, we identified 126 olive varieties with distinct genotypes in China. These varieties could be divided into two clusters, and it was revealed that the grouping of the varieties has a certain relationship with their origin. Herein, it is concluded that these single-copy orthologous nuclear sequences could be used to construct a universal fingerprint database of olives across different laboratories and platforms inexpensively. Based on such a database, variety identification can be performed easily by any laboratory, which would further facilitate olive breeding and variety exchange globally.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01434-9.

Transposable element evolution in plant genome ecosystems

The relationship of transposable elements (TEs) with their host genomes has usually been seen as an arms race between TEs and their host genomes. Consequently, TEs are supposed to amplify by bursts of transposition, when the TE escapes host surveillance, followed by long periods of TE quiescence and efficient host control. Recent data obtained from an increasing number of assembled plant genomes and resequencing population datasets show that TE dynamics is more complex and varies among TE families and their host genomes. This variation ranges from large genomes that accommodate large TE populations to genomes that are very active in TE elimination, and from inconspicuous elements with very low activity to elements with high transposition and elimination rates. The dynamics of each TE family results from a long history of interaction with the host in a genome populated by many other TE families, very much like an evolving ecosystem.

Whole genome scanning of a Mediterranean basin hotspot collection provides new insights into olive tree biodiversity and biology

Olive tree (Olea europaea L. subsp. europaea var. europaea) is one of the most important species of the Mediterranean region and one of the most ancient species domesticated. The availability of whole genome assemblies and annotations of olive tree cultivars and oleaster (O. europaea subsp. europaea var. sylvestris) has contributed to a better understanding of genetic and genomic differences between olive tree cultivars. However, compared to other plant species there is still a lack of genomic resources for olive tree populations that span the entire Mediterranean region. In the present study we developed the most complete genomic variation map and the most comprehensive catalog/resource of molecular variation to date for 89 olive tree genotypes originating from the entire Mediterranean basin, revealing the genetic diversity of this commercially significant crop tree and explaining the divergence/similarity among different variants. Additionally, the monumental ancient tree 'Throuba Naxos' was studied to characterize the potential origin or routes of olive tree domestication. Several candidate genes known to be associated with key agronomic traits, including olive oil quality and fruit yield, were uncovered by a selective sweep scan to be under selection pressure on all olive tree chromosomes. To further exploit the genomic and phenotypic resources obtained from the current work, genome-wide association analyses were performed for 23 morphological and two agronomic traits. Significant associations were detected for eight traits that provide valuable candidates for fruit tree breeding and for deeper understanding of olive tree biology.

Combined Transcriptomic and Metabolomic Approach Revealed a Relationship between Light Control, Photoprotective Pigments, and Lipid Biosynthesis in Olives

Olive possesses excellent nutritional and economic values for its main healthy products. Among them, a high content of antioxidant compounds, balanced during the ripening process, are produced under genetic and environmental control, resulting in high variability among cultivars. The genes involved in these complex pathways are mainly known, but despite many studies which indicated the key role of light quality and quantity for the synthesis of many metabolites in plants, limited information on these topics is available in olive. We carried out a targeted gene expression profiling in three olive cultivars, Cellina di Nardò, Ruveia, and Salella, which were selected for their contrasting oleic acid and phenolic content. The -omics combined approach revealed a direct correlation between a higher expression of the main flavonoid genes and the high content of these metabolites in 'Cellina di Nardò'. Furthermore, it confirmed the key role of FAD2-2 in the linoleic acid biosynthesis. More interestingly, in all the comparisons, a co-regulation of genes involved in photoperception and circadian clock machinery suggests a key role of light in orchestrating the regulation of these pathways in olive. Therefore, the identified genes in our analyses might represent a useful tool to support olive breeding, although further investigations are needed.

Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast

BACKGROUND

Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors.

RESULTS

We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide consistent estimates of [Formula: see text]50 non-reference TE insertions per strain and that Ty2 has the highest number of non-reference TE insertions in a species-wide panel of [Formula: see text]1000 yeast genomes. Finally, we show that best-in-class predictors for yeast applied to resequencing data have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences revealed previously for experimentally-induced Ty1 insertions to spontaneous insertions for other copia-superfamily retrotransposons in yeast.

CONCLUSION

McClintock ( https://github.com/bergmanlab/mcclintock/ ) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors in other species.

OeBAS and CYP716C67 catalyze the biosynthesis of health-beneficial triterpenoids in olive (Olea europaea) fruits

The bioactive properties of olive (Olea europaea) fruits and olive oil are largely attributed to terpenoid compounds, including diverse triterpenoids such as oleanolic, maslinic and ursolic acids, erythrodiol, and uvaol. They have applications in the agri-food, cosmetics, and pharmaceutical industries. Some key steps involved in the biosynthesis of these compounds are still unknown. Genome mining, biochemical analysis, and trait association studies have been used to identify major gene candidates controlling triterpenoid content of olive fruits. Here, we identify and functionally characterize an oxidosqualene cyclase (OeBAS) required for the production of the major triterpene scaffold β-amyrin, the precursor of erythrodiol, oleanolic and maslinic acids, and a cytochrome P450 (CYP716C67) that mediates 2α oxidation of the oleanane- and ursane-type triterpene scaffolds to produce maslinic and corosolic acids, respectively. To confirm the enzymatic functions of the entire pathway, we have reconstituted the olive biosynthetic pathway for oleanane- and ursane-type triterpenoids in the heterologous host, Nicotiana benthamiana. Finally, we have identified genetic markers associated with oleanolic and maslinic acid fruit content on the chromosomes carrying the OeBAS and CYP716C67 genes. Our results shed light on the biosynthesis of olive triterpenoids and provide new gene targets for germplasm screening and breeding for high triterpenoid content.

Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast

BACKGROUND

Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors.

RESULTS

We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide a consistent and biologically meaningful view of non-reference TE insertions in a species-wide panel of ∼1000 yeast genomes, as evaluated by coverage-based abundance estimates and expected patterns of tRNA promoter targeting. Finally, we show that best-in-class predictors for yeast have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences first revealed experimentally for Ty1 to natural insertions and related copia-superfamily retrotransposons in yeast.

CONCLUSION

McClintock (https://github.com/bergmanlab/mcclintock/) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors for other species.

OliveAtlas: A Gene Expression Atlas Tool for Olea europaea

The olive (Olea europaea L.) is an ancient crop of great importance in the Mediterranean basin due to the production of olive oil and table olives, which are important sources of fat and have benefits for human health. This crop is expanding and increasing its production worldwide and five olive genomes have recently been sequenced, representing a wild olive and important cultivars in terms of olive oil production, intensive agriculture, and adaptation to the East Asian climate. However, few bioinformatic and genomic resources are available to assist olive research and breeding, and there are no platforms to query olive gene expression data. Here, we present OliveAtlas, an interactive gene expression atlas for olive with multiple bioinformatics tools and visualization methods, enabling multiple gene comparison, replicate inspection, gene set enrichment, and data downloading. It contains 70 RNA-seq experiments, organized in 10 data sets representing the main olive plant organs, the pollen germination and pollen tube elongation process, and the response to a collection of biotic and abiotic stresses, among other experimental conditions. OliveAtlas is a web tool based on easyGDB with expression data based on the 'Picual' genome reference and gene annotation.

PEATmoss: A Gene Expression Atlas for Bryophytes

PEATmoss is an interactive gene expression atlas for bryophytes, which originally unified Physcomitrium patens RNA-seq and microarray expression data from multiple gene annotation versions. This atlas includes more than 100 experiments of P. patens, is expanding to host Anthoceros agrestis and Marchantia polymorpha, and aims to host data from more species in the future. PEATmoss has multiple visualization methods and tools for data downloading and is connected to the Physcomitrium patens Gene Model Lookup DB (PpGML DB), which links P. patens genes to annotations and resources from several databases and contains tools for gene version lookup and sequence and annotation extraction. Among the new features available in PEATmoss are dataset privacy control, multispecies menu, interactive color scale, co-expression network visualization, and replicate data downloading.

Identification of Genetic Markers and Genes Putatively Involved in Determining Olive Fruit Weight

The fruit size of a cultivated olive tree is consistently larger than its corresponding wild relatives because fruit size is one of the main traits associated with olive tree domestication. Additionally, large fruit size is one of the main objectives of modern olive breeding programs. However, as the long juvenile period is one main hindrance in classic breeding approaches, obtaining genetic markers associated with this trait is a highly desirable tool. For this reason, GWAS analysis of both genetic markers and the genes associated with fruit size determination, measured as fruit weight, was herein carried out in 50 genotypes, of which 40 corresponded to cultivated and 10 to wild olive trees. As a result, 113 genetic markers were identified, which showed a very high statistically significant correlation with fruit weight variability, p < 10−10. These genetic markers corresponded to 39 clusters of genes in linkage disequilibrium. The analysis of a segregating progeny of the cross of “Frantoio” and “Picual” cultivars allowed us to confirm 10 of the 18 analyzed clusters. The annotation of the genes in each cluster and the expression pattern of the samples taken throughout fruit development by RNAseq enabled us to suggest that some studied genes are involved in olive fruit weight determination.

A review of strategies used to identify transposition events in plant genomes

Transposable elements (TEs) were initially considered redundant and dubbed 'junk DNA'. However, more recently they were recognized as an essential element of genome plasticity. In nature, they frequently become active upon exposition of the host to stress conditions. Even though most transposition events are neutral or even deleterious, occasionally they may happen to be beneficial, resulting in genetic novelty providing better fitness to the host. Hence, TE mobilization may promote adaptability and, in the long run, act as a significant evolutionary force. There are many examples of TE insertions resulting in increased tolerance to stresses or in novel features of crops which are appealing to the consumer. Possibly, TE-driven de novo variability could be utilized for crop improvement. However, in order to systematically study the mechanisms of TE/host interactions, it is necessary to have suitable tools to globally monitor any ongoing TE mobilization. With the development of novel potent technologies, new high-throughput strategies for studying TE dynamics are emerging. Here, we present currently available methods applied to monitor the activity of TEs in plants. We divide them on the basis of their operational principles, the position of target molecules in the process of transposition and their ability to capture real cases of actively transposing elements. Their possible theoretical and practical drawbacks are also discussed. Finally, conceivable strategies and combinations of methods resulting in an improved performance are proposed.

Physiological and Structural Responses of Olive Leaves Related to Tolerance/Susceptibility to Verticillium dahliae

Verticillium wilt of olive (VWO), caused by the soil borne fungus Verticillium dahliae, is one of the most relevant diseases affecting this crop worldwide. One of the best VWO management strategies is the use of tolerant cultivars. Scarce information is available about physiological and structural responses in the leaves of olive cultivars displaying different levels of tolerance to VWO. To identify links between this phenotype and variations in functional characteristics of the leaves, this study examined the structural and physiological traits and the correlations among them in different olive varieties. This evaluation was conducted in the presence/absence of V. dahliae. On the one hand, no leaf trait but the area was related to VWO tolerance in the absence of the pathogen. On the other hand, after inoculation, susceptible cultivars showed lower leaf area and higher leaf mass per area and dry matter content. Furthermore, at the physiological level, these plants showed severe symptoms resembling water stress. Analyzing the relationships among physiological and structural traits revealed differences between tolerant and susceptible cultivars both in the absence and in the presence of V. dahliae. These results showed that olive leaves of VWO-tolerant and VWO-susceptible cultivars adopt different strategies to cope with the pathogen.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

EasyGDB: a low-maintenance and highly customizable system to develop genomics portals

SUMMARY

EasyGDB is an easy-to-implement low-maintenance tool developed to create genomic data management web platforms. It can be used for any species, group of species, or multiple genome or annotation versions. EasyGDB provides a framework to develop a web portal that includes the general information about species, projects and members, and bioinformatics tools such as file downloads, BLAST, genome browser, annotation search, gene expression visualization, annotation and sequence download, and gene ids and orthologs lookup. The code of EasyGDB facilitates data maintenance and update for non-experienced bioinformaticians, using BLAST databases to store and retrieve sequence data in gene annotation pages and bioinformatics tools, and JSON files to customize metadata. EasyGDB is a highly customizable tool. Any section and tool can be enabled or disabled like a switch through a single configuration file. This tool aims to simplify the development of genomics portals in non-model species, providing a modern web style with embedded interactive bioinformatics tools to cover all the common needs derived from genomics projects.

AVAILABILITY AND IMPLEMENTATION

The code and manual to use EasyGDB can be found at https://github.com/noefp/easy_gdb.

Amplification of LTRs of extrachromosomal linear DNAs (ALE-seq) identifies two active Oryco LTR retrotransposons in the rice cultivar Dongjin

Long terminal repeat retrotransposons (LTR-RTs) make up a considerable portion of plant genomes. New insertions of these active LTR-RTs modify gene structures and functions and play an important role in genome evolution. Therefore, identifying active forms of LTR-RTs could uncover the effects of these elements in plants. Extrachromosomal linear DNA (eclDNA) forms during LTR-RT replication; therefore, amplification LTRs of eclDNAs followed by sequencing (ALE-seq) uncover the current transpositional potential of the LTR-RTs. The ALE-seq protocol was validated by identification of Tos17 in callus of Nipponbare cultivar. Here, we identified two active LTR-RTs belonging to the Oryco family on chromosomes 6 and 9 in rice cultivar Dongjin callus based on the ALE-seq technology. Each Oryco family member has paired LTRs with identical sequences and internal domain regions. Comparison of the two LTR-RTs revealed 97% sequence identity in their internal domains and 65% sequence identity in their LTRs. These two putatively active Oryco LTR-RT family members could be used to expand our knowledge of retrotransposition mechanisms and the effects of LTR-RTs on the rice genome.

The Singular Evolution of Olea Genome Structure

The current view of plant genome evolution proposes that genome size has mainly been determined by polyploidisation and amplification/loss of transposons, with a minor role played by other repeated sequences, such as tandem repeats. In cultivated olive (Olea europaea subsp. europaea var. europaea), available data suggest a singular model of genome evolution, in which a massive expansion of tandem-repeated sequences accompanied changes in nuclear architecture. This peculiar scenario highlights the importance of focusing on Olea genus evolution, to shed light on mechanisms that led to its present genomic structure. Next-generation sequencing technologies, bioinformatics and in situ hybridisation were applied to study the genomic structure of five related Olea taxa, which originated at different times from their last common ancestor. On average, repetitive DNA in the Olea taxa ranged from ~59% to ~73% of the total genome, showing remarkable differences in terms of composition. Among repeats, we identified 11 major families of tandem repeats, with different abundances in the analysed taxa, five of which were novel discoveries. Interestingly, overall tandem repeat abundance was inversely correlated to that of retrotransposons. This trend might imply a competition in the proliferation of these repeat classes. Indeed, O. paniculata, the species closest to the Olea common ancestor, showed very few tandem-repeated sequences, while it was rich in long terminal repeat retrotransposons, suggesting that the amplification of tandem repeats occurred after its divergence from the Olea ancestor. Furthermore, some tandem repeats were physically localised in closely related O. europaea subspecies (i.e., cultivated olive and O. europaea subsp. cuspidata), which showed a significant difference in tandem repeats abundance. For 4 tandem repeats families, a similar number of hybridisation signals were observed in both subspecies, apparently indicating that, after their dissemination throughout the olive genome, these tandem repeats families differentially amplified maintaining the same positions in each genome. Overall, our research identified the temporal dynamics shaping genome structure during Olea speciation, which represented a singular model of genome evolution in higher plants.

Utility of EST-SNP Markers for Improving Management and Use of Olive Genetic Resources: A Case Study at the Worldwide Olive Germplasm Bank of Córdoba

Olive, the emblematic Mediterranean fruit crop, owns a great varietal diversity, which is maintained in ex situ field collections, such as the World Olive Germplasm Bank of Córdoba (WOGBC), Spain. Accurate identification of WOGBC, one of the world’s largest collections, is essential for efficient management and use of olive germplasm. The present study is the first report of the use of a core set of 96 EST-SNP markers for the fingerprinting of 1273 accessions from 29 countries, including both field and new acquired accessions. The EST-SNP fingerprinting made possible the accurate identification of 668 different genotypes, including 148 detected among the new acquired accessions. Despite the overall high genetic diversity found at WOGBC, the EST-SNPs also revealed the presence of remarkable redundant germplasm mostly represented by synonymy cases within and between countries. This finding, together with the presence of homonymy cases, may reflect a continuous interchange of olive cultivars, as well as a common and general approach for their naming. The structure analysis revealed a certain geographic clustering of the analysed germplasm. The EST-SNP panel under study provides a powerful and accurate genotyping tool, allowing for the foundation of a common strategy for efficient safeguarding and management of olive genetic resources.

Genome-wide exploration of oil biosynthesis genes in cultivated olive tree varieties (Olea europaea): insights into regulation of oil biosynthesis

Genome-wide oil biosynthesis was explored by de novo sequencing two cultivated olive tree (Olea europaea) varieties (cv. Ayvalik and Picual). This is the first report of the former variety sequencing. As outgroups, raw reads of cv. Leccino and scaffold-level assembly of cv. Farga were also retrieved. Each of these four cultivars was chromosome-scale assembled into 23 pseudochromosomes, with 1.31 Gbp (Farga), 0.93 Gbp (Ayvalik), 0.7 Gbp (Picual), and 0.54 Gbp (Leccino) in size. Ab initio gene finding was performed on these assemblies, using wild olive tree (oleaster)-trained parameters. High numbers of gene models were predicted and anchored to the pseudochromosomes: 69,028 (Ayvalik), 55,073 (Picual), 63,785 (Farga), and 40,449 (Leccino). Using previously reported oil biosynthesis genes from wild olive tree genome project, the following homologous sequences were identified: 1,355 (Ayvalik), 1,269 (Farga), 812 (Leccino), and 774 (Picual). Of these, 358 sequences were commonly shared by all cultivars. Besides, some sequences were cultivar unique: Ayvalik (126), Farga (118), Leccino (46), and Picual (52). These putative sequences were assigned to various GO terms, ranging from lipid metabolism to stress tolerance, from signal transactions to development, and to many others, implicating that oil biosynthesis is synergistically regulated with involvement of various other pathways.

Several Isoforms for Each Subunit Shared by RNA Polymerases are Differentially Expressed in the Cultivated Olive Tree (Olea europaea L.)

Plants contain five nuclear RNA polymerases, with RNA pols IV and V in addition to conserved eukaryotic RNA pols I, II, and III. These transcriptional complexes share five common subunits, which have been extensively analyzed only in yeasts. By taking advantage of the recently published olive tree cultivar (Olea europaea L. cv. Picual) genome, we performed a genome-wide analysis of the genomic composition corresponding to subunits common to RNA pols. The cultivated olive tree genome is quite complex and contains many genes with several copies. We also investigated, for the first time, gene expression patterns for subunits common to RNA pols using RNA-Seq under different economically and biologically relevant conditions for the cultivar "Picual": tissues/organs, biotic and abiotic stresses, and early development from seeds. Our results demonstrated the existence of a multigene family of subunits common to RNA pols, and a variable number of paralogs for each subunit in the olive cultivar "Picual." Furthermore, these isoforms display specific and differentiated expression profiles depending on the isoform and growth conditions, which may be relevant for their role in olive tree biology.

Characterization of Differentially Expressed Genes under Salt Stress in Olive

Climate change, currently taking place worldwide and also in the Mediterranean area, is leading to a reduction in water availability and to groundwater salinization. Olive represents one of the most efficient tree crops to face these scenarios, thanks to its natural ability to tolerate moderate salinity and drought. In the present work, four olive cultivars (Koroneiki, Picual, Royal de Cazorla and Fadak86) were exposed to high salt stress conditions (200 mM of NaCl) in greenhouse, in order to evaluate their tolerance level and to identify key genes involved in salt stress response. Molecular and physiological parameters, as well as plant growth and leaves’ ions Na⁺ and K⁺ content were measured. Results of the physiological measurements showed Royal de Cazorla as the most tolerant cultivar, and Fadak86 and Picual as the most susceptible ones. Ten candidate genes were analyzed and their complete genomic, CDS and protein sequences were identified. The expression analysis of their transcripts through reverse transcriptase quantitative PCR (RT-qPCR) demonstrated that only OeNHX7, OeP5CS, OeRD19A and OePetD were upregulated in tolerant cultivars, thus suggesting their key role in the activation of a salt tolerance mechanism.

Comparative Analysis of Genotyping by Sequencing and Whole-Genome Sequencing Methods in Diversity Studies of Olea europaea L

Olive, Olea europaea L., is a tree of great economic and cultural importance in the Mediterranean basin. Thousands of cultivars have been described, of which around 1200 are conserved in the different olive germplasm banks. The genetic characterisation of these cultivars can be performed in different ways. Whole-genome sequencing (WGS) provides more information than the reduced representation methods such as genotype by sequencing (GBS), but at a much higher cost. This may change as the cost of sequencing continues to drop, but, currently, genotyping hundreds of cultivars using WGS is not a realistic goal for most research groups. Our aim is to systematically compare both methodologies applied to olive genotyping and summarise any possible recommendations for the geneticists and molecular breeders of the olive scientific community. In this work, we used a selection of 24 cultivars from an olive core collection from the World Olive Germplasm Collection of the Andalusian Institute of Agricultural and Fisheries Research and Training (WOGBC), which represent the most of the cultivars present in cultivated fields over the world. Our results show that both methodologies deliver similar results in the context of phylogenetic analysis and popular population genetic analysis methods such as clustering. Furthermore, WGS and GBS datasets from different experiments can be merged in a single dataset to perform these analytical methodologies with proper filtering. We also tested the influence of the different olive reference genomes in this type of analysis, finding that they have almost no effect when estimating genetic relationships. This work represents the first comparative study between both sequencing techniques in olive. Our results demonstrate that the use of GBS is a perfectly viable option for replacing WGS and reducing research costs when the goal of the experiment is to characterise the genetic relationship between different accessions. Besides this, we show that it is possible to combine variants from GBS and WGS datasets, allowing the reuse of publicly available data.

Elucidation of the Origin of the Monumental Olive Tree of Vouves in Crete, Greece

The olive tree of Vouves in Crete, is considered the oldest producing olive tree in the world with an estimated age exceeding 4000 years. In the present study, we sequenced two samples (from the bottom and the top of the tree) to elucidate the genetic relation of this ancient tree with other olive cvs as well as to gain some insights about its origin. Our results showed that both samples have different genetic origins, proving that this ancient tree has been grafted at least one time. On the basis of whole genome sequences the sample from the top of the Vouves tree showed relation of the same order than half-siblings to one accession corresponding to the present-day Greek cv ‘Mastoidis’. Nevertheless, in the framework of a microsatellite analysis it was found to cluster with the ‘Mastoidis’ samples. The Vouves rootstock (bottom sample) showed a clear grouping with the oleaster samples in a similar way to that of ‘Megaritiki’ Greek cv although it does not show any signal of introgression from them. The genomic analyses did not show a strong relation of this sample with the present-day Greek cvs analyzed in this study so it cannot be proved that it has been used as a source for cultivated olive tree populations represented by available genome sequences. Nevertheless, on the basis of microsatellite analyses, the Vouves rootstock showed affinity with two present-day Greek cvs, one “ancient” rootstock from continental Greece as well as monumental trees from Cyprus. The analysis of the impact of the variants in the gene space revealed an enrichment of genes associated to pathways related with carbohydrate and amino acid metabolism. This is in agreement with what has been found before in the sweep regions related with the process of domestication. The absence of oleaster gene flow, its old age and its variant profile, similar to other cultivated populations, makes it an excellent reference point for domestication studies.

K-mer counting and curated libraries drive efficient annotation of repeats in plant genomes

The annotation of repetitive sequences within plant genomes can help in the interpretation of observed phenotypes. Moreover, repeat masking is required for tasks such as whole-genome alignment, promoter analysis, or pangenome exploration. Although homology-based annotation methods are computationally expensive, k-mer strategies for masking are orders of magnitude faster. Here, we benchmarked a two-step approach, where repeats were first called by k-mer counting and then annotated by comparison to curated libraries. This hybrid protocol was tested on 20 plant genomes from Ensembl, with the k-mer-based Repeat Detector (Red) and two repeat libraries (REdat, last updated in 2013, and nrTEplants, curated for this work). Custom libraries produced by RepeatModeler were also tested. We obtained repeated genome fractions that matched those reported in the literature but with shorter repeated elements than those produced directly by sequence homology. Inspection of the masked regions that overlapped genes revealed no preference for specific protein domains. Most Red-masked sequences could be successfully classified by sequence similarity, with the complete protocol taking less than 2 h on a desktop Linux box. A guide to curating your own repeat libraries and the scripts for masking and annotating plant genomes can be obtained at https://github.com/Ensembl/plant-scripts.

How to Choose a Good Marker to Analyze the Olive Germplasm (Olea europaea L.) and Derived Products

The olive tree (Olea europaea L.) is one of the most cultivated crops in the Mediterranean basin. Its economic importance is mainly due to the intense production of table olives and oil. Cultivated varieties are characterized by high morphological and genetic variability and present a large number of synonyms and homonyms. This necessitates the introduction of a rapid and accurate system for varietal identification. In the past, the recognition of olive cultivars was based solely on analysis of the morphological traits, however, these are highly influenced by environmental conditions. Therefore, over the years, several methods based on DNA analysis were developed, allowing a more accurate and reliable varietal identification. This review aims to investigate the evolving history of olive tree characterization approaches, starting from the earlier morphological methods to the latest technologies based on molecular markers, focusing on the main applications of each approach. Furthermore, we discuss the impact of the advent of next generation sequencing and the recent sequencing of the olive genome on the strategies used for the development of new molecular markers.

Coupling the endophytic microbiome with the host transcriptome in olive roots

The connection between olive genetic responses to environmental and agro-climatic conditions and the composition, structure and functioning of host-associated, belowground microbiota has never been studied under the holobiont conceptual framework. Two groups of cultivars growing under the same environmental, pedological and agronomic conditions, and showing highest (AH) and lowest (AL) Actinophytocola relative abundances, were earlier identified. We aimed now to: i) compare the root transcriptome profiles of these two groups harboring significantly different relative abundances in the above-mentioned bacterial genus; ii) examine their rhizosphere and root-endosphere microbiota co-occurrence networks; and iii) connect the root host transcriptome pattern to the composition of the root microbial communities by correlation and co-occurrence network analyses. Significant differences in olive gene expression were found between the two groups. Co-occurrence networks of the root endosphere microbiota were clearly different as well. Pearson's correlation analysis enabled a first portray of the interaction occurring between the root host transcriptome and the endophytic community. To further identify keystone operational taxonomic units (OTUs) and genes, subsequent co-occurrence network analysis showed significant interactions between 32 differentially expressed genes (DEGs) and 19 OTUs. Overall, negative correlation was detected between all upregulated genes in the AH group and all OTUs except of Actinophytocola. While two groups of olive cultivars grown under the same conditions showed significantly different microbial profiles, the most remarkable finding was to unveil a strong correlation between these profiles and the differential gene expression pattern of each group. In conclusion, this study shows a holistic view of the plant-microbiome communication.

Genome Wide MeDIP-Seq Profiling of Wild and Cultivated Olives Trees Suggests DNA Methylation Fingerprint on the Sensory Quality of Olive Oil

Secondary metabolites are particularly important to humans due to their pharmaceutical properties. Moreover, secondary metabolites are key compounds in climate change adaptation in long-living trees. Recently, it has been described that the domestication of Olea subspecies had no major selection signature on coding variants and was mainly related to changes in gene expression. In addition, the phenotypic plasticity in Olea subspecies was linked to the activation of transposable elements in the genes neighboring. Here, we investigated the imprint of DNA methylation in the unassigned fraction of the phenotypic plasticity of the Olea subspecies, using methylated DNA immuno-precipitation sequencing (MeDIP-seq) for a high-resolution genome-wide DNA methylation profiling of leaves and fruits during fruit development in wild and cultivated olives from Turkey. Notably, the methylation profiling showed a differential DNA methylation in secondary metabolism responsible for the sensory quality of olive oil. Here, we highlight for the first time the imprint of DNA methylation in modulating the activity of the Linoleate 9S lipoxygenase in the biosynthesis of volatile aromatic compounds. Unprecedently, the current study reveals the methylation status of the olive genome during fruit ripening.

Olive Fruit Development and Ripening: Break on through to the "-Omics" Side

The olive tree (Olea europaea L. subsp. europaea) is the most important perennial crop in the Mediterranean region, producing table olives and oil, both appreciated for their nutraceutical value. Although olive oil quality traits have been extensively studied, much less attention has been paid to olive drupe. Olive drupe ripening is an extremely complex process involving numerous physiological and molecular changes that are unique in this fruit crop species. This review underlines the contribution of "-omics" techniques and of the recent advances in bioinformatics and analytical tools, notably next-generation sequencing and mass spectrometry, for the characterization of the olive ripening syndrome. The usage of high-dimensional datasets, such as transcriptomics, proteomics, and metabolomics, will provide a systematical description of the molecular-specific processes regulating olive fruit development and ripening. However, the incomplete sequence of the O. europaea L. reference genome has largely hampered the utilization of omics tools towards olive drupe research. Due to this disadvantage, the most reported -omics studies on fruit trees concern metabolomics and only a few transcriptomics and proteomics. In this review, up-to-date applications of -omics technologies towards olive drupe biology are addressed, and future perspectives in olive fruit research are highlighted.

Functional Traits of Olive Varieties and Their Relationship with the Tolerance Level towards Verticillium Wilt

Verticillium wilt of olive (VWO), caused by the soil-borne pathogen Verticillium dahliae, is considered one of the most important diseases affecting this tree crop. One of the best VWO management measures is the use of tolerant cultivars. Remarkably, no information is available about olive functional traits and their potential relationship with tolerance to V. dahliae. Twenty-five selected functional traits (for leaf, stem, root and whole plant) were evaluated in six olive varieties differing in their VWO tolerance level to identify possible links between this phenotype and functional traits’ variation. High intervarietal diversity was found among cultivars and several functional traits were related with VWO tolerance. Tolerant varieties showed higher leaf area, dry matter content (leaf, stem and plant) and mass fraction for stems, but lower for leaves. Significant differences were also detected for root functional traits, tolerant cultivars displaying larger fine root diameter and lignin content but smaller specific length and area of thick and fine roots. Correlations were found among functional traits both within varieties and between levels of tolerance/susceptibility to VWO. Associations were observed between biomass allocation, dry matter content and VWO tolerance. The most relevant difference between tolerant and susceptible cultivars was related to root system architecture.

A Complex Gene Network Mediated by Ethylene Signal Transduction TFs Defines the Flower Induction and Differentiation in Olea europaea L

The olive tree (Olea europaea L.) is a typical Mediterranean crop, important for olive and oil production. The high tendency to bear fruits in an uneven manner, defined as irregular or alternate bearing, results in a significant economic impact for the high losses in olives and oil production. Buds from heavy loaded (‘ON’) and unloaded (‘OFF’) branches of a unique olive tree were collected in July and the next March to compare the transcriptomic profiles and get deep insight into the molecular mechanisms regulating floral induction and differentiation. A wide set of DEGs related to ethylene TFs and to hormonal, sugar, and phenylpropanoid pathways was identified in buds collected from ‘OFF’ branches. These genes could directly and indirectly modulate different pathways, suggesting their key role during the lateral bud transition to flowering stage. Interestingly, several genes related to the flowering process appeared as over-expressed in buds from March ‘OFF’ branches and they could address the buds towards flower differentiation. By this approach, interesting candidate genes related to the switch from vegetative to reproductive stages were detected and analyzed. The functional analysis of these genes will provide tools for developing breeding programs to obtain olive trees characterized by more constant productivity over the years.

Verticillium wilt resistant and susceptible olive cultivars express a very different basal set of genes in roots

BACKGROUND

Olive orchards are threatened by a wide range of pathogens. Of these, Verticillium dahliae has been in the spotlight for its high incidence, the difficulty to control it and the few cultivars that has increased tolerance to the pathogen. Disease resistance not only depends on detection of pathogen invasion and induction of responses by the plant, but also on barriers to avoid the invasion and active resistance mechanisms constitutively expressed in the absence of the pathogen. In a previous work we found that two healthy non-infected plants from cultivars that differ in V. dahliae resistance such as 'Frantoio' (resistant) and 'Picual' (susceptible) had a different root morphology and gene expression pattern. In this work, we have addressed the issue of basal differences in the roots between Resistant and Susceptible cultivars.

RESULTS

The gene expression pattern of roots from 29 olive cultivars with different degree of resistance/susceptibility to V. dahliae was analyzed by RNA-Seq. However, only the Highly Resistant and Extremely Susceptible cultivars showed significant differences in gene expression among various groups of cultivars. A set of 421 genes showing an inverse differential expression level between the Highly Resistant to Extremely Susceptible cultivars was found and analyzed. The main differences involved higher expression of a series of transcription factors and genes involved in processes of molecules importation to nucleus, plant defense genes and lower expression of root growth and development genes in Highly Resistant cultivars, while a reverse pattern in Moderately Susceptible and more pronounced in Extremely Susceptible cultivars were observed.

CONCLUSION

According to the different gene expression patterns, it seems that the roots of the Extremely Susceptible cultivars focus more on growth and development, while some other functions, such as defense against pathogens, have a higher expression level in roots of Highly Resistant cultivars. Therefore, it seems that there are constitutive differences in the roots between Resistant and Susceptible cultivars, and that susceptible roots seem to provide a more suitable environment for the pathogen than the resistant ones.

Genomic evidence for recurrent genetic admixture during the domestication of Mediterranean olive trees (Olea europaea L.)

BACKGROUND

Olive tree (Olea europaea L. subsp. europaea, Oleaceae) has been the most emblematic perennial crop for Mediterranean countries since its domestication around 6000 years ago in the Levant. Two taxonomic varieties are currently recognized: cultivated (var. europaea) and wild (var. sylvestris) trees. However, it remains unclear whether olive cultivars derive from a single initial domestication event followed by secondary diversification, or whether cultivated lineages are the result of more than a single, independent primary domestication event. To shed light into the recent evolution and domestication of the olive tree, here we analyze a group of newly sequenced and available genomes using a phylogenomics and population genomics framework.

RESULTS

We improved the assembly and annotation of the reference genome, newly sequenced the genomes of twelve individuals: ten var. europaea, one var. sylvestris, and one outgroup taxon (subsp. cuspidata)-and assembled a dataset comprising whole genome data from 46 var. europaea and 10 var. sylvestris. Phylogenomic and population structure analyses support a continuous process of olive tree domestication, involving a major domestication event, followed by recurrent independent genetic admixture events with wild populations across the Mediterranean Basin. Cultivated olives exhibit only slightly lower levels of genetic diversity than wild forms, which can be partially explained by the occurrence of a mild population bottleneck 3000-14,000 years ago during the primary domestication period, followed by recurrent introgression from wild populations. Genes associated with stress response and developmental processes were positively selected in cultivars, but we did not find evidence that genes involved in fruit size or oil content were under positive selection. This suggests that complex selective processes other than directional selection of a few genes are in place.

CONCLUSIONS

Altogether, our results suggest that a primary domestication area in the eastern Mediterranean basin was followed by numerous secondary events across most countries of southern Europe and northern Africa, often involving genetic admixture with genetically rich wild populations, particularly from the western Mediterranean Basin.

Genetic Resources of Olea europaea L. in the Garda Trentino Olive Groves Revealed by Ancient Trees Genotyping and Parentage Analysis of Drupe Embryos

The area of the Garda Lake within the Trentino province (north of Italy) is the northernmost part of Europe where the Mediterranean species Olea europaea L. is traditionally cultivated. ‘Casaliva’ is claimed as the main variety traditionally grown in the Garda Trentino area (GT) from which a world renowned niche extra virgin olive oil is produced. Since a dominant presence of ‘Casaliva’ would link the fruit set success and yield to a self-pollination compatibility system, a deep genetic survey of the olive tree population in the GT has been performed with the aim of establishing the actual varietal composition and of understanding from which pollen donor the ‘Casaliva’ olives originate. Forty-four different genetic profiles were observed among the 205 leaf samples collected from 106 ancient trees through the analysis of 20 nuclear microsatellite markers. The varietal composition in modern orchards was also explored and the vast majority of the additional 151 trees analyzed showed the same genotype as the ancient accessions of ‘Casaliva’. The results support the long historical link of ‘Casaliva’ with the GT and, besides a high varietal homogeneity, they also revealed the presence of olive genetic resources essential to fruit production. In fact, the parentage analysis of 550 embryos from drupes of ‘Casaliva’ evidenced that a cross-fertilization system is favored and a list of candidate cultivars most suitable as local pollinizers of ‘Casaliva’ was identified.

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

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.

Genome Wild Analysis and Molecular Understanding of the Aquaporin Diversity in Olive Trees (Olea Europaea L.)

Cellular aquaporin water channels (AQPs) constitute a large family of transmembrane proteins present throughout all kingdoms of life, playing important roles in the uptake of water and many solutes across the membranes. In olive trees, AQP diversity, protein features and their biological functions are still largely unknown. This study focuses on the structure and functional and evolution diversity of AQP subfamilies in two olive trees, the wild species Olea europaea var. sylvestris (OeuAQPs) and the domesticated species Olea europaea cv. Picual (OleurAQPs), and describes their involvement in different physiological processes of early plantlet development and in biotic and abiotic stress tolerance in the domesticated species. A scan of genomes from the wild and domesticated olive species revealed the presence of 52 and 79 genes encoding full-length AQP sequences, respectively. Cross-genera phylogenetic analysis with orthologous clustered OleaAQPs into five established subfamilies: PIP, TIP, NIP, SIP, and XIP. Subsequently, gene structures, protein motifs, substrate specificities and cellular localizations of the full length OleaAQPs were predicted. Functional prediction based on the NPA motif, ar/R selectivity filter, Froger's and specificity-determining positions suggested differences in substrate specificities of Olea AQPs. Expression analysis of the OleurAQP genes indicates that some genes are tissue-specific, whereas few others show differential expressions at different developmental stages and in response to various biotic and abiotic stresses. The current study presents the first detailed genome-wide analysis of the AQP gene family in olive trees and it provides valuable information for further functional analysis to infer the role of AQP in the adaptation of olive trees in diverse environmental conditions in order to help the genetic improvement of domesticated olive trees.

Gene Expression Pattern in Olive Tree Organs (Olea europaea L.)

The olive tree (Olea europaea L.) was one of the first plant species in history to be domesticated. Throughout olive domestication, gene expression has undergone drastic changes that may affect tissue/organ-specific genes. This is an RNA-seq study of the transcriptomic activity of different tissues/organs from adult olive tree cv. "Picual" under field conditions. This analysis unveiled 53,456 genes with expression in at least one tissue, 32,030 of which were expressed in all organs and 19,575 were found to be potential housekeeping genes. In addition, the specific expression pattern in each plant part was studied. The flower was clearly the organ with the most exclusively expressed genes, 3529, many of which were involved in reproduction. Many of these organ-specific genes are generally involved in regulatory activities and have a nuclear protein localization, except for leaves, where there are also many genes with a plastid localization. This was also observed in stems to a lesser extent. Moreover, pathogen defense and immunity pathways were highly represented in roots. These data show a complex pattern of gene expression in different organs, and provide relevant data about housekeeping and organ-specific genes in cultivated olive.