Like mother like son: Transgenerational memory and cross-tolerance from drought to heat stress are identified in chloroplast proteome and seed provisioning in Pinus radiata.

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.

Multiomics analyses reveal the central role of the nucleolus and its machinery during heat stress acclimation in Pinus radiata.

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.

Allele-Specific Reprogramming of Cancer Metabolism by the Long Non-coding RNA CCAT2.

Altered energy metabolism is a cancer hallmark as malignant cells tailor their metabolic pathways to meet their energy requirements. Glucose and glutamine are the major nutrients that fuel cellular metabolism, and the pathways utilizing these nutrients are often altered in cancer. Here, we show that the long ncRNA CCAT2, located at the 8q24 amplicon on cancer risk-associated rs6983267 SNP, regulates cancer metabolism in vitro and in vivo in an allele-specific manner by binding the Cleavage Factor I (CFIm) complex with distinct affinities for the two subunits (CFIm25 and CFIm68). The CCAT2 interaction with the CFIm complex fine-tunes the alternative splicing of Glutaminase (GLS) by selecting the poly(A) site in intron 14 of the precursor mRNA. These findings uncover a complex, allele-specific regulatory mechanism of cancer metabolism orchestrated by the two alleles of a long ncRNA.

Integrative analysis in Pinus revealed long-term heat stress splicing memory.

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.

Biomarker metabolite mating of viable frozen-thawed in vitro-produced bovine embryos with pregnancy-competent recipients leads to improved birth rates.

Selection of competent recipients before embryo transfer (ET) is indispensable for improving pregnancy and birth rates in cattle. However, pregnancy prediction can fail when the competence of the embryo is ignored. We hypothesized that the pregnancy potential of biomarkers could improve with information on embryonic competence. In vitro-produced embryos cultured singly for 24 h (from d 6 to 7) were transferred to d 7 synchronized recipients as fresh or after freezing and thawing. Recipient blood was collected on d 0 (estrus; n = 108) and d 7 (4-6 h before ET; n = 107) and plasma was analyzed by nuclear magnetic resonance (1H+NMR). Spent embryo culture medium (CM) was collected and analyzed by ultra-high-performance liquid chromatography tandem mass spectrometry in a subset of n = 70 samples. Concentrations of metabolites quantified in plasma (n = 35) were statistically analyzed as a function of pregnancy diagnosed on d 40, d 62 and birth. Univariate analysis with plasma metabolites consisted of a block study with controllable fixed factors (i.e., embryo cryopreservation, recipient breed, and day of blood collection; Wilcoxon test and t-test). Metabolite concentrations in recipients and embryos were independently analyzed by iterations that reclassified embryos or recipients using the support vector machine. Iterations identified some competent embryos, but mostly competent recipients that had a pregnancy incompetent partner embryo. Misclassified recipients that could be classified as competent were reanalyzed in a new iteration to improve the predictive model. After subsequent iterations, the predictive potential of recipient biomarkers was recalculated. On d 0, creatine, acetone and l-phenylalanine were the most relevant biomarkers at d 40, d 62, and birth, and on d 7, l-glutamine, l-lysine, and ornithine. Creatine was the most representative biomarker within blocks (n = 20), with a uniform distribution over pregnancy endpoints and type of embryos. Biomarkers showed higher abundance on d 7 than d 0, were more predictive for d 40 and d 62 than at birth, and the pregnancy predictive ability was lower with frozen-thawed (F-T) embryos. Six metabolic pathways differed between d 40 pregnant recipients for fresh and F-T embryos. Within F-T embryos, more recipients were misclassified, probably due to pregnancy losses, but were accurately identified when combined with embryonic metabolite signals. After recalculation, 12 biomarkers increased receiver operator characteristic-area under the curve (>0.65) at birth, highlighting creatine (receiver operator characteristic-area under the curve = 0.851), and 5 new biomarkers were identified. Combining metabolic information of recipient and embryos improves the confidence and accuracy of single biomarkers.

