Population Genetic Diversity of Quercus ilex subsp. ballota (Desf.) Samp. Reveals Divergence in Recent and Evolutionary Migration Rates in the Spanish Dehesas

A novel framework of smart monitoring to face the challenges of tree management in historic gardens

Historic gardens are green spaces characterised by tree stands with several veteran specimens of high artistic and cultural value. Such valuable plant components have to cope with biotic and abiotic stress factors as well as ongoing senescence processes. Maintaining tree health is therefore crucial to preserve their ecosystem services, but also to protect the monument and visitor health. In this context, finding smart, fast and cost-effective management solutions to monitor health and detect critical conditions for both stands and individual veteran trees can promote garden conservation. For this reason, we developed a novel framework based on Sentinel2 imagery, LiDAR sources and automatic cameras to identify risk spots regarding trees in historic gardens. The pilot study area consists of two closed Italian gardens from the 16th century, which were analysed as a unique Historic Garden System (HGS). The tree health status at stand level was assessed using a criterion based on the Normalized Difference Vegetation Index weighed on tree volume (NDVIt) and validated by a visual crown defoliation assessment. At the tree level, the health status of four veteran trees defined by the NDVIt was also evaluated using green chromatic coordinates (GCC) obtained from digital images acquired by cameras at daily intervals during one growing season. The 33% of the tree population was classified as being in poor health, i.e. "at risk". Veteran trees classified as "at risk" showed an anticipation of phenological phases and a lower GCC compared to reference trees. Despite variability determined by Sentinel medium resolution, the proposed framework showed good accuracy (0.74) for monitoring historical gardens. The semi-automatic risk point mapping system tested here proved to be effective in facilitating the management of historic gardens, which in turn could be applied in the wider context of urban greening.

Proteomic and Metabolomic Analysis of the Quercus ilex-Phytophthora cinnamomi Pathosystem Reveals a Population-Specific Response, Independent of Co-Occurrence of Drought

Holm oak (Quercus ilex) is considered to be one of the major structural elements of Mediterranean forests and the agrosilvopastoral Spanish "dehesa", making it an outstanding example of ecological and socioeconomic sustainability in forest ecosystems. The exotic Phytophthora cinnamomi is one of the most aggressive pathogens of woody species and, together with drought, is considered to be one of the main drivers of holm oak decline. The effect of and response to P. cinnamomi inoculation were studied in the offspring of mother trees from two Andalusian populations, Cordoba and Huelva. At the two locations, acorns collected from both symptomatic (damaged) and asymptomatic (apparently healthy) trees were sampled. Damage symptoms, mortality, and chlorophyll fluorescence were evaluated in seedlings inoculated under humid and drought conditions. The effect and response depended on the population and were more apparent in Huelva than in Cordoba. An integrated proteomic and metabolomic analysis revealed the involvement of different metabolic pathways in response to the pathogen in both populations, including amino acid metabolism pathways in Huelva, and terpenoid and flavonoid biosynthesis in Cordoba. However, no differential response was observed between seedlings inoculated under humid and drought conditions. A protective mechanism of the photosynthetic apparatus was activated in response to defective photosynthetic activity in inoculated plants, which seemed to be more efficient in the Cordoba population. In addition, enzymes and metabolites of the phenylpropanoid and flavonoid biosynthesis pathways may have conferred higher resistance in the Cordoba population. Some enzymes are proposed as markers of resilience, among which glyoxalase I, glutathione reductase, thioredoxin reductase, and cinnamyl alcohol dehydrogenase are candidates.

Quercus ilex L. dieback is genetically determined: Evidence provided by dendrochronology, δ13C and SSR genotyping

