Next Article in Journal
Impaired Intracellular Ca2+ Dynamics, M-Band and Sarcomere Fragility in Skeletal Muscles of Obscurin KO Mice
Next Article in Special Issue
Full-Length Transcriptome and RNA-Seq Analyses Reveal the Mechanisms Underlying Waterlogging Tolerance in Kiwifruit (Actinidia valvata)
Previous Article in Journal
The Versatility in the Applications of Dithiocarbamates
Previous Article in Special Issue
Molecular Basis to Integrate Microgravity Signals into the Photoperiodic Flowering Pathway in Arabidopsis thaliana under Spaceflight Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Changing Temperature Conditions during Somatic Embryo Maturation Result in Pinus pinaster Plants with Altered Response to Heat Stress

1
Agrarian and Environmental Sciences Department, Institute for Research on Environmental Sciences (IUCA), High Polytechnic School, University of Zaragoza, Ctra. Cuarte s/n, 22071 Huesca, Spain
2
Plant Biology Department, Biotechnology and Biomedicine (BiotecMed) Institute, Universitat de València, Vicent Andrés Estellés s/n, 46100 Burjassot, Spain
3
Department of Life Sciences, Universidade de Coimbra, 3000-456 Coimbra, Portugal
4
Department of Forestry Science, NEIKER, 01192 Arkaute, Spain
5
Departamento de Producción Vegetal, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
6
Laboratory of Growth Regulators, Faculty of Science of Palacký University, Institute of Experimental Botany of the Czech Academy of Sciences, Šlechtitelů 27, 78371 Olomouc, Czech Republic
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(3), 1318; https://doi.org/10.3390/ijms23031318
Submission received: 29 December 2021 / Revised: 22 January 2022 / Accepted: 23 January 2022 / Published: 24 January 2022
(This article belongs to the Special Issue Biotic and Abiotic Stress Effects on Plant Structure and Physiology)

Abstract

:
Under the global warming scenario, obtaining plant material with improved tolerance to abiotic stresses is a challenge for afforestation programs. In this work, maritime pine (Pinus pinaster Aiton) plants were produced from somatic embryos matured at different temperatures (18, 23, or 28 °C, named after M18, M23, and M28, respectively) and after 2 years in the greenhouse a heat stress treatment (45 °C for 3 h/day for 10 days) was applied. Temperature variation during embryo development resulted in altered phenotypes (leaf histology, proline content, photosynthetic rates, and hormone profile) before and after stress. The thickness of chlorenchyma was initially larger in M28 plants, but was significantly reduced after heat stress, while increased in M18 plants. Irrespective of their origin, when these plants were subjected to a heat treatment, relative water content (RWC) and photosynthetic carbon assimilation rates were not significantly affected, although M18 plants increased net photosynthesis rate after 10 days recovery (tR). M18 plants showed proline contents that increased dramatically (2.4-fold) when subjected to heat stress, while proline contents remained unaffected in M23 and M28 plants. Heat stress significantly increased abscisic acid (ABA) content in the needles of maritime pine plants (1.4-, 3.6- and 1.9-fold in M18, M23, and M28 plants, respectively), while indole-3-acetic acid content only increased in needles from M23 plants. After the heat treatment, the total cytokinin contents of needles decreased significantly, particularly in M18 and M28 plants, although levels of active forms (cytokinin bases) did not change in M18 plants. In conclusion, our results suggest that maturation of maritime pine somatic embryos at lower temperature resulted in plants with better performance when subjected to subsequent high temperature stress, as demonstrated by faster and higher proline increase, lower increases in ABA levels, no reduction in active cytokinin, and a better net photosynthesis rate recovery.

1. Introduction

Maritime pine (Pinus pinaster Aiton) is the most abundant conifer in the Mediterranean basin, showing a surprisingly ample ecological niche that covers more than 4 million ha. This species grows from mild environments in the Atlantic to harsh, dry Mediterranean climates, and from sea level to more than 2000 m in the southern regions [1]. As a consequence, large genetic variation in adaptive traits has been found between provenances for this plant species, which has an important ecological and economic value [2,3,4]. Accordingly, maritime pine is used as a model species for the study of conifer adaptation to abiotic stress [5,6]. For instance, maritime pine is the most drought-tolerant Pinus species in Southern Europe forests, where a combination of high temperatures and droughts has been linked to a decrease in tree growth driving forest mortality [7,8]. These adverse conditions are expected to be aggravated in the near future under the global warming scenario, therefore obtaining plant material with improved tolerance to heat and hydric stresses is a challenge for afforestation programs in this region.
Breeding of forest species by conventional techniques is hampered by the large size of the individuals and their long-generation time. In addition, traits related to abiotic stress tolerance are regulated by multiple loci that are widespread across the whole genome [9]. Biotechnological tools, such as somatic embryogenesis, are alternative approaches for producing superior individuals in conifers [10]. Somatic embryogenesis (SE) is a vegetative propagation technology that can be used to induce multiple embryos from somatic cells, and that allows the implementation of deployment strategies for improved reforestation materials [11]. One of these strategies, known as “priming,” is the induction of adaptive responses by exposing plant material at early embryogenic stages to stressful conditions, which results in an epigenetic memory that makes plants more tolerant to future stress exposure [12]. This memory is mediated by mechanisms such as sustained alterations in gene expression, changes in hormonal profiles, and accumulation of signaling proteins, and transcription factors [13]. Evidence of partly epigenetic regulation of some traits, such as the timing of bud set and active growth cessation, that were affected by environmental conditions (temperature and photoperiod) during embryo development [14], provided a way to modify cell lines by adjusting the temperature during SE processes. This induced variation is added to the unintentional selection among cell lines due to genotypical differences in the ability to perform through the process that operates in all SE stages (initiation, proliferation, and maturation), as well as in germination, plant development, and acclimatization steps, affecting overall plant yield [15].
Efficient SE protocols have been developed for commercial pine species [11], and in particular to maritime pine materials from different Spanish provenances [6,16,17]. In addition, we have recently reported that a short-time high-temperature priming applied to megagametophytes induced heat tolerance in SE-derived maritime pine plants [18]. The high temperatures used for priming in these experiments (up to 50 °C) affected the efficiency of the SE process, as occurred when milder temperatures were employed along the different steps of the protocol. Thus, increasing the temperature from 23 to 28 °C during the induction and proliferation phases of P. pinaster somatic embryogenesis resulted in higher rates of embryo production [6]. In Pinus radiata, temperature conditions at the induction step also affected the production of somatic embryos [19], whereas priming with high temperatures during the maturation phase in P. halepensis and P. radiata did not affect the number of embryos obtained; however, the regenerated plants showed better performance under drought stress [20].
Heat stress in tree species not only affects plant morphology, but also alters leaf anatomy, photosynthesis, dark respiration, stomatal conductance, and transpiration, among others [21]. At cellular and subcellular levels, the reduction in water content caused by heat negatively influences cell division and growth, as well as the functioning of photosystem II and rubisco activity. Other deleterious effects are protein denaturation, enzyme inactivation, or membrane dysfunction [22]. In response to stress, plants could alter their metabolisms, particularly by producing compatible solutes, such as proline, that are able to maintain cell turgor by osmotic adjustment. Furthermore, proline has additional roles, such as stabilizing proteins, membranes, and subcellular structures, protecting cellular functions by scavenging ROS, and acting as a signalling molecule [23]. Thus, proline plays an important role in maintaining the metabolism and growth of plants under abiotic stress conditions and many reports indicate a positive relationship between proline accumulation and tolerance of plants to various abiotic stresses [24].
Plant hormones involved in heat stress response include ABA, brassinosteroids, CKs, salicylic acid, jasmonic acid, and ethylene. Abscisic acid (ABA) is the main regulator of abiotic stress tolerance in plants. ABA signaling can be very fast without involving transcriptional activity. A good example is the control of stomatal aperture by this hormone through the biochemical regulation of ion and water transport processes. In addition, ABA regulates biomass allocation in woody plants, mainly affecting the leaf area [25]. Heat stress studies of many other plants have recently focused on physiological and molecular changes, while little is known about how heat stress affects anatomical structures [26].
The effects of various hormones on plant stress tolerance have been characterized in many plant species, but the results are often inconsistent. This could be due to variations of plant species and cultivars, growth stages, stress intensity and duration, and possible interactions among hormones [27]. Priming for improving thermotolerance in Pinus spp. by short heat treatments during SE induction affected DNA methylation and changed the expression pattern of stress-related genes in Pinus halepensis [28] and in P. radiata [29], and also altered gene expression patterns and the hormonal profile in P. pinaster plants [18]. Temperature during induction also had significant effects on the cytokinin (CK) profile of embryogenic lines of Pinus halepensis [30] and P. radiata [31]. Furthermore, CKs were involved in the drought tolerance of radiata pine plants derived from primed somatic embryos [32]. Temperature during maturation also affected the CK and auxin profiles of radiata pine somatic embryos [33], and thermotolerance acquisition was related to hormone profile, photosynthetic pigments, and carbohydrate contents [34,35]. Kolb and Robberecht [36] reported that the survival of seedlings of Pinus ponderosa in high temperatures and droughts was mediated by higher stomatal conductance and transpiration rates.
In this work, we studied the response to heat stress of maritime pine plants derived from somatic embryos induced and proliferated at 28 °C, but matured at either 18, 23, or 28 °C. Our aim was to characterize histological and physiological changes induced in SE- derived maritime pine plants by exposure to low or high temperatures during the maturation phase of the propagation protocol, and to analyze whether these changes play a role in plant adaptation to high temperature conditions.