Nucleus and chloroplast: A necessary understanding to overcome heat stress in Pinus radiata.

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.

Impaired condensin complex and Aurora B kinase underlie mitotic and chromosomal defects in hyperdiploid B-cell ALL.

B-cell acute lymphoblastic leukemia (ALL; B-ALL) is the most common pediatric cancer, and high hyperdiploidy (HyperD) identifies the most common subtype of pediatric B-ALL. Despite HyperD being an initiating oncogenic event affiliated with childhood B-ALL, the mitotic and chromosomal defects associated with HyperD B-ALL (HyperD-ALL) remain poorly characterized. Here, we have used 54 primary pediatric B-ALL samples to characterize the cellular-molecular mechanisms underlying the mitotic/chromosome defects predicated to be early pathogenic contributors in HyperD-ALL. We report that HyperD-ALL blasts are low proliferative and show a delay in early mitosis at prometaphase, associated with chromosome-alignment defects at the metaphase plate leading to robust chromosome-segregation defects and nonmodal karyotypes. Mechanistically, biochemical, functional, and mass-spectrometry assays revealed that condensin complex is impaired in HyperD-ALL cells, leading to chromosome hypocondensation, loss of centromere stiffness, and mislocalization of the chromosome passenger complex proteins Aurora B kinase (AURKB) and Survivin in early mitosis. HyperD-ALL cells show chromatid cohesion defects and an impaired spindle assembly checkpoint (SAC), thus undergoing mitotic slippage due to defective AURKB and impaired SAC activity, downstream of condensin complex defects. Chromosome structure/condensation defects and hyperdiploidy were reproduced in healthy CD34+ stem/progenitor cells upon inhibition of AURKB and/or SAC. Collectively, hyperdiploid B-ALL is associated with a defective condensin complex, AURKB, and SAC.

Non-Invasive Identification of Sex in Cultured Bovine Embryos by UHPLC-MS/MS Metabolomics.

INTRODUCTION: Different gene expression between male and female bovine embryos leads to metabolic differences. OBJECTIVE: We used UHPLC-MS/MS to identify sex metabolite biomarkers in embryo culture medium (CM). METHODS: Embryos were produced in vitro under highly variable conditions, i.e., fertilized with 7 bulls, two breeds, and cultured with BSA or BSA + serum until Day-6. On Day-6, embryos were cultured individually for 24 h. CM of Day-7 embryos (86 female and 81 male) was collected, and Day-6 and Day-7 embryonic stages recorded. RESULTS: A study by sample subsets with fixed factors (culture, bull breed, and Day-6 and Day-7 stages) tentatively identified 31 differentially accumulated metabolites through 182 subsets. Day-6 and Day-7 stage together affected 13 and 11 metabolites respectively, while 19 metabolites were affected by one or another stage and/or day. Culture supplements and individual bull changed 19 and 15 metabolites, respectively. Single bull exerted the highest influence (20 metabolites with the significantly highest p values). Lipid (93 subsets; 11 metabolites) and amino acid (55 subsets; 13 metabolites) were the most relevant classes for sex identification. CONCLUSIONS: Single biomarker led to inefficient sex diagnosis, while metabolite combinations accurately identified sex. Our study is a first in non-invasive sex identification in cattle by overcoming factors that induce metabolic variation.

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.

The rainbow protocol: A sequential method for quantifying pigments, sugars, free amino acids, phenolics, flavonoids and MDA from a small amount of sample.

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.

Comparative proteomics of Pinus-Fusarium circinatum interactions reveal metabolic clues to biotic stress resistance.