Quercus ilex L. dieback has been reported in several Mediterranean forests, revealing different degree of crown damages even in close sites, as observed in two Q. ilex forest stands in southern Tuscany (IT). In this work, we applied a novel approach combining dendrochronological, tree-ring δ13C and genetic analysis to test the hypothesis that different damage levels observed in a declining (D) and non-declining (ND) Q. ilex stands are connected to population features linked to distinct response to drought. Furthermore, we investigated the impact of two major drought events (2012 and 2017), that occurred in the last fifteen years in central Italy, on Q. ilex growth and intrinsic water use efficiency (WUEi). Overall, Q. ilex showed slightly different ring-width patterns between the two stands, suggesting a lower responsiveness to seasonal climatic variations for trees at D stand, while Q. ilex at ND stand showed changes in the relationship between climatic parameters and growth across time. The strong divergence in δ13C signals between the two stands suggested a more conservative use of water for Q. ilex at ND compared to D stand that may be genetically driven. Q. ilex at ND resulted more resilient to drought compared to trees at D, probably thanks to its safer water strategy. Genotyping analysis based on simple-sequence repeat (SSR) markers revealed the presence of different Q. ilex populations at D and ND stands. Our study shows intraspecific variations in drought response among trees grown in close. In addition, it highlights the potential of combining tree-ring δ13C data with SSR genotyping for the selection of seed-bearing genotypes aimed to preserve Mediterranean holm oak ecosystem and improve its forest management.

A first draft genome of holm oak (Quercus ilex subsp. ballota), the most representative species of the Mediterranean forest and the Spanish agrosylvopastoral ecosystem "dehesa"

The holm oak (Quercus ilex subsp. ballota) is the most representative species of the Mediterranean Basin and the agrosylvopastoral Spanish "dehesa" ecosystem. Being part of our life, culture, and subsistence since ancient times, it has significant environmental and economic importance. More recently, there has been a renewed interest in using the Q. ilex acorn as a functional food due to its nutritional and nutraceutical properties. However, the holm oak and its related ecosystems are threatened by different factors, with oak decline syndrome and climate change being the most worrying in the short and medium term. Breeding programs informed by the selection of elite genotypes seem to be the most plausible biotechnological solution to rescue populations under threat. To achieve this and other downstream analyses, we need a high-quality and well-annotated Q. ilex reference genome. Here, we introduce the first draft genome assembly of Q. ilex using long-read sequencing (PacBio). The assembled nuclear haploid genome had 530 contigs totaling 842.2 Mbp (N50 = 3.3 Mbp), of which 448.7 Mb (53%) were repetitive sequences. We annotated 39,443 protein-coding genes of which 94.80% were complete and single-copy genes. Phylogenetic analyses showed no evidence of a recent whole-genome duplication, and high synteny of the 12 chromosomes between Q. ilex and Quercus lobata and between Q. ilex and Quercus robur. The chloroplast genome size was 142.3 Kbp with 149 protein-coding genes successfully annotated. This first draft should allow for the validation of omics data as well as the identification and functional annotation of genes related to phenotypes of interest such as those associated with resilience against oak decline syndrome and climate change and higher acorn productivity and nutraceutical value.

Population structure and genetic diversity of Tamarix chinensis as revealed with microsatellite markers in two estuarine flats

Background

Tamarix chinensis Lour. is a 3-6-meter-tall small tree with high salt- and alkali- tolerance and aggressive invasiveness, mainly distributed in the eastern part of China in warm-temperate and subtropical climate zones, yet there is little information available regarding genetic diversity and population structure.

Methods

A total of 204 individuals of nine T. chinensis populations were investigated for genetic diversity and population structure using a set of 12 highly polymorphic microsatellite markers.

Results

The total number of alleles detected was 162, the average number of effective allele was 4.607, the average polymorphism information content (PIC) value of the 12 loci was 0.685, and the mean observed heterozygosity (Ho) and the mean expected heterozygosity (He) was 0.653 and 0.711, respectively. Analysis of molecular variance (AMOVA) showed a 5.32% genetic variation among T. chinensis populations. Despite a low population differentiation, Bayesian clustering analysis, discriminant analysis of principal components (DAPC) and the unweighted pair group method with arithmetic mean (UPGMA) clearly identified three genetic clusters correlated to the populations' geographic origin: the northern populations including those from Yellow River Delta, the Fangshan (FS) population from Beijing, the Changyi (CY) population from Bohai Bay, the Huanjiabu (HHJ) population from Hangzhou Bay, and the remaining two populations from Hangzhou Bay. There was a significant relationship between the genetic distance and geographical distance of the paired populations. Gene flow (Nm) was 4.254 estimated from FST.

Conclusion

