In-depth analysis of the Quercus suber metabolome under drought stress and recovery reveals potential key metabolic players

Legacy effects of premature defoliation in response to an extreme drought event modulate phytochemical profiles with subtle consequences for leaf herbivory in European beech

Extreme droughts can have long-lasting effects on forest community dynamics and species interactions. Yet, our understanding of how drought legacy modulates ecological relationships is just unfolding. We tested the hypothesis that leaf chemistry and herbivory show long-term responses to premature defoliation caused by an extreme drought event in European beech (Fagus sylvatica L.). For two consecutive years after the extreme European summer drought in 2018, we collected leaves from the upper and lower canopy of adjacently growing drought-stressed and unstressed trees. Leaf chemistry was analyzed and leaf damage by different herbivore-feeding guilds was quantified. We found that drought had lasting impacts on leaf nutrients and on specialized metabolomic profiles. However, drought did not affect the primary metabolome. Drought-related phytochemical changes affected damage of leaf-chewing herbivores whereas damage caused by other herbivore-feeding guilds was largely unaffected. Drought legacy effects on phytochemistry and herbivory were often weaker than between-year or between-canopy strata variability. Our findings suggest that a single extreme drought event bears the potential to long-lastingly affect tree-herbivore interactions. Drought legacy effects likely become more important in modulating tree-herbivore interactions since drought frequency and severity are projected to globally increase in the coming decades.

MBROLE3: improved functional enrichment of chemical compounds for metabolomics data analysis

MBROLE (Metabolites Biological Role) facilitates the biological interpretation of metabolomics experiments. It performs enrichment analysis of a set of chemical compounds through statistical analysis of annotations from several databases. The original MBROLE server was released in 2011 and, since then, different groups worldwide have used it to analyze metabolomics experiments from a variety of organisms. Here we present the latest version of the system, MBROLE3, accessible at http://csbg.cnb.csic.es/mbrole3. This new version contains updated annotations from previously included databases as well as a wide variety of new functional annotations, such as additional pathway databases and Gene Ontology terms. Of special relevance is the inclusion of a new category of annotations, 'indirect annotations', extracted from the scientific literature and from curated chemical-protein associations. The latter allows to analyze enriched annotations of the proteins known to interact with the set of chemical compounds of interest. Results are provided in the form of interactive tables, formatted data to download, and graphical plots.

Laser Microdissection of Woody and Suberized Plant Tissues for RNA-Seq Analysis

An accurate profile of gene expression at a cellular level can contribute to a better understanding of biological processes and complexities involved in regulatory mechanism of woody plants. Laser microdissection is one technique that allows isolation of specific, target cells or tissue from a heterogeneous cell population. This technique entails microscopic visualization of the selected tissue and use a laser beam to separate the desired cells from surrounding tissue. Initial identification of these cells is made based on morphology and/or histological staining. Some works have been made in several tissues and plant models. However, there are few studies of laser microdissection application in woody species, particularly, lignified and suberized cells. Moreover, the presence of high level of suberin in cell walls can be a big challenge for the application of this approach. In our study it was developed a technique for tissue isolation, using laser microdissection of four different plant cell types (phellogen, lenticels, cortex and xylem) from woody tissues of cork oak (Quercus suber), followed by RNA extraction and RNA-Seq. We tested several methodologies regarding laser microdissection, cryostat equipments, fixation treatments, duration of single-cells collection and number of isolated cells by laser microdissection and RNA extraction procedures. A simple and efficient protocol for tissue isolation by laser microdissection and RNA purification was obtained, with a final method validation of RNA-Seq analysis. The optimized methodology combining RNA-Seq for expression analysis will contribute to elucidate the molecular pathways associated with different development processes of the xylem and phellem in oaks, including the lenticular channels formation.

Understanding the Seasonal Effect of Metabolite Production in Terminalia catappa L. Leaves through a Concatenated MS- and NMR-Based Metabolomics Approach

Terminalia catappa L. (Combretaceae) is a medicinal plant that is part of the Brazilian biodiversity; this plant is popularly used for the treatment of a wide range of diseases. To better understand the chemical composition of T. catappa in different seasons, we conducted a thorough study using LC-MS and NMR data analysis techniques. The study helped obtain a chemical profile of the plant ethanolic extracts in different seasons of the year (spring, summer, autumn, and winter). The dereplication of LC-HRMS data allowed the annotation of 90 compounds in the extracts of T. catappa (hydrolyzable tannins, ellagic acid derivatives, and glycosylated flavonoids). Triterpenes and C-glycosyl flavones were the compounds that significantly contributed to differences observed between T. catappa plant samples harvested in autumn/winter and spring, respectively. The variations observed in the compound composition of the plant leaves may be related to processes induced by environmental stress and leaf development. Data fusion applied in the metabolomic profiling study allowed us to identify metabolites with greater confidence, and provided a better understanding regarding the production of specialized metabolites in T. catappa leaves under different environmental conditions, which may be useful to establish appropriate quality criteria for the standardization of this medicinal plant.