Fusarium circinatum, causing pine pitch canker (PPC), affects conifers productivity and health worldwide. Selection and breeding for resistance arises as the most promising approach to fight PPC. Therefore, it is crucial to explore the response of hosts with varying levels of susceptibility to PPC to unveil the genes/pathways behind these phenotypes. We evaluated the dynamics of the needle proteome of a susceptible (Pinus radiata) and a relatively resistant (Pinus pinea) species upon F. circinatum inoculation by GeLC-MS/MS. Integration with physiological data and validation of key genes by qPCR allowed to identify core pathways regulating these contrasting responses. In P. radiata, the pathogen may target both the secondary metabolism to negatively regulate immune response and chloroplast redox proteins to increase energy-producing pathways for amino acid production in its favour. In contrast, chloroplast redox regulation may assure redox homeostasis in P. pinea, as well as nonenzymatic antioxidants. The presence of membrane trafficking-related proteins exclusively in P. pinea likely explains its defence response against F. circinatum. A crosstalk between abscisic acid and epigenetic regulation of gene expression is also proposed in PPC response. These results are useful to support breeding programs aiming to achieve PPC resistance.

Chromatin regulation by Histone H4 acetylation at Lysine 16 during cell death and differentiation in the myeloid compartment.

Histone H4 acetylation at Lysine 16 (H4K16ac) is a key epigenetic mark involved in gene regulation, DNA repair and chromatin remodeling, and though it is known to be essential for embryonic development, its role during adult life is still poorly understood. Here we show that this lysine is massively hyperacetylated in peripheral neutrophils. Genome-wide mapping of H4K16ac in terminally differentiated blood cells, along with functional experiments, supported a role for this histone post-translational modification in the regulation of cell differentiation and apoptosis in the hematopoietic system. Furthermore, in neutrophils, H4K16ac was enriched at specific DNA repeats. These DNA regions presented an accessible chromatin conformation and were associated with the cleavage sites that generate the 50 kb DNA fragments during the first stages of programmed cell death. Our results thus suggest that H4K16ac plays a dual role in myeloid cells as it not only regulates differentiation and apoptosis, but it also exhibits a non-canonical structural role in poising chromatin for cleavage at an early stage of neutrophil cell death.

Dual RNA-Sequencing Analysis of Resistant (Pinus pinea) and Susceptible (Pinus radiata) Hosts during Fusarium circinatum Challenge.

Fusarium circinatum causes one of the most important diseases of conifers worldwide, the pine pitch canker (PPC). However, no effective field intervention measures aiming to control or eradicate PPC are available. Due to the variation in host genetic resistance, the development of resistant varieties is postulated as a viable and promising strategy. By using an integrated approach, this study aimed to identify differences in the molecular responses and physiological traits of the highly susceptible Pinus radiata and the highly resistant Pinus pinea to F. circinatum at an early stage of infection. Dual RNA-Seq analysis also allowed to evaluate pathogen behavior when infecting each pine species. No significant changes in the physiological analysis were found upon pathogen infection, although transcriptional reprogramming was observed mainly in the resistant species. The transcriptome profiling of P. pinea revealed an early perception of the pathogen infection together with a strong and coordinated defense activation through the reinforcement and lignification of the cell wall, the antioxidant activity, the induction of PR genes, and the biosynthesis of defense hormones. On the contrary, P. radiata had a weaker response, possibly due to impaired perception of the fungal infection that led to a reduced downstream defense signaling. Fusarium circinatum showed a different transcriptomic profile depending on the pine species being infected. While in P. pinea, the pathogen focused on the degradation of plant cell walls, active uptake of the plant nutrients was showed in P. radiata. These findings present useful knowledge for the development of breeding programs to manage PPC.

Integrative analysis of the nuclear proteome in Pinus radiata reveals thermopriming coupled to epigenetic regulation.

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.

Integrated Physiological, Proteomic, and Metabolomic Analysis of Ultra Violet (UV) Stress Responses and Adaptation Mechanisms in Pinus radiata.