T. chinensis possessed high genetic diversity comparable to tree species, and although the population differentiation is shallow, our results classified the sampled populations according to sampling localities, suggesting the different origins of the study populations.

Predicting potential reforestation areas by Quercus ilex (L.) species using machine learning algorithms: case of upper Ziz, southeastern Morocco

The selection of appropriate areas for reforestation remains a complex task because of influence by several factors, which requires the use of new techniques. Based on the accurate outcomes obtained through machine learning in prior investigations, the current study evaluates the capacities of six machine learning techniques (MLT) for delineating optimal areas for reforestation purposes specifically targeting Quercus ilex, an important local species to protect soil and water in upper Ziz, southeast Morocco. In the initial phase, the remaining stands of Q. ilex were identified, and at each site, measurements were taken for a set of 12 geo-environmental parameters including slope, aspect, elevation, geology, distance to stream, rainfall, slope length, plan curvature, profile curvature, erodibility, soil erosion, and land use/land cover. Subsequently, six machine learning algorithms were applied to model optimal areas for reforestation. In terms of models' performance, the results were compared, and the best were obtained by Bagging (area under the curve (AUC) = 0.98) and Naive Bayes (AUC = 0.97). Extremely favorable areas represent 8% and 17% of the study area according to Bagging and NB respectively, located to the west where geological unit of Bathonian-Bajocian with low erodibility index (K) and where rainfall varies between 250 and 300 mm/year. This work provides a roadmap for decision-makers to increase the chances of successful reforestation at lower cost and in less time.

Dramatic impact of future climate change on the genetic diversity and distribution of ecologically relevant Western Mediterranean Carex (Cyperaceae)

Anticipating the evolutionary responses of species to ongoing climate change is essential to propose effective management and conservation measures. The Western Mediterranean Basin constitutes one of the hotspots of biodiversity where the effects of climate change are expected to be more dramatic. Plant species with ecological relevance constitute ideal models to evaluate and predict the impact of climate change on ecosystems. Here we investigate these impacts through the spatio-temporal comparison of genetic diversity/structure (AFLPs), potential distribution under different future scenarios of climate change, and ecological space in two Western Mediterranean sister species of genus Carex. Both species are ecologically key in their riparian habitats, but display contrasting distribution patterns, with one widespread in the Iberian Peninsula and North Africa (C. reuteriana), while the other (C. panormitana) is a restricted, probably endangered, Central Mediterranean endemic. At present, we found a strong genetic structure driven by geography in both species, and lower values of genetic diversity and a narrower ecological space in C. panormitana than in C. reuteriana, while the allelic rarity was higher in the former than in C. reuteriana subspecies. Future projections predict an overall dramatic reduction of suitable areas for both species under all climate change scenarios, which could be almost total for C. panormitana. In addition, gene diversity was inferred to decrease in all taxa, with genetic structure reinforcing in C. reuteriana by the loss of admixture among populations. Our findings stress the need for a reassessment of C. panormitana conservation status under IUCN Red List criteria and the implementation of conservation measures.

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

The effect and responses to drought stress were analyzed in Quercus ilex L. seedlings using a nontargeted metabolomic approach, implementing the approaches of previous studies in which other -omics platforms, transcriptomics, and proteomics were employed. This work aimed to characterize the Q. ilex leaf metabolome, determining possible mechanisms and molecular markers of drought tolerance and identifying putative bioactive compounds. Six-month-old seedling leaves subjected to drought stress imposed by water withholding under high-temperature and irradiance conditions were collected when leaf fluorescence decreased by 20% (day 17) and 45% (day 24) relative to irrigated seedlings. A total of 3934 compounds were resolved, with 616 being variable and 342 identified, which belonged to five chemical families. Out of the identified compounds, 33 were variable, mostly corresponding to amino acids, carboxylic acids, benzenoids, flavonoids and isoprenoids. Epigallocatechin, ellagic acid, pulegone, indole-3-acrylic acid and dihydrozeatin-O-glucoside were up-accumulated under drought conditions at both sampling times. An integrated multi-omics analysis of phenolic compounds and related enzymes was performed, revealing that some enzymes involved in the flavonoid pathways (chalcone synthase, anthocyanidin synthase and anthocyanidin reductase) were up-accumulated at day 24 in non-irrigated seedlings. Some putative markers of tolerance to drought in Q. ilex are proposed for assisting breeding programs based on the selection of elite genotypes.