Multiomic Data Integration in the Analysis of Drought-Responsive Mechanisms in Quercus ilex Seedlings

The integrated analysis of different omic layers can provide new knowledge not provided by their individual analysis. This approach is also necessary to validate data and reveal post-transcriptional and post-translational mechanisms of gene expression regulation. In this work, we validated the possibility of applying this approach to non-model species such as Quercus ilex. Transcriptomics, proteomics, and metabolomics from Q. ilex seedlings subjected to drought-like conditions under the typical summer conditions in southern Spain were integrated using a non-targeted approach. Two integrative approaches, PCA and DIABLO, were used and compared. Both approaches seek to reduce dimensionality, preserving the maximum information. DIABLO also allows one to infer interconnections between the different omic layers. For easy visualization and analysis, these interconnections were analyzed using functional and statistical networks. We were able to validate results obtained by analyzing the omic layers separately. We identified the importance of protein homeostasis with numerous protease and chaperones in the networks. We also discovered new key processes, such as transcriptional control, and identified the key function of transcription factors, such as DREB2A, WRKY65, and CONSTANS, in the early response to drought.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, "dehesa." It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors' group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

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.

Insights on Fructans and Resistance of Plants to Drought Stress

Drought, one of the major abiotic stresses affecting plants, is characterized by a decrease of water availability, resulting in a decrease of the water potential (Ψ) of the cells. One of the strategies of plants in resisting to this low Ψ and related stresses is regulating their water-plant relation and the interplay between Ψsolutes and the turgor pressure (Ψp). This regulation avoids the dehydration induced by low Ψ and is resulting from the accumulation of specific molecules which induce higher tolerance to water deficit and also other mechanisms that prevent or repair cell damages. In plants, fructans, the non-structural carbohydrates (NSC), have other physiological functions than carbon reserve. Among these roles, fructans have been implicated in protecting plants against water deficit caused by drought. As an efficient strategy to survive to this abiotic stress, plants synthesize fructans in response to osmotic pressure in order to osmoregulate the cellular flux, therefore, protecting the membrane damage and maintaining Ψp. Although different studies have been conducted to elucidate the mechanisms behind this strategy, still the concept itself is not well-understood and many points remain unclear and need to be elucidated in order to understand the causal relation between water deficit and fructans accumulation during water scarcity. This understanding will be a key tool in developing strategies to enhance crop tolerance to stressful dry conditions, particularly under the changing climate prediction. This review aims to give new insights on the roles of fructans in the response and resistance of plants to water deficit and their fate under this severe environmental condition.

Metabolic differences of two constructive species in saline-alkali grassland in China

BACKGROUND

Salinization of soil is an urgent problem that restricts agroforestry production and environmental protection. Substantial accumulation of metal ions or highly alkaline soil alters plant metabolites and may even cause plant death. To explore the differences in the response strategies between Suaeda salsa (S. salsa) and Puccinellia tenuiflora (P. tenuiflora), two main constructive species that survive in saline-alkali soil, their metabolic differences were characterized.

RESULT

Metabolomics was conducted to study the role of metabolic differences between S. salsa and P. tenuiflora under saline-alkali stress. A total of 68 significantly different metabolites were identified by GC-MS, including 9 sugars, 13 amino acids, 8 alcohols, and 34 acids. A more detailed analysis indicated that P. tenuiflora utilizes sugars more effectively and may be saline-alkali tolerant via sugar consumption, while S. salsa utilizes mainly amino acids, alcohols, and acids to resist saline-alkali stress. Measurement of phenolic compounds showed that more C6C3C6-compounds accumulated in P. tenuiflora, while more C6C1-compounds, phenolic compounds that can be used as signalling molecules to defend against stress, accumulated in S. salsa.

CONCLUSIONS

Our observations suggest that S. salsa resists the toxicity of saline-alkali stress using aboveground organs and that P. tenuiflora eliminates this toxicity via roots. S. salsa has a stronger habitat transformation ability and can provide better habitat for other plants.

UHPLC-(ESI)-HRMS and NMR-Based Metabolomics Approach to Access the Seasonality of Byrsonima intermedia and Serjania marginata From Brazilian Cerrado Flora Diversity