Globally expected changes in environmental conditions, especially the increase of UV irradiation, necessitate extending our knowledge of the mechanisms mediating tree species adaptation to this stress. This is crucial for designing new strategies to maintain future forest productivity. Studies focused on environmentally realistic dosages of UV irradiation in forest species are scarce. Pinus spp. are commercially relevant trees and not much is known about their adaptation to UV. In this work, UV treatment and recovery of Pinus radiata plants with dosages mimicking future scenarios, based on current models of UV radiation, were performed in a time-dependent manner. The combined metabolome and proteome analysis were complemented with measurements of + physiological parameters and gene expression. Sparse PLS analysis revealed complex molecular interaction networks of molecular and physiological data. Early responses prevented phototoxicity by reducing photosystem activity and the electron transfer chain together with the accumulation of photoprotectors and photorespiration. Apart from the reduction in photosynthesis as consequence of the direct UV damage on the photosystems, the primary metabolism was rearranged to deal with the oxidative stress while minimizing ROS production. New protein kinases and proteases related to signaling, coordination, and regulation of UV stress responses were revealed. All these processes demonstrate a complex molecular interaction network extending the current knowledge on UV-stress adaptation in pine.

Proteomics and comparative genomics of Nitrososphaera viennensis reveal the core genome and adaptations of archaeal ammonia oxidizers.

Ammonia-oxidizing archaea (AOA) are among the most abundant microorganisms and key players in the global nitrogen and carbon cycles. They share a common energy metabolism but represent a heterogeneous group with respect to their environmental distribution and adaptions, growth requirements, and genome contents. We report here the genome and proteome of Nitrososphaera viennensis EN76, the type species of the archaeal class Nitrososphaeria of the phylum Thaumarchaeota encompassing all known AOA. N. viennensis is a soil organism with a 2.52-Mb genome and 3,123 predicted protein-coding genes. Proteomic analysis revealed that nearly 50% of the predicted genes were translated under standard laboratory growth conditions. Comparison with genomes of closely related species of the predominantly terrestrial Nitrososphaerales as well as the more streamlined marine Nitrosopumilales [Candidatus (Ca.) order] and the acidophile "Ca. Nitrosotalea devanaterra" revealed a core genome of AOA comprising 860 genes, which allowed for the reconstruction of central metabolic pathways common to all known AOA and expressed in the N. viennensis and "Ca Nitrosopelagicus brevis" proteomes. Concomitantly, we were able to identify candidate proteins for as yet unidentified crucial steps in central metabolisms. In addition to unraveling aspects of core AOA metabolism, we identified specific metabolic innovations associated with the Nitrososphaerales mediating growth and survival in the soil milieu, including the capacity for biofilm formation, cell surface modifications and cell adhesion, and carbohydrate conversions as well as detoxification of aromatic compounds and drugs.

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.

Natural variation of DNA methylation and gene expression may determine local adaptations of Scots pine populations.

Long-lived conifers are vulnerable to climate change because classical evolutionary processes are slow in developing adaptive responses. Therefore, the capacity of a genotype to adopt different phenotypes is important. Gene expression is the primary mechanism that converts genome-encoded information into phenotypes, and DNA methylation is employed in the epigenetic regulation of gene expression. We investigated variations in global DNA methylation and gene expression between three Scots pine (Pinus sylvestris L.) populations located in northern and southern Finland using mature seeds. Gene expression levels were studied in six DNA methyltransferase (DNMT) genes, which were characterized in this study, and in 19 circadian clock genes regulating adaptive traits. In embryos, expression diversity was found for three DNMT genes, which maintain DNA methylation. The expression of two DNMT genes was strongly correlated with climate variables, which suggests a role for DNA methylation in local adaptation. For adaptation-related genes, expression levels showed between-population variation in 11 genes in megagametophytes and in eight genes in embryos, and many of these genes were linked to climate factors. Altogether, our results suggest that differential DNA methylation and gene expression contribute to local adaptation in Scots pine populations and may enhance the fitness of trees under rapidly changing climatic conditions.

System-wide analysis of short-term response to high temperature in Pinus radiata.

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.