Antagonistic Potential of Native Trichoderma spp. against Phytophthora cinnamomi in the Control of Holm Oak Decline in Dehesas Ecosystems

Phytophthora root rot caused by the pathogen Phytophthora cinnamomi is one of the main causes of oak mortality in Mediterranean open woodlands, the so-called dehesas. Disease control is challenging; therefore, new alternative measures are needed. This study focused on searching for natural biocontrol agents with the aim of developing integrated pest management (IPM) strategies in dehesas as a part of adaptive forest management (AFM) strategies. Native Trichoderma spp. were selectively isolated from healthy trees growing in damaged areas by P. cinnamomi root rot, using Rose Bengal selective medium. All Trichoderma (n = 95) isolates were evaluated against P. cinnamomi by mycelial growth inhibition (MGI). Forty-three isolates presented an MGI higher than 60%. Twenty-one isolates belonging to the highest categories of MGI were molecularly identified as T. gamsii, T. viridarium, T. hamatum, T. olivascens, T. virens, T. paraviridescens, T. linzhiense, T. hirsutum, T. samuelsii, and T. harzianum. Amongst the identified strains, 10 outstanding Trichoderma isolates were tested for mycoparasitism, showing values on a scale ranging from 3 to 4. As far as we know, this is the first report referring to the antagonistic activity of native Trichoderma spp. over P. cinnamomi strains cohabiting in the same infected dehesas. The analysis of the tree health status and MGI suggest that the presence of Trichoderma spp. might diminish or even avoid the development of P. cinnamomi, protecting trees from the worst effects of P. cinnamomi root rot.

Unraveling the genetic diversity and structure of Quercus liaotungensis population through analysis of microsatellite markers

Background

Quercus liaotungensis Koidz. is an ecologically and economically important tree species widely distributed in Northern China. However, the effective assessment, utilization, and protection of Q. liaotungensis resources remain unexplored.

Methods

In total, 120 samples obtained from 12 Q. liaotungensis populations of Northern China were investigated for genetic diversity and structure using 19 simple sequence repeat (SSR) primer pairs.

Results

The total number of alleles detected was 293, the average number of effective allele (Ne) was 6.084, the genetic differentiation coefficient (Fst) was 0.033, and the mean observed heterozygosity (Ho) and expected heterozygosity (He) were 0.690 and 0.801, respectively. Moreover, analysis of molecular variance (AMOVA) showed a 5.5% genetic variation among 12 Q. liaotungensis populations, indicating that a high level of genetic diversity and a low degree of genetic differentiation among Q. liaotungensis populations. STRUCTURE and cluster analysis divided the 12 Q. liaotungensis populations into the following three subpopulations: Bashang Plateau subpopulation (SH), Liaodong Peninsula subpopulation (NC), and Loess Plateau subpopulation (other 10 populations). The cluster analysis based on 19 climatic factors was consistent with the genetic structure. A positive correlation was found between genetic distance and geographical distance (r = 0.638, p = 0.028) by the Mantel test, and two boundaries were found among the 12 Q. liaotungensis populations by the Barrier analysis, indicating that Q. liaotungensis populations existed isolated by geographical distance and physical barrier.

Conclusion

This study suggests that geographical isolation, physical barrier, climatic types, and natural hybridization promote the formation of genetic structures, which can contribute to future protection and genetic improvement of Q. liaotungensis.

Proteomics Data Analysis for the Identification of Proteins and Derived Proteotypic Peptides of Potential Use as Putative Drought Tolerance Markers for Quercus ilex