2. Results

2.1. Plant Height and Histological Determinations in Maritime Pine Needles

Variation in maturation temperature during maritime pine somatic embryo production affected plant development, in such a way that M23 plants were on average higher than M18 and M28 plants (11.9 ± 2.2 cm front to 8.8 ± 2.6 and 9.1 ± 2.4 cm, respectively; p < 0.001).
Histological analysis also demonstrated that needles from M28 plants where thicker than those from M23 and M18 plants (Table 1 and Figure 1a–c). Nevertheless, after heat treatment (t10), needle thickness was reduced in M28 but increased in M18 and M23 plants (Table 1 and Figure 1d–f). When focusing on specific tissues, we found that the epidermis cell layer in needles from M18 plants (Figure 1a) was significantly thicker (p = 0.018) than in M23 (Figure 1b) and M28 plants (Figure 1c). After heat stress, the thickness of this tissue was significantly reduced in M18 plants, as well as in needles from M28 plants (Figure 1d,f; Table 1; p = 0.045 and p < 0.001, respectively). In contrast, needles from plants derived from the SE process performed in standard temperature conditions (M23 plants) did not show this histological variation when submitted to heat stress (Figure 1b,e).
Significant differences in thickness of chlorenchyma were also observed in transversal sections of needles sampled in maritime pine plants derived from somatic embryos matured at different temperatures. In M28 plants, this leaf tissue was initially larger than in M18 and M23 plants (Table 1; Figure 1a–c). At the end of the heat stress treatment, the thickness of chlorenchyma was significantly reduced in M28 plants, but increased in M18 and M23 plants (Table 1 and Figure 1d–f).

2.2. Water Relations and Osmotic Adjustment

The relative water content (RWC) of P. pinaster needles was similar in plants derived from SE matured at different temperatures (p = 0.152) and was not affected by the 10-day heat treatment (p = 0.260), with a mean value of 74.0% ± 5.0% (Supplementary Figure S1).
The proline content of needles varied among groups of maritime pine plants (Figure 2), and along the heat stress experiment. The basal levels of this amino acid were significantly lower in M18 than in M23 and M28 plants (p = 0.019), but at the end of the high temperature treatment (t10), the plants derived from somatic embryos matured at 18 °C showed significantly higher levels of this amino acid than did the M23 and M28 plants (0.76 ± 0.03 front to 0.29 ± 0.03 and 0.39 ± 0.11µg/g FW, respectively; p < 0.001). The great increase (2.4-fold) observed in M18 plants was maintained after the recovery time (p < 0.001), while in M23 and M28 plants the proline content did not change after the stress, and finally increased at tR in M23 plants (1.5-fold, p = 0.002) but decreased significantly in M28 plants (p = 0.014).

2.3. Photosynthesiselated Parameters

Several parameters related to photosynthesis (AN, net photosynthetic CO2 assimilation, gs, stomatal conductance, and Ci substomatal CO2) were measured at three point-times during the experiment: before the heat treatment (t0), at the end of the heat treatment (t10), and after the recovery period (tR).
Maritime pine M18 plants showed, on average, lower photosynthetic rates than M23 and M28 plants (p = 0.008); in particular, the initial photosynthetic rates in these plants were significantly lower than those determined in plants obtained under standard temperature conditions (7.4 ± 2.4 front to 10.9 ± 3.8 μmol CO2 m−2 s−1, p = 0.037). Although this parameter was slightly reduced after stress (t10 samples) in the three groups of plants, heat treatment did not have a significant effect on net photosynthetic CO2 assimilation, and initial levels were regained in tR samples (Table 2). Similarly, heat stress reduced stomatal conductance, but the effect was significant only in M28 plants (p = 0.032). Levels of substomatal CO2 concentration (Ci) were significantly reduced in tR samples from M18 (p = 0.034) and M28 (p = 0.024) plants, but not in needles from M23 plants (Table 2). Higher photosynthetic rates and Ci concentration were significantly correlated to higher stomatal conductance (Spearman coefficient of correlation ρ = 0.466, and ρ = 0.698, respectively, p < 0.001). Furthermore, Ci decreased with stomatal conductance by the effect of heat stress, suggesting stomatal limitations to net photosynthetic rates.

2.4. Hormone Content

The levels of three hormones, ABA, IAA, and CKs, were also determined in maritime pine needles, before and after the heat stress treatment.
Needles sampled in maritime pine plants derived from SE matured at 28 °C (M28) showed basal levels of ABA (t0) significantly higher than those from M18 and M23 plants (90.1 ± 19.2 pmol/g FW front to 55.5 ± 5.7 and 41.0 ± 8.7 pmol/g FW, respectively). Heat stress treatment significantly increased ABA levels in the three groups of plants (p < 0.001), particularly in those derived from SE matured in standard conditions (M23), which showed a 3.6-fold increase as compared to the basal ABA level (Figure 3a). After the recovery period, M18 plants showed ABA levels higher than those observed at the end of the stress, while ABA content decreased significantly in M23 and M28 plants. It is interesting to note that basal ABA levels were not recovered in M23 plants, while in M28 plants they were lower at tR than at t0. In contrast, levels of indole-3-acetic acid (IAA) (Figure 3b) did not vary throughout the heat stress experiment in M18 (p = 0.097) nor in M28 plants (p = 0.307), while the content of this hormone significantly increased in M23 plants (p = 0.004) at the end of the stress (t10).
Analyses of maritime pine needles determined 22 cytokinins (Supplementary Table S1) that varied among groups of plants derived from somatic embryos matured at different temperatures, and also throughout the heat stress experiment, as demonstrated by the PCA analysis (Figure 4). Basal levels of reversible CK metabolites, O-glucosides of trans-/cis-zeatin (tZ/cZ) and trans-/cis-zeatin ribosides, were higher in M28 plants than in M18 and M23 plants. Heat stress induced the accumulation of CKs ribosides in M23 and M28 plants, while needles of M18 plants accumulated dihydro-zeatin and the irreversible metabolite dihydrozeatin-9-glucoside. After recovery, M23 plants showed O-glucosides levels similar to those observed in M28 plants before the stress treatment, while M18 plants showed a different pattern of CKs.
In M18 plants, the total level of cZ- and tZ-type CKs were also significantly reduced after the stress treatment (t10), and raised again after the recovery period, while the total content of DZ types decreased significantly at this time (Table 3). These results contrasted with those observed in M23 plants, in which the total content of DZ-type CKs was significantly lower at t0 than in M18 and in M28 plants, and increased significantly in tR samples, as occurred with levels of cZ- and tZ-type CKs. In M28 plants, basal levels of cZ- and tZ-type cytokinins were significantly higher than in M18 and M23 plants (6.81 ± 1.14 pmol/g FW front to 2.65 ± 0.16 and 2.26 ± 0.06 pmol/g FW, p = 0.001, and 3.73 ± 0.62 pmol/g FW front to 1.45 ± 0.08 and 1.38 ± 0.05 pmol/g FW, p = 0.002; respectively). Both groups of compounds were significantly reduced by the heat stress treatment (t10), while the contents were partially recovered in tR samples (Table 3).
Table 4 summarizes the effect of the treatments on CKs grouped by biological form. Active CK levels were, on average, significantly higher in M18 than in M23 plants (0.17 ± 0.01 front to 0.14 ± 0.03 pmol/g FW, p = 0.040), although they did not differ significantly from the mean content (0.16 ± 0.06 pmol/g FW) determined in M28 plants. Furthermore, this higher level of CK bases in needles of M18 plants did not change during the heat stress experiment (p = 0.479), while it increased significantly after the recovery period in M23 plants (p = 0.007) and dropped drastically in M28 plants (p = 0.002). Major CK forms in the three groups of plants were O-glucosides (up to 96% in M28 plants after recovery); therefore, this group of metabolites drove the pattern of variation of the total CKs. In plants derived from SE performed in standard conditions (M23), heat stress slightly reduced the total CK contents, which were significantly increased after the recovery period (p = 0.002). In contrast, in M18 and M28 plants, the total CK levels decreased significantly after heat stress (p = 0.002) and these values were not fully restored after the recovery period.

