ABSTRACT
Recent studies have documented plant responses to climate change extensively, particularly to single-stress exposures. However, critical factors for stress survival, such as sexual differentiation, are not often considered. The dioicous Marchantia polymorpha stands as an evolutionary milestone, potentially preserving ancestral traits from the early colonizers. In this study, we employed proteomic analyses complemented with physiological monitoring to investigate combined heat and drought responses in Tak-1 (male) and Tak-2 (female) accessions of this liverwort. Additionally, targeted transcriptomics was conducted using different natural populations from contrasting environments. Our findings revealed sex-biased dynamics among natural accessions, particularly evident under control conditions and during early stress responses. Although Tak-2 exhibited greater diversity than Tak-1 under control conditions, male accession demonstrated distinct and more rapid stress sensing and signaling. These differences in stress response appeared to be strongly related to sex-specific plasticity influenced by geoclimatic origin. Furthermore, we established distinct protein gene ages and genomic distribution trends, underscoring the importance of protein diversification over time. This study provides an evolutionary perspective on sexual divergence and stress emergence employing a systems biology approach, which allowed for the establishment of global and sex-specific interaction networks in the stress response.
ABSTRACT
Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.
Subject(s)
Arabidopsis , Oryza , Genomics , Arabidopsis/genetics , Oryza/genetics , Zea mays/genetics , Plants/genetics , Chromatin , Evolution, Molecular , Genome, PlantABSTRACT
Global warming is causing rapid changes in mean annual temperature and more severe drought periods. These are major contributors of forest dieback, which is becoming more frequent and widespread. In this work, we investigated how the transcriptome of Pinus radiata changed during initial heat stress response and acclimation. To this end, we generated a high-density dataset employing Illumina technology. This approach allowed us to reconstruct a needle transcriptome, defining 12 164 and 13 590 transcripts as down- and up-regulated, respectively, during a time course stress acclimation experiment. Additionally, the combination of transcriptome data with other available omics layers allowed us to determine the complex inter-related processes involved in the heat stress response from the molecular to the physiological level. Nucleolus and nucleoid activities seem to be a central core in the acclimating process, producing specific RNA isoforms and other essential elements for anterograde-retrograde stress signaling such as NAC proteins (Pra_vml_051671_1 and Pra_vml_055001_5) or helicase RVB. These mechanisms are connected by elements already known in heat stress response (redox, heat-shock proteins, or abscisic acid-related) and with others whose involvement is not so well defined such as shikimate-related, brassinosteriods, or proline proteases together with their potential regulatory elements. This work provides a first in-depth overview about molecular mechanisms underlying the heat stress response and acclimation in P. radiata.
Subject(s)
Pinus , Pinus/metabolism , Multiomics , Hot Temperature , Acclimatization/genetics , Heat-Shock Response/geneticsABSTRACT
How different stressors impact plant health and memory when they are imposed in different generations in wild ecosystems is still scarce. Here, we address how different environments shape heritable memory for the next generation in seeds and seedlings of Pinus radiata, a long-lived species with economic interest. The performance of the seedlings belonging to two wild clonal subpopulations (optimal fertirrigation vs. slightly stressful conditions) was tested under heat stress through physiological profiling and comparative time-series chloroplast proteomics. In addition, we explored the seeds conducting a physiological characterization and targeted transcriptomic profiling in both subpopulations. Our results showed differential responses between them, evidencing a cross-stress transgenerational memory. Seedlings belonging to the stressed subpopulation retained key proteins related to Photosystem II, chloroplast-to-nucleus signalling and osmoprotection which helped to overcome the applied heat stress. The seeds also showed a differential gene expression profile for targeted genes and microRNAs, as well as an increased content of starch and secondary metabolites, molecules which showed potential interest as biomarkers for early selection of primed plants. Thus, these finds not only delve into transgenerational cross-stress memory in trees, but also provide a new biotechnological tool for forest design.
Subject(s)
Ecosystem , Pinus , Female , Humans , Proteome/metabolism , Pinus/genetics , Droughts , Mothers , Nuclear Family , Seedlings/physiology , Heat-Shock Response , Seeds/genetics , Chloroplasts , Stress, PhysiologicalABSTRACT
BACKGROUND: Promoter hypermethylation of tumour suppressor genes is frequently observed during the malignant transformation of colorectal cancer (CRC). However, whether this epigenetic mechanism is functional in cancer or is a mere consequence of the carcinogenic process remains to be elucidated. RESULTS: In this work, we performed an integrative multi-omic approach to identify gene candidates with strong correlations between DNA methylation and gene expression in human CRC samples and a set of 8 colon cancer cell lines. As a proof of concept, we combined recent CRISPR-Cas9 epigenome editing tools (dCas9-TET1, dCas9-TET-IM) with a customized arrayed gRNA library to modulate the DNA methylation status of 56 promoters previously linked with strong epigenetic repression in CRC, and we monitored the potential functional consequences of this DNA methylation loss by means of a high-content cell proliferation screen. Overall, the epigenetic modulation of most of these DNA methylated regions had a mild impact on the reactivation of gene expression and on the viability of cancer cells. Interestingly, we found that epigenetic reactivation of RSPO2 in the tumour context was associated with a significant impairment in cell proliferation in p53-/- cancer cell lines, and further validation with human samples demonstrated that the epigenetic silencing of RSPO2 is a mid-late event in the adenoma to carcinoma sequence. CONCLUSIONS: These results highlight the potential role of DNA methylation as a driver mechanism of CRC and paves the way for the identification of novel therapeutic windows based on the epigenetic reactivation of certain tumour suppressor genes.
Subject(s)
Colonic Neoplasms , DNA Methylation , Humans , DNA Demethylation , Epigenesis, Genetic , Carcinogenesis , Mixed Function Oxygenases , Proto-Oncogene ProteinsABSTRACT
Due to the current climate change, many studies have described main drivers in abiotic stress. Recent findings suggest that alternative splicing (AS) has a critical role in controlling plant responses to high temperature. AS is a mechanism that allows organisms to create an assortment of RNA transcripts and proteins using a single gene. However, the most important roles of AS in stress could not be rigorously addressed because research has been focused on model species, covering only a narrow phylogenetic and lifecycle spectrum. Thus, AS degree of diversification among more dissimilar taxa in heat response is still largely unknown. To fill this gap, the present study employs a systems biology approach to examine how the AS landscape responds to and 'remembers' heat stress in conifers, a group which has received little attention even though their position can solve key evolutionary questions. Contrary to angiosperms, we found that potential intron retention may not be the most prevalent type of AS. Furthermore, our integrative analysis with metabolome and proteome data places splicing as the main source of variation during the response. Finally, we evaluated possible acquired long-term splicing memory in a diverse subset of events, and although this mechanism seems to be conserved in seed plants, AS dynamics are divergent. These discoveries reveal the particular way of remembering past temperature changes in long-lived plants and open the door to include species with unique features to determine the extent of conservation in gene expression regulation.
Subject(s)
Gene Expression Regulation, Plant , Pinus , Gene Expression Regulation, Plant/genetics , Pinus/genetics , Phylogeny , RNA Splicing , Heat-Shock Response/genetics , Plants/geneticsABSTRACT
The recovery and maintenance of plant homeostasis under stressful environments are complex processes involving organelle crosstalk for a coordinated cellular response. Here, we revealed through nuclear and chloroplast subcellular proteomics, biochemical cell profiles and targeted transcriptomics how chloroplasts and nuclei developed their responses under increased temperatures in a long-lived species (Pinus radiata). Parallel to photosynthetic impairment and reactive oxygen species production in the chloroplast, a DNA damage response was triggered in the nucleus followed by an altered chromatin conformation. In addition, in the nuclei, we found several proteins, such as HEMERA or WHIRLY, which change their locations from the chloroplasts to the nuclei carrying the stress message. Additionally, our data showed a deep rearrangement of RNA metabolism in both organelles, revealing microRNAs and AGO1 as potential regulators of the acclimation mechanisms. Altogether, our study highlights the synchronisation among the different stages required for thermotolerance acquisition in P. radiata, pointing out the role of chromatin conformation and posttranscriptional gene regulation in overcoming heat stress and assuring plant survival for the following years.
Subject(s)
Cell Nucleus/physiology , Chloroplasts/physiology , Heat-Shock Response , Pinus/physiology , Plant Proteins/physiology , Proteome/physiology , MicroRNAs/metabolism , RNA, Plant/metabolism , Signal TransductionABSTRACT
The elucidation of plant health status requires quantifying multiple molecular metabolism markers. Until now, the extraction of these biomarkers is performed independently, with different extractions and protocols. This approach is inefficient, since it increases laboratory time, amount of sample, and could introduce biases or difficulties when comparing data. To limit these drawbacks, we introduce a versatile protocol for quantifying seven of the most commonly analysed biomarkers (photosynthetic pigments, free amino acids, soluble sugars, starch, phenolic compounds, flavonoids and malondialdehyde) covering substantial parts of plant metabolism, requiring only a minimum sample amount and common laboratory instrumentation. The procedures of this protocol rely on classic methods that have been updated to allow their sequential use, increasing reproducibility, sensibility and easiness to obtain quantitative results. Our method has been tested and validated over an extended diversity of organisms (Arabidopsis thaliana, Solanum lycopersicum, Olea europaea, Quercus ilex, Pinus pinaster and Chlamydomonas reinhardtii), tissues (leaves, roots and seeds) and stresses (cold, drought, heat, ultraviolet B and nutrient deficiency). Its application will allow increasing the number of parameters that can be monitored at once while decreasing sample handling and consequently, increasing the capacity of the laboratory.
Subject(s)
Amino Acids/analysis , Coloring Agents/analysis , Flavonoids/analysis , Metabolomics/methods , Sugars/analysis , Chemical Fractionation/methods , Chlamydomonas reinhardtii/metabolism , Solanum lycopersicum/metabolism , Malondialdehyde/analysis , Olea/metabolism , Phenols/analysis , Quercus/metabolism , Reproducibility of ResultsABSTRACT
Bud maturation is a physiological process that implies a set of morphophysiological changes that lead to the transition of growth patterns from young to mature. This transition defines tree growth and architecture, and in consequence traits such as biomass production and wood quality. In Pinus pinaster Aiton, a conifer of great timber value, bud maturation is closely related to polycyclism (multiple growth periods per year). This process causes a lack of apical dominance, and consequently increased branching that reduces its timber quality and value. However, despite its importance, little is known about bud maturation. In this work, proteomics and metabolomics were employed to study apical and basal sections of young and mature buds in P. pinaster. Proteins and metabolites in samples were described and quantified using (n)UPLC-LTQ-Orbitrap. The datasets were analyzed employing an integrative statistical approach, which allowed the determination of the interactions between proteins and metabolites and the different bud sections and ages. Specific dynamics of proteins and metabolites such as histones H3 and H4, ribosomal proteins L15 and L12, chaperonin TCP1, 14-3-3 protein gamma, gibberellins A1, A3 and A8, strigolactones and abscisic acid, involved in epigenetic regulation, proteome remodeling, hormonal signaling and abiotic stress pathways showed their potential role during bud maturation. Candidates and pathways were validated employing interaction databases and targeted transcriptomics. These results increase our understanding of the molecular processes behind bud maturation, a key step towards improving timber production and natural pine forests management in a future scenario of climate change. However, further studies are necessary using different P. pinaster populations that show contrasting wood quality and stress tolerance in order to generalize the results.
Subject(s)
Pinus , Climate Change , Epigenesis, Genetic , Pinus/genetics , Trees , WoodABSTRACT
Based on the hypothesis that embryo development is a crucial stage for the formation of stable epigenetic marks that could modulate the behaviour of the resulting plants, in this study, radiata pine somatic embryogenesis was induced at high temperatures (23 °C, eight weeks, control; 40 °C, 4 h; 60 °C, 5 min) and the global methylation and hydroxymethylation levels of emerging embryonal masses and somatic plants were analysed using LC-ESI-MS/ MS-MRM. In this context, the expression pattern of six genes previously described as stress-mediators was studied throughout the embryogenic process until plant level to assess whether the observed epigenetic changes could have provoked a sustained alteration of the transcriptome. Results indicated that the highest temperatures led to hypomethylation of both embryonal masses and somatic plants. Moreover, we detected for the first time in a pine species the presence of 5-hydroxymethylcytosine, and revealed its tissue specificity and potential involvement in heat-stress responses. Additionally, a heat shock protein-coding gene showed a down-regulation tendency along the process, with a special emphasis given to embryonal masses at first subculture and ex vitro somatic plants. Likewise, the transcripts of several proteins related with translation, oxidative stress response, and drought resilience were differentially expressed.
ABSTRACT
Cork oak (Quercus suber L.) is a species of ecological, social and economic importance in the Mediterranean region. Given its xerophytic adaptability, the study of cork oak's response to drought stress conditions may provide important data in the global scenario of climate change. The mechanisms behind cork oak's adaptation to drought conditions can inform the design and development of tools to better manage this species under the changing climate patterns. Metabolomics is one of the most promising omics layers to capture a snapshot of a particular physiological state and to identify putative biomarkers of stress tolerance. Drastic changes were observed in the leaf metabolome of Q. suber between the different experimental conditions, namely at the beginning of the drought stress treatment, after one month under drought and post rehydration. All experimental treatments were analyzed through sPLS to inspect for global changes and stress and rehydration responses were analyzed independently for specific alterations. This allowed a more in-depth study and a search for biomarkers specific to a given hydric treatment. The metabolome analyses showed changes in both primary and secondary metabolism, but highlighted the role of secondary metabolism. In addition, a compound-specific response was observed in stress and rehydration. Key compounds such as L-phenylalanine and epigallocatechin 3-gallate were identified in relation to early drought response, terpenoid leonuridine and the flavonoid glycoside (-)-epicatechin-3'-O-glucuronide in long-term drought response, and flavone isoscoparine was identified in relation to the recovery process. The results here obtained provide novel insights into the biology of cork oak, highlighting pathways and metabolites potentially involved in the response of this species during drought and recovery that may be essential for its adaptation to long periods of drought. It is expected that this knowledge can encourage further functional studies in order to validate potential biomarkers of drought and recovery that maybe used to support decision-making in cork oak breeding programs.
Subject(s)
Droughts , Metabolome/physiology , Plant Proteins/metabolism , Quercus/physiology , Plant Leaves/physiology , Stress, PhysiologicalABSTRACT
BACKGROUND: The exposure of microalgae and plants to low UV-C radiation dosages can improve their biomass composition and stress tolerance. Despite UV-C sharing these effects with UV-A/B but at much lower dosages, UV-C sensing and signal mechanisms are still mostly unknown. Thus, we have described and integrated the proteometabolomic and physiological changes occurring in Chlamydomonas reinhardtii-a simple Plantae model-into the first 24 h after a short and low-intensity UV-C irradiation in order to reconstruct the microalgae response system to this stress. RESULTS: The microalgae response was characterized by increased redox homeostasis, ROS scavenging and protein damage repair/avoidance elements. These processes were upregulated along with others related to the modulation of photosynthetic electron flux, carbon fixation and C/N metabolism. These changes, attributed to either direct UV-C-, ROS- or redox unbalances-associated damage, trigger a response process involving novel signaling intermediaries and effectors such as the translation modulator FAP204, a PP2A-like protein and a novel DYRK kinase. These elements were found linked to the modulation of Chlamydomonas biomass composition (starch accumulation) and proliferation, within an UV-C response probably modulated by different epigenetic factors. CONCLUSION: Chosen multiomics integration approach was able to describe many fast changes, including biomass composition and ROS stress tolerance, as a response to a low-intensity UV-C stress. Moreover, the employed omics and systems biology approach placed many previously unidentified protein and metabolites at the center of these changes. These elements would be promising targets for the characterization of this stress response in microalgae and plants and the engineering of more productive microalgae strains.
ABSTRACT
Microalgae are gaining attention in industry for their high value-added biomolecules and biomass production and for studying fundamental processes in biology. The introduction of novel approaches for understanding and modeling molecular networks at different omic levels is paramount for increasing the productivity of these organisms. However, the construction of these networks requires high quality datasets with, if possible, perfectly overlapping datasets. The employ of different materials for different biomolecule isolation protocols, even if they come from the same homogenate, is one of the commonest issues affecting quality. Hence, a new method has been developed, allowing for the combined extraction of different levels including total metabolites, or their pigments or lipid fractions along nucleic acids (DNA and RNA) and/or proteins from the same sample reducing biological and time variation between levels data.
Subject(s)
High-Throughput Screening Assays/methods , Mass Spectrometry/methods , Microalgae/chemistry , Biological Factors/chemistry , Biomass , DNA/chemistry , Lipids/chemistry , Pigments, Biological/chemistry , Proteins/chemistry , RNA/chemistryABSTRACT
The evolution of next-generation sequencing and high-throughput technologies has created new opportunities and challenges in data science. Currently, a classic proteomics analysis can be complemented by going a step beyond the individual analysis of the proteome by using integrative approaches. These integrations can be focused either on inferring relationships among proteins themselves, with other molecular levels, phenotype, or even environmental data, giving the researcher new tools to extract and determine the most relevant information in biological terms. Furthermore, it is also important the employ of visualization methods that allow a correct and deep interpretation of data.To carry out these analyses, several bioinformatics and biostatistical tools are required. In this chapter, different workflows that enable the creation of interaction networks are proposed. Resulting networks reduce the complexity of original datasets, depicting complex statistical relationships (through PLS analysis and variants), functional networks (STRING, shinyGO), and a combination of both approaches. Recently developed methods for integrating different omics levels, such as coinertial analyses or DIABLO, are also described. Finally, the use of Cytoscape or Gephi was described for the representation and mining of the different networks.This approach constitutes a new way of acquiring a deeper knowledge of the function of proteins, such as the search for specific connections of each group to identify differentially connected modules, which may reflect involved protein complexes and key pathways.
Subject(s)
High-Throughput Screening Assays/methods , Protein Interaction Maps/physiology , Proteins/metabolism , Proteomics/methods , Amino Acid Sequence , Computational Biology/methods , Humans , Proteome/metabolism , SoftwareABSTRACT
The complexity of the plant cell proteome, exhibiting thousands of proteins whose abundance varies in several orders of magnitude, makes impossible to cover most of the plant proteins using standard shotgun-based approaches. Despite this general description of plant proteomes, the complexity is not a big issue (current protocols and instrumentation allow for the identification of several thousand proteins per injection), low or medium abundant proteins cannot be detected most of times, being necessary to fraction or perform targeted analyses in order to detect and quantify them. Among fractioning choices, cell fractioning in its different organelles is a good strategy for gaining not only a deeper coverage of the proteome but also the basis for understanding organelle function, protein dynamics, and trafficking within the cell, as nuclear and chloroplast communication. This approach is used routinely in many labs working with model species; however, the available protocols focusing on tree species are scarce. In this chapter, we provide a simple but robust protocol for isolating nuclei and chloroplasts in pine needles that is fully compatible with later mass spectrometry-based proteome analysis.
Subject(s)
Cell Nucleus/metabolism , Chloroplasts/metabolism , Proteome/metabolism , Tracheophyta/metabolism , Cell Fractionation/methods , Mass Spectrometry/methods , Plant Cells/metabolism , Plant Proteins/metabolism , Proteomics/methodsABSTRACT
Despite it being an important issue in the context of climate change, for most plant species it is not currently known how abiotic stresses affect nuclear proteomes and mediate memory effects. This study examines how Pinus radiata nuclei respond, adapt, 'remember', and 'learn' from heat stress. Seedlings were heat-stressed at 45 °C for 10 d and then allowed to recover. Nuclear proteins were isolated and quantified by nLC-MS/MS, the dynamics of tissue DNA methylation were examined, and the potential acquired memory was analysed in recovered plants. In an additional experiment, the expression of key gene genes was also quantified. Specific nuclear heat-responsive proteins were identified, and their biological roles were evaluated using a systems biology approach. In addition to heat-shock proteins, several clusters involved in regulation processes were discovered, such as epigenomic-driven gene regulation, some transcription factors, and a variety of RNA-associated functions. Nuclei exhibited differential proteome profiles across the phases of the experiment, with histone H2A and methyl cycle enzymes in particular being accumulated in the recovery step. A thermopriming effect was possibly linked to H2A abundance and over-accumulation of spliceosome elements in recovered P. radiata plants. The results suggest that epigenetic mechanisms play a key role in heat-stress tolerance and priming mechanisms.
Subject(s)
Pinus , Proteome , Epigenesis, Genetic , Gene Expression Regulation, Plant , Pinus/genetics , Pinus/metabolism , Plant Proteins/genetics , Plant Proteins/metabolism , Stress, Physiological , Tandem Mass SpectrometryABSTRACT
Seed germination and early seedling development have been studied in the recalcitrant species Quercus ilex using targeted transcriptional, hormonal, and sugar analysis. Embryos and seedlings were collected at eight morphologically defined developmental stages, S0-S7. A typical triphasic water uptake curve was observed throughout development, accompanied by a decrease in sucrose and an increase in glucose and fructose. Low levels of abscisic acid (ABA) and high levels of gibberellins (GAs) were observed in mature seeds. Post-germination, indole-3-acetic acid (IAA), increased, whereas GA remained high, a pattern commonly observed during growth and development. The abundance of transcripts from ABA-related genes was positively correlated with the changes in the content of the phytohormone. Transcripts of the drought-related genes Dhn3 and GolS were more abundant at S0, then decreased in parallel with increasing water content. Transcripts for Gapdh, and Nadh6 were abundant at S0, supporting the occurrence of an active metabolism in recalcitrant seeds at the time of shedding. The importance of ROS during germination is manifest in the high transcript levels for Sod and Gst, found in mature seeds. The results presented herein help distinguish recalcitrant (e.g., Q. ilex) seeds from their orthodox counterparts. Our results indicate that recalcitrance is established during seed development but not manifest until germination (S1-S3). Post-germination the patterns are quite similar for both orthodox and recalcitrant seeds.
ABSTRACT
The predicted accentuation of drought events highlights the importance of optimize plants capacity to tolerate drought, but also the capacity to recovery from it, especially in species, as olive tree (Olea europaea L.), that grows in particularly susceptible regions. Three different concentrations (10, 100 and 1000 µM) of salicylic acid (SA), a stress signaling phytohormone, was sprayed on 3-year-old potted olive trees subjected to three successive drought and rewatering events. Trees responses to SA application are concentration dependent, being 100 µM the most effective concentration to improve drought tolerance and recovery capacity. During drought events, this effectiveness was achieved by osmolytes accumulation, leaf water status maintenance, reduced photosynthetic systems drought-associated damages, and by optimizing shoot/root ratio. The better plant fitness during drought allowed a fast recovery of the physiological functions upon rewatering and reduced the necessity to invest in extra repair damages, allowing the regrowth. The intense abscisic acid (ABA) signal close to upper epidermis in stressed controls suggests a "memory" of the worst water status displayed by those plants. SA attenuated the limitation of total biomass accumulation imposed by drought, mainly in root system, increased water use efficiency and lead to a higher intense signal of indoleacetic acid (IAA) in leaves during recovery period. In summary, in a suitable concentration, SA demonstrate to be a promising tool to increase drought adaptability of olive trees.
Subject(s)
Olea/growth & development , Plant Growth Regulators/pharmacology , Salicylic Acid/pharmacology , Dehydration , Dose-Response Relationship, Drug , Olea/drug effects , Olea/physiology , Plant Leaves/metabolism , Water/metabolismABSTRACT
The integrative omics approach is crucial to identify the molecular mechanisms underlying high-temperature response in non-model species. Based on future scenarios of heat increase, Pinus radiata plants were exposed to a temperature of 40°C for a period of 5 days, including recovered plants (30 days after last exposure to 40°C) in the analysis. The analysis of the metabolome using complementary mass spectrometry techniques (GC-MS and LC-Orbitrap-MS) allowed the reliable quantification of 2,287 metabolites. The analysis of identified metabolites and highlighter metabolic pathways across heat time exposure reveal the dynamism of the metabolome in relation to high-temperature response in P. radiata, identifying the existence of a turning point (on day 3) at which P. radiata plants changed from an initial stress response program (shorter-term response) to an acclimation one (longer-term response). Furthermore, the integration of metabolome and physiological measurements, which cover from the photosynthetic state to hormonal profile, suggests a complex metabolic pathway interaction network related to heat-stress response. Cytokinins (CKs), fatty acid metabolism and flavonoid and terpenoid biosynthesis were revealed as the most important pathways involved in heat-stress response in P. radiata, with zeatin riboside (ZR) and isopentenyl adenosine (iPA) as the key hormones coordinating these multiple and complex interactions. On the other hand, the integrative approach allowed elucidation of crucial metabolic mechanisms involved in heat response in P. radiata, as well as the identification of thermotolerance metabolic biomarkers (L-phenylalanine, hexadecanoic acid, and dihydromyricetin), crucial metabolites which can reschedule the metabolic strategy to adapt to high temperature.
ABSTRACT
Pinus radiata seedlings, the most widely planted pine species in the world, were exposed to temperatures within a range mimicking future scenarios based on current models of heat increase. The short-term heat response in P. radiata was studied in detail by exploring the metabolome, proteome and targeted transcriptome. The use of complementary mass spectrometry techniques, GC-MS and LC-Orbitrap-MS, together with novel bioinformatics tools allowed the reliable quantification of 2,075 metabolites and 901 protein groups. Integrative analyses of different functional levels and plant physiological status revealed a complex molecular interaction network of positive and negative correlations between proteins and metabolites involved in short-term heat response, including three main physiological functions as: 1) A hormone subnetwork, where fatty acids, flavonoids and hormones presented a key role; 2) An oxidoreductase subnetwork, including several dehydrogenase and peroxidase proteins; and 3) A heat shock protein subnetwork, with numerous proteins that contain a HSP20 domain, all of which were overexpressed at the transcriptional level. Integrated analysis pinpointed the basic mechanisms underlying the short-term physiological reaction of P. radiata during heat response. This approach was feasible in forest species and unmasked two novel candidate biomarkers of heat resistance, PHO1 and TRANSCRIPTION FACTOR APFI, and a MITOCHONDRIAL SMALL HEAT SHOCK PROTEIN, for use in future breeding programs.