Drought is one of the main causes of mortality in holm oak (Quercus ilex) seedlings used in reforestation programs. Although this species shows high adaptability to the extreme climate conditions prevailing in Southern Spain, its intrinsic genetic variability may play a role in the differential response of some populations and individuals. The aim of this work was to identify proteins and derived proteotypic peptides potentially useful as putative markers for drought tolerance in holm oak by using a targeted post-acquisition proteomics approach. For this purpose, we used a set of proteins identified by shotgun (LC-MSMS) analysis in a drought experiment on Q. ilex seedlings from four different provenances (viz. the Andalusian provinces Granada, Huelva, Cadiz and Seville). A double strategy involving the quantification of proteins and target peptides by shotgun analysis and post-acquisition data analysis based on proteotypic peptides was used. To this end, an initial list of proteotypic peptides from proteins highly represented under drought conditions was compiled that was used in combination with the raw files from the shotgun experiment to quantify the relative abundance of the fragment’s ion peaks with the software Skyline. The most abundant peptides under drought conditions in at least two populations were selected as putative markers of drought tolerance. A total of 30 proteins and 46 derived peptides belonging to the redox, stress-related, synthesis,-folding and degradation, and primary and secondary metabolism functional groups were thus identified. Two proteins (viz., subtilisin and chaperone GrpE protein) were found at increased levels in three populations, which make them especially interesting for validation drought tolerance markers in subsequent experiments.

Molecular Research on Stress Responses in Quercus spp.: From Classical Biochemistry to Systems Biology through Omics Analysis

The genus Quercus (oak), family Fagaceae, comprises around 500 species, being one of the most important and dominant woody angiosperms in the Northern Hemisphere. Nowadays, it is threatened by environmental cues, which are either of biotic or abiotic origin. This causes tree decline, dieback, and deforestation, which can worsen in a climate change scenario. In the 21st century, biotechnology should take a pivotal role in facing this problem and proposing sustainable management and conservation strategies for forests. As a non-domesticated, long-lived species, the only plausible approach for tree breeding is exploiting the natural diversity present in this species and the selection of elite, more resilient genotypes, based on molecular markers. In this direction, it is important to investigate the molecular mechanisms of the tolerance or resistance to stresses, and the identification of genes, gene products, and metabolites related to this phenotype. This research is being performed by using classical biochemistry or the most recent omics (genomics, epigenomics, transcriptomics, proteomics, and metabolomics) approaches, which should be integrated with other physiological and morphological techniques in the Systems Biology direction. This review is focused on the current state-of-the-art of such approaches for describing and integrating the latest knowledge on biotic and abiotic stress responses in Quercus spp., with special reference to Quercus ilex, the system on which the authors have been working for the last 15 years. While biotic stress factors mainly include fungi and insects such as Phytophthora cinnamomi, Cerambyx welensii, and Operophtera brumata, abiotic stress factors include salinity, drought, waterlogging, soil pollutants, cold, heat, carbon dioxide, ozone, and ultraviolet radiation. The review is structured following the Central Dogma of Molecular Biology and the omic cascade, from DNA (genomics, epigenomics, and DNA-based markers) to metabolites (metabolomics), through mRNA (transcriptomics) and proteins (proteomics). An integrated view of the different approaches, challenges, and future directions is critically discussed.

Phytochemical composition and variability in Quercus ilex acorn morphotypes as determined by NIRS and MS-based approaches

Quercus acorns have been used as a staple food since archeological times; currently, there is a renewed interest in the direct use of the acorns as nuts or as a source of consumable flour. In this regard, the phytochemical composition and variability of various morphotypes of acorns from Quercus ilex have been determined by using two platforms, near infrared spectroscopy (NIRS) and mass spectrometry (MS)-based metabolomics. According to NIRS, starch was the most abundant compound (50-60%). UHPLC-QToF analysis was used for untargeted metabolomic analysis and 192 metabolites were annotated. Principal component analysis (PCA) discriminated the morphotypes; 50 compounds out of 192 identified had the highest load over the first two PCA components (explaining 67.2% variability) and can be potential markers of variability. The data presented here support the value of acorns as a source of nutritionally rich compounds thus justifying their use a good alternative to other nuts.

Responses and Differences in Tolerance to Water Shortage under Climatic Dryness Conditions in Seedlings from Quercus spp. and Andalusian Q. ilex Populations