3. Discussion

Variation in temperature conditions during the maturation step of the maritime pine SE-mediated propagation protocol resulted in modifications in the phenotype of the produced plants in terms of height, needle anatomy, and physiological and hormonal profiles, before and after a short-time heat stress treatment.
The most adequate conditions for maritime pine embryo maturation (23 °C) reported previously [6] produced higher plants after 2 years growing in a greenhouse than those derived from embryos matured at lower (18 °C) or higher temperature (28 °C); although 28 °C had been found to be the optimal temperature for the induction and proliferation stages [6]. The thicker epidermis found in M18 plants could be associated with a better adaptation to abiotic stress, as described in other species [37,38]. Exposure to high temperatures might produce foliar damage that is visible, depending on the species, the duration of exposure, and the time of year [39]. In our experiment, maritime pine plants were able to tolerate 45 °C for 3 h/day for 10 days, without any sign of visible damage at first sight. In contrast, a 30 min exposure to 48 °C produced 50% foliar damage in seedlings from Picea glauca [39]. Regarding anatomy, in our study, heat stress decreased chlorenchyma in M28 whereas it increased in M18 especially, which might indicate a better adaptation of these leaves to high temperatures. Similarly, palisade mesophyll was increased by heat in brassica heat-tolerant WS-1 but not in the sensitive WS-6 genotype [40]; in addition, a correlation between mesophyll cells damaged by heat stress and tolerance has been reported in salvia [41] and rhododendron [26].
Heat stress significantly affects the photosynthetic machinery and the chlorophyll and pigments biosynthetic pathways, thus disturbing the overall morpho-physiology of higher plants [42]. Heat sensitivity is intrinsically linked to drought avoidance; plants can either cool through transpiration, thus risking drought stress, or speedily close stomata to avoid drought stress while risking heat damage [43,44]. Plant responses to abnormal or extreme temperature changes are primarily mediated by hormones that might activate ROS production [45]. Heat tolerance is acquired when, in response to the excessive production of ROS, plants initiate heat stress responsive pathways for ROS-scavenging, such as the accumulation of proline and other osmolytes. ABA, the hormone regulating stomatal aperture under heat stress, has been suggested to be responsible for the initiation of other physiological responses, such as the accumulation of proline in stressed plants [24,27]. Proline has diverse roles, such as stabilizing proteins, membranes, and subcellular structures, protecting cellular functions by scavenging ROS, and acting as a signaling molecule [23]. Accumulation of such solutes may enhance stress tolerance via maintaining the cellular redox potential by lowering the generation of ROS at values suitable for metabolism when stomata are closed during osmotic stress, which limits carbon uptake and subsequently decreases NADPH consumption by the Calvin cycle [23]. This could explain why, in our abiotic stress experiment, neither RWC (Supplementary Figure S1) nor the photosynthesis-related parameters (Table 2) were significantly affected, since we only observed a reduction in stomatal conductance in needles from M28 plants sampled after the recovery period, in which a significant reduction in proline content was also detected. Needles sampled in M18 plants at t0 contained significantly lower levels of proline than M23 and M28 plants. Moreover, basal photosynthesis rates were significantly lower in M18 plants than in M23 plants and this parameter was significantly correlated with stomatal conductance (ρ = 0.493, p = 0.038), and inversely correlated with substomatal CO2 concentration (ρ = −0.474, p = 0.047). Heat stress treatment significantly increased the proline contents in M18 plants, and these higher contents were maintained after the recovery period, which might explain the increased net photosynthesis rate at tR in these M18 plants. This protective effect of proline was also suggested by [46] to explain the stability of the relative electron transfer rate under heat stress in seedlings of Pinus halepensis. Proline was also thought to mediate in alleviating the adverse effects of heat stress on photosynthesis in wheat [47] and pepper [48], as well as in the antioxidant response under heat conditions of Quercus suber [49] and Olea europaea [50]. De Diego et al. [51] also reported high levels of proline in radiata pine after recovery of a drought stress treatment, suggesting the possible implication of this metabolite in improving tolerance against future stress situations. In a previous work, we demonstrated that heat priming (50 °C) during SE induction improved performance under heat stress in Pinus pinaster, since primed plants showed better osmotic adjustment and higher increases in CKs, chlorophyll, soluble sugars, and starch contents, and their photosynthesis rates were less affected [18].
The hormone profile of SE-derived maritime pine plants was affected by variation in temperature during embryo maturation. Thus, needles of M28 plants showed significantly higher contents of ABA and total CK, mainly of CK bases and glucosides, and of cZ- and tZ-types, while M18 plants accumulated mainly iP-types, and both M18 and M28 plants showed higher contents of DZ-types. Higher basal levels of ABA, IAA, and CKs were also observed in primed maritime pine plants studied in our previous work, in which short temperature treatments at 37 or 50 °C were applied during the induction phase of the SE propagation protocol [18]. Variation in the hormone profiles of plants derived from somatic embryos matured at different temperatures resulted in significant coefficients of correlation between basal levels of ABA and total CKs (ρ = 0.804, p = 0.002); in addition, the total CK contents correlated with the content of active CK bases (ρ = 0.719, p = 0.001). Hormone profiles were associated with variation in needles histology, since the thicker epidermis observed in M18 plants resulted in a significant negative coefficient of correlation between epidermis thickness and IAA content (ρ = −0.729, p = 0.026), and the thicker chlorenchyma observed in M28 plants resulted in significant coefficients of correlation between this parameter and ABA and the total contents of CKs (ρ = 0.683, p = 0.020 and ρ = 0.613, p = 0.026, respectively).
As mentioned, M28 maritime pine plants, which derived from somatic embryos matured at a warmer temperature than that of the standard protocol (28 front to 23 °C), accumulated higher contents of ABA (Figure 3a) but did not increase proline levels after heat stress, while M18 plants, obtained under colder conditions, showed a high increase in the concentration of this osmolyte in response to heat stress (Figure 2). Therefore, the prolonged exposure (12 weeks) of embryogenic lines at the maturation step to higher or lower temperatures also altered the performance under heat stress of the regenerated maritime pine plants. Heat or cold priming-induced cross-tolerance is very common in plants, and often results from the synergistic co-activation of multiple stress signaling pathways, involving plant hormones, as well as ROS and transcription factors [52]. For instance, heat priming during SE induced drought tolerance in Pinus radiata [32] and drought priming induced heat tolerance in Festuca arundinacea [53]. Our results also suggest that SE maturation at low temperatures produces plants with increased resilience to further heat stress.
Cytokinins have also a positive role in plant response to high temperature. Earlier studies found that the CK contents are generally reduced by heat stress, but recent studies indicate that initially active CK contents increase and are rapidly depleted, suggesting that CKs could serve as primary modulators in temperature-sensing responses. In this sense, heat acclimation and inhibition of CK degradation positively affected heat stress tolerance of Arabidopsis [54]. In our heat stress experiment, proline and the contents of CK bases were positively correlated in t10 samples (ρ = 0.703, p = 0.035), but inversely correlated with total CK contents (ρ = −0.850, p = 0.004). After the recovery period, higher proline contents were also significantly associated with higher levels of CK bases and ABA (ρ = 0.833, p = 0.005 and ρ = 0.917, p = 0.001, respectively), suggesting heat stress adaptation.
Needles sampled after heat stress in M18 and M28 plants showed a significant reduction in total CK levels, although basal contents were regained after the recovery period (Table 4). In contrast, heat treatment did not reduce total CK content in needles from M23 plants, while these rates were significantly increased after the recovery time.
Taken together, our results show the main role of proline, ABA, and cytokinin bases in the heat-stress response of maritime pine plants. In Pinus radiata, heat-induced thermotolerance is mediated by ABA and SA at the first stage [34], probably due to the plant’s urgent need to regulate stomatal closure and counteract the increase in oxidative membrane damage, while proline, total sugars, IAA, and CKs seem to be more relevant in longer exposures and recovery.
In conclusion, temperature conditions during somatic embryo maturation affected the regenerated maritime pine (Pinus pinaster) plants at anatomical and physiological levels, resulting in altered performance under heat stress conditions. Our results suggest that those plants derived from the lower temperature (18 °C) treatment were more adapted to subsequent high temperature stress, since these plants showed faster and higher proline accumulation, increased photosynthesis after recovery, a lower increase in ABA levels, and no reduction in active CK contents. Further research is required to study other counterplayers involved in this phenotype and whether it is maintained for a long time after transferring plants to field conditions. In the meantime, we state that in vitro propagation techniques might allow the generation of plants with better adaptation to extreme conditions to broaden the pathway for obtaining maritime pine plants that could be more able to cope with the global warming challenge.