Seasonality is one of the major environmental factors that exert influence over the synthesis and accumulation of secondary metabolites in medicinal plants. The application of the metabolomics approach for quality control of plant extracts is essentially important because it helps one to establish a standard metabolite profile and to analyze factors that affect the effectiveness of the medicinal plants. The Brazilian Cerrado flora is characterized by a rich diversity of native plant species, and a number of these plant species have been found to have suitable medicinal properties. Some of these plant species include Byrsonima intermedia and Serjania marginata. To better understand the chemical composition of these plant species, we conducted a study using the state-of-the-art techniques including the HPLC system coupled to an Exactive-Orbitrap high resolution mass spectrometer with electrospray ionization interface UHPLC-(ESI)-HRMS and by NMR being performed 2D J-resolved and proton NMR spectroscopy. For the analysis, samples were harvested bimonthly during two consecutive years. UHPLC-(ESI)-HRMS data were preprocessed and the output data uploaded into an in-house Excel macro for peak dereplication. MS and NMR data were concatenated using the data fusion method and submitted to multivariate statistical analysis. The dereplication of LC-HRMS data helped in the annotation of the major compounds present in the extracts of the three plant species investigated allowing the annotation of 68 compounds in the extracts of B. intermedia (cinnamic acids, phenolic acids derived from galloyl quinic and shikimic acid, proanthocyanidins, glycosylated flavonoids, triterpenes and other phenols) and 81 compounds in the extracts of S. marginata (phenolic acids, saponins, proanthocyanidins, glycosylated flavonoids among other compounds). For a better assessment of the great number of responses, the significance of the chemical variables for the differentiation and correlation of the seasons was determined using the variable importance on projection (VIP) technique and through the application of the false discovery rate (FDR) estimation. The statistical data obtained showed that seasonal factors played an important role on the production of metabolites in each plant species. Temperature conditions, drought and solar radiation were found to be the main factors that affected the variability of phenolic compounds in each species.

Metabolomics and transcriptomics to decipher molecular mechanisms underlying ectomycorrhizal root colonization of an oak tree

Mycorrhizas are known to have a positive impact on plant growth and ability to resist major biotic and abiotic stresses. However, the metabolic alterations underlying mycorrhizal symbiosis are still understudied. By using metabolomics and transcriptomics approaches, cork oak roots colonized by the ectomycorrhizal fungus Pisolithus tinctorius were compared with non-colonized roots. Results show that compounds putatively corresponding to carbohydrates, organic acids, tannins, long-chain fatty acids and monoacylglycerols, were depleted in ectomycorrhizal cork oak colonized roots. Conversely, non-proteogenic amino acids, such as gamma-aminobutyric acid (GABA), and several putative defense-related compounds, including oxylipin-family compounds, terpenoids and B6 vitamers were induced in mycorrhizal roots. Transcriptomic analysis suggests the involvement of GABA in ectomycorrhizal symbiosis through increased synthesis and inhibition of degradation in mycorrhizal roots. Results from this global metabolomics analysis suggest decreases in root metabolites which are common components of exudates, and in compounds related to root external protective layers which could facilitate plant-fungal contact and enhance symbiosis. Root metabolic pathways involved in defense against stress were induced in ectomycorrhizal roots that could be involved in a plant mechanism to avoid uncontrolled growth of the fungal symbiont in the root apoplast. Several of the identified symbiosis-specific metabolites, such as GABA, may help to understand how ectomycorrhizal fungi such as P. tinctorius benefit their host plants.

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.

Metabolomics Reveal Induction of ROS Production and Glycosylation Events in Wheat Upon Exposure to the Green Leaf Volatile Z-3-Hexenyl Acetate

The activation and priming of plant defense upon perception of green leaf volatiles (GLVs) have often been reported. However, information as to which metabolic pathways in plants are affected by GLVs remains elusive. We report the production of reactive oxygen species in the tip of young wheat leaves followed by activation of antioxidant-related enzyme activity. In this study, we aimed to uncover metabolic signatures upon exposure to the GLV Z-3-hexenyl acetate (Z-3-HAC). By using an untargeted metabolomics approach, we observed changes in the phenylpropanoid pathways which yield metabolites that are involved in many anti-oxidative processes. Furthermore, exposure to GLV, followed by infection with Fusarium graminearum (Fg), induced significantly greater changes in the phenylpropanoid pathway compared to a sole Z-3-HAC treatment. Fragmentation of a selection of metabolites, which are significantly more upregulated in the Z-3-HAC + Fg treatment, showed D-glucose to be present as a substructure. This suggests that Z-3-HAC induces early glycosylation processes in plants. Additionally, we identified the presence of hexenyl diglycosides, which indicates that aerial Z-3-HAC is metabolized in the leaves by glycosyltransferases. Together these data indicate that GLV Z-3-HAC is taken up by leaves and incites oxidative stress. This subsequently results in the modulation of the phenylpropanoid pathway and an induction of glycosylation processes.