Analyzing differences in tolerance to drought in Quercus spp., and the characterization of these responses at the species and individual population level, are imperative for the selection of resilient elite genotypes in reforestation programs. The main objective of this work was to evaluate differences in the response and tolerance to water shortage under in five Quercus spp. and five Andalusian Q. ilex populations at the inter- and intraspecies level. Six-month-old seedlings grown in perlite were subjected to drought treatments by withholding water for 28 days under mean 37 °C temperature, 28 W m−2 solar irradiance, and 41% humidity. The use of perlite as the substrate enabled the establishment of severe drought stress with reduction in water availability from 73% (field capacity) to 28% (dryness), corresponding to matric potentials of 0 and −30 kPa. Damage symptoms, mortality rate, leaf water content, photosynthetic, and biochemical parameters (amino acids, sugars, phenolics, and pigments) were determined. At the phenotypic level, based on damage symptoms and mortality, Q. ilex behaved as the most drought tolerant species. Drought caused a significant decrease in leaf fluorescence, photosynthesis rate, and stomatal conductance in all Quercus spp. analyzed, being less pronounced in Q. ilex. There were not differences between irrigated and non-irrigated Q. ilex seedlings in the content of sugar and photosynthetic pigments, while the total amino acid and phenolic content significantly increased under drought conditions. As a response to drought, living Q. ilex seedlings adjust stomata opening and gas exchange, and keep hydrated, photosynthetically active, and metabolically competent. At the population level, based on damage symptoms, mortality, and physiological parameters, the eastern Andalusian populations were more tolerant than the western ones. These observations inform the basis for the selection of resilient genotypes to be used in breeding and reforestation programs.

The Dehesas of Extremadura, Spain: A Potential for Socio-economic Development Based on Agritourism Activities

Dehesas form an agroforestry system which is highly developed in Extremadura as they occupy over 1,000,000 hectares. This is an extensive production system based on exploitation by stockbreeding with certain complements of agricultural and extractive products. It is however underexploited for activities other than agrolivestock production. As a result this research starts from the hypothesis that greater socio-economic development is possible in areas near this type of ecosystem if we consider that it is extremely attractive to tourism, at least in some of its forms. Complementing agrolivestock exploitation with tourist exploitation means that it is necessary to be familiar with the situation of the current availability of accommodation. In order to do so we analyse the main parameters characterising rural accommodation in Extremadura, Spain, the results of which reflect low occupancy for a large part of the year and consequently seasonal variation is high. This fact contrasts with the huge potential of dehesas for agritourism. In order to carry out this study we resort to the use of spatial statistics, in particular the grouping analysis. In its configuration we consider the location of the rural accommodation and its proximity to areas of dehesas made up of holm oaks (Quercus ilex) and cork oaks (Quecus suber) together with protected natural spaces, basing ourselves also on proximity to livestock trails and natural swimming pools. The results lead to the creation of 5 homogeneous groups of which 3 correspond to accommodation establishments capable of setting up agritourism and agriecotourism initiatives; this affects 45% of the establishments. However, no current initiatives of this type exist despite the fact that the National Reference Centre for Agritourism is located in Extremadura. From this can be inferred the need for implementing tourist policies to encourage the generating of tourist products specifically intended to exploit the potential of the dehesa as a complement to the current availability of rural accommodation.

Proteomics, Holm Oak (Quercus ilex L.) and Other Recalcitrant and Orphan Forest Tree Species: How do They See Each Other?

Proteomics has had a big impact on plant biology, considered as a valuable tool for several forest species, such as Quercus, Pines, Poplars, and Eucalyptus. This review assesses the potential and limitations of the proteomics approaches and is focused on Quercus ilex as a model species and other forest tree species. Proteomics has been used with Q. ilex since 2003 with the main aim of examining natural variability, developmental processes, and responses to biotic and abiotic stresses as in other species of the genus Quercus or Pinus. As with the progress in techniques in proteomics in other plant species, the research in Q. ilex moved from 2-DE based strategy to the latest gel-free shotgun workflows. Experimental design, protein extraction, mass spectrometric analysis, confidence levels of qualitative and quantitative proteomics data, and their interpretation are a true challenge with relation to forest tree species due to their extreme orphan and recalcitrant (non-orthodox) nature. Implementing a systems biology approach, it is time to validate proteomics data using complementary techniques and integrate it with the -omics and classical approaches. The full potential of the protein field in plant research is quite far from being entirely exploited. However, despite the methodological limitations present in proteomics, there is no doubt that this discipline has contributed to deeper knowledge of plant biology and, currently, is increasingly employed for translational purposes.