4. Materials and Methods

4.1. Plant Material

Embryogenic lines were derived from three OP mother trees (1007, 1046, and 1058) belonging to the Galician Tree Breeding Program (Conselleria do Medio Rural, Xunta de Galicia, Spain) [6]. Maritime pine plants were obtained by an SE-mediated propagation protocol described by Cano et al. [16], except that the induction and proliferation phases were performed at 28 °C, and the maturation step (12 weeks) was carried out at 18, 23, or 28 °C as described by Arrillaga et al. [6]. Plants derived from somatic embryos developed in these different conditions (M18, M23, and M28 plants) were acclimated to ex vitro conditions and cultivated in 0.5 L pots in the greenhouse (SCSIE, Central Service for Experimental Research, University of Valencia, Spain) for 2 years with a 16 h photoperiod and 200–300 W·m−2 irradiance. The plants were watered and fertilized weekly. Twelve plants from each group were selected, randomly distributed, and used for the heat stress experiment. Plants were watered at field capacity and subjected to an increasing linear gradient (±0.1 °C/min) of air temperature that was set from 22 to 45 °C for 4 h, using an aerothermal heather, then maintained for 3 h at this temperature before decreasing to the initial conditions (22 ± 3 °C). This gradient was repeated for 10 days. Mature needles from the maritime pine plants were sampled before the stress experiment (t0), at the end of the stress experiment (t10), and after a further recovery period of 10 days (tR). Samples were frozen in liquid nitrogen and stored at −80 °C until analysis.

4.2. Osmotic Adjustment Determination

Relative water content (RWC) was determined in the manner described by Escandón et al. [35], using three needle fragments (1 cm long) in each of six replicates per group of plants and sampling point. Samples’ fresh weight (FW) was registered and the needles were maintained with deionized water for 24 h in the dark at 4 °C, after which the turgid weight (TW) was recorded. Then, needles were dried at 80 °C for 72 h, and dry weight (DW) was registered. RWC was calculated by using the following equation: RWC (%) = (FW − DW)/(TW − DW) × 100.
Proline was quantified in needles as described by Ábrahám et al. [55], with some modifications, in six replicates per treatment. Approximately 150–300 mg of frozen needles were homogenized in 3% sulfosalicylic acid (5 μL/mg FW), and centrifuged at 13,000 rpm for 5 min. A mixture with 100 μL of 3% sulfosalicylic acid, 200 μL of glacial acetic acid, and 200 μL of acidic ninhydrin was added to 100 μL of the supernatant of the extract, and the resulting mixture was vortexed and incubated at 96 °C for 1 h. After that, the reaction was finished on ice for 10 min. Samples were extracted in 0.5 mL of toluene and vortexed for 20 s, and the formation of the two phases was observed. The absorbance of the chromophore-containing toluene phase was read at 520 nm (Eppendorf BioSpectrometer® basic, Hamburg, Germany) using toluene as the blank reagent, and proline concentration was determined from an L-proline (Sigma-Aldrich®, Merck KGaA, Darmstadt, Germany) standard curve with six points (0–150 μg/mL).

4.3. Photosynthesis-related Parameters

Photosynthesis-related parameters were determined in needles of M18, M23, and M28 maritime pine plants sampled before (t0) and after (t10) the heat stress treatment, as well as after a recovery period of a further 10 days (tR). Determinations of net photosynthetic CO2 assimilation rates (μmol CO2 m−2 s−1), stomatal conductance (mol H2O m−2 s−1), and substomatal CO2 concentration (μmol CO2 mol−1) were performed in the morning using six plants from each group and time, by using a portable photosynthesis system LICOR 6400 (LI-COR, Lincoln, NE, USA). The system was combined with a chamber for conifers (LI 6400-05). The conditions in the measuring chamber were controlled at a flow rate of 500 mol s−1, a saturating PAR of 1200 μmol m−2 s−1, 400 ppm CO2, and 60–70% relative humidity. We estimated the percentage of variation for each parameter at t10 and tR, as compared to t0 values.

4.4. Hormone Content Determinations

The endogenous contents of ABA, IAA, and CKs were determined in Pinus pinaster needles sampled from 2-year-old M18, M23, and M28 plants, before (t0) and after (t10) the heat treatment, and after the recovery period (tR). The extraction, purification, and quantification of CK metabolites were performed as described by Svačinová et al. [56] using multi-StageTip technology based on C18, SDB-RPS, and Cation-SR sorbents (Affinisep, AttractSPETM, Petit Couronne, France). For each group of plants and time, two samples were extracted in 1 mL of modified Bieleski solution containing stable isotope-labelled internal standards (0.2 pmol of CK bases, ribosides, and N9- and N7-glucosides; 0.5 pmol of CK O-glucosides and nucleotides). After extraction, from each sample, three technical replicates of 300 μL were purified using microSPE columns and then eluates were evaporated to dryness.
The concentration levels of IAA and ABA were determined according to the modified method described by Šimura et al. [57]. Briefly, samples (four replicates from each group of plants and time) containing 10 mg FW were extracted in aqueous solution of 50% acetonitrile (v/v). Crude extracts were loaded onto conditioned Oasis HLB columns (30 mg/1 mL, Waters) and washed with an aqueous solution of 30% acetonitrile (v/v). Flow-through fractions containing purified analytes were collected and evaporated to dryness in a vacuum evaporator.
All samples were dissolved in 30 μL of mobile phase and then analyzed using an Acquity I-Class system (Waters, Milford, MA, USA) combined with a triple quadrupole mass spectrometer (XevoTMTQ-S MS, Waters, MS Technologies, Manchester, UK). A mixture of stable isotope-labelled standards of hormones was added to validate the LC-MS/MS method, and concentration levels were calculated using the isotope dilution method. All data were processed with MassLynx V4.2 software (Waters, Milford, MA, USA).

4.5. Histological Analyses

Maritime pine needles sampled before and after the heat stress treatment (24 replicates for each group of plants), were fixed in Karnovsky solution and stained with toluidine blue. Transversal sections were photographed under a microscope, and we determined the mean thickness of the epidermis (8 points × 10 needles per group of plants), and of the chlorenchyma (2 measures × 24 needles), using ImageJ software.

4.6. Statistical Analyses

Data recorded in the different experiments were subjected to analysis of variance using SPSS software (IBM Statistics, Armonk, NY, USA). When they did not adjust to a normal distribution (Kolmogorov–Smirnoff test), significant differences were assessed using the Kruskal–Wallis analysis of variance test. Correlation among parameters was also estimated by the Spearman ρ coefficient [58].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23031318/s1.

Author Contributions

Conceptualization, I.A., J.S., and E.S.; methodology, all authors; formal analysis, E.S., E.C., S.G.N., I.P., O.N., C.P., and M.A.P.-O.; investigation, all authors.; writing—review and editing, all authors; supervision, E.S., I.A.; funding acquisition, I.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MICINN and by the European Regional Development Fund Project (AGL2016-76143-C4-01-R). This work was also supported by the Czech Science Foundation (20-22875S) and the Ministry of Education, Youth and Sports of the Czech Republic, from the European Regional Development Fund (CZ.02.1.01/0.0/0.0/16_019/0000827).

Acknowledgments

