1. Introduction
In 2004, Prof. Jesús V. Jorrín Novo (Agroforestry and Plant Biochemistry, Proteomics and Systems Biology Research Group) and Prof. Rafael M. Navarro-Cerrillo (Evaluation and Restoration of Agricultural and Forestry Systems Research Group), both at the University of Cordoba (Spain), started a pioneering joint project devoted to studying the molecular biology of the holm oak (Quercus ilex L.), as well as other forest tree species such as Pinus. This collaboration was born out of a double aim. Firstly, to understand different aspects of their biology, focusing on developmental processes (seed maturation and germination) and responses to environmental biotic and abiotic stresses using a molecular language throughout the backward flow of the central dogma of molecular biology, from phenotype towards genotype. Secondly, to decipher the molecular basis of the enormous phenotypic variability observed in Q. ilex, as well as the genetic structure of different populations correlated to edaphoclimatic conditions. However, we soon became aware of the project’s ambition, considering that the latest molecular techniques (e.g., modern omics techniques), used in biomedical research and, to a lesser extent, in plant model experimental systems such as Arabidopsis thaliana, were not being applied to this species. In general, the biological characteristics of forest trees have made them a challenging experimental system to work with, for example, due to their longevity, long regeneration period, allogamy, high phenotypic variability, and seed recalcitrance. For this reason, field experiments are difficult to carry out, making it necessary to include many individuals to reach consistent and reliable conclusions with biological significance [1]. The plant material selected for the experiments is young seedlings grown under greenhouse-controlled conditions, preferentially sampling leaves or roots. The experimental design must be scheduled on a yearly basis as, in the case of Q. ilex, seeds are collected in November–December. The recalcitrant character of the seeds makes it necessary to germinate them immediately after collection and perform the experiments by June–July without the possibility of repeating them until the next acorn harvest, one year later. Moreover, from our experience, standard and routine protocols and commercial kits successfully used in other plant systems, e.g., A. thaliana, do not work with Q. ilex, which requires an extra effort to optimize each step of the methodology (e.g., Echevarria-Zomeño et al. [2]).
The first pioneers’ biochemical papers on Q. ilex were published from the mid-1990s onwards [3,4,5]. From then until the 2020s, there were hardly any papers published on Q. ilex’s biochemistry and molecular biology, and the number devoted to the genus Quercus was quite low, which means that Q. ilex is still, like most Quercus spp., an orphan in molecular studies. Table 1 shows the number of publications found in the Web of Science (WOS) database at the time of writing this review, reporting on the use of omics approaches in Quercus, Q. ilex, and other organisms including several plant species and forest trees. The numbers illustrate that molecular research on forest tree species, especially Q. ilex, is minimal compared to that carried out with model and herbaceous plants, particularly compared to the number of papers on humans. Some reviews on forest tree species and Quercus have previously been published by the authors’ research group [6,7,8].
2. Quercus ilex, the Dominant Forest Tree Species in the Mediterranean Area (Why)
The genus Quercus (family Fagaceae, common name oak) emerged during global warming 55 million years ago [9], and is one of the largest genera of flowering woody plants, evergreen or deciduous, that produces acorn fruits; it dominates so-called oak forests and woodlands, natural or modified ecosystems, and constitutes a source of biodiversity [10]. The genus Quercus is one of the most relevant clades in the Northern Hemisphere in terms of number of species and ecological dominance [11]; in fact, the Plants of the World Online database [12] contains 460 Quercus spp. distributed throughout the Northern Hemisphere to Malaysia and Colombia.
Oaks have played an important role in human life since the prehistoric period and throughout different civilizations, being a source of a wide variety of goods and services, from wood for fuel and construction to food for human and animal consumption [13]. Despite being one of the most important woody species in terms of diversity, ecological dominance, and economic value [11], their molecular biology is scarcely known, and oaks can be considered, in general terms, an orphan of molecular studies, with some exceptions. Among those, it is worth mentioning Quercus robur L., considered the model within the genus, Quercus lobata Née, Quercus suber L. and Q. ilex, with this last species constituting our object of study for the past 20 years [8].
Within the genus, Q. ilex (Figure 1) is the dominant tree species in the Mediterranean forest and the agrosilvopastoral ecosystem “dehesa” (Spain) or “montado” (Portugal) (Figure S1) [14]. Hence, in Spain, out of the 15 Mha of forested area (50 Mha total Spanish land surface), 30% (4.5 Mha) corresponds to Q. ilex, mostly in “dehesas” forests (3,1 Mha; [15,16,17]). The most profitable product of Q. ilex is the acorns, which constitute the food of the Iberian pig, from whom high-quality meat products are created [18,19]. However, the true relevance of the Q. ilex plantations, not translatable into economic parameters, is related to environmental aspects such as carbon balance, water resources, biodiversity and, ultimately, desertification prevention [20].
Currently, there is a renewed interest in using orphan crops, including Q. ilex and other forest trees, for dietary diversification, sustainable food production and, ultimately, to improve nutrition, which is a challenge, especially in marginal areas [21,22,23]. Q. ilex fruits are currently undervalued and underexploited compared to other equivalent nuts such as the chestnut (Castanea sativa Mill.) and walnut (Juglans regia L.). Therefore, alternative economic uses are being considered for acorns, including as food [23], beverages [24], biofuel [25], pharmaceuticals, and in other industrial applications [26].
3. Biotechnology and Molecular-Assisted Breeding Programs for Resilience and Nutritional Purposes (What for)
Beyond the economic benefits, Q. ilex, like other forest tree species and their associated natural or transformed ecosystems, represent intangible wealth, a legacy we owe it to future generations to preserve. Consequently, forest management and conservation nowadays constitute a worldwide priority that should be reflected in international and national policies, integrating economic and environmental issues to generate richness and welfare [27]. The reality, however, is far from the theory, and a worrisome loss of forest spaces and increasing tree mortality has been occurring in recent decades [28], a problem that may worsen in a climate change scenario [29]. Q. ilex is a good example considering the increase in mortality observed over the last decade [30]. This is associated with anthropogenic factors and a combination of biotic and abiotic stresses. So-called oak decline syndrome is a result of combined stresses. In Q. ilex, the combination of drought episodes with the presence of the oomycete Phytophthora cinnamomi Rands is the main cause of the decline [31,32,33].
Faced with this concerning reality, it is imperative to orchestrate effective measures to ensure Q. ilex conservation, sustainable management, and successful reforestation to ultimately safeguard its ecosystem. Modern biotechnology, based on molecular knowledge, provides solutions that can be integrated with other chemical and biological alternatives [8]. The biological characteristics of Q. ilex preclude the biotechnological strategies successfully used in crops, such as genetic crosses, genetic engineering, and the search for agrochemicals. Thus, a unique feasible approach is based on the exploitation of biodiversity to select elite genotypes that will be used for clonal propagation later [34,35,36]. In fact, as discussed by Sork [37], the enormous phenotypic variability found in Q. ilex populations could be the result of its adaptation to environmental conditions, favoring the appearance of resilient phenotypes. It is, therefore, necessary to conduct the morphological, physiological, and molecular characterization of such phenotypic variability using molecular marker techniques. This will allow for the selection of elite genotypes in terms of higher germination rates and acorn production, desirable quality traits related to nutritional values, or adaptation and resilience to adverse biotic and abiotic stresses [7,8]. To date, the number of scientific works dealing with the molecular aspects of Q. ilex is rather limited due to its intrinsic difficulty as an experimental system (Table 1; [8]). This includes the high variability found in natural populations of Q. ilex, while Populus and Eucalyptus studies can use clones with sequenced genomes.
Despite the challenges it entails, we are addressing the molecular characterization of biodiversity in Q. ilex following the central dogma of molecular biology and using the omics cascade, which explains the flow of genetic information from DNA (genomics, epigenomics, and DNA-based markers) to metabolites (metabolomics), through mRNA (transcriptomics) and proteins (proteomics).
The use of state-of-the-art genomics and functional genomic approaches in model plants (e.g., A. thaliana) and crops (e.g., Oryza sativa L. and Triticum aestivum L.) has contributed significantly to our knowledge of the genetic underpinnings of phenotypes of interest, such as those related to acclimatation and adaptation to a changing environment, as well as to the selection of candidate marker genes to be used in breeding programs. In this direction, the research carried out by Prof. V. Sork on Q. lobata [9] and by Prof. C. Plomión in Q. robur [38] offers excellent roadmaps to follow in other Quercus spp., identifying geographic patterns of phenotypic and genotypic variations and performing complementary studies on pollen and seed movement, genomics, and epigenomics.
4. The “Omics” and Their Integration with Other Approaches (How)
During their 20 years of Q. ilex research, the authors have optimized many techniques and protocols, using both in silico and wet lab experiments, from classical biochemistry to the latest omics techniques. These customized procedures have been utilized in the characterization of mechanisms and the identification of genes and gene products associated with elite genotypes, more resilient to adverse environmental conditions (e.g., decline syndrome) and with higher nutritional value. Omics molecular data have been integrated with morphometric, physiological, and classical biochemical results, and with those of edaphoclimatic conditions, to correlate gene and gene products with specific phenotypes and environmental conditions, following a systems biology approach. In the case of acorns, metabolomics and proteomics data have been correlated with those of classical biochemistry (starch, sugars, amino acids, phenolics and flavonoids determination) and NIRS (near-infrared spectroscopy) [34,36,39]. The UHPLC-QToF (ultra-high performance liquid chromatography–quadrupole time-of-flight mass spectrometry) technique was used for untargeted metabolomic analysis. A total of 192 metabolites were annotated. Several acorn morphotypes were discriminated by principal component analysis (PCA), identifying 50 compounds out of 192 annotated with the highest load over the first two PCA components (explaining 67.2% of variability). These compounds could be considered potential markers of variability in Q. ilex acorns [36]. When determining the alimentary use of acorns, as nuts or for flour-derived products, sweetness is a desirable property [40] which is related to the water content and the chemical composition. The starch and sugar content is correlated to the sweetness, while tannins and other phenolics contribute to the astringent or bitter characteristics [41,42]. To date, we have not found a correlation between acorn size and chemical composition. Regarding seeds, they are of particular interest for studies aimed at unravelling the molecular bases of acorn recalcitrance and the differences between viable and non-viable seeds. In this matter, protein signatures for the different stages of the maturation and germination phases have been defined [43,44,45]. The interpretation of the omics data concerning the study of seedlings’ response to individual (drought or P. cinnamomi) or combined (drought and P. cinnamomi) stresses must take into consideration the physiological state of the seedlings sampled: mortality and damaged seedlings, water content and hydric regime, and photosynthesis-related parameters, as determined by a FluorPen FP100 portable leaf fluorimeter from Photon Systems Instruments (Drasiv, Czech Republic), a portable IR CO2 gas analyzer (IRGA) equipped with a light source integrated in a leaf chamber fluorometer and CO2 injector system (LiCor Li6400XT, LiCor, Inc.; Lincoln, NE, USA), and by pigment quantitation [33,35,46,47,48].
The following subsections are devoted to different DNA markers and omics techniques.
4.1. DNA Molecular Markers
A genetic marker is a sequence of DNA, present in specific genes, genotypes, populations, or species, that can be used in the identification of biological samples through hybridization, restriction mapping, or PCR techniques. Its high level of polymorphism, simple inheritance, and occurrence along the entire genome makes it preferable to morphological or biochemical markers. Practical applications in the plant field include studies of traceability, botanical systematics, genetic structure, divergence, phylogenetic relationships, and evolution. There are a wide variety of genetic markers that can be classified according to the technique they use and are known by the following acronyms: RFLP (restriction fragment length polymorphism), RAPD (random amplification of polymorphic DNA), AFLP (amplified fragment length polymorphism), SCAR (sequence characterized amplified region), SNP (single-nucleotide polymorphism), and SSR (single sequence repeat). Scotti-Saintagne et al. [49] compared different molecular markers (isozymes, AFLPs, SCARs, SSRs, and SNPs) to evaluate their usefulness for differentiating between two closely related oak species (Q. robur L. and Quercus petraea (Matt.) Liebl.).
DNA-based markers were employed in a study of a large number of species within the genus Quercus, as reviewed by Simeone et al. [50], Backs and Ashley [51], Lazic et al. [51], and Liang et al. [52], and are summarized in Table S1. In general, data from DNA marker techniques clearly show a high level of genetic variation within the genus, especially in Q. ilex [53]. Several studies have reported on Q. ilex × Q. suber introgression events, more frequently in the direction ilex (female) × suber (male), which favors the immigration of Q. suber beyond its main range, in regions already colonized by Q. ilex [54]. Lopez de Heredia et al. [55]) analyzed, by ddRAD sequencing, adult hybrids (sensu lato) in a mixed stand and their open-pollinated progenies, estimating introgression levels in adults and seedlings, and concluded that the contribution of Q. suber to the genome is higher in progenies than in adults, which suggests preferential backcrossing with this parental species. Despite intraspecific ploidy variation within the genus Quercus being uncommon, Dzialuk et al. [56] have reported on the use of microsatellites to verify triploidy in Q. petraea and Q. robur. As an anecdotal curiosity, SSRs have been used in forensic botany, allowing for the identification of a leaf sample taken from a suspect’s vehicle as Quercus and placing the individual at or near a crime scene [57].
Among the different types of DNA markers available, our research group has used SSR markers for genetic structure and diversity studies in several populations of Q. ilex [58]. This type of marker is defined as having short tandemly repeated sequences, from one to six bases, distributed along the nuclear, mitochondrial, and chloroplast genomes [59]. Unlike other primer-based marker techniques, SSRs exploit the high level of homology occurring in the flanking regions of the SSRs among taxonomically related species. This feature enables cross-transferability—namely, markers developed for one species can be used for its relatives, which is especially important in species with little-studied genomes [60]. In fact, using SSRs, pure individuals from hybrids between Q. ilex and other species of the genus Quercus have been identified, such as Q. ilex–Q. suber [61,62,63] and Q. ilex–Quercus coccifera L. [64], and the analysis of the genetic structure and diversity of Spanish populations has been carried out [58,65]. We analyzed the genetic structure of several Q. ilex Andalusian populations using a battery of 10 chloroplast SSRs to examine maternal inheritance, and 10 nuclear SSRs to examine chromosomal recombination [58]. High levels of spatial differentiation were found for chloroplast DNA, indicating little seed dispersal between populations. As for nuclear SSRs, these exhibited much weaker spatial organization (differentiation between populations) than chloroplast DNA. The Betic Cordillera (Cádiz) population appeared to be consistently well separated from the northern Sierra Morena populations, suggesting that the Guadalquivir Valley has played an important role in determining population divergence (Figure S1). Currently, this battery of SSRs is being successfully used to characterize symptomatic and asymptomatic Q. ilex and Q. suber individuals located in areas of Spain affected by oak decline (Figure 2). So far, we have observed a possible differentiation between symptomatic and asymptomatic individuals analyzed in those areas (unpublished data). In a similar study developed in Q. ilex, PCR-RFLP and SSR markers show a clear geographic structure of genetic diversity [66].
4.2. Genomics
In the broadest and most general sense, genomics encompasses the structural and functional analysis of the complete genome of an organism, including mapping, sequencing, and gene expression profiling. At the time of writing, the number of land plant (Embryophyta) species whose reference genome has been made publicly available on GenBank was 897, with the figures for the Magnoliopsida (flowering plants) phylum, Fagales order, Fagaceae family, and Quercus spp. being 855, 20, 12, and 8, respectively. The value for forest tree species remains anecdotal compared with that of herbaceous models and crop species. By 2018, the first draft genomes of Quercus were released—those of Q. robur [38,67], and Q. lobata [9,68]— followed a few years later by Q. suber [69], Q. mongolica Fisch. ex Ledeb. [70], Quercus aquifolioides Rehder and E.H.Wilson, Quercus wislizeni A. DC., Quercus gilva Blume, and Quercus glauca (Thunb.) (Table 2).
As early as 1998, the genome sizes of four deciduous and three evergreen Quercus were reported, with values ranging from 1.88 to 2.00 pg/2C, revealing an average GC homogeneous content of 39.9%. In the case of Q. petraea, the DNA content differed between two populations, as confirmed by cytogenetic observations, displaying extra chromosomes in some individuals [71]. The characteristics and assembly statistics of the Quercus spp. whose genome has been sequenced are summarized in Table 2.
The achievement of genome sequencing and the analysis of gene sequences of forest species open the way for debate, speculation, and the formulation of hypotheses about the origin of mutations and evolutionary changes in long-lived organisms, how these contribute to coping with changes, and how to predict dynamics in a climate change scenario [38,72]. In this context, Sork et al. [9], while performing a high-quality de novo genome assembly of Q. lobata, traced its evolutionary history; interestingly, the analysis of 39,373 mapped protein-coding genes revealed duplications consistent with the genetic and phenotypic diversity underlying oak evolutionary success. Fitzek et al. [73] performed novel work with DNA barcoding to study introgression. By using restriction site-associated DNA sequencing, an SNP barcoding system was developed to characterize sympatric eastern North American white oak species and hybrids.
Our group is working on the first draft of the Q. ilex genome, using DNA extracted from leaves of an adult individual located in Aldea de Cuenca, Fuente Obejuna, Córdoba (coordinates: 38°3′95″ N; 5°55′40″ E) (Figure 1). For this purpose, the single-molecule real-time (SMRT) technology of the PacBio platform was used in collaboration with the University of Delaware. The genome dataset has been submitted to the NCBI BioProject repository (ID: 687489) and will be available soon. Previously, the genome size of Q. ilex was estimated by flow cytometry at 1.94 Gb/2C, and its 12 homologous chromosome pairs were identified using fluorescence microscopy (2n = 2x = 24) [6]. The preliminary assembled genome has a total length of 842 Mb.
4.3. Epigenomics
Epigenomics studies the stable inheritance of phenotypes resulting from changes in gene expression without alterations in the DNA sequence. Such changes are involved in the regulation of many biological processes and are mediated by different epigenetic mechanisms (DNA methylation, histone modification and micro-RNAs) [74]. To approach their study, cutting-edge platforms, techniques, and protocols have been continuously improved over the last decade [75]. Epigenomics is an emerging field of research in forest species that will increase our understanding of phenotypic plasticity and adaptative responses to environmental cues, including those related to climate change [76]. Despite scientific evidence supporting the involvement of epigenetic modifications in adaptation and responses to the environment, it is unknown to what extent genetic and epigenetic factors determine important traits related to fitness under these circumstances [77]. Such knowledge can be translated into biotechnology to address the worrisome increase in tree mortality observed in recent decades [78].
When we conducted a search on WOS with the keywords Quercus and epigenetics, 16 results were found (Table S2). It is worth noting the work led by V. Sork’s group on Q. lobata, which can be considered a model system for the genus Quercus in epigenomics. They investigated how adaptive evolution occurs not only through genetic mutations but also through epigenetic mechanisms, such as DNA methylation [79]. Both genetic and epigenetic patterns showed local adaptation suggesting that CpG methyl polymorphisms are involved in the adaptation.
To our knowledge, only one paper, by Rico et al. [80], deals with the study of methylation marks in Q. ilex using the MSAP (methylation-sensitive amplification polymorphism) technique. The authors investigated the epigenetic modifications that may contribute to drought acclimation in natural forest trees, resulting in an increase in hypermethylated loci and a decrease in fully methylated plants exposed to drought. As a first approach to tackling epigenomics in Q. ilex, we are analyzing the global DNA methylation patterns at three developmental stages (leaves of adult trees and seedlings as well as embryos) from the Q. ilex individual used in whole-genome sequencing (Figure 1). We are employing a RAD-Seq variant and the MSAP-Seq (methylation-sensitive amplified polymorphism sequencing) technique, which has previously been described for Hordeum vulgare [81] and Populus alba [82] (Figure 3). Preliminarily, a higher degree of cytosine methylation was detected in adults (58%), followed by seedlings (33%) and embryos (9%) (unpublished data). Our results display a similar pattern to that reported for Q. ilex by Rico et al. [80], with 190 loci: 69% unmethylated, 22% full methylated, and 8% hemimethylated. Multivariate analyses conducted to explore the methylation status at MSAP loci also showed a clear distinction between developmental stages and organs. Regardless of the low number of loci commonly observed between developmental stages, some of them were present in trees and seedlings. To identify those putative differentially methylated genes expressed at the three developmental stages and two organs, the MSAP pre-PCR products obtained are being sequenced by the Illumina Inc. (San Diego, CA, USA) platform.
4.4. Transcriptomics
Transcriptomics deals with the study of the transcriptome, namely, the total set of RNA transcripts present in any cell or biological system at specific developmental and physiological stages and under specific environmental conditions. This definition reflects the dynamic character of the transcriptome, as changes in the RNA profiles determine the functional and phenotypic characteristics of the biological systems. The first reference transcriptomes for the genus Quercus were those of Q. robur and Q. petraea, obtained by Sanger and pyrosequencing technology [83]. Currently, there are more than 300,000 plant RNA sequencing data from organs, tissues, cells, developmental stages, and treatments available in different databases (reviewed by Lim et al. [84]). The platforms and methodologies used to obtain these data have evolved over time, from microarrays to what is now most common: RNA-seq [85,86]. In any case, these untargeted omics data must be validated by quantitative real-time PCR. However, data validation by quantitative real-time PCR data must consider the occurrence of alternative transcripts [87]. One of the key methodological issues at this point is the selection of suitable reference genes whose expression profile does not vary in response to the process under study [88]. In this context, we selected and evaluated 12 candidate reference genes not differentially expressed in Q. ilex seedlings in response to drought stress [89]. From the data obtained, and using the Bestkeeper, GeNorm, NormFinder and ΔCt comparative method algorithms, actin, GAPDH, and β-tubulin were selected as the most steady and reliable candidate reference genes (unpublished data).
When we conducted a search on WOS with the keywords of plants and transcriptomics, it resulted in 6680 items, with around 4% (238 items) corresponding to forest species. When the search was restricted to Quercus, a total of 28 items were displayed, with Q. robur, Q. suber, Q. lobata, and Q. ilex being the most represented species, agreeing with what was stated in the previous section and shown in Table S3.
The priority in current research into Quercus spp. is the response to stresses and climate change conditions, and this has been approached using a comparative transcriptomic approach in different species [90,91,92,93]. A comparison of the least drought-tolerant deciduous Q. robur, quite tolerant deciduous Q. pubescens, and the most tolerant evergreen Q. ilex revealed common and species-specific responses [94]. Oney-Birol et al. [95] carried out a comparative transcriptome analysis in response to drought in different Quercus spp., identifying candidate genes that can contribute to adaptative divergence among hybridizing species. They found 398,042 variants across loci and 4532 variants in 139 candidate genes.
Transcriptomic analyses have also been very useful in shedding light on beneficial symbiotic and cooperative relationships, such as the interactions of Quercus–Mycorrhizas and Quercus–gall-forming insects, respectively [96,97]. The latter is an interesting example of cooperation between organisms, in which the gall-forming insects modify the plant’s genome to redirect the host plant tissue; there have been changes reported in 28% of the whole transcriptome between galled and ungalled tissue [98]. One of the changes observed in the inner larval capsule is the suppression of genes of the plant immune system by utilizing the host RNA interference mechanism. Furthermore, transcriptomics has been applied, together with DNA marker techniques (see Section 4.1), in the study of introgression between different Quercus spp. [95] and to unravel evolutionary aspects. Cokus et al. [99] developed a de novo assembled and annotated reference transcriptome for Q. lobata and Quercus garryana Douglas ex Hook, and identified SNPs, hypothesizing on the divergent selection for pathogen resistance election, which might explain ecological differences between species. Other aspects of marked practical interest studied through transcriptomics are Q. suber phellogen related to the quality of cork [100,101] and the molecular basis of seed recalcitrance [102].
The transcriptomic research carried out by our group in Q. ilex has mainly focused on the generation of RNA-Seq datasets to explore gene expression in several biological contexts, such as those related to environmental stresses and plant–pathogen interactions. The availability of such RNA-seq datasets has allowed for significant advances in the analysis of the stress responses in Q. ilex. The first transcriptome was generated from a mixture of tissues (leaf, root, and embryo) using the Illumina sequencing platform [103] and was further improved by a second one using the Ion Torrent platform [104]. A high number of transcripts—5405—were associated with stress response, of which 46 were involved in drought stress (Figure 4). Recently, the group has published a paper focused on the transcriptome (Illumina platform) and proteome (LC-MS/MS) profiles of six-month-old seedlings subjected to severe drought condition [89]. The quality and confidence of the mRNA and protein identifications and quantifications were assessed, identifying a total of 25,169 transcripts and 3312 proteins. Despite the poor correlation observed between mRNA and protein, four gene products (viz., FtSH6, CLPB1, CLPB3, and HSP22) were upregulated at both omics levels and, hence, were proposed as potential drought-tolerance markers to be used in the selection of elite, resilient genotypes and breeding programs.
Beyond the global RNA-seq analysis, the first transcriptomic analysis reported by our group dealt with the recalcitrant nature of Q. ilex acorns using a targeted transcriptomic approach combined with a hormonal and sugar analysis in mature seeds, germinated seeds, and seedlings, in an attempt to understand the physiological changes that take place during seed germination and seedling growth in Q. ilex [105]. Eleven genes were selected for the comparative analysis profile encoding proteins involved in desiccation tolerance (DHN3 and GOLS), the regulation of ABA (abscisic acid) signaling (OCP3, SKP1 and SDIR1), metabolism (FDH, GADPH, RBLC and NADH6), and oxidative stress (SOD1 and GST). The transcript abundance of ABA-related genes (OCP3, SKP1, and SDIR1) was positively correlated with changes in phytohormone content (low levels of ABA were observed in mature seeds). Transcripts of drought-related genes (DHN3 and GOLS), metabolism-related genes (GAPDH and NADH6) and oxidative stress-related genes (SOD1 and GST) were more abundant in mature seeds. From our data it can be concluded that recalcitrance is established during seed development, which distinguishes between orthodox and recalcitrant seeds. In contrast, post-germination molecular events are quite similar for both types of seeds.
4.5. Proteomics
Proteomics is the scientific discipline whose object of study is the proteome, understood as the total set of proteins present in a biological unit (subcellular fraction, cell, tissue, organ, individual, ecosystem, or purified extracts) at a certain point in development, and under specific environmental conditions. Over time, an increasing number of proteomics branches have emerged, attempting to identify, catalog, and characterize a specific structural or functional group of proteins (proteases, phosphoproteins, membrane proteins, etc.) [106]. The field of plant proteomics was extensively reviewed from a conceptual, historical, methodological, and biological point of view by the authors’ group from 2003 [107] to 2021 [106]. Several of these studies discuss, in depth, the contribution of proteomics to plant biology [8,106,108,109,110,111], forest trees [6] or Quercus [7,8,112] knowledge. The field of proteomics, in relation to its application to the study of Mediterranean forest trees, was also reviewed by Pinheiro et al. [113]. Only three Quercus spp. have been examined using proteomics; Q. robur (1 publication by Sergeant et al. [114]), Q. suber (7 publications), and Q. ilex (20 publications).
The first report on the proteome in Quercus analyzed the protein profiles of Q. suber somatic and gametic in vitro culture-derived embryos using differential gel electrophoresis coupled with MALDI-TOF/TOF [115]. Subsequent papers dealt with somatic embryogenesis [116], responses to the oomycete P. cinnamomi [117], and cork-forming tissues, framed in a project aimed at studying cork quality [118,119].
The review by Rey et al. [7] discussed “How proteomics sees Q. ilex.” To some extent, the proteomes of individual organs (pollen, seeds, leaves, and roots) at different stages (mature trees and seedlings) under optimal and stress conditions have been characterized using classic 2-DE-based MS (number of identified proteins in a single experiment in the range of 400–600) and shotgun approaches (number of identified proteins in a single experiment in the range of 2000–4000 for data-dependent acquisition, and up to 6000 for data-independent acquisition). While the number of proteins confidently identified and quantified depends on the proteomic platform and database employed, our results indicate that the different platforms are complementary, providing new information, and hence allowing for deeper proteome coverage [8,120,121]. In this regard, the creation of a Q. ilex protein database based on available DNA/RNA sequences increased the quality (scores) and number of identified proteins in our experiments by almost three-fold [89,122].
By the mid-2000s, we had published the first proteome reference map of Q. ilex leaves, revealing the great plasticity of the organ, tree variability, and differences between developmental stages [123,124]. With the objective of characterizing Q. ilex populations, the high biological variability found in leaves prompted us to study acorns, which present a more stable proteome. Acorn 2-DE protein profile data allowed for the identification of group populations, correlating them with geographical location and climate conditions [34]. The same objective was pursued by analyzing the pollen proteome profile by 2-DE-MS and shotgun (LC-MSMS), corroborating the former data in terms of variability by provenance [125]. The effect and responses to stresses, drought and P. cinnamomi among individuals and populations were evaluated by analyzing the leaf proteome of 8–10-month-old seedlings. Echevarría-Zomeño et al. [46] reported that changes under drought and water-recovering conditions affected the proteins of photosynthesis, carbohydrate and nitrogen metabolism, as well as stress-related functional categories. A comparative proteomic analysis with two Q. ilex populations, differing in drought tolerance, revealed a general decrease in protein abundance, caused by drought, more pronounced in the less tolerant population and mainly affecting those involved in ATP synthesis and photosynthesis [48]. In parallel, these responses were also analyzed in the roots of three-month-old seedlings [126,127], with results pointing to metabolic adjustment (sugar mobilization and phenolic metabolism) and the activation of the defense system. Recently, we have gone one step further by proposing proteins and derived proteotypic peptides as markers of tolerance to drought to be used for the selection of elite tolerant genotypes in a breeding program. A panel of 30 proteins and 46 derived peptides was proposed, with two proteins (viz., subtilisin and chaperone GrpE protein) being upregulated under water stress in the three highly drought-tolerant populations studied [128].
Q. ilex responses to P. cinnamomi have been studied in one-year-old seedlings from two Andalusian provenances by analyzing changes in the leaf protein profile [47]. The general trend observed was a decrease in protein abundance upon inoculation, albeit slight in the least-susceptible population. The changes observed mostly affect chloroplast proteins involved in photosynthesis, the Calvin cycle, and carbohydrate metabolism, which agrees with what has been previously described in plants subject to drought. We have recently switched to studying combined stress responses in Q. ilex, recreating what occurs in nature and focusing on drought and P. cinnamomi, the two main causes of decline syndrome in Spain. San Eufrasio et al. [33] studied the effect and response of eight-month-old seedlings from three Q. ilex Andalusian populations—differing in their stress tolerance phenotype—to individual and combined stresses (P. cinnamomi and drought) using morphological, physiological, biochemical, and proteomic analyses. Photosynthetic proteins tended to be downregulated while those of stress-responsive proteins were upregulated. Although no treatment-specific response was observed for any functional group, protein changes corresponding to quantitative differences were mainly observed in response to the combined stresses. As a result, aldehyde dehydrogenase, glucose-6-phosphate isomerase, 50S ribosomal protein L5, and α-1,4-glucan-protein synthase (UDP-forming) are proposed as putative markers for resilience in Q. ilex.
Seed maturation and germination have also been studied in Q. ilex using a proteomic approach to understand the molecular basis of acorn recalcitrance. The protein profiles of the different parts of the seed, embryo, cotyledon, and teguments revealed a compartmentalization of metabolic pathways and a division of metabolic tasks [45]. Proteomic results and those from complementary approaches (gel-based and gel-free analyses) demonstrated that recalcitrance is established during seed development, while post-germination events were similar for both orthodox and recalcitrant seeds [44,121]. Our current research is evolving from a holistic to a targeted perspective. For example, Escandón et al. [129] have characterized several proteases and protease inhibitors in Q. ilex seeds by using proteomics combined with an in silico analysis of its transcriptome and in vitro and in gel activity assays. As mentioned above, proteases play a key role in stress responses and seed germination in Q. ilex, justifying the research focus on these enzyme families. Finally, proteomics generates knowledge with a clear translational potential in relation to the identification of allergens and bioactive peptides. Thus, Pedrosa et al. [130] characterized the first allergen reported in the pollen of Q. ilex, Que i 1. In a recent work carried out by our group (unpublished data), a total of 55 bioactive peptides encrypted in 2468 proteins were identified in the species-specific Q. ilex transcriptome database [89]. Of them, 12 bioactive peptides encrypted in 22 proteins were found in acorn protein extract by a shotgun analysis. The bioactive peptide NLAG was the most represented, being encrypted in seven tryptic peptides from four proteins, namely, enolase (qilexprot_67398), ketol–acid reductoisomerase (qilexprot_50910), 26S proteasome non-ATPase regulatory subunit 2 (qilexprot_34420) and Lysyl-tRNA synthetase (qilexprot_15726). It has been reported to be an inhibitor of the angiotensin-converting enzyme (ACE) and, therefore, to have antihypertensive activity [131].
4.6. Metabolomics
Metabolomics, the last omics to appear on the biological research scene, studies the metabolome, the total set of small molecules (Mr lower than 1500 Da) of a biological system—from subcellular fractions to whole organisms—as the product or intermediate of different metabolic pathways. From six papers published by 2001, the number increased to 4562 by 2021 (search on WOS). Of the 14,464 items found on WOS, 4302 corresponded to plants, 95 to forest species, 13 to Quercus, and 9 to Q. ilex. The lowest analytical potential of metabolomic techniques, resulting in hundreds of metabolites identified in a single experiment, compared to tens of thousands of transcripts and thousands of proteins identified in equivalent transcriptomics and proteomic experiments, is counterbalanced by being the closest to the phenotype in the flux of genetic information.
The metabolome is highly complex and dynamic and includes a wide range of physicochemical heterogeneous molecules from primary and secondary metabolic pathways, representing a major challenge for the analytical tools employed. These include classical biochemistry procedures based on colorimetric reactions and non-destructive NIRS, along with mass spectrometry and NMR (nuclear magnetic resonance) platforms for untargeted global analysis. As with the other omics, metabolomics requires the use of algorithms and databases for metabolite annotation, quantitation, and statistical analysis [132].
The field of plant metabolomics, and specifically forest tree metabolomics, has been subject to periodic review [133,134,135,136,137]. Metabolomics, either alone or in combination with other omics techniques, has generated new knowledge on the molecular basis of biological processes in systems biology, such as responses to stresses in herbaceous plants and crops and, to a lesser extent in forest trees [138,139]. Within the genus Quercus, metabolomic studies have been reported in Q. robur, Q. suber, Quercus pubescens Mills, Q, petraea, Quercus aliena Blume, Quercus agrifolia Née, and Q. ilex. Differences in the metabolite profile between species, populations, and individuals have been studied to understand the ecological and evolutionary role of phytochemicals [140]. Buche et al. [141], using untargeted LC-MS, described relevant compounds of Q. robur and Q. petraea wood samples, proposing the use of metabolite signatures as markers of wood quality. In the same direction, the phenolic composition of phellem and its relationship with wood quality has been evaluated in Q. suber [119]. In Q. robur, differences in primary and secondary metabolites between insect-resistant and -susceptible individuals were found. While tolerant trees were enriched in secondary metabolites (tannins, flavonoids, and terpenoids), susceptible ones were enriched in carbohydrate and amino acid derivatives, supporting the hypothesis of a constitutive defense against insects [142,143].
A few studies have dealt with responses to drought and climate change conditions. In a four-year Q. pubescens drought trial, GC-MS metabolomics revealed that tolerance implicates the tricarboxylic acid cycle shunt through the glyoxylate pathway, and the galactose metabolism by reducing carbohydrate storage and increasing proteolytic activity [92]. Changes in the phenolic pattern in Q. suber during drought and recovery treatments have been reported [144], and the process of metabolite resorption from plant leaves during senescence has been evaluated under drought and warming conditions in Q. rubra [145].
Concerning metabolomic studies in response to pathogens, seasonal variations in the phloem phenolic profile and in response to Phytophthora ramorum Werres, De Cock and Man in ‘t Veld were evaluated in Q. agrifolia [146], finding a correlation between the phenolic profile and resistance level. A topic of great interest, but scarcely represented in the literature, is the truffle–Quercus interaction. In this regard, Li et al. [147] investigated the Tuber indicum P. Micheli ex F. H. Wigg–Q. aliena symbiosis at the metabolomic level.
With respect to Q. ilex, Prof. J. Peñuelas’ group evaluated seasonal-, climate change-, and drought-induced metabolic changes through field experiments, with a special emphasis on volatiles [148,149]. Observed changes under moderate drought conditions included a higher content of sugars and phenolics, with these changes being related to folivory [150]. Drought tolerance is linked to signaling and osmoregulation by hexoses, and to antioxidant protectors such as polyols and amino acids [151]. The work carried out by Rivas-Ubach et al. [152] revealed differences in the metabolite profile depending on the soil characteristics and tree size.
Our group approached Q. ilex metabolomics with different objectives using complementary platforms, including classical biochemical colorimetric assays, NIRS, and untargeted metabolomics by GC-MS and LC-MS platforms. The first published work in this field [34] reported on natural variability in Q. ilex Andalusian populations based on acorn morphometry and chemical composition determined by NIRS analysis. Two main clusters were detected, corresponding to the northern and southern provenances, with the former having higher protein, sugar, linoleic, and oleic acid content. Because of the renewed interest in the alimentary use of acorns, either as nuts or for flour-derived products, the phytochemical composition variability among morphotypes was determined by LC-MS [36]. Of the 3259 (positive mode) and 1841 (negative mode) resolved compounds, 192 metabolites were annotated, with differences mostly in the carbohydrates, amino acids, and secondary metabolites. Some of the identified compounds had bioactivity (antimicrobial, insecticide, allelochemical, antioxidant, antitumor, proangiogenic, etc.), proving the nutraceutical value of the acorns. The most recent work published by our group focused on the response to drought stress in Q. ilex seedlings by using a non-targeted metabolomic approach [153]. The leaf metabolome has been described to determine possible mechanisms and molecular markers of drought tolerance, as well as identify putative bioactive compounds. A total of 3934 compounds were identified: 616 variable and 342 annotated (Figure 5). Of the annotated compounds, 33 were variable, mostly corresponding to amino acids, carboxylic acids, benzenoids, flavonoids, and isoprenoids.
By using a multiomics approach, including transcriptomics (NGS-Illumina), proteomics (shotgun LC-MS/MS), and metabolomics (GC-MS), Q. ilex metabolism has been partially reconstructed [154]. These analyses resulted in the identification of 62,629 transcripts, 2380 protein species, and 62 metabolites. Of the 127 metabolic pathways reported in the KEGG pathway database, 123 can be visualized using the described methodology. The TCA cycle was the pathway most represented, with 5 out of 10 metabolites, 6 out of 8 protein enzymes, and 8 out of 8 enzyme transcripts.
Ongoing research is addressing changes in the metabolite profile in response to individual or combined (drought and P. cinnamomi) stresses, and the search for resilience metabolite markers; furthermore, we are evaluating metabolic changes during seed germination.
5. Conclusions and Future Directions
In 2004, our group started a research project on molecular biology with forest tree species, mainly focused on Q. ilex. The general objective was to understand the biology of the species, follow the central dogma of molecular science using the omics cascade, and integrate the omics data obtained with those from morphology, physiology, and classical biochemistry in the direction of systems biology. For this, our specific objectives were to characterize variability, genes, and gene products implicated in the response to biotic and abiotic stresses related to decline syndrome and climate change, later evolving towards studies of seed recalcitrance and ultimately to acorn phytochemical analysis, related to its nutritional and nutraceutical value as seeds and flour-derived products. The knowledge generated is intended to be translated into breeding programs, assisted by molecular markers and based on the selection of resilient and productive elite genotypes, and into seed conservation strategies to aid the search for bioactive compounds. Nowadays, we are embracing new objectives, such as the analysis of plant biostimulants [155,156] and the priming phenomenon [157] for protecting Q. ilex against stresses related to decline syndrome. Preliminary unpublished results show that seeds treated with BTH (benzothiadiazole), a chemical analogue of the defense-related molecule salicylic acid, display an increased and accelerated germination rate, reducing seedling infection by P. cinnamomi.
This project is proving challenging due to the recalcitrance of the species as an experimental system, the enormous intra- and interpopulation variability found, and characteristics such as its long life cycle and allogamy. Field data are difficult to interpret in biological terms, and greenhouse-controlled experiments are limited to a short period of time each year. Beyond doubt, the use of clonal individuals generated by vegetative propagation or somatic embryogenesis will favor molecular research and data analysis [158,159]. The allogamous character of Q. ilex makes it extremely difficult to draw conclusions about inheritance, since acorns from individual trees, and the seedlings and plants derived from them, are highly variable due to the unknown origin of the pollen. Reliable interpretations of the results would only be possible by including a high number of biological replicates in the experiment, which is incompatible with omics studies. In this regard, it would be convenient to conduct the experiments under controlled pollination conditions, although no optimized protocols for this have been reported in the literature so far. As an example, in a recent experiment to assess the response to the combined stress of drought and P. cinnamomi in the offspring of symptomatic and asymptomatic Q. ilex trees surveyed in areas of decline, all possible combinations, with symptomatic and asymptomatic seedlings coming from symptomatic and asymptomatic mother trees, were seen (unpublished data). This clearly illustrates the complexity of the resilient phenotype in which genetic, epigenetic, and environmental factors are involved, and which remains poorly investigated in healthy or diseased individuals of Q. ilex, in particular, and of the genus Quercus in general.
However, our current research with Q. ilex, together with that reported on other Quercus spp., mainly Q. lobata, Q. robur, and Q. suber, confirm that the most recent and cutting-edge molecular biology techniques successfully used with model plants and crops are suitable for this species. We intend, following the path of other Quercus spp. research (e.g.,
The sequencing of the Q. ilex genome opens many possibilities for future research, of which we will prioritize the implementation of state-of-the-art approaches such as a genome-wide association study (GWAS) [161], as well as the identification of SNP variants [162]. Epigenetics is a priority in our current and future research [76]. We have just started to address this field by measuring the global methylation level of cytosines at different development stages and organs using the MSAP technique, and have concluded that it is not a very informative technique for this purpose. These results have determined the switch to bisulfite sequencing to characterize DNA methylation markers, in addition to expanding our studies of epigenetic mechanisms to the modification of histones and small RNAs, mostly linked to priming and stress memory [163]. The availability of genomic, transcriptomic, and proteomic information in our database will greatly facilitate the identification of pivotal enzymes in epigenetic mechanisms.
Finally, it is worth mentioning our intention to open up a new line of research, that of functional or reverse genomics [164], to validate those genes and gene products proposed as candidate markers of resilience [33,35,89]. This task is unachievable within Q. ilex, since the genetic resources and protocols for transformation and regeneration are not available. For this, it is necessary the use heterologous functional validation systems, such as A. thaliana or Nicotiana benthamiana, as has already been reported for Q. suber genes related to floral organ identity and drought response [165,166].
In summary, the objectives that we pursued with our research on Q. ilex were previously proposed by Plomión et al. [167] and achieved to a certain extent with other forest tree species [168].
Conceptualization, M.Á.C., M.-D.R. and J.V.J.-N.; methodology, all authors; writing—original draft preparation, M.-D.R. and J.V.J.-N.; writing—review and editing, all authors; funding acquisition, J.V.J.-N. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
Not applicable.
C.L.-H. is grateful for the contract “Severo Ochoa Program (BP17-112, Government of Principado de Asturias, Spain)”. M.Á.C. and M.-D.R. are grateful for the awards of Ramón y Cajal (RYC-2017-23706), and the Juan de la Cierva-Incorporación (IJC2018-035272-I) contracts from the Spanish Ministry of Science, Innovation and Universities, respectively.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Tree used to obtain the first draft genome of Q. ilex, located in Aldea de Cuenca, Fuente Obejuna, Cordoba, Andalusia, Spain (38°19′46″ N, 5°33′15″ W).
Enlarge this image.Figure 2. Electropherogram of Udt4 (chloroplast) and QrZAG15 (nuclear) SSR locus in Q. ilex (A) and Q. suber (B).
Enlarge this image.Figure 3. Examples of MSAP polymorphisms in electropherograms corresponding to Q. ilex adult tree leaves (A), embryos (B) and seedling leaves (C) using MspI + EcoRI digestion and the primer combination EcoRI + ACG and MspI + TTA.
![View Image - Figure 4. Pie charts of GO classification of assembled Q. ilex transcripts including the first 15 transcripts assigned to biological processes (A), molecular functions (B) and cellular components (C). The number in parentheses indicates the number of transcripts. The remaining transcripts assigned to each GO category are shown in Guerrero-Sánchez et al. [104].](https://media.proquest.com/media/hms/THUM/8/0bL3O?_a=ChgyMDI1MDcxNDAyMDM0NDAxODo1OTY0ODUSBzEzODMzMjMaCk9ORV9TRUFSQ0giDDIxNi43My4yMTcuNCoHMjAzMjM0MTIKMjcxMTM2ODA5MDoLSW5saW5lSW1hZ2VCATBSBk9ubGluZVoDSUxJYgRUSFVNagoyMDIyLzA5LzAxcgoyMDIyLzA5LzE0egCCASBQLTEwMDAwODQtMjAxMzk1LVNUQUZGLW51bGwtbnVsbJIBBk9ubGluZcoBZ01vemlsbGEvNS4wIEFwcGxlV2ViS2l0LzUzNy4zNiAoS0hUTUwsIGxpa2UgR2Vja287IGNvbXBhdGlibGU7IENsYXVkZUJvdC8xLjA7ICtjbGF1ZGVib3RAYW50aHJvcGljLmNvbSnSARJTY2hvbGFybHkgSm91cm5hbHOaAgdQcmVQYWlkqgIbT1M6RU1TLUZ1bGxUZXh0LWdldEZ1bGxUZXh0ygIPQXJ0aWNsZXxGZWF0dXJl4gIBTvICAPoCAVmCAwNXZWKKAxxDSUQ6MjAyNTA3MTQwMjAzNDQwMTg6MjkwMDA2&_s=slkHCsG2F4vz6vFzlFWr6aKetFY%3D "View Image - Figure 4. Pie charts of GO classification of assembled Q. ilex transcripts including the first 15 transcripts assigned to biological processes (A), molecular functions (B) and cellular components (C). The number in parentheses indicates the number of transcripts. The remaining transcripts assigned to each GO category are shown in Guerrero-Sánchez et al. [104].")
Enlarge this image.Figure 4. Pie charts of GO classification of assembled Q. ilex transcripts including the first 15 transcripts assigned to biological processes (A), molecular functions (B) and cellular components (C). The number in parentheses indicates the number of transcripts. The remaining transcripts assigned to each GO category are shown in Guerrero-Sánchez et al. [104].
![View Image - Figure 5. Pie chart of the chemical families identified in Q. ilex seedlings subjected to drought stress. The number in parentheses indicates the annotated compounds [148].](https://media.proquest.com/media/hms/THUM/10/0bL3O?_a=ChgyMDI1MDcxNDAyMDM0NDAxODo1OTY0ODUSBzEzODMzMjMaCk9ORV9TRUFSQ0giDDIxNi43My4yMTcuNCoHMjAzMjM0MTIKMjcxMTM2ODA5MDoLSW5saW5lSW1hZ2VCATBSBk9ubGluZVoDSUxJYgRUSFVNagoyMDIyLzA5LzAxcgoyMDIyLzA5LzE0egCCASBQLTEwMDAwODQtMjAxMzk1LVNUQUZGLW51bGwtbnVsbJIBBk9ubGluZcoBZ01vemlsbGEvNS4wIEFwcGxlV2ViS2l0LzUzNy4zNiAoS0hUTUwsIGxpa2UgR2Vja287IGNvbXBhdGlibGU7IENsYXVkZUJvdC8xLjA7ICtjbGF1ZGVib3RAYW50aHJvcGljLmNvbSnSARJTY2hvbGFybHkgSm91cm5hbHOaAgdQcmVQYWlkqgIbT1M6RU1TLUZ1bGxUZXh0LWdldEZ1bGxUZXh0ygIPQXJ0aWNsZXxGZWF0dXJl4gIBTvICAPoCAVmCAwNXZWKKAxxDSUQ6MjAyNTA3MTQwMjAzNDQwMTg6MjkwMDA2&_s=aWml63rph7swzcGpLRHfNiNQ5Nk%3D "View Image - Figure 5. Pie chart of the chemical families identified in Q. ilex seedlings subjected to drought stress. The number in parentheses indicates the annotated compounds [148].")
Enlarge this image.Figure 5. Pie chart of the chemical families identified in Q. ilex seedlings subjected to drought stress. The number in parentheses indicates the annotated compounds [148].
Number of papers found in the Web of Science (WOS) database (28 July 2022) from a search using the name of the organism (first column) and the omics approach employed (second to fifth columns) as keywords (included in the abstract).
| Organisms | Genomics | Transcriptomics | Proteomics | Metabolomics |
|---|---|---|---|---|
| Quercus | 27 | 9 | 20 | 9 |
| Quercus ilex | 0 | 3 | 15 | 4 |
| Populus | 161 | 29 | 34 | 28 |
| Pinus | 70 | 11 | 28 | 12 |
| Eucalyptus | 55 | 8 | 26 | 14 |
| Arabidopsis | 1197 | 293 | 723 | 344 |
| Rice | 1137 | 153 | 452 | 312 |
| Saccharomyces | 548 | 101 | 387 | 192 |
| Humans | 11,064 | 2842 | 11,458 | 4975 |
Assembly statistics of the whole-genome sequencing of Quercus spp.
| Quercus spp. | Genome Size (Mb) | Contig N50 (Mb) | Number of Scaffolds | GC Content (%) | Repetitive |
Protein-Coding Genes | NCBI Genbank b |
|---|---|---|---|---|---|---|---|
| Q. robur a | 814.3 | 1.35 | 1409 | 35.5 | 53.30 | 25,808 | GCA_932294415 [ |
| Q. lobata | 845.9 | 0.97 | 2002 | 35.0 | 54.00 | 39,373 | GCA_001633185 [ |
| Q. suber | 953.3 | 0.81 | 23,347 | 36.0 | 11.96 | 33,658 | GCA_002906115 [ |
| Q. mongolica | 809.8 | 2.40 | 321 | 35.5 | 53.75 | 36,553 | GCA_011696235 [ |
| Q. aquifoloides | 926.5 | 1.40 | 212 | 36.5 | NA c | NA | GCA_019022515 |
| Q. wislizeni | 724.8 | 9.00 | 358 | 35.0 | NA | NA | GCA_023055345 |
| Q. gilva | 889.8 | 28.30 | 515 | 35.5 | NA | NA | GCA_023621385 |
| Q. glauca | 903.1 | 7.60 | 415 | 35.5 | NA | NA | GCA_023736055 |
a The haploid genome of Q. robur. b NCBI Genbank (National Center for Biotechnology Information Genbank database). c NA indicates data have not been published yet. This information has been obtained from the NCBI Genbank.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Valledor, L.; Jorrín-Novo, J.V. Back to the basics: Maximizing the information obtained by quantitative two-dimensional gel electrophoresis analyses by an appropriate experimental design and statistical analyses. J. Proteom.; 2011; 74, pp. 1-18. [DOI: https://dx.doi.org/10.1016/j.jprot.2010.07.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20656082]
2. Echevarría-Zomeño, S.; Abril, N.; Ruiz-Laguna, J.; Jorrín-Novo, J.; Maldonado-Alconada, A.M. Simple, rapid and reliable methods to obtain high quality RNA and genomic DNA from Quercus ilex L. leaves suitable for molecular biology studies. Acta Physiol. Plant.; 2012; 34, pp. 793-805. [DOI: https://dx.doi.org/10.1007/s11738-011-0880-z]
3. Loreto, F.; Ciccioli, P.; Cecinato, A.; Brancaleoni, E.; Frattoni, M.; Fabozzi, C.; Tricoli, D. Evidence of the photosynthetic origin of monoterpenes emitted by Quercus ilex L. leaves by 13C labeling. Plant Physiol.; 1996; 110, pp. 1317-1322. [DOI: https://dx.doi.org/10.1104/pp.110.4.1317] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12226263]
4. Fischbach, R.J.; Zimmer, I.; Steinbrecher, R.; Pfichner, A.; Schnitzler, J.P. Monoterpene synthase activities in leaves of Picea abies (L.) Karst. and Quercus ilex L. Phytochemistry; 2000; 54, pp. 257-265. [DOI: https://dx.doi.org/10.1016/S0031-9422(00)00119-9]
5. Cherbuy, B.; Joffre, R.; Gillon, D.; Rambal, S. Internal remobilization of carbohydrates, lipids, nitrogen and phosphorus in the Mediterranean evergreen oak Quercus ilex. Tree Physiol.; 2001; 21, pp. 9-17. [DOI: https://dx.doi.org/10.1093/treephys/21.1.9]
6. Abril, N.; Gion, J.M.; Kerner, R.; Müller-Starck, G.; Navarro-Cerrillo, R.M.; Plomion, C.; Renaut, J.; Valledor, L.; Jorrín-Novo, J.V. Proteomics research on forest trees, the most recalcitrant and orphan plant species. Phytochemistry; 2011; 72, pp. 1219-1242. [DOI: https://dx.doi.org/10.1016/j.phytochem.2011.01.005]
7. Rey, M.-D.; Castillejo, M.Á.; Sánchez-Lucas, R.; Guerrero-Sanchez, V.M.; López-Hidalgo, C.; Romero-Rodríguez, C.; Valero-Galván, J.; Sghaier-Hammami, B.; Simova-Stoilova, L.; Echevarría-Zomeño, S. et al. Proteomics, holm oak (Quercus ilex L.) and other recalcitrant and orphan forest tree species: How do they see each other?. Int. J. Mol. Sci.; 2019; 20, 692. [DOI: https://dx.doi.org/10.3390/ijms20030692]
8. Escandón, M.; Castillejo, M.Á.; Jorrín-Novo, J.V.; Rey, M.D. Molecular research on stress responses in Quercus spp. : From classical biochemistry to Systems Biology through Omics analysis. Forests; 2021; 12, 364.
9. Sork, V.L.; Cokus, S.J.; Fitz-Gibbon, S.T.; Zimin, A.V.; Puiu, D.; Garcia, J.A.; Gugger, G.F.; Henriquez, C.L.; Zhen, Y.; Lohmueller, K.E. et al. High-quality genome and methylomes illustrate features underlying evolutionary success of oaks. Nat. Commun.; 2022; 13, 2047. [DOI: https://dx.doi.org/10.1038/s41467-022-29584-y]
10. Nixon, K. Infrageneric classification of Quercus (Fagaceae) and typification of sectional names. Ann. For. Sci.; 1993; 50, pp. 25s-34s. [DOI: https://dx.doi.org/10.1051/forest:19930701]
11. Backs, J.R.; Ashley, M.V. Quercus Genetics: Insights into the past, present, and future of oaks. Forests; 2021; 12, 1628. [DOI: https://dx.doi.org/10.3390/f12121628]
12. Kew Royal Botanic Gardens. Available online: https://powo.science.kew.org/taxon/urn:lsid:ipni.org:names:325819-2#children (accessed on 15 June 2022).
13. Leroy, T.; Plomion, C.; Kremer, A. Oak symbolism in the light of genomics. New Phytol.; 2020; 226, pp. 1012-1017. [DOI: https://dx.doi.org/10.1111/nph.15987]
14. de Rigo, D.; Caudullo, G. Quercus ilex in Europe: Distribution, habitat, usage and threats. European Atlas of Forest Tree Species; San-Miguel-Ayanz, J.; de Rigo, D.; Caudullo, G.; Houston Durrant, T.; Mauri, A. Publications Office of the European Union: Luxembourg, 2016; pp. 152-153.
15. Pulido, F.J.; Diaz, M.; Hidalgo de Trucios, S.J. Size structure and regeneration of Spanish holm oak Quercus ilex forests and dehesas: Effects of agroforestry use in their long-term sustainability. For. Ecol. Manag.; 2016; 146, pp. 1-13. [DOI: https://dx.doi.org/10.1016/S0378-1127(00)00443-6]
16. Plieninger, T.; Pulido, F.J.; Konold, W. Effects of land-use history on size structure of holm oak stands in Spanish dehesas: Implications for conservation and restoration. Environ. Conserv.; 2003; 30, pp. 61-70. [DOI: https://dx.doi.org/10.1017/S0376892903000055]
17. Plieninger, T.; Flinzberger, L.; Hetman, M.; Horstmannshoff, I.; Reinhard-Kolempas, M.; Topp, E.; Moreno, G.; Huntsinger, L. Dehesas as high nature value farming systems: A social-ecological synthesis of drivers, pressures, state, impacts, and responses. Ecol. Soc.; 2021; 26, 23. [DOI: https://dx.doi.org/10.5751/ES-12647-260323]
18. Gea-Izquierdo, G.; Natalini, F.; Cardillo, E. Holm oak death is accelerated but not sudden and expresses drought legacies. Sci. Total Environ.; 2021; 754, 141793. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2020.141793]
19. Campos, P.; Álvarez, A.; Mesa, B.; Oviedo, J.L.; Ovando, P.; Caparrós, A. Total income and ecosystem service sustainability index: Accounting applications to holm oak dehesa case study in Andalusia-Spain. Land Use Policy; 2020; 97, 104692. [DOI: https://dx.doi.org/10.1016/j.landusepol.2020.104692]
20. Natalini, F.; Alejano, R.; Vázquez-Piqué, J.; Cañellas, I.; Gea-Izquierdo, G. The role of climate change in the widespread mortality of holm oak in open woodlands of Southwestern Spain. Dendrochronologia; 2016; 38, pp. 51-60. [DOI: https://dx.doi.org/10.1016/j.dendro.2016.03.003]
21. Talabi, A.O.; Vikram, P.; Thushar, S.; Rahman, H.; Ahmadzai, H.; Nhamo, N.; Shahid, M.; Singh, R.K. Orphan crops: A best fit for dietary enrichment and diversification in highly deteriorated marginal environments. Front. Plant Sci.; 2022; 13, 839704. [DOI: https://dx.doi.org/10.3389/fpls.2022.839704]
22. Kumar, A.; Anju, T.; Kumar, S.; Chhapekar, S.S.; Sreedharan, S.; Singh, S.; Choi, S.R.; Ramchiary, N.; Lim, Y.P. Integrating Omics and Gene Editing Tools for Rapid Improvement of Traditional Food Plants for Diversified and Sustainable Food Security. Int. J. Mol. Sci.; 2021; 22, 8093. [DOI: https://dx.doi.org/10.3390/ijms22158093]
23. Godeto, Y.G.; Bachheti, A.; Husen, A.; Pandey, D.P.; Bachheti, R.K. Forest-based edible seeds and nuts for health care and disease control. Non-Timber Forest Products; Springer: Cham, Switzerland, 2021; pp. 145-174.
24. Sardão, R.; Amaral, R.A.; Alexandre, E.M.; Saraiva, J.A.; Pintado, M. Effect of high-pressure processing to improve the safety and quality of a Quercus acorn beverage. LWT; 2021; 149, 111858. [DOI: https://dx.doi.org/10.1016/j.lwt.2021.111858]
25. Chao, B.; Liu, R.; Zhang, X.; Zhang, X.; Tan, T. Tannin extraction pretreatment and very high gravity fermentation of acorn starch for bioethanol production. Bioresour. Technol.; 2017; 241, pp. 900-907. [DOI: https://dx.doi.org/10.1016/j.biortech.2017.06.026] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28629106]
26. Vinha, A.F.; Costa, A.S.G.; Barreira, J.C.M.; Pacheco, R.; Oliveira, M.B.P.P. Chemical and antioxidant profiles of acorn tissues from Quercus spp.: Potential as new industrial raw materials. Ind. Crop. Prod.; 2016; 94, pp. 143-151. [DOI: https://dx.doi.org/10.1016/j.indcrop.2016.08.027]
27. Sotirov, M.; Pokorny, B.; Kleinschmit, D.; Kanowski, P. International forest governance and policy: Institutional architecture and pathways of influence in global sustainability. Sustainability; 2020; 12, 7010. [DOI: https://dx.doi.org/10.3390/su12177010]
28. Senf, C.; Sebald, J.; Seidl, R. Increasing canopy mortality affects the future demographic structure of Europe’s forests. One Earth; 2021; 4, pp. 749-755. [DOI: https://dx.doi.org/10.1016/j.oneear.2021.04.008]
29. Petit-Cailleux, C.; Davi, H.; Lefèvre, F.; Verkerk, P.J.; Fady, B.; Lindner, M.; Oddou-Muratorio, S. Tree mortality risks under climate change in Europe: Assessment of silviculture practices and genetic conservation networks. Front. Ecol. Evol.; 2022; 9, 706414. [DOI: https://dx.doi.org/10.3389/fevo.2021.706414]
30. Ogaya, R.; Liu, D.; Barbeta, A.; Peñuelas, J. Stem mortality and forest dieback in a 20-years experimental drought in a Mediterranean holm oak forest. Front. For. Glob. Change; 2020; 2, 89. [DOI: https://dx.doi.org/10.3389/ffgc.2019.00089]
31. Brasier, C. Phytophthora cinnamomi and oak decline in southern Europe. Environmental constraints including climate change. Ann. For. Sci.; 1996; 53, pp. 347-358.
32. Ruiz-Gómez, F.; Pérez-de-Luque, A.; Sánchez-Cuesta, R.; Quero, J.; Cerrillo, R.M.N. Differences in the response to acute drought and Phytophthora cinnamomi rands infection in Quercus ilex L. seedlings. Forests; 2018; 9, 634. [DOI: https://dx.doi.org/10.3390/f9100634]
33. San-Eufrasio, B.; Castillejo, M.Á.; Labella-Ortega, M.; Ruiz-Gómez, F.J.; Navarro-Cerrillo, R.M.; Tienda-Parrilla, M.; Jorrín-Novo, J.V.; Rey, M.D. Effect and response of Quercus ilex subsp. ballota [Desf.] Samp. seedlings from three contrasting Andalusian populations to individual and combined Phytophthora cinnamomi and drought stresses. Front. Plant Sci.; 2021; 12, 722802.
34. Valero-Galván, J.; Valledor, L.; Navarro-Cerrillo, R.M.; Gil Pelegrín, E.; Jorrín-Novo, J.V. Studies of variability in Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.) through acorn protein profile analysis. J. Proteom.; 2011; 74, pp. 1244-1255.
35. San-Eufrasio, B.; Sánchez-Lucas, R.; López-Hidalgo, C.; Guerrero-Sanchez, V.M.; Castillejo, M.; Maldonado-Alconada, A.M.; Jorrín-Novo, J.V.; Rey, M.-D. Responses and differences in tolerance to water shortage under climatic dryness conditions in seedlings from Quercus spp. and Andalusian Q. ilex populations. Forests; 2020; 11, 707. [DOI: https://dx.doi.org/10.3390/f11060707]
36. López-Hidalgo, C.; Trigueros, M.; Menéndez, M.; Jorrin-Novo, J.V. Phytochemical composition and variability in Quercus ilex acorn morphotypes as determined by NIRS and MS-based approaches. Food Chem.; 2021; 338, 127803. [DOI: https://dx.doi.org/10.1016/j.foodchem.2020.127803]
37. Sork, V.L. Genomic studies of local adaptation in natural plant populations. J. Hered.; 2017; 109, pp. 3-15. [DOI: https://dx.doi.org/10.1093/jhered/esx091]
38. Plomion, C.; Aury, J.; Amselem, J.; Leroy, T.; Murat, F.; Duplessis, S.; Faye, S.; Francillonne, N.; Labadie, K.; Provost, G.L. et al. Oak genome reveals facets of long lifespan. Nat. Plants; 2018; 4, pp. 440-452. [DOI: https://dx.doi.org/10.1038/s41477-018-0172-3]
39. Valero Galván, J.; Jorrín Novo, J.V.; Cabrera, A.G.; Ariza, D.; García-Olmo, J.; Navarro-Cerrillo, R.M. Population variability based on the morphometry and chemical composition of the acorn in Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.). Eur. J. For. Res.; 2012; 131, pp. 893-904. [DOI: https://dx.doi.org/10.1007/s10342-011-0563-8]
40. Pignone, D.; Laghetti, G. On sweet acorn (Quercus spp.) cake tradition in Italian cultural and ethnic islands. Genet. Resour. Crop Evol.; 2010; 57, pp. 1261-1266. [DOI: https://dx.doi.org/10.1007/s10722-010-9625-x]
41. Stiti, B.; Mezni, F.; Zarroug, Y.; Slama, A.; Fkiri, S. Agri-food valuation of oak fruits in human consumption: Formulation of desserts. J. Food Chem. Nanotechnol.; 2022; 8, pp. 38-42.
42. Vinha, A.F.; Barreira, J.C.; Costa, A.S.; Oliveira, M.B.P. A new age for Quercus spp. fruits: Review on nutritional and phytochemical composition and related biological activities of acorns. Compr. Rev. Food Sci.; 2016; 15, pp. 947-981. [DOI: https://dx.doi.org/10.1111/1541-4337.12220]
43. Sghaier-Hammami, B.; Castillejo, M.Á.; Baazaoui, N.; Jorrín-Novo, J.V.; Escandón, M. GeLC-Orbitrap/MS and 2-DE-MALDITOF/TOF comparative proteomics analysis of seed cotyledons from the non-orthodox Quercus ilex tree species. J. Proteom.; 2021; 233, 104087. [DOI: https://dx.doi.org/10.1016/j.jprot.2020.104087]
44. Sghaier-Hammami, B.; Hammami, S.; Baazaoui, N.; Gómez-Díaz, C.; Jorrín-Novo, J.V. Dissecting the seed maturation and germination processes in the non-orthodox Quercus ilex species based on protein signatures as revealed by 2-DE coupled to MALDI-TOF/TOF proteomics strategy. Int. J. Mol. Sci.; 2020; 21, 4870. [DOI: https://dx.doi.org/10.3390/ijms21144870]
45. Sghaier-Hammami, B.; Redondo-López, I.; Valero-Galván, J.; Jorrín-Novo, J.V. Protein profile of cotyledon, tegument, and embryonic axis of mature acorns from a non-orthodox plant species: Quercus ilex. Planta; 2016; 243, pp. 369-396. [DOI: https://dx.doi.org/10.1007/s00425-015-2404-3]
46. Echevarría-Zomeño, S.; Ariza, D.; Jorge, I.; Lenz, C.; del Campo, A.D.; Jorrín, J.; Navarro, R.M. Changes in the protein profile of Quercus ilex leaves in response to drought stress and recovery. J. Plant Physiol.; 2009; 166, pp. 233-245. [DOI: https://dx.doi.org/10.1016/j.jplph.2008.05.008]
47. Sghaier-Hammami, B.; Valero-Galván, J.; Romero-Rodríguez, M.C.; Navarro-Cerrillo, R.M.; Abdelly, C.; Jorrín-Novo, J.V. Physiological and proteomics analyses of Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.) responses to Phytophthora cinnamomi. Plant Physiol. Biochem.; 2013; 71, pp. 191-202.
48. Valero-Galván, J.; González-Fernández, R.; Navarro-Cerrillo, R.; Gil-Pelegrin, E.; Jorrín-Novo, J. Physiological and proteomic analyses of drought stress response in Holm oak provenances. J. Proteome Res.; 2013; 12, pp. 5110-5123. [DOI: https://dx.doi.org/10.1021/pr400591n]
49. Scotti-Saintagne, C.; Mariette, S.; Porth, I.; Goicoechea, P.G.; Barreneche, T.; Bodénes, C.; Burg, K.; Kremer, A. Genome scanning for interspecific differentiation between two closely related oak species [Quercus robur L. and Q. petraea (Matt.) Liebl.]. Genetics; 2004; 168, pp. 1615-1626. [DOI: https://dx.doi.org/10.1534/genetics.104.026849]
50. Simeone, M.C.; Piredda, R.; Papini, A.; Vessella, F.; Schirone, B. Application of plastid and nuclear markers to DNA barcoding of Euro-Mediterranean oaks (Quercus, Fagaceae): Problems, prospects and phylogenetic implications. Bot. J. Linn.; 2013; 172, pp. 478-499. [DOI: https://dx.doi.org/10.1111/boj.12059]
51. Lazic, D.; Hipp, A.L.; Carlson, J.E.; Gailing, O. Use of genomic resources to assess adaptive divergence and introgression in oaks. Forests; 2021; 12, 690. [DOI: https://dx.doi.org/10.3390/f12060690]
52. Liang, Y.Y.; Shi, Y.; Yuan, S.; Zhou, B.F.; Chen, X.Y.; An, Q.Q.; Ingvarsson, P.K.; Plomión, C.; Wang, B. Linked selection shapes the landscape of genomic variation in three oak species. New Phytol.; 2022; 233, pp. 555-568. [DOI: https://dx.doi.org/10.1111/nph.17793]
53. Coelho, A.C.; Lima, M.B.; Neves, D.; Cravador, A. Genetic diversity of two evergreen oaks [Quercus suber (L.) and Quercus ilex subsp rotundifolia (Lam.)] in Portugal using AFLP markers. Silvae Genet.; 2006; 55, pp. 105-118. [DOI: https://dx.doi.org/10.1515/sg-2006-0016]
54. Belahbib, N.; Pemonge, M.H.; Ouassou, A.; Sbay, H.; Kremer, A.; Petit, R.J. Frequent cytoplasmic exchanges between oak species that are not closely related: Quercus suber and Q. ilex in Morocco. Mol. Ecol.; 2001; 10, pp. 2003-2012. [DOI: https://dx.doi.org/10.1046/j.0962-1083.2001.01330.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11555243]
55. Lopez de Heredia, U.; Mora-Márquez, F.; Goicoechea, P.G.; Guillardín-Calvo, L.; Simeone, M.C.; Soto, Á. ddRAD sequencing-based identification of genomic boundaries and permeability in Quercus ilex and Q. suber hybrids. Front. Plant Sci.; 2020; 11, 564414. [DOI: https://dx.doi.org/10.3389/fpls.2020.564414]
56. Dzialuk, A.; Chybicki, I.; Welc, M.; Śliwińska, E.; Burczyk, J. Presence of triploids among oak species. Ann. Bot.; 2007; 99, pp. 959-964. [DOI: https://dx.doi.org/10.1093/aob/mcm043] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17409100]
57. Craft, K.J.; Owens, J.D.; Ashley, M.V. Application of plant DNA markers in forensic botany: Genetic comparison of Quercus evidence leaves to crime scene trees using microsatellites. Forensic Sci. Int.; 2007; 165, pp. 64-70. [DOI: https://dx.doi.org/10.1016/j.forsciint.2006.03.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16632287]
58. Fernández i Marti, A.; Romero-Rodríguez, C.; Navarro-Cerrillo, R.; Abril, N.; Jorrín-Novo, J.; Dodd, R. Population genetic diversity of Quercus ilex subsp. ballota (Desf.) Samp. reveals divergence in recent and evolutionary migration rates in the Spanish dehesas. Forests; 2018; 9, 337.
59. Grover, A.; Sharma, P.C. Development and use of molecular markers: Past and present. Crit. Rev. Biotechnol.; 2016; 36, pp. 290-302. [DOI: https://dx.doi.org/10.3109/07388551.2014.959891]
60. Shekhar, C.; Rawat, A.; Bhandari, M.S.; Barthwal, S.; Ginwal, H.S.; Meena, R.K. Cross-transferability-based identification and validation of simple sequence repeat (SSR) markers in oaks of western Himalayas. Silvae Genet.; 2021; 70, pp. 108-116. [DOI: https://dx.doi.org/10.2478/sg-2021-0009]
61. Soto, A.; Lorenzo, Z.; Gil, L. Nuclear microsatellite markers for the identification of Quercus ilex L. and Q. suber L. hybrids. Silvae Genet.; 2003; 52, pp. 63-66.
62. Soto, A.; Lorenzo, Z.; Gil, L. Differences in fine-scale genetic structure and dispersal in Quercus ilex L. and Q. suber L.: Consequences for regeneration of Mediterranean open woods. Heredity; 2007; 99, pp. 601-607. [DOI: https://dx.doi.org/10.1038/sj.hdy.6801007]
63. Lopez De Heredia, U.; Sanchez, H.; Soto, A. Molecular evidence of bidirectional introgression between Quercus suber and Quercus ilex. IForest; 2018; 11, 338. [DOI: https://dx.doi.org/10.3832/ifor2570-011]
64. Ortego, J.; Bonal, R. Natural hybridisation between kermes (Quercus coccifera L.) and holm oaks (Q. ilex L.) revealed by microsatellite markers. Plant Biol.; 2010; 12, pp. 234-238. [DOI: https://dx.doi.org/10.1111/j.1438-8677.2009.00244.x]
65. Guzmán, B.; Rodríguez López, C.M.; Forrest, A.; Cano, E.; Vargas, P. Protected areas of Spain preserve the neutral genetic diversity of Quercus ilex L. irrespective of glacial refugia. Tree Genet. Genomes; 2015; 11, 124. [DOI: https://dx.doi.org/10.1007/s11295-015-0950-2]
66. Fineschi, S.; Cozzolino, S.; Migliaccio, M.; Musacchio, A.; Innocenti, M.; Vendramin, G.G. Sicily represents the Italian reservoir of chloroplast DNA diversity of Quercus ilex L. (Fagaceae). Ann. For. Sci.; 2005; 62, pp. 79-84. [DOI: https://dx.doi.org/10.1051/forest:2004085]
67. Plomion, C.; Aury, J.; Amselem, J.; Alaeitabar, T.; Barbe, V.; Belser, C.; Bergès, H.; Bodénès, C.; Boudet, N.; Boury, C. et al. Decoding the oak genome: Public release of sequence data, assembly, annotation and publication strategies. Mol. Ecol. Resour.; 2016; 16, pp. 254-265. [DOI: https://dx.doi.org/10.1111/1755-0998.12425]
68. Sork, V.; Fitz-Gibbon, S.; Puiu, D.; Crepeau, M.; Gugger, P.F.; Sherman, R.M.; Stevens, K.; Langley, C.; Pellegrini, M.; Salzberg, S. First draft assembly and annotation of the genome of a California Endemic oak Quercus lobata Née (Fagaceae). G3 Genes|Genomes|Genet.; 2016; 6, pp. 3485-3495. [DOI: https://dx.doi.org/10.1534/g3.116.030411]
69. Ramos, A.M.; Usie, A.; Barbosa, P.; Barros, P.M.; Capote, T.; Chaves, I.; Simoes, F.; Abreu, I.; Carrasquinho, I.; Faro, C. et al. The draft genome sequence of cork oak. Sci. Data; 2018; 5, 180069. [DOI: https://dx.doi.org/10.1038/sdata.2018.69]
70. Ai, W.; Liu, Y.; Mei, M.; Zhang, X.; Tan, E.; Liu, H.; Han, X.; Zhan, H.; Lu, X. A chromosome-scale genome assembly of the Mongolian oak (Quercus mongolica). Mol. Ecol. Resour.; 2022; 22, pp. 2396-2410. [DOI: https://dx.doi.org/10.1111/1755-0998.13616]
71. Zoldoš, V.; Papeš, D.; Brown, S.; Panaud, O.; Šiljak-Yakovlev, S. Genome size and base composition of seven Quercus species: Inter- and intra-population variation. Genome; 1998; 41, pp. 162-168. [DOI: https://dx.doi.org/10.1139/g98-006]
72. Saleh, D.; Chen, J.; Leplé, J.C.; Leroy, T.; Truffaut, L.; Dencausse, B.; Lalanne, C.; Labadie, K.; Lesur, I.; Bert, D. et al. Genome-wide evolutionary response of European oaks during the Anthropocene. Evol. Lett.; 2022; 6, pp. 4-20. [DOI: https://dx.doi.org/10.1002/evl3.269]
73. Fitzek, E.; Delcamp, A.; Guichoux, E.; Hahn, M.; Lobdell, M.; Hipp, A.L. A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting. Ecol. Evol.; 2018; 8, pp. 5837-5851. [DOI: https://dx.doi.org/10.1002/ece3.4122]
74. Kakoulidou, I.; Avramidou, E.V.; Baránek, M.; Brunel-Muguet, S.; Farrona, S.; Johannes, F.; Kaiserli, E.; Lieberman-Lazarovich, M.; Martinelli, F.; Mladenov, V. et al. Epigenetics for crop improvement in times of global change. Biology; 2021; 10, 766. [DOI: https://dx.doi.org/10.3390/biology10080766] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34439998]
75. Spillane, C.; McKeown, P. Plant Epigenetics and Epigenomics Methods and Protocols; Methods in Molecular Biology Humana: New York, NY, USA, 2020; Volume 2093.
76. Amaral, J.; Ribeyre, Z.; Vigneaud, J.; Sow, M.D.; Fichot, R.; Messier, C.; Pinto, G.; Nolet, P.; Maury, S. Advances and promises of epigenetics for forest trees. Forests; 2020; 11, 976. [DOI: https://dx.doi.org/10.3390/f11090976]
77. Browne, L.; MacDonald, B.; Fitz-Gibbon, S.; Wright, J.W.; Sork, V.L. Genome-wide variation in DNA methylation predicts variation in leaf traits in an ecosystem-foundational oak species. Forests; 2021; 12, 569. [DOI: https://dx.doi.org/10.3390/f12050569]
78. García-García, I.; Méndez-Cea, B.; Martín-Gálvez, D.; Seco, J.I.; Gallego, F.J.; Linares, J.C. Challenges and perspectives in the epigenetics of climate change-induced forests decline. Front. Plant Sci.; 2022; 12, 797958. [DOI: https://dx.doi.org/10.3389/fpls.2021.797958]
79. Platt, A.; Gugger, P.F.; Pellegrini, M.; Sork, V.L. Genome-wide signature of local adaptation linked to variable CpG methylation in oak populations. Mol. Ecol.; 2015; 24, pp. 3823-3830. [DOI: https://dx.doi.org/10.1111/mec.13230]
80. Rico, L.; Ogaya, R.; Barbeta, A.; Penuelas, J. Changes in DNA methylation fingerprint of Quercus ilex trees in response to experimental field drought simulating projected climate change. Plant Biol.; 2014; 16, pp. 419-427. [DOI: https://dx.doi.org/10.1111/plb.12049]
81. Chwialkowska, K.; Korotko, U.; Kosinska, J.; Szarejko, I.; Kwasniewski, M. Methylation sensitive amplification polymorphism sequencing (MSAP-Seq)—A method for high-throughput analysis of differentially methylated CCGG sites in plants with large genomes. Front. Plant Sci.; 2017; 8, 2056. [DOI: https://dx.doi.org/10.3389/fpls.2017.02056]
82. Guarino, F.; Heinze, B.; Castiglione, S.; Cicatelli, A. Epigenetic analysis through MSAP-NGS coupled technology: The case study of white poplar monoclonal populations/stands. Int. J. Mol. Sci.; 2020; 21, 7393. [DOI: https://dx.doi.org/10.3390/ijms21197393]
83. Ueno, S.; Le Provost, G.; Léger, V.; Klopp, C.; Noirot, C.; Frigerio, J.-M.; Salin, F.; Salse, J.; Abrouk, M.; Murat, F. et al. Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: Oak. BMC Genom.; 2010; 11, 650. [DOI: https://dx.doi.org/10.1186/1471-2164-11-650]
84. Lim, P.K.; Zheng, X.; Goh, J.C.; Mutwil, M. Exploiting plant transcriptomic databases: Resources, tools, and approaches. Plant Comm.; 2022; 3, 100323. [DOI: https://dx.doi.org/10.1016/j.xplc.2022.100323]
85. Meyers, B.C.; Galbraith, D.W.; Nelson, T.; Agrawal, V. Methods for transcriptional profiling in plants. Be fruitful and replicate. Plant Physiol.; 2004; 135, pp. 637-652. [DOI: https://dx.doi.org/10.1104/pp.104.040840]
86. Lowe, R.; Shirley, N.; Bleackley, M.; Dolan, S.; Shafee, T. Transcriptomics technologies. PLoS Comput. Biol.; 2017; 13, e1005457. [DOI: https://dx.doi.org/10.1371/journal.pcbi.1005457]
87. Ling, Y.; Alshareef, S.; Butt, H.; Lozano-Juste, J.; Li, L.; Galal, A.A.; Moustafa, A.; Momin, A.A.; Tashkandi, M.; Richardson, D.N. et al. Pre-mRNA splicing repression triggers abiotic stress signaling in plants. Plant J.; 2017; 89, pp. 291-309. [DOI: https://dx.doi.org/10.1111/tpj.13383]
88. Marum, L.; Miguel, A.; Ricardo, C.P.; Miguel, C. Reference gene selection for quantitative real-time PCR normalization in Quercus suber. PLoS ONE; 2012; 7, e35113. [DOI: https://dx.doi.org/10.1371/annotation/13c5a136-9db4-43a9-aad3-f73acb064d0a]
89. Guerrero-Sánchez, V.M.; Castillejo, M.Á.; López-Hidalgo, C.; Maldonado-Alconada, A.M.; Jorrín-Novo, J.V.; Rey, M.D. Changes in the transcript and protein profiles of Quercus ilex seedlings in response to drought stress. J. Proteom.; 2021; 243, 104263. [DOI: https://dx.doi.org/10.1016/j.jprot.2021.104263]
90. Gugger, P.F.; Peñaloza-Ramírez, J.M.; Wright, J.W.; Sork, V.L. Whole-transcriptome response to water stress in a California endemic oak, Quercus lobata. Tree Physiol.; 2017; 37, pp. 632-644. [DOI: https://dx.doi.org/10.1093/treephys/tpw122]
91. Natali, L.; Vangelisti, A.; Guidi, L.; Remorini, D.; Cotrozzi, L.; Lorenzini, G.; Nali, C.; Pellegrini, E.; Trivellini, A.; Vernieri, P. et al. How Quercus ilex L. saplings face combined salt and ozone stress: A transcriptome analysis. BMC Genom.; 2018; 19, 872. [DOI: https://dx.doi.org/10.1186/s12864-018-5260-2]
92. Mevy, J.P.; Loriod, B.; Liu, X.; Corre, E.; Torres, M.; Büttner, M.; Haguenauer, A.; Reiter, I.M.; Fernandez, C.; Gauquelin, T. Response of Downy Oak (Quercus pubescens Willd.) to climate change: Transcriptome assembly, differential gene analysis and targeted metabolomics. Plants; 2020; 9, 1149. [DOI: https://dx.doi.org/10.3390/plants9091149] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32899727]
93. Soltani, N.; Best, T.; Grace, D.; Nelms, C.; Shumaker, K.; Romero-Severson, J.; Moses, D.; Schuster, S.; Staton, M.; Carlson, J. et al. Transcriptome profiles of Quercus rubra responding to increased O3 stress. BMC Genom.; 2020; 21, 160. [DOI: https://dx.doi.org/10.1186/s12864-020-6549-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32059640]
94. Madritsch, S.; Wischnitzki, E.; Kotrade, P.; Ashoub, A.; Burg, A.; Fluch, S.; Brüggemann, W.; Sehr, E.M. Elucidating drought stress tolerance in European oaks through cross-species transcriptomics. G3 Genes|Genomes|Genet.; 2019; 9, pp. 3181-3199. [DOI: https://dx.doi.org/10.1534/g3.119.400456] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31395652]
95. Oney-Birol, S.; Fitz-Gibbon, S.; Chen, J.M.; Gugger, P.F.; Sork, V.L. Assessment of shared alleles in drought-associated candidate genes among southern California white oak species (Quercus sect. Quercus). BMC Genet.; 2018; 19, 88. [DOI: https://dx.doi.org/10.1186/s12863-018-0677-9]
96. Bouffaud, M.L.; Herrmann, S.; Tarkka, M.T.; Bönn, M.; Feldhahn, L.; Buscot, F. Oak displays common local but specific distant gene regulation responses to different mycorrhizal fungi. BMC Genom.; 2020; 21, 399. [DOI: https://dx.doi.org/10.1186/s12864-020-06806-5]
97. Sebastiana, M.; Gargall-Garriga, A.; Sardans, J.; Pérez-Trujillo, M.; Monteiro, F.; Figueiredo, A.; Maia, M.; Nascimento, R.; Silva, M.S.; Ferreira, A.N. et al. Metabolomics and transcriptomics to decipher molecular mechanisms underlying ectomycorrhizal root colonization of an oak tree. Sci Rep.; 2021; 11, 8576. [DOI: https://dx.doi.org/10.1038/s41598-021-87886-5]
98. Martinson, E.O.; Werren, J.H.; Egan, S.P. Tissue-specific gene expression shows a cynipid wasp repurposes oak host gene networks to create a complex and novel parasite-specific organ. Mol. Ecol.; 2022; 31, pp. 3228-3240. [DOI: https://dx.doi.org/10.1111/mec.16159]
99. Cokus, S.J.; Gugger, P.F.; Sork, V.L. Evolutionary insights from de novo transcriptome assembly and SNP discovery in California white oaks. BMC Genom.; 2015; 16, 552. [DOI: https://dx.doi.org/10.1186/s12864-015-1761-4]
100. Lopes, S.T.; Sobral, D.; Costa, B.; Perdiguero, P.; Chaves, I.; Costa, A.; Miguel, C.M. Phellem versus xylem: Genome-wide transcriptomic analysis reveals novel regulators of cork formation in cork oak. Tree Physiol.; 2020; 40, pp. 129-141. [DOI: https://dx.doi.org/10.1093/treephys/tpz118]
101. Leal, A.R.; Sapeta, H.; Beeckman, T.; Barros, P.M.; Oliveira, M.M. Spatiotemporal development of suberized barriers in cork oak taproots. Tree Physiol.; 2022; 42, pp. 1269-1285. [DOI: https://dx.doi.org/10.1093/treephys/tpab176]
102. Li, D.; Li, Y.; Qian, J.; Liu, X.; Xu, H.; Zhang, G.; Ren, J.; Wang, L.; Zhang, L.; Yu, H. Comparative transcriptome analysis revealed candidate genes potentially related to desiccation sensitivity of recalcitrant Quercus variabilis seeds. Front Plant Sci.; 2021; 12, 717563. [DOI: https://dx.doi.org/10.3389/fpls.2021.717563]
103. Guerrero-Sanchez, V.M.; Maldonado-Alconada, A.M.; Amil-Ruiz, F.; Jorrín-Novo, J.V. Holm oak (Quercus ilex) transcriptome. De novo sequencing and assembly analysis. Front. Mol. Biosci.; 2017; 4, 70.
104. Guerrero-Sanchez, V.M.; Maldonado-Alconada, A.M.; Amil-Ruiz, F.; Verardi, A.; Jorrín-Novo, J.V.; Rey, M.-D. Ion Torrent and lllumina, two complementary RNA-seq platforms for constructing the Holm oak (Quercus ilex) transcriptome. PLoS ONE; 2019; 14, e0210356. [DOI: https://dx.doi.org/10.1371/journal.pone.0210356]
105. Romero-Rodríguez, M.C.; Abril, N.; Sánchez-Lucas, R.; Jorrín-Novo, J.V. Multiplex staining of 2-DE gels for an initial phosphoproteome analysis of germinating seeds and early grown seedlings from a non-orthodox specie: Quercus ilex L. subsp. ballota [Desf.]. Samp. Front. Plant Sci.; 2015; 6, 620. [DOI: https://dx.doi.org/10.3389/fpls.2015.00620]
106. Jorrín-Novo, J.V. Proteomics and plant biology: Contributions to date and a look towards the next decade. Expert Rev. Proteom.; 2021; 18, pp. 93-103. [DOI: https://dx.doi.org/10.1080/14789450.2021.1910028]
107. Canovas, F.M.; Dumas-Gaudot, E.; Recorbet, G.; Jorrón-Novo, J.V.; Mock, H.P.; Rossignol, M. Plant proteome analysis. Proteomics; 2004; 4, pp. 285-298. [DOI: https://dx.doi.org/10.1002/pmic.200300602]
108. Jorrín-Novo, J.V.; Maldonado, A.M.; Echevarria-Zomeno, S.; Valledor, L.; Castillejo, M.A.; Curto, M.; Valero-Galván, J.; Sghaier-Hammami, B.; Donoso, G.; Redondo, I. Plant proteomics update (2007–2008): Second-generation proteomic techniques, an appropriate experimental design, and data analysis to fulfill MIAPE standards, increase plant proteome coverage and expand biological knowledge. J. Proteom.; 2009; 72, pp. 285-314. [DOI: https://dx.doi.org/10.1016/j.jprot.2009.01.026] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19367730]
109. Jorrín-Novo, J.V.; Pascual, J.; Sánchez-Lucas, R.; Romero-Rodríguez, M.C.; Rodríguez-Ortega, M.J.; Lenz, C.; Valledor, L. Fourteen years of plant proteomics reflected in Proteomics: Moving from model species and 2DE-based approaches to orphan species and gel-free platforms. Proteomics; 2015; 15, pp. 1089-1112. [DOI: https://dx.doi.org/10.1002/pmic.201400349] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25487722]
110. Jorrín-Novo, J.V. What is new in (plant) proteomics methods and protocols: The 2015–2019 Quinquennium. Plant Proteomics; Humana: New York, NY, USA, 2020; pp. 1-10.
111. Rey, M.D.; Valledor, L.; Castillejo, M.Á.; Sánchez-Lucas, R.; López-Hidalgo, C.; Guerrero-Sánchez, V.M.; Colina, F.J.; Escandón, M.; Maldonado-Alconada, A.M.; Jorrín-Novo, J.V. Recent advances in MS-based plant proteomics: Proteomics data validation through integration with other classic and -omics approaches. Progress in Botany; Cánovas, F.; Lüttge, U.; Leuschner, C.; Risueño, M. Springer: Cham, Switzerland, 2019; Volume 81, ISBN 978-3-030-36327-7
112. Jorrín-Novo, J.V.; Navarro-Cerrillo, R.M. Variabilidad y respuesta a distintos estreses en poblaciones de encina (Quercus ilex L.) en Andalucía mediante una aproximación proteómica. Ecosistemas; 2014; 23, pp. 99-107.
113. Pinheiro, C.; Guerra-Guimaraes, L.; David, T.S.; Vieira, A. Proteomics: State of the art to study Mediterranean woody species under stress. Environ. Exp. Bot.; 2014; 103, pp. 117-127. [DOI: https://dx.doi.org/10.1016/j.envexpbot.2014.01.010]
114. Sergeant, K.; Spiess, N.; Renaut, J.; Wilhelm, E.; Hausman, J.F. One dry summer: A leaf proteome study on the response of oak to drought exposure. J. Proteom.; 2011; 74, pp. 1385-1395. [DOI: https://dx.doi.org/10.1016/j.jprot.2011.03.011]
115. Gómez, A.; López, J.A.; Pintos, B.; Camafeita, E.; Bueno, M.A. Proteomic analysis from haploid and diploid embryos of Quercus suber L. identifies qualitative and quantitative differential expression patterns. Proteomics; 2009; 9, pp. 4355-4367. [DOI: https://dx.doi.org/10.1002/pmic.200900179]
116. Gomez-Garay, A.; Lopez, J.A.; Camafeita, E.; Bueno, M.A.; Pintos, B. Proteomic perspective of Quercus suber somatic embryogenesis. J. Proteom.; 2013; 93, pp. 314-325. [DOI: https://dx.doi.org/10.1016/j.jprot.2013.06.006]
117. Coelho, A.C.; Pires, R.; Schütz, G.; Santa, C.; Manadas, B.; Pinto, P. Disclosing proteins in the leaves of cork oak plants associated with the immune response to Phytophthora cinnamomi inoculation in the roots: A long-term proteomics approach. PLoS ONE; 2021; 16, e0245148. [DOI: https://dx.doi.org/10.1371/journal.pone.0245148]
118. Ricardo, C.P.; Martins, I.; Francisco, R.; Sergeant, K.; Pinheiro, C.; Campos, A.; Renaut, J.; Fevereiro, P. Proteins associated with cork formation in Quercus suber L. stem tissues. J. Proteom.; 2011; 74, pp. 1266-1278. [DOI: https://dx.doi.org/10.1016/j.jprot.2011.02.003]
119. Pinheiro, C.; Wienkoop, S.; de Almeida, J.F.; Brunetti, C.; Zarrouk, O.; Planchon, S.; Gori, A.; Tattini, M.; Ricardo, C.P.; Renaut, J. et al. Phellem cell-wall components are discriminants of cork quality in Quercus suber. Front. Plant Sci.; 2019; 10, 944. [DOI: https://dx.doi.org/10.3389/fpls.2019.00944]
120. Escandón, M.; Jorrín-Novo, J.V.; Castillejo, M.Á. Application and optimization of label-free shotgun approaches in the study of Quercus ilex. J. Proteom.; 2020; 233, 104082. [DOI: https://dx.doi.org/10.1016/j.jprot.2020.104082]
121. Romero-Rodríguez, M.C.; Jorrín-Novo, J.V.; Castillejo, M.Á. Towards characterizing seed germination and seedling establishment in the non- orthodox forest tree species Quercus ilex through complementary gel and gel-free proteomic approaches. J. Proteom.; 2019; 197, pp. 60-70. [DOI: https://dx.doi.org/10.1016/j.jprot.2018.11.003]
122. Romero-Rodríguez, M.C.; Pascual, J.; Valledor, L.; Jorrín-Novo, J. Improving the quality of protein identification in non-model species. Characterization of Quercus ilex seed and Pinus radiata needle proteomes by using SEQUEST and custom databases. J. Proteom.; 2014; 105, pp. 85-91. [DOI: https://dx.doi.org/10.1016/j.jprot.2014.01.027]
123. Jorge, I.; Navarro, R.M.; Lenz, C.; Ariza, D.; Porras, C.; Jorrín-Novo, J. The Holm oak leaf proteome: Analytical and biological variability in the protein expression level assessed by 2-DE and protein identification tandem mass spectrometry de novo sequencing and sequence similarity searching. Proteomics; 2005; 5, pp. 222-234. [DOI: https://dx.doi.org/10.1002/pmic.200400893]
124. Jorge, I.; Navarro, R.M.; Lenz, C.; Ariza, D.; Jorrín-Novo, J.V. Variation in the holm oak leaf proteome at different plant developmental stages, between provenances and in response to drought stress. Proteomics; 2006; 6, pp. S207-S214. [DOI: https://dx.doi.org/10.1002/pmic.200500364]
125. Valero-Galván, J.; Valledor, L.; Fernandez, R.G.; Cerrillo, R.M.N.; Jorrín-Novo, J.V. Proteomic analysis of Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.) pollen. J. Proteom.; 2012; 75, pp. 2736-2744.
126. Simova-Stoilova, L.P.; Romero-Rodríguez, M.C.; Sánchez-Lucas, R.; Navarro-Cerrillo, R.M.; Medina-Aunon, J.A.; Jorrín-Novo, J.V. 2-DE proteomics analysis of drought treated seedlings of Quercus ilex supports a root active strategy for metabolic adaptation in response to water shortage. Front. Plant Sci.; 2015; 6, 627. [DOI: https://dx.doi.org/10.3389/fpls.2015.00627]
127. Simova-Stoilova, L.P.; López-Hidalgo, C.; Sanchez-Lucas, R.; Valero-Galvan, J.; Romero-Rodríguez, C.; Jorrin-Novo, J.V. Holm oak proteomic response to water limitation at seedling establishment stage reveals specific changes in different plant parts as well as interaction between roots and cotyledons. Plant Sci.; 2018; 276, pp. 1-13. [DOI: https://dx.doi.org/10.1016/j.plantsci.2018.07.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30348307]
128. San-Eufrasio, B.; Bigatton, E.D.; Guerrero-Sánchez, V.M.; Chaturvedi, P.; Jorrín-Novo, J.V.; Rey, M.D.; Castillejo, M.Á. Proteomics data analysis for the identification of proteins and derived proteotypic peptides of potential use as putative drought tolerance markers for Quercus ilex. Int. J. Mol. Sci.; 2021; 22, 3191. [DOI: https://dx.doi.org/10.3390/ijms22063191] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33800973]
129. Escandón, M.; Bigatton, E.D.; Guerrero-Sánchez, V.M.; Hernández-Lao, T.; Rey, M.D.; Jorrín-Novo, J.V.; Castillejo, M.A. Identification of proteases and protease inhibitors in seeds of the recalcitrant forest tree species Quercus ilex. Front. Plant Sci.; 2022; 13, 907042. [DOI: https://dx.doi.org/10.3389/fpls.2022.907042] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35832232]
130. Pedrosa, M.; Guerrero-Sanchez, V.M.; Canales-Bueno, N.; Loli-Ausejo, D.; Castillejo, M.Á.; Quirce, S.; Jorrin-Novo, J.V.; Rodriguez-Perez, R. Quercus ilex pollen allergen, Que i 1, responsible for pollen food allergy syndrome caused by fruits in Spanish allergic patients. Clin. Exp. Allergy; 2020; 50, pp. 815-823. [DOI: https://dx.doi.org/10.1111/cea.13679]
131. Saadi, S.; Saari, N.; Anwar, F.; Abdul Hamid, A.; Ghazali, H.M. Recent advances in food biopeptides: Production, biological functionalities, and therapeutic applications. Biotechnol. Adv.; 2015; 33, pp. 80-116. [DOI: https://dx.doi.org/10.1016/j.biotechadv.2014.12.003]
132. Misra, B.B. New software tools, databases, and resources in metabolomics: Updates from 2020. Metabolomics; 2021; 17, 49. [DOI: https://dx.doi.org/10.1007/s11306-021-01796-1]
133. Hong, J.; Yang, L.; Zhang, D.; Shi, J. Plant metabolomics: An indispensable Systems Biology tool for plant science. Int. J. Mol. Sci.; 2016; 17, 767. [DOI: https://dx.doi.org/10.3390/ijms17060767]
134. Kumar, R.; Bohra, A.; Pandey, A.K.; Pandey, M.K.; Kumar, A. Metabolomics for plant improvement: Status and prospects. Front. Plant Sci.; 2017; 8, 1302. [DOI: https://dx.doi.org/10.3389/fpls.2017.01302]
135. Rodrigues, A.M.; Ribeiro-Barros, A.I.; António, C. Experimental design and sample preparation in forest tree metabolomics. Metabolites; 2019; 9, 285. [DOI: https://dx.doi.org/10.3390/metabo9120285]
136. Tsugawa, H.; Rai, A.; Saito, K.; Nakabayashi, R. Metabolomics and complementary techniques to investigate the plant phytochemical cosmos. Nat. Prod. Rep.; 2021; 38, pp. 1729-1759. [DOI: https://dx.doi.org/10.1039/D1NP00014D]
137. Patel, M.K.; Pandey, S.; Kumar, M.; Haque, M.I.; Pal, S.; Yadav, N.S. Plants metabolome study: Emerging tools and techniques. Plants; 2021; 20, 2409. [DOI: https://dx.doi.org/10.3390/plants10112409]
138. Misra, B.B.; Chen, S.X. Advances in understanding CO2 responsive plant metabolomes in the era of climate change. Metabolomics; 2015; 11, pp. 1478-1491. [DOI: https://dx.doi.org/10.1007/s11306-015-0825-4]
139. Rodrigues, A.M.; Miguel, C.; Chaves, I.; António, C. Mass spectrometry-based forest tree metabolomics. Mass Spectrom. Rev.; 2021; 40, pp. 126-157. [DOI: https://dx.doi.org/10.1002/mas.21603]
140. Wetzel, W.C.; Whitehead, S.R. The many dimensions of phytochemical diversity: Linking theory to practice. Ecol. Lett.; 2020; 23, pp. 16-32. [DOI: https://dx.doi.org/10.1111/ele.13422]
141. Buche, G.; Colas, C.; Fougère, L.; Giordanengo, T.; Destandau, E. Untargeted UHPLC-Q-TOF-HRMS based determination of discrimating compounds for oak species Quercus robur L. and Quercus petraea Liebl. identification. Phytochem. Anal.; 2021; 32, pp. 660-671. [DOI: https://dx.doi.org/10.1002/pca.3013]
142. Kersten, B.; Ghirardo, A.; Schnitzler, J.P.; Kanawati, B.; Schmitt-Kopplin, P.; Fladung, M.; Schroeder, H. Integrated transcriptomics and metabolomics decipher differences in the resistance of pedunculate oak to the herbivore Tortrix viridana L. BMC Genom.; 2013; 14, 737. [DOI: https://dx.doi.org/10.1186/1471-2164-14-737]
143. Bertić, M.; Schroeder, H.; Kersten, B.; Fladung, M.; Orgel, F.; Buegger, F.; Schnitzler, J.P.; Ghirardo, A. European oak chemical diversity—from ecotypes to herbivore resistance. New Phytol.; 2021; 232, pp. 818-834. [DOI: https://dx.doi.org/10.1111/nph.17608]
144. Almeida, T.; Pinto, G.; Correia, B.; Gonçalves, S.; Meijón, M.; Escandón, M. In-depth analysis of the Quercus suber metabolome under drought stress and recovery reveals potential key metabolic players. Plant Sci.; 2020; 299, 110606. [DOI: https://dx.doi.org/10.1016/j.plantsci.2020.110606]
145. Suseela, V.; Tharayil, N.; Xing, B.; Dukes, J.S. Warming and drought differentially influence the production and resorption of elemental and metabolic nitrogen pools in Quercus rubra. Glob. Change Biol.; 2015; 21, pp. 4177-4195. [DOI: https://dx.doi.org/10.1111/gcb.13033]
146. Conrad, A.O.; McPherson, B.A.; Wood, D.L.; Madden, L.V.; Bonello, P. Constitutive phenolic biomarkers identify naïve Quercus agrifolia resistant to Phytophthora ramorum, the causal agent of sudden oak death. Tree Physiol.; 2017; 37, pp. 1686-1696. [DOI: https://dx.doi.org/10.1093/treephys/tpx116] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29036534]
147. Li, Q.; Yan, L.; Ye, L.; Zhou, J.; Zhang, B.; Peng, W.; Zhang, X.; Li, X. Chinese black truffle (Tuber indicum) alters the ectomycorrhizosphere and endoectomycosphere microbiome and metabolic profiles of the host tree Quercus aliena. Front. Microbiol.; 2018; 9, 2202. [DOI: https://dx.doi.org/10.3389/fmicb.2018.02202] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30283422]
148. Rivas-Ubach, A.; Pérez-Trujillo, M.; Sardans, J.; Gargallo-Garriga, A.; Parella, T.; Penuelas, J. Ecometabolomics: Optimized NMR-based method. Methods Ecol. Evol.; 2013; 4, pp. 464-473. [DOI: https://dx.doi.org/10.1111/2041-210X.12028]
149. Sardans, J.; Gargallo-Garriga, A.; Pérez-Trujillo, M.; Parella, T.J.; Seco, R.; Filella, I.; Penuelas, J. Metabolic responses of Quercus ilex seedlings to wounding analysed with nuclear magnetic resonance profiling. Plant Biol.; 2014; 16, pp. 395-403. [DOI: https://dx.doi.org/10.1111/plb.12032] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23590498]
150. Rivas-Ubach, A.; Gargallo-Garriga, A.; Sardans, J.; Oravec, M.; Mateu-Castell, L.; Pérez-Trujillo, M.; Ogaya, R.; Urban, P.; Peñuelas, J. Drought enhances folivory by shifting foliar metabolomes in Quercus ilex trees. New Phytol.; 2014; 202, pp. 874-885. [DOI: https://dx.doi.org/10.1111/nph.12687]
151. Rodríguez-Calcerrada, J.; Rodrigues, A.M.; Perdiguero, P.; António, C.; Atkin, O.K.; Li, M.; Collada, C.; Gil, L. A molecular approach to drought-induced reduction in leaf CO2 exchange in drought-resistant Quercus ilex. Physiol. Plant.; 2018; 162, pp. 394-408. [DOI: https://dx.doi.org/10.1111/ppl.12649]
152. Rivas-Ubach, A.; Barbeta, A.; Sardans, J.; Guenther, A.; Ogaya, R.; Oravec, M.; Urban, O.; Penuelas, J. Topsoil depth substantially influences the responses to drought of the foliar metabolomes of Mediterranean forests. Perspect. Plant Ecol. Evol. Syst.; 2016; 21, pp. 41-54. [DOI: https://dx.doi.org/10.1016/j.ppees.2016.06.001]
153. Tienda-Parrilla, M.; López-Hidalgo, C.; Guerrero-Sanchez, V.M.; Infantes-González, Á.; Valderrama-Fernández, R.; Castillejo, M.Á.; Jorrín-Novo, J.V.; Rey, M.D. Untargeted MS-based metabolomics analysis of the responses to drought stress in Quercus ilex L. leaf seedlings and the identification of putative compounds related to tolerance. Forests; 2022; 13, 551. [DOI: https://dx.doi.org/10.3390/f13040551]
154. López-Hidalgo, C.; Guerrero-Sanchez, V.M.; Gomez-Galvez, I.; Sánchez-Lucas, R.; Castillejo-Sánchez, M.A.; Maldonado-Alconada, A.M.; Valledor, L.; Jorrín-Novo, J.V. A multi-omics analysis pipeline for the metabolic pathway reconstruction in the orphan species Quercus ilex. Front. Plant Sci.; 2018; 9, 935. [DOI: https://dx.doi.org/10.3389/fpls.2018.00935]
155. García-García, A.L.; García-Machado, F.J.; Borges, A.A.; Morales-Sierra, S.; Boto, A.; Jiménez-Arias, D. Pure organic active compounds against abiotic stress: A biostimulant overview. Front. Plant Sci.; 2020; 11, 575829. [DOI: https://dx.doi.org/10.3389/fpls.2020.575829]
156. Marthandan, V.; Geetha, R.; Kumutha, K.; Renganathan, V.G.; Karthikeyan, A.; Ramalingam, J. Seed priming: A feasible strategy to enhance drought tolerance in crop plants. Int. J. Mol. Sci.; 2020; 21, 8258. [DOI: https://dx.doi.org/10.3390/ijms21218258]
157. Mauch-Mani, B.; Baccelli, I.; Luna, E.; Flors, V. Defense Priming: An Adaptive Part of Induced Resistance. Annu. Rev. Plant Biol.; 2017; 68, pp. 485-512. [DOI: https://dx.doi.org/10.1146/annurev-arplant-042916-041132] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28226238]
158. Martínez, M.T.; San-José, M.D.C.; Arrillaga, I.; Cano, V.; Morcillo, M.; Cernadas, M.J.; Corredoira, E. Holm oak somatic embryogenesis: Current status and future perspectives. Front. Plant Sci.; 2019; 10, 239. [DOI: https://dx.doi.org/10.3389/fpls.2019.00239]
159. Martínez, M.T.; Arrillaga, I.; Sales, E.; Pérez-Oliver, M.A.; González-Mas, M.d.C.; Corredoira, E. Micropropagation, Characterization, and conservation of Phytophthora cinnamomi-tolerant holm oak mature trees. Forests; 2021; 12, 1634. [DOI: https://dx.doi.org/10.3390/f12121634]
160. Singh, A.; Shannon, C.P.; Gautier, B.; Rohart, F.; Vacher, M.; Tebbutt, S.J.; Lê Cao, K.A. DIABLO: An integrative approach for identifying key molecular drivers from multi-omics assays. Bioinformatics; 2019; 35, pp. 3055-3062. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty1054]
161. Stocks, J.J.; Metheringham, C.L.; Plumb, W.J.; Lee, S.J.; Kelly, L.J.; Nichols, R.A.; Buggs, R.J.A. Genomic basis of European ash tree resistance to ash dieback fungus. Nat. Ecol. Evol.; 2019; 3, pp. 1686-1696. [DOI: https://dx.doi.org/10.1038/s41559-019-1036-6]
162. Du, Q.; Lu, W.; Quan, M.; Xiao, L.; Song, F.; Li, P.; Zhou, D.; Xie, J.; Wang, L.; Zhang, D. Genome-wide association studies to improve wood properties: Challenges and prospects. Front. Plant Sci.; 2018; 9, 1912. [DOI: https://dx.doi.org/10.3389/fpls.2018.01912]
163. Espinas, N.A.; Saze, H.; Saijo, Y. Epigenetic control of defense signaling and priming in plants. Front. Plant Sci.; 2016; 7, 1201. [DOI: https://dx.doi.org/10.3389/fpls.2016.01201]
164. Kumar, S.; Fladung, M. Forest tree transgenesis and functional genomics: From fast forward to reverse genetics. Silvae Genet.; 2003; 52, pp. 229-232.
165. Sobral, R.; Costa, M.M.R. Role of floral organ identity genes in the development of unisexual flowers of Quercus suber L. Sci. Rep.; 2017; 7, 10368. [DOI: https://dx.doi.org/10.1038/s41598-017-10732-0]
166. Magalhães, A.P.; Verde, N.; Reis, F.; Martins, I.; Costa, D.; Lino-Neto, T.; Castro, P.H.; Tavares, R.M.; Azevedo, H. RNA-seq and gene network analysis uncover activation of an ABA-dependent signalosome during the cork oak root response to drought. Front. Plant Sci.; 2016; 6, 1195. [DOI: https://dx.doi.org/10.3389/fpls.2015.01195]
167. Plomión, C.; Bastien, C.; Bogeat-Triboulot, M.B.; Bouffier, L.; Déjardin, A.; Duplessis, S.; Fady, B.; Heuertz, M.; Le Gac, A.-L.; Le Provost, G. et al. Forest tree genomics: 10 achievements from the past 10 years and future prospects. Ann. For. Sci.; 2016; 73, pp. 77-103. [DOI: https://dx.doi.org/10.1007/s13595-015-0488-3]
168. Tuskan, G.A.; Groover, A.T.; Schmutz, J.; DiFazio, S.P.; Myburg, A.; Grattapaglia, D.; Smart, L.B.; Yin, T.; Aury, J.-M.; Kremer, A. et al. Hardwood Tree Genomics: Unlocking Woody Plant Biology. Front. Plant Sci.; 2018; 9, 1799. [DOI: https://dx.doi.org/10.3389/fpls.2018.01799] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30619389]
169. Beck, P.; Caudullo, G.; Mauri, A.; De Rigo, D.; Durrant, T.; San-Miguel-Ayanz, J. Tree Species Distribution Data and Maps for Europe; EUR 30148 EN Publications Office of the European Union: Luxembourg, 2020.
170. Herrera-Arroyo, M.L.; Sork, V.L.; González-Rodríguez, A.; Rocha-Ramírez, V.; Vega, E.; Oyama, K. Seed-mediated connectivity among fragmented populations of Quercus castanea (Fagaceae) in a Mexican landscape. Am. J. Bot.; 2013; 100, pp. 1663-1671. [DOI: https://dx.doi.org/10.3732/ajb.1200396]
171. Ueno, S.; Taguchi, Y.; Tsumura, Y. Microsatellite markers derived from Quercus mongolica var. crispula (Fagaceae) inner bark expressed sequence tags. Genes Genet. Syst.; 2008; 83, pp. 179-187.
172. Guichoux, E.; Lagache, L.; Wagner, S.; Léger, P.; Petit, R.J. Two highly validated multiplexes (12-plex and 8-plex) for species delimitation and parentage analysis in oaks (Quercus spp.). Mol. Ecol. Resour.; 2011; 11, pp. 578-585. [DOI: https://dx.doi.org/10.1111/j.1755-0998.2011.02983.x]
173. Di Pietro, R.; Conte, A.L.; Di Marzio, P.; Fortini, P.; Farris, E.; Gianguzzi, L.; Müller, M.; Rosati, L.; Spampinato, G.; Gailing, O. Does the genetic diversity among pubescent white oaks in southern Italy, Sicily and Sardinia islands support the current taxonomic classification?. Eur. J. For. Res.; 2021; 140, pp. 355-371. [DOI: https://dx.doi.org/10.1007/s10342-020-01334-z]
174. Gugerli, F.; Walser, J.C.; Dounavi, K.; Holderegger, R.; Finkeldey, R. Coincidence of small-scale spatial discontinuities in leaf morphology and nuclear microsatellite variation of Quercus petraea and Q. robur in a mixed forest. Ann. Bot.; 2007; 99, pp. 713-722. [DOI: https://dx.doi.org/10.1093/aob/mcm006]
175. Beatty, G.E.; Montgomery, W.I.; Spaans, F.; Tosh, D.G.; Provan, J. Pure species in a continuum of genetic and morphological variation: Sympatric oaks at the edge of their range. Ann. Bot.; 2016; 117, pp. 541-549. [DOI: https://dx.doi.org/10.1093/aob/mcw002]
176. Gugerli, F.; Brodbeck, S.; Holderegger, R. Utility of multilocus genotypes for taxon assignment in stands of closely related European white Oaks from Switzerland. Ann. Bot.; 2008; 102, pp. 855-863. [DOI: https://dx.doi.org/10.1093/aob/mcn164]
177. Lopes, T.; Pinto, G.; Loureiro, J.; Costa, A.; Santos, C. Determination of genetic stability in long-term somatic embryogenic cultures and derived plantlets of cork oak using microsatellite markers. Tree Physiol.; 2006; 26, pp. 1145-1152. [DOI: https://dx.doi.org/10.1093/treephys/26.9.1145]
178. Castillo-Mendoza, E.; Salinas-Sánchez, D.; Valencia-Cuevas, L.; Zamilpa, A.; Tovar-Sánchez, E. Natural hybridisation among Quercus glabrescens, Q. rugosa and Q. obtusata (Fagaceae): Microsatellites and secondary metabolites markers. Plant Biol.; 2019; 21, pp. 110-121. [DOI: https://dx.doi.org/10.1111/plb.12899]
179. Peñaloza-Ramírez, J.M.; González-Rodríguez, A.; Mendoza-Cuenca, L.; Caron, H.; Kremer, A.; Oyama, K. Interspecific gene flow in a multispecies oak hybrid zone in the Sierra Tarahumara of Mexico. Ann. Bot.; 2010; 105, pp. 389-399. [DOI: https://dx.doi.org/10.1093/aob/mcp301]
180. Craft, K.J.; Ashley, M.V. Pollen-mediated gene flow in isolated and continuous stands of bur oak, Quercus macrocarpa (Fagaceae). Am. J. Bot.; 2010; 97, pp. 1999-2006. [DOI: https://dx.doi.org/10.3732/ajb.0900390]
181. Ramírez-Valiente, J.A.; Etterson, J.R.; Deacon, N.J.; Cavender-Bares, J. Evolutionary potential varies across populations and traits in the neotropical oak Quercus oleoides. Tree Physiol.; 2019; 39, pp. 427-439. [DOI: https://dx.doi.org/10.1093/treephys/tpy108]
182. Backs, J.R.; Ashley, M.V. Evolutionary history and gene flow of an endemic island oak: Quercus pacifica. Am. J. Bot.; 2016; 103, pp. 2115-2125. [DOI: https://dx.doi.org/10.3732/ajb.1600259]
183. Li, Y.; Zhang, X.; Wang, L.; Sork, V.L.; Mao, L.; Fang, Y. Influence of Pliocene and Pleistocene climates on hybridization patterns between two closely related oak species in China. Ann. Bot.; 2022; 129, pp. 231-245. [DOI: https://dx.doi.org/10.1093/aob/mcab140]
184. Li, Y.; Zhang, X.; Fang, Y. Landscape features and climatic forces shape the genetic structure and evolutionary history of an oak species (Quercus chenii) in East China. Front. Plant Sci.; 2019; 10, 1060. [DOI: https://dx.doi.org/10.3389/fpls.2019.01060]
185. Tamaki, I.; Obora, T.; Ohsawa, T.; Matsumoto, A.; Saito, Y.; Ide, Y. Different population size change and migration histories created genetic diversity of three oaks in Tokai region, central Japan. J. Plant Res.; 2021; 134, pp. 933-946. [DOI: https://dx.doi.org/10.1007/s10265-021-01323-2]
186. Degen, B.; Yanbaev, Y.; Mader, M.; Ianbaev, R.; Bakhtina, S.; Schroeder, H.; Blanc-Jolivet, C. Impact of Gene Flow and Introgression on the Range Wide Genetic Structure of Quercus robur (L.) in Europe. Forests.; 2021; 12, 1425. [DOI: https://dx.doi.org/10.3390/f12101425]
187. Cavender-Bares, J.; González-Rodríguez, A.; Eaton, D.A.; Hipp, A.A.; Beulke, A.; Manos, P.S. Phylogeny and biogeography of the American live oaks (Quercus subsection Virentes): A genomic and population genetics approach. Mol. Ecol.; 2015; 24, pp. 3668-3687. [DOI: https://dx.doi.org/10.1111/mec.13269]
188. Ainsworth, E.A.; Tranel, P.J.; Drake, B.G.; Long, S.P. The clonal structure of Quercus geminata revealed by conserved microsatellite loci. Mol. Ecol.; 2003; 12, pp. 527-532. [DOI: https://dx.doi.org/10.1046/j.1365-294X.2003.01749.x]
189. Guo, B.; Hao, X.; Han, L.; Zhai, Y.; Zhou, S.; Chen, S.; Ren, D.; An, X. Unraveling the genetic diversity and structure of Quercus liaotungensis population through analysis of microsatellite markers. PeerJ; 2021; 9, e10922. [DOI: https://dx.doi.org/10.7717/peerj.10922] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33954024]
190. Cao, Y.; Zhang, D.Y.; Zeng, Y.F.; Bai, W.N. Recent demographic histories of temperate deciduous trees inferred from microsatellite markers. BMC Ecol. Evol.; 2021; 21, 88. [DOI: https://dx.doi.org/10.1186/s12862-021-01805-w]
191. Chatwin, W.B.; Carpenter, K.K.; Jimenez, F.R.; Elzinga, D.B.; Johnson, L.A.; Maughan, P.J. Microsatellite primer development for post oak, Quercus stellata (Fagaceae). Appl. Plant Sci.; 2014; 2, 1400070. [DOI: https://dx.doi.org/10.3732/apps.1400070]
192. Wang, X.; Li, J.; Li, Y. Isolation and characterization of microsatellite markers for an endemic tree in East Asia, Quercus variabilis (Fagaceae). Appl. Plant Sci.; 2015; 3, 1500032. [DOI: https://dx.doi.org/10.3732/apps.1500032]
193. Degen, B.; Streiff, R.; Ziegenhagen, B. Comparative study of genetic variation and differentiation of two pedunculate oak (Quercus robur) stands using microsatellite and allozyme loci. Heredity; 1999; 83, pp. 597-603. [DOI: https://dx.doi.org/10.1038/sj.hdy.6886220] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10620033]
194. Zhang, Y.Y.; Fang, Y.M.; Yu, M.K.; Li, X.X.; Xia, T. Molecular characterization and genetic structure of Quercus acutissima germplasm in China using microsatellites. Mol. Biol. Rep.; 2013; 40, pp. 4083-4090. [DOI: https://dx.doi.org/10.1007/s11033-013-2486-6]
195. Browne, L.; Mead, A.; Horn, C.; Chang, K.; Celikkol, Z.A.; Henriquez, C.L.; Ma, F.; Beraut, E.; Meyer, R.S.; Sork, V.L. Experimental DNA demethylation associates with changes in growth and gene expression of oak tree seedlings. G3 Genes|Genomes|Genet.; 2020; 10, pp. 1019-1028. [DOI: https://dx.doi.org/10.1534/g3.119.400770] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31941723]
196. Rodríguez-Sanz, H.; Manzanera, J.A.; Solís, M.T.; Gómez-Garay, A.; Pintos, B.; Risueño, M.C.; Testillano, P.S. Early markers are present in both embryogenesis pathways from microspores and immature zygotic embryos in cork oak, Quercus suber L. BMC Plant Biol.; 2014; 14, 224. [DOI: https://dx.doi.org/10.1186/s12870-014-0224-4]
197. Pérez, M.; Viejo, M.; LaCuesta, M.; Toorop, P.; Cañal, M.J. Epigenetic and hormonal profile during maturation of Quercus suber L. somatic embryos. J. Plant Physiol.; 2015; 173, pp. 51-61. [DOI: https://dx.doi.org/10.1016/j.jplph.2014.07.028] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25462078]
198. Correia, B.; Valledor, L.; Meijón, M.; Rodriguez, J.L.; Dias, M.C.; Santos, C.; Cañal, M.J.; Rodriguez, R.; Pinto, G. Is the interplay between epigenetic markers related to the acclimation of cork oak plants to high temperatures?. PLoS ONE; 2013; 8, e53543.
199. Ribeiro, T.; Viegas, W.; Morais-Cecílio, L. Epigenetic marks in the mature pollen of Quercus suber L.(Fagaceae). Sex. Plant Reprod.; 2009; 22, pp. 1-7. [DOI: https://dx.doi.org/10.1007/s00497-008-0083-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20033450]
200. Vičić, V.; Barišić, D.; Horvat, T.; Biruš, I.; Zoldos, V. Epigenetic characterization of chromatin in cycling cells of pedunculate oak, Quercus robur L. Tree Genet. Genomes.; 2013; 9, pp. 1247-1256. [DOI: https://dx.doi.org/10.1007/s11295-013-0632-x]
201. Inácio, V.; Barros, P.M.; Costa, A.; Roussado, C.; Gonçalves, E.; Costa, R.; Graça, J.; Oliveira, M.M.; Morais-Cecílio, L. Differential DNA methylation patterns are related to phellogen origin and quality of Quercus suber cork. PLoS ONE; 2017; 12, e0169018. [DOI: https://dx.doi.org/10.1371/journal.pone.0169018]
202. Gugger, P.F.; Fitz-Gibbon, S.; PellEgrini, M.; Sork, V.L. Species-wide patterns of DNA methylation variation in Quercus lobata and their association with climate gradients. Mol. Ecol.; 2016; 25, pp. 1665-1680. [DOI: https://dx.doi.org/10.1111/mec.13563]
203. Ramos, M.; Rocheta, M.; Carvalho, L.; Graça, J.; Morais-Cecilio, L. Expression analysis of DNA methyltransferase and co-repressor genes in Quercus suber phellogen: An attempt to correlate with cork quality. BMC Proc.; 2011; 5, P169. [DOI: https://dx.doi.org/10.1186/1753-6561-5-S7-P169]
204. Rossi, F.; Crnjar, A.; Comitani, F.; Feliciano, R.; Jahn, L.; Malim, G.; Southgate, L.; Kay, E.; Oakey, R.; Buggs, R. et al. Extraction and high-throughput sequencing of oak heartwood DNA: Assessing the feasibility of genome-wide DNA methylation profiling. PLoS ONE; 2021; 16, e0254971. [DOI: https://dx.doi.org/10.1371/journal.pone.0254971]
205. Tarkka, M.T.; Herrmann, S.; Wubet, T.; Feldhahn, L.; Recht, S.; Kurth, F.; Mailänder, S.; Bönn, M.; Neef, M.; Angay, O. et al. OakContig DF 159.1, a reference library for studying differential gene expression in Quercus robur during controlled biotic interactions: Use for quantitative transcriptomic profiling of oak roots in ectomycorrhizal symbiosis. New Phytol.; 2013; 199, pp. 529-540. [DOI: https://dx.doi.org/10.1111/nph.12317]
206. Maboreke, H.R.; Feldhahn, L.; Bönn, M.; Tarkka, M.T.; Buscot, F.; Herrmann, S.; Menzel, R.; Ruess, L. Transcriptome analysis in oak uncovers a strong impact of endogenous rhythmic growth on the interaction with plant-parasitic nematodes. BMC Genom.; 2016; 17, 627. [DOI: https://dx.doi.org/10.1186/s12864-016-2992-8]
207. Graf, M.; Bönn, M.; Feldhahn, L.; Kurth, F.; Grams, T.E.E.; Herrmann, S.; Tarkka, M.; Buscot, F.; Scheu, S. Collembola interact with mycorrhizal fungi in modifying oak morphology, C and N incorporation and transcriptomics. R. Soc. Open Sci.; 2019; 6, 181869. [DOI: https://dx.doi.org/10.1098/rsos.181869]
208. Lesur, I.; Le Provost, G.; Bento, P.; Da Silva, C.; Leplé, J.-C.; Murat, F.; Ueno, S.; Bartholomé, J.; Lalanne, C.; Ehrenmann, F. et al. The oak gene expression atlas: Insights into Fagaceae genome evolution and the discovery of genes regulated during bud dormancy release. BMC Genom.; 2015; 16, 112. [DOI: https://dx.doi.org/10.1186/s12864-015-1331-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25765701]
209. Pereira-Leal, J.B.; Abreu, I.A.; Alabaça, C.S.; Almeida, M.H.; Almeida, P.; Almeida, T.; Amorim, M.I.; Araújo, S.; Azevedo, H.; Badia, A. et al. A comprehensive assessment of the transcriptome of cork oak (Quercus suber) through EST sequencing. BMC Genom.; 2014; 15, 371. [DOI: https://dx.doi.org/10.1186/1471-2164-15-371] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24885229]
210. Rocheta, M.; Sobral, R.; Magalhães, J.; Amorim, M.I.; Ribeiro, T.; Pinheiro, M.; Egas, C.; Morais-Cecãlio, L.; Costa, M.M.R. Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber. Front. Plant Sci.; 2014; 5, 599. [DOI: https://dx.doi.org/10.3389/fpls.2014.00599] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25414713]
211. Sebastiana, M.; Vieira, B.; Lino-Neto, T.; Monteiro, F.; Figueiredo, A.; Sousa, L.; Pais, M.S.; Tavares, R.; Paulo, O.S. Oak root response to ectomycorrhizal symbiosis establishment: RNA-Seq derived transcript identification and expression profiling. PLoS ONE; 2014; 9, e98376.
212. Fernández-Piñán, S.; Boher, P.; Soler, M.; Figueras, M.; Serra, O. Transcriptomic analysis of cork during seasonal growth highlights regulatory and developmental processes from phellogen to phellem formation. Sci. Rep.; 2021; 11, 12053. [DOI: https://dx.doi.org/10.1038/s41598-021-90938-5]
213. Meireles, B.; Usié, A.; Barbosa, P.; Fortes, A.M.; Folgado, A.; Chaves, I.; Carrasquinho, I.; Costa, R.L.; Gonçalves, S.; Teixeira, R.T. et al. Characterization of the cork formation and production transcriptome in Quercus cerris× suber hybrids. Physiol. Mol. Biol. Plants; 2018; 24, pp. 535-549. [DOI: https://dx.doi.org/10.1007/s12298-018-0526-3]
214. Guerrero-Sanchez, V.M.; Maldonado-Alconada, A.M.; Sánchez-Lucas, R.; Rey, M.D. Specific protein database creation from transcriptomics data in nonmodel species: Holm oak (Quercus ilex L.). Plant Proteom.; 2020; 4, pp. 57-68.
215. Ueno, S.; Klopp, C.; Leplé, J.C.; Derory, J.; Noirot, C.; Léger, V.; Prince, E.; Kremer, A.; Plomion, C.; Le Provost, G. Transcriptional profiling of bud dormancy induction and release in oak by next-generation sequencing. BMC Genom.; 2013; 14, 236. [DOI: https://dx.doi.org/10.1186/1471-2164-14-236]
216. Torre, S.; Tattini, M.; Brunetti, C.; Fineschi, S.; Fini, A.; Ferrini, F.; Sebastiani, F. RNA-seq analysis of Quercus pubescens leaves: De novo transcriptome assembly, annotation and functional markers development. PLoS ONE; 2014; 9, e112487. [DOI: https://dx.doi.org/10.1371/journal.pone.0112487]
217. Dong, X.; Huang, L.; Chen, Q.; Lv, Y.; Sun, H.; Liang, Z. Physiological and anatomical differences and differentially expressed genes reveal yellow leaf coloration in shumard oak. Plants; 2020; 9, 169. [DOI: https://dx.doi.org/10.3390/plants9020169]
218. Teixeira, R.T.; Fortes, A.M.; Bai, H.; Pinheiro, C.; Pereira, H. Transcriptional profiling of cork oak phellogenic cells isolated by laser microdissection. Planta; 2018; 247, pp. 317-338. [DOI: https://dx.doi.org/10.1007/s00425-017-2786-5]
219. Kurth, F.; Feldhahn, L.; Bönn, M.; Herrmann, S.; Buscot, F.; Tarkka, M.T. Large scale transcriptome analysis reveals interplay between development of forest trees and a beneficial mycorrhiza helper bacterium. BMC Genom.; 2015; 16, 658. [DOI: https://dx.doi.org/10.1186/s12864-015-1856-y]
220. Teixeira, R.T.; Fortes, A.M.; Pinheiro, C.; Pereira, H. Comparison of good-and bad-quality cork: Application of high-throughput sequencing of phellogenic tissue. J. Exp. Bot.; 2014; 65, pp. 4887-4905. [DOI: https://dx.doi.org/10.1093/jxb/eru252] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24958897]
221. Miguel, A.; de Vega-Bartol, J.; Marum, L.; Chaves, I.; Santo, T.; Leitão, J.; Varela, M.C.; Miguel, C.M. Characterization of the cork oak transcriptome dynamics during acorn development. BMC Plant Biol.; 2015; 15, 158. [DOI: https://dx.doi.org/10.1186/s12870-015-0534-1]
222. Boher, P.; Soler, M.; Sánchez, A.; Hoede, C.; Noirot, C.; Paiva, J.A.P.; Serra, O.; Figueras, M. A comparative transcriptomic approach to understanding the formation of cork. Plant Mol. Biol.; 2018; 96, pp. 103-118. [DOI: https://dx.doi.org/10.1007/s11103-017-0682-9]
223. Busta, L.; Serra, O.; Kim, O.T.; Molinas, M.; Peré-Fossoul, I.; Figueras, M.; Jetter, R. Oxidosqualene cyclases involved in the biosynthesis of triterpenoids in Quercus suber cork. Sci. Rep.; 2020; 10, 8011. [DOI: https://dx.doi.org/10.1038/s41598-020-64913-5]
224. Derory, J.; Léger, P.; Garcia, V.; Schaeffer, J.; Hauser, M.T.; Salin, F.; Luschnig, C.; Plomion, C.; Glössl, J.; Kremer, A. Transcriptome analysis of bud burst in sessile oak (Quercus petraea). New Phytol.; 2006; 170, pp. 723-738. [DOI: https://dx.doi.org/10.1111/j.1469-8137.2006.01721.x]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; María-Dolores Rey 1
; Labella-Ortega, Mónica 1
; Tienda-Parrilla, Marta 1 ; Hernández-Lao, Tamara 1 ; Honrubia-Gómez, Irene 1 ; Ramírez-García, Javier 1 ; Guerrero-Sanchez, Víctor M 2
; López-Hidalgo, Cristina 3
; Valledor, Luis 4
; Navarro-Cerrillo, Rafael M 5
; Jorrin-Novo, Jesús V 1
1 Agroforestry and Plant Biochemistry, Proteomics and Systems Biology, Department of Biochemistry and Molecular Biology, University of Cordoba, UCO-CeiA3, 14014 Cordoba, Spain
2 Agroforestry and Plant Biochemistry, Proteomics and Systems Biology, Department of Biochemistry and Molecular Biology, University of Cordoba, UCO-CeiA3, 14014 Cordoba, Spain; Cardiovascular Proteomics Laboratory, Centro Nacional de Investigaciones Cardiovasculares (CNIC), 28029 Madrid, Spain
3 Agroforestry and Plant Biochemistry, Proteomics and Systems Biology, Department of Biochemistry and Molecular Biology, University of Cordoba, UCO-CeiA3, 14014 Cordoba, Spain; Plant Physiology, Department of Organisms and Systems Biology, University Institute of Biotechnology of Asturias (IUBA), University of Oviedo, 33006 Asturias, Spain
4 Plant Physiology, Department of Organisms and Systems Biology, University Institute of Biotechnology of Asturias (IUBA), University of Oviedo, 33006 Asturias, Spain
5 Evaluation and Restoration of Agronomic and Forest Systems ERSAF, Department of Forest Engineering, University of Córdoba, 14014 Cordoba, Spain