We thank Alex Alborch and María Cano for their technical assistance in tissue culture and greenhouse maintenance; Hana Martínková and Ivan Petřík for their technical assistance with plant hormone profiling; and the SCSIE (Central Service for Experimental Research, University of Valencia) for greenhouse maintenance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Zas, R.; Sampedro, L.; Solla, A.; Vivas, M.; Lombardero, M.J.; Alía, R.; Rozas, V. Dendroecology in common gardens: Population differentiation and plasticity in resistance, recovery and resilience to extreme drought events in Pinus pinaster. Agric. For. Meteorol. 2020, 291, 108060. [Google Scholar] [CrossRef]
  2. Alía, R.; Martín, S. EUFORGEN Technical Guidelines for Genetic Conservation and Use for Maritime Pine (Pinus pinaster); International Plant Genetic Resources Institute; Bioversity International: Rome, Italy, 2003; p. 6. [Google Scholar]
  3. Nguyen-Queyrens, A.; Bouchet-Lannat, F. Osmotic adjustment in three-year-old seedlings of five provenances of maritime pine (Pinus pinaster) in response to drought. Tree Physiol. 2003, 23, 397–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Menéndez-Gutiérrez, M.; Alonso, M.; Toval, G.; Díaz, R. Variation in pinewood nematode susceptibility among Pinus pinaster Ait. provenances from the Iberian Peninsula and France. Ann. For. Sci. 2017, 74, 76. [Google Scholar] [CrossRef] [Green Version]
  5. Arrillaga, I.; Guevara, M.A.; Muñoz-Bertomeu, J.; Lázaro-Gimeno, D.; Sáez-Laguna, E.; Díaz, L.M.; Torralba, L.; Mendoza-Poudereux, I.; Segura, J.; Cervera, M.T. Selection of haploid cell lines from megagametophyte cultures of maritime pine as a DNA source for massive sequencing of the species. Plant Cell Tiss. Org. Cult. 2014, 118, 147–155. [Google Scholar] [CrossRef]
  6. Arrillaga, I.; Morcillo, M.; Zanón, I.; Lario, F.; Segura, J.; Sales, E. New approaches to optimize somatic embryogenesis in maritime pine. Front. Plant Sci. 2019, 10, 138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Andivia, E.; Ruiz-Benito, P.; Díaz-Martínez, P.; Carro-Martínez, N.; Zavala, M.A.; Madrigal-González, J. Inter-specific tolerance to recurrent droughts of pine species revealed in saplings rather than adult trees. For. Ecol. Manag. 2020, 459, 117858. [Google Scholar] [CrossRef]
  8. Neumann, M.; Mues, V.; Moreno, A.; Hasenauer, H.; Seidl, R. Climate variability drives recent tree mortality in Europe. Glob. Chang. Biol. 2017, 23, 4788–4797. [Google Scholar] [CrossRef]
  9. Lebedev, V.G.; Lebedeva, T.N.; Chernodubov, A.I.; Shestibratov, K.A. Genomic selection for forest tree improvement: Methods, achievements and perspectives. Forests 2020, 11, 1190. [Google Scholar] [CrossRef]
  10. Lelu-Walter, M.A.; Klimaszewska, K.; Miguel, C.; Aronen, T.; Hargreaves, C.; Teyssier, C.; Trontin, J.F. Somatic embryogenesis for more effective breeding and deployment of improved varieties in Pinus spp.: Bottlenecks and recent advances. In Somatic Embryogenesis: Fundamental Aspects and Applications; Loyola-Vargas, V.M., Ochoa-Alejo, N., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 319–365. [Google Scholar]
  11. Lelu-Walter, M.A.; Thompson, D.; Harvengt, L.; Sanchez, L.; Toribio, M.; Pâques, L.E. Somatic embryogenesis in forestry with a focus on Europe: State-of-the-art, benefits, challenges and future direction. Tree Genet. Genomes 2013, 9, 883–899. [Google Scholar] [CrossRef]
  12. Turgut-Kara, N.; Arikan, B.; Celik, H. Epigenetic memory and priming in plants. Genetica 2020, 148, 47–54. [Google Scholar] [CrossRef]
  13. Galviz, Y.C.F.; Ribeiro, R.V.; Souza, G.M. Yes, plants do have memory. Theor. Exp. Plant Physiol. 2020, 32, 195–202. [Google Scholar] [CrossRef]
  14. Kvaalen, H.; Johnsen, Ø. Timing of bud set in Picea abies is regulated by a memory of temperature during zygotic and somatic embryogenesis. New Phytol. 2008, 177, 49–59. [Google Scholar] [CrossRef]
  15. Egertsdotter, U. Plant physiological and genetical aspects of the somatic embryogenesis process in conifers. Scand. J. For. Res. 2019, 34, 360–369. [Google Scholar] [CrossRef] [Green Version]
  16. Cano, M.; Morcillo, A.; Humánez, A.; Mendoza-Poudereux, I.; Alborch, A.; Segura, J.; Arrillaga, I. Maritime pine (Pinus pinaster Aiton). In Step Wise Protocols for Somatic Embryogenesis of Important Woody Plants, 2nd ed.; Jain, S.M., Gupta, P., Eds.; Springer International Publishing: Cham, Switzerland, 2018; Volume 1, pp. 167–179. [Google Scholar]
  17. Humánez, A.; Blasco, M.; Brisa, C.; Segura, J.; Arrillaga, I. Somatic embryogenesis from different tissues of Spanish populations of maritime pine. Plant Cell Tissue Organ Cult. 2012, 111, 373–383. [Google Scholar] [CrossRef]
  18. Pérez-Oliver, M.A.; Haro, J.G.; Pavlović, I.; Novák, O.; Segura, J.; Sales, E.; Arrillaga, I. Priming maritime pine megagametophytes during somatic embryogenesis improved plant adaptation to heat stress. Plants 2021, 10, 446. [Google Scholar] [CrossRef]
  19. García-Mendiguren, O.; Montalbán, I.A.; Goicoa, T.; Ugarte, M.D.; Moncaleán, P. Environmental conditions at the initial stages of Pinus radiata somatic embryogenesis affect the production of somatic embryos. Trees 2016, 30, 949–958. [Google Scholar] [CrossRef]
  20. Marques do Nascimento, A.M.; Alves-Barroso, P.; Ferreira do Nascimento, N.F.; Goicoa, T.; Ugarte, M.D.; Montalbán, I.A.; Moncaleán, P. Pinus spp. somatic embryo conversion under high temperature: Effect on the morphological and physiological characteristics of plantlets. Forests 2020, 11, 1181. [Google Scholar] [CrossRef]
  21. Chen, W.L.; Yang, W.J.; Lo, H.F.; Yeh, D.M. Physiology, anatomy, and cell membrane thermostability selection ofleafy radish (Raphanus sativus var. Oleiformis Pers.) with different tolerance under heat stress. Sci. Hortic. 2014, 179, 367–375. [Google Scholar] [CrossRef]
  22. Estravis-Barcala, M.; Mattera, M.G.; Soliani, C.; Bellora, N.; Opgenoorth, L.; Heer, K.; Arana, M.V. Molecular bases of responses to abiotic stress in trees. J. Exp. Bot. 2020, 71, 3765–3779. [Google Scholar] [CrossRef]
  23. Kaur, G.; Asthir, B. Proline: A key player in plant abiotic stress tolerance. Biol. Plant. 2015, 59, 609–619. [Google Scholar] [CrossRef]
  24. Ghosh, U.K.; Islam, M.N.; Siddiqui, M.N.; Cao, X.; Khan, M.A.R. Proline, a multifaceted signalling molecule in plant responses to abiotic stress: Understanding the physiological mechanisms. Plant Biol. 2021. [Google Scholar] [CrossRef]
  25. Yu, D.; Wildhagen, H.; Tylewicz, S.; Miskolczi, P.C.; Bhalerao, R.P.; Polle, A. Abscisic acid signalling mediates biomass trade-off and allocation in poplar. New Phytol. 2019, 223, 1192–1203. [Google Scholar] [CrossRef]
  26. Shen, H.F.; Zhao, B.; Xu, J.J.; Liang, W.; Huang, W.M.; Li, H.H. Effects of heat stress on changes in physiology and anatomy in two cultivars of Rhododendron. S. Afr. J. Bot. 2017, 112, 338–345. [Google Scholar] [CrossRef]
  27. Li, M.; Jannasch, A.H.; Jiang, Y. Growth and Hormone Alterations in Response to Heat Stress in Perennial Ryegrass Accessions Differing in Heat Tolerance. J. Plant Growth Regul. 2020, 39, 1022–1102. [Google Scholar] [CrossRef]
  28. Pereira, C.; Castander-Olarieta, A.; Sales, E.; Montalbán, I.A.; Canhoto, J.; Moncaleán, P. Heat stress in Pinus halepensis somatic embryogenesis induction: Effect in DNA methylation and differential expression of stress-related genes. Plants 2021, 10, 2333. [Google Scholar] [CrossRef]
  29. Castander-Olarieta, A.; Pereira, C.; Sales, E.; Meijón, M.; Arrillaga, I.; Cañal, M.J.; Goicoa, T.; Ugarte, M.D.; Moncaleán, P.; Montalbán, I.A. Induction of radiata pine somatic embryogenesis at high temperatures provokes a long-term decrease in DNA methylation/hydroxymethylation and differential expression of stress-related genes. Plants 2020, 9, 1762. [Google Scholar] [CrossRef]
  30. Pereira, C.; Castander-Olarieta, A.; Montalbán, I.A.; Pěnčík, A.; Petřík, I.; Pavlovic, I.; De Medeiros Oliveira, E.; Fraga, H.; Guerra, M.P.; Novak, O.; et al. Embryonal masses induced at high temperatures in Aleppo pine: Cytokinin profile and cytological characterization. Forests 2020, 11, 807. [Google Scholar] [CrossRef]
  31. Castander-Olarieta, A.; Pereira, C.; Montalbán, I.A.; Pěnčík, A.; Petřík, I.; Pavlović, I.; Novák, O.; Strnad, M.; Moncaleán, P. Quantification of endogenous aromatic cytokinins in Pinus radiata embryonal masses after application of heat stress during initiation of somatic embryogenesis. Trees 2020, 35, 1075–1080. [Google Scholar] [CrossRef]
  32. Castander-Olarieta, A.; Moncaleán, P.; Pereira, C.; Pěnčík, A.; Petřík, I.; Pavlovic, I.; Novák, O.; Strnad, M.; Goicoa, T.; Ugarte, M.D.; et al. Cytokinins are involved in drought tolerance of Pinus radiata plants originating from embryonal masses induced at high temperatures. Tree Physiol. 2021, 41, 912–926. [Google Scholar] [CrossRef]
  33. Moncaleán, P.; García-Mendiguren, O.; Novak, O.; Strnad, M.; Goicoa, T.; Ugarte, M.D.; Montalbán, I.A. Temperature and water availability during maturation affect the cytokinins and auxins profile of radiata pine somatic embryos. Front. Plant Sci. 2018, 9, 1898. [Google Scholar] [CrossRef] [Green Version]
  34. Escandón, M.; Cañal, M.; Pascual, J.; Pinto, G.; Correia, B.; Amaral, J.; Meijón, M. Integrated physiological and hormonal profile of heat-induced thermotolerance in Pinus radiata. Tree Physiol. 2016, 36, 63–77. [Google Scholar] [CrossRef] [PubMed]
  35. Escandón, M.; Valledor, L.; Pascual, J.; Pinto, G.; Cañal, M.J.; Meijón, M. System-wide analysis of short-term response to high temperature in Pinus radiata. J. Exp. Bot. 2017, 68, 3629–3641. [Google Scholar] [CrossRef] [PubMed]
  36. Kolb, P.F.; Robberecht, R. High temperature and drought stress effects on survival of Pinus ponderosa seedlings. Tree Physiol. 1996, 16, 665–672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Turunen, M.; Latola, K. UV-B radiation and acclimation in timberline plants. Environ. Pollut. 2005, 137, 390–403. [Google Scholar] [CrossRef]
  38. Li, W.D.; Biswas, D.K.; Xu, H.; Xu, C.Q.; Wang, X.Z.; Liu, J.K.; Jiang, G.M. Photosynthetic responses to chromosome doubling in relation to leaf anatomy in Lonicera japonica subjected to water stress. Funct. Plant. Biol. 2009, 36, 783–792. [Google Scholar] [CrossRef]
  39. Teskey, R.; Wertin, T.; Bauweraerts, I.; Ameye, M.; McGuire, M.; Steppe, K. Review: Responses of tree species to heat waves and extreme heat events. Plant Cell Environ. 2015, 38, 1699–1712. [Google Scholar] [CrossRef]
  40. Yuan, L.; Tang, L.; Zhu, S.; Hou, J.; Chen, G.; Liu, F.; Liu, S.; Wang, C. Influence of heat stress on leaf morphology and nitrogen–carbohydrate metabolisms in two wucai (Brassica campestris L.) genotypes. Acta Soc. Bot. Pol. 2017, 86, 3554. [Google Scholar] [CrossRef]
  41. Natarajan, S.; Kuehny, J.S. Morphological, Physiological, and Anatomical Characteristics Associated with Heat Preconditioning and Heat Tolerance in Salvia splendens. J. Amer. Soc. Hort. Sci. 2008, 133, 527–534. [Google Scholar] [CrossRef] [Green Version]
  42. Muhammad, I.; Shalmani, A.; Ali, M.; Yang, Q.H.; Ahmad, H.; Li, F.B. Mechanisms regulating the dynamics of photosynthesis under abiotic stresses. Front. Plant Sci. 2021, 11, 2310. [Google Scholar] [CrossRef]
  43. Konôpková, A.; Kurjak, D.; Kmeť, J.; Klumpp, R.; Longauer, R.; Ditmarová, Ľ.; Gömöry, D. Differences in photochemistry and response to heat stress between silver fir (Abies alba Mill.) provenances. Trees 2018, 32, 73–86. [Google Scholar] [CrossRef]
  44. Kunert, N.; Hajek, P.; Hietz, P.; Morris, H.; Rosner, S.; Tholen, D. Summer temperatures reach the thermal tolerance threshold of photosynthetic decline in temperate conifers. Plant Biol. 2021. [Google Scholar] [CrossRef] [PubMed]
  45. Devireddy, A.R.; Tschaplinski, T.J.; Tuskan, G.A.; Muchero, W.; Chen, J.G. Role of reactive oxygen species and hormones in plant responses to temperature changes. Int. J. Mol. Sci. 2021, 22, 8843. [Google Scholar] [CrossRef] [PubMed]
  46. Birami, B.; Gattmann, M.; Heyer, A.G.; Grote, R.; Arneth, A.; Ruehr, N.K. Heat waves alter carbon allocation and increase mortality of Aleppo pine under dry conditions. Front. For. Glob. Chang. 2018, 1, 8. [Google Scholar] [CrossRef]
  47. Khan, M.I.R.; Iqbal, N.; Masood, A.; Per, T.S.; Khan, N.A. Salicylic acid alleviates adverse effects of heat stress on photosynthesis through changes in proline production and ethylene formation. Plant Signal. Behav. 2013, 8, e26374. [Google Scholar] [CrossRef] [Green Version]
  48. Rajametov, S.N.; Yang, E.Y.; Cho, M.C.; Chae, S.Y.; Jeong, H.B.; Chae, W.B. Heat-tolerant hot pepper exhibits constant photosynthesis via increased transpiration rate, high proline content and fast recovery in heat stress condition. Sci. Rep. 2021, 11, 14328. [Google Scholar] [CrossRef]
  49. Correia, B.; Rodriguez, J.L.; Valledor, L.; Almeida, T.; Santos, C.; Cañal, M.J.; Pinto, G. Analysis of the expression of putative heat-stress related genes in relation to thermotolerance of cork oak. J. Plant Physiol. 2014, 171, 399–406. [Google Scholar] [CrossRef]
  50. Dias, M.C.; Santos, C.; Silva, S.; Pinto, D.C.G.A.; Silva, A.M.S. Physiological and Metabolite Reconfiguration of Olea europaea to cope and recover from a heat or high UV-B shock. J. Agric. Food Chem. 2020, 68, 11339–11349. [Google Scholar] [CrossRef]
  51. De Diego, N.; Saiz-Fernández, I.; Rodríguez, J.L.; Pérez-Alfocea, P.; Sampedro, M.C.; Barrio, R.J.; Lacuesta, M.; Moncaleán, P. Metabolites and hormones are involved in the intraspecific variability of drought hardening in radiata pine. J. Plant Physiol. 2015, 188, 64–71. [Google Scholar] [CrossRef]
  52. Hossain, M.A.; Li, Z.; Hoque, T.S.; Burritt, D.J.; Fujita, M.; Munné-Bosch, S. Heat or cold priming-induced cross-tolerance to abiotic stresses in plants: Key regulators and posible mechanisms. Protoplasma 2018, 255, 399–412. [Google Scholar] [CrossRef]
  53. Zhang, X.; Wang, X.; Zhuang, L.; Gao, Y.; Huang, B. Abscisic acid mediation of drought priming-enhanced heat tolerance in tall fescue (Festuca arundinacea) and Arabidopsis. Physiol. Plant. 2019, 167, 488–501. [Google Scholar] [CrossRef]
  54. Prerostova, S.; Dobrev, P.I.; Kramna, B.; Gaudinova, A.; Knirsch, V.; Spichal, L.; Zatloukal, M.; Vankova, R. Heat Acclimation and Inhibition of Cytokinin Degradation Positively Affect Heat Stress Tolerance of Arabidopsis. Front. Plant Sci. 2020, 11, 87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ábrahám, E.; Hourton-cabassa, C.; Erdei, L.; Szabados, L. Methods for determination of proline in plants. Methods Mol. Biol. 2010, 639, 317–331. [Google Scholar] [CrossRef] [PubMed]
  56. Svačinová, J.; Novák, O.; Plačková, L.; Lenobel, R.; Holík, J.; Strnad, M.; Doležal, K. A new approach for cytokinin isolation from Arabidopsis tissues using miniaturized purification: Pipette tip solid-phase extraction. Plant Methods 2012, 8, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Šimura, J.; Antoniadi, I.; Široká, J.; Tarkowská, D.; Strnad, M.; Ljung, K.; Novák, O. Plant hormonomics: Multiple phyto-hormone profiling by targeted metabolomics. Plant Physiol. 2018, 177, 476–489. [Google Scholar] [CrossRef] [Green Version]
  58. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017. [Google Scholar]
Figure 1. Transversal sections of needles from maritime pine plants derived from somatic embryos matured at 18 °C (a), 23 °C (b), or 28 °C (c) before (a, b, and c, respectively) or after (d, e and f, respectively) a heat treatment for 10 days at 45 °C for 3 h/day. Note epidermis and chlorenchymatic tissues. Bar, 50 µm.
Figure 1. Transversal sections of needles from maritime pine plants derived from somatic embryos matured at 18 °C (a), 23 °C (b), or 28 °C (c) before (a, b, and c, respectively) or after (d, e and f, respectively) a heat treatment for 10 days at 45 °C for 3 h/day. Note epidermis and chlorenchymatic tissues. Bar, 50 µm.
Ijms 23 01318 g001
Figure 2. Proline content in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 10 days at 45 °C for 3 h/day. Plants were sampled at the beginning (t0) and at the end (t10) of the stress treatment, and after a recovery period of a further 10 days (tR). Data are the mean ± SD of three replicates. Within each group of plants, values followed by the same lowercase letter did not differ significantly according to the Tukey-b test, while capital letters indicate mean separation among groups of plants sampled at the same time.
Figure 2. Proline content in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 10 days at 45 °C for 3 h/day. Plants were sampled at the beginning (t0) and at the end (t10) of the stress treatment, and after a recovery period of a further 10 days (tR). Data are the mean ± SD of three replicates. Within each group of plants, values followed by the same lowercase letter did not differ significantly according to the Tukey-b test, while capital letters indicate mean separation among groups of plants sampled at the same time.
Ijms 23 01318 g002
Figure 3. Abscisic acid (a) and indole-acetic acid (b) contents in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 45 °C for 3 h/day, for 10 days. Plants were sampled at the beginning of the stress treatment (t0), at the end of the stress treatment (t10), and after a recovery period of a further 10 days (tR). Data are the mean ± SD of four replicates. Within each group of plants, values followed by the same lowercase letter did not differ significantly according to the Tukey-b test, while capital letters indicate mean separation among groups of plants sampled at the same time.
Figure 3. Abscisic acid (a) and indole-acetic acid (b) contents in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 45 °C for 3 h/day, for 10 days. Plants were sampled at the beginning of the stress treatment (t0), at the end of the stress treatment (t10), and after a recovery period of a further 10 days (tR). Data are the mean ± SD of four replicates. Within each group of plants, values followed by the same lowercase letter did not differ significantly according to the Tukey-b test, while capital letters indicate mean separation among groups of plants sampled at the same time.
Ijms 23 01318 g003
Figure 4. Results of the analysis when determined metabolites were grouped by CK types. Total iP-type CK content varied significantly throughout the heat stress experiment (p < 0.001), and this variation depended on the group of plants, as total iP levels decreased in M18 plants at the end of the stress (from 0.34 ± 0.04 to 0.22 ± 0.03 pmol/g FW), while it increased in M23 and M28 plants (from 0.26 ± 0.04 to 0.45 ± 0.06 pmol/g FW, and from 0.28 ± 0.06 to 0.50 ± 0.12 pmol/g FW, respectively). After the recovery period, total iP contents in M18 and M28 plants returned to levels similar to those at t0, while remaining unchanged in M23 plants.
Figure 4. Results of the analysis when determined metabolites were grouped by CK types. Total iP-type CK content varied significantly throughout the heat stress experiment (p < 0.001), and this variation depended on the group of plants, as total iP levels decreased in M18 plants at the end of the stress (from 0.34 ± 0.04 to 0.22 ± 0.03 pmol/g FW), while it increased in M23 and M28 plants (from 0.26 ± 0.04 to 0.45 ± 0.06 pmol/g FW, and from 0.28 ± 0.06 to 0.50 ± 0.12 pmol/g FW, respectively). After the recovery period, total iP contents in M18 and M28 plants returned to levels similar to those at t0, while remaining unchanged in M23 plants.
Ijms 23 01318 g004
Table 1. Morphology and histology of needles from 2-year-old maritime pine plants derived from somatic embryos matured at 18, 23, or 28 °C (M18, M23, and M28 plants, respectively) that were sampled before (t0) and after (t10) a heat stress treatment for 10 days at 45 °C for 3 h/day. Data are the mean ± SD of at least 10 needles; values followed by the same letter were not significantly different according to Tukey-b (x) or Tamhane (y) tests.
Table 1. Morphology and histology of needles from 2-year-old maritime pine plants derived from somatic embryos matured at 18, 23, or 28 °C (M18, M23, and M28 plants, respectively) that were sampled before (t0) and after (t10) a heat stress treatment for 10 days at 45 °C for 3 h/day. Data are the mean ± SD of at least 10 needles; values followed by the same letter were not significantly different according to Tukey-b (x) or Tamhane (y) tests.
PlantsNeedles Thickness x
(µm)
Epidermis Thickness y
(µm)
Chlorenchyma Thickness x
(µm)
t0t10t0t10t0t10
M18553 ± 30 b651 ± 32 a20.3 ± 0.9 a18.6 ± 2.3 ab143.0 ± 11.8 c184.8 ± 13.6 a
M23542 ± 42 b657 ± 57 a18.8 ± 1.4 a18.9 ± 1.5 a139.8 ± 19.0 c159.0 ± 18.3 b
M28635 ± 24 a557 ± 31 b19.1 ± 1.0 a16.6 ± 0.8 b166.8 ± 9.8 b136.4 ± 8.1 c
Table 2. Photosynthesis-related parameters in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C), before (t0) and after (t10) a heat stress treatment (for 10 days at 45 °C for 3 h/day), as well as after a recovery period of further 10 days (tR). Six samples were employed for determining net photosynthetic CO2 assimilation rates (AN; µmol CO2 m−2 s−1), stomatal conductance (gs; mol H2O m−2 s−1), and substomatal CO2 concentration (Ci; µmol CO2 mol−1).
Table 2. Photosynthesis-related parameters in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C), before (t0) and after (t10) a heat stress treatment (for 10 days at 45 °C for 3 h/day), as well as after a recovery period of further 10 days (tR). Six samples were employed for determining net photosynthetic CO2 assimilation rates (AN; µmol CO2 m−2 s−1), stomatal conductance (gs; mol H2O m−2 s−1), and substomatal CO2 concentration (Ci; µmol CO2 mol−1).
ParameterPlantst0t10tR
AN xM187.4 ± 2.4 b7.0 ± 4.5 (−5.2%) a9.0 ± 2.7 (22.4%) a
M2310.9 ± 3.8 a8.7 ± 2.9 (−20.5%) a10.6 ± 2.0 (−3.2%) a
M289.9 ± 2.0 a9.1 ± 0.7 (−8.4%) a9.5 ± 1.0 (−4.0%) a
gs yM180.24 ± 0.06 a0.17 ± 0.06 (−29.2%) a0.17 ± 0.12 (−29.2%) a
M230.24 ± 0.11 a0.16 ± 0.07 (−33.3%) a0.17 ± 0.09 (−29.2%) a
M280.25 ± 0.12 a0.19 ± 0.09 (−24.0%) ab0.13 ± 0.01 (−48.0%) b
M18313 ± 31 a294 ± 54 (−6.1%) ab264 ± 39 (−15.4%) b
(Ci) zM23282 ± 12 a270 ± 36 (−4.3%) a255 ± 35 (−9.6%) a
M28292 ± 21 a271 ± 52 (−7.2%) ab 235 ± 39 (−19.4%) b
x, mean separation within each sampling time according to Tukey-b test. y, mean separation within each group of plants according to Kruskal-Wallis ANOVA test with Bonferroni correction for multiple comparisons. z, mean separation within each group of plants according to the Tukey-b test.
Table 3. Cytokinin content in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 10 days at 45 °C. Plants were sampled at the beginning of the stress treatment (t0), at the end of the stress treatment (t10), and after a recovery period of a further 10 days (tR). Data are the mean ± SE of six replicates and 22 cytokinins grouped by type: cis-zeatin (cZ), trans-zeatin (tZ), dihydro-zeatin (DZ), and isopentenyladenine (iP). Within each group of plants, values followed by the same letter did not differ significantly (p < 0.05), according to Kruskal-Wallis one-way ANOVA.
Table 3. Cytokinin content in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 10 days at 45 °C. Plants were sampled at the beginning of the stress treatment (t0), at the end of the stress treatment (t10), and after a recovery period of a further 10 days (tR). Data are the mean ± SE of six replicates and 22 cytokinins grouped by type: cis-zeatin (cZ), trans-zeatin (tZ), dihydro-zeatin (DZ), and isopentenyladenine (iP). Within each group of plants, values followed by the same letter did not differ significantly (p < 0.05), according to Kruskal-Wallis one-way ANOVA.
PlantsTimeTotal iP-TypesTotal tZ-TypesTotal DZ-TypesTotal cZ-Types
M18t00.34 ± 0.04 a1.45 ± 0.08 a2.30 ± 0.07 a2.65 ± 0.16 a
t100.22 ± 0.03 b0.34 ± 0.05 b0.86 ± 0.12 ab0.87 ± 0.08 b
tR0.35 ± 0.04 a0.47 ± 0.05 ab0.62 ± 0.03 b1.87 ± 0.04 a
M23t00.26 ± 0.04 b1.38 ± 0.05 ab1.56 ± 0.09 ab2.26 ± 0.06 b
t100.45 ± 0.06 a0.95 ± 0.11 b0.87 ± 0.15 b2.67 ± 0.22 ab
tR0.45 ± 0.08 a2.72 ± 0.52 a4.12 ± 0.50 a4.22 ± 0.50 a
M28t00.28 ± 0.06 b3.73 ± 0.62 a2.44 ± 0.34 a6.81 ± 1.14 a
t100.50 ± 0.12 a0.82 ± 0.05 b0.64 ± 0.09 b2.90 ± 0.28 b
tR0.28 ± 0.03 b2.56 ± 0.14 a2.19 ± 0.16 a4.60 ± 0.32 ab
Table 4. Cytokinin content in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 10 days at 45 °C for 3 h/day. Plants were sampled at the beginning of the stress treatment (t0), at the end of the stress treatment (t10), and after a recovery period of a further 10 days (tR). Data are the mean ± SE of six replicates and 22 compounds grouped by conjugated structure, and the percentages refer to total CK content.
Table 4. Cytokinin content in needles of maritime pine plants derived from somatic embryos matured at different temperatures (18, 23, or 28 °C) that were employed in a heat stress experiment for 10 days at 45 °C for 3 h/day. Plants were sampled at the beginning of the stress treatment (t0), at the end of the stress treatment (t10), and after a recovery period of a further 10 days (tR). Data are the mean ± SE of six replicates and 22 compounds grouped by conjugated structure, and the percentages refer to total CK content.
PlantsTimeActive FormsTransport FormsPrecursorsReversible MetabolitesIrreversible MetabolitesTotal CK y
CKs Bases xCK Ribosides yCK Nucleotides yCK O-glucosides yCK N-glucosides y
M18t00.16 ± 0.03 (2.4%) a0.06 ± 0.01 (0.9%) b0.20 ± 0.02 (3.0%) a6.28 ± 0.17 (93.3%) a0.03 ± 0.00 (0.4%) a6.73 ± 0.18 a
t100.16 ± 0.03 (7.1%) a0.13 ± 0.02 (5.9%) ab0.09 ± 0.01 (3.0%) b1.89 ± 0.21 (82.6) b0.01 ± 0.00 (0.5%) ab2.28 ± 0.26 b
tR0.18 ± 0.03 (5.3%) a0.24 ± 0.03 (7.1%) a0.16 ± 0.01 (4.8%) ab2.74 ± 0.06 (82.6%) ab0.01 ± 0.00 (0.2%) b3.31 ± 0.11 ab
M23t00.12 ± 0.04 (2.2%) b0.05 ± 0.01 (0.9%) b0.20 ± 0.04 (3.6%) b5.07 ± 0.09 (92.8%) ab0.03 ± 0.02 (0.5%) ab5.46 ± 0.13 ab
t100.12 ± 0.01 (2.4%) b0.29 ± 0.03 (5.9%) a0.59 ± 0.09 (12.0%) a3.93 ± 0.37 (79.6%) b0.01 ± 0.00 (0.1%) b4.94 ± 0.47 b
tR0.17 ± 0.02 (1.5%) a0.09 ± 0.00 (0.8%) ab0.29 ± 0.07 (2.5%) ab10.92 ± 1.48 (94.9%) a0.04 ± 0.01 (0.3%) a11.50± 1.55 a
M28t00.21 ± 0.06 (1.6%) a0.04 ± 0.02 (0.3%) b0.14 ± 0.04 (1.1%) b12.82 ± 2.04 (96.7%) a0.04 ± 0.01 (0.3%) a13.25 ± 2.15 a
t100.17 ± 0.04 (3.5%) a0.33 ± 0.02 (6.9%) a0.41 ± 0.08 (8.4%) a3.94 ± 0.42 (81.1%) b0.01 ± 0.02 (0.1%) b4.86 ± 0.48 b
tR0.09 ± 0.02 (1.0%) b0.07 ± 0.00 (0.7%) ab0.20 ± 0.03 (2.0%) ab9.24 ± 0.62 (96.0%) ab0.03 ± 0.00 (0.3%) ab9.63 ± 0.60 ab
x, mean separation within each group of plants according to the Tukey-b test. y, mean separation within each group of plants according to the Kruskal-Wallis test with Bonferroni correction for multiple comparisons.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sales, E.; Cañizares, E.; Pereira, C.; Pérez-Oliver, M.A.; Nebauer, S.G.; Pavlović, I.; Novák, O.; Segura, J.; Arrillaga, I. Changing Temperature Conditions during Somatic Embryo Maturation Result in Pinus pinaster Plants with Altered Response to Heat Stress. Int. J. Mol. Sci. 2022, 23, 1318. https://doi.org/10.3390/ijms23031318

AMA Style

Sales E, Cañizares E, Pereira C, Pérez-Oliver MA, Nebauer SG, Pavlović I, Novák O, Segura J, Arrillaga I. Changing Temperature Conditions during Somatic Embryo Maturation Result in Pinus pinaster Plants with Altered Response to Heat Stress. International Journal of Molecular Sciences. 2022; 23(3):1318. https://doi.org/10.3390/ijms23031318

Chicago/Turabian Style

Sales, Ester, Eva Cañizares, Catia Pereira, María Amparo Pérez-Oliver, Sergio G. Nebauer, Iva Pavlović, Ondřej Novák, Juan Segura, and Isabel Arrillaga. 2022. "Changing Temperature Conditions during Somatic Embryo Maturation Result in Pinus pinaster Plants with Altered Response to Heat Stress" International Journal of Molecular Sciences 23, no. 3: 1318. https://doi.org/10.3390/ijms23031318

APA Style

Sales, E., Cañizares, E., Pereira, C., Pérez-Oliver, M. A., Nebauer, S. G., Pavlović, I., Novák, O., Segura, J., & Arrillaga, I. (2022). Changing Temperature Conditions during Somatic Embryo Maturation Result in Pinus pinaster Plants with Altered Response to Heat Stress. International Journal of Molecular Sciences, 23(3), 1318. https://doi.org/10.3390/ijms23031318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop