Next Article in Journal
Anti-Inflammatory Effects of Huberia peruviana Cogn. Methanol Extract by Inhibiting Src Activity in the NF-κB Pathway
Previous Article in Journal
Preharvest Spray Hexanal Formulation Enhances Postharvest Quality in ‘Honeycrisp’ Apples by Regulating Phospholipase D and Calcium Sensor Proteins Genes
Previous Article in Special Issue
Comparative Transcriptomics of Non-Embryogenic and Embryogenic Callus in Semi-Recalcitrant and Non-Recalcitrant Upland Cotton Lines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heat Stress in Pinus halepensis Somatic Embryogenesis Induction: Effect in DNA Methylation and Differential Expression of Stress-Related Genes

by
Cátia Pereira
1,2,
Ander Castander-Olarieta
2,
Ester Sales
3,
Itziar A. Montalbán
2,
Jorge Canhoto
1,* and
Paloma Moncaleán
2,*
1
Center for Functional Ecology, Department of Life Sciences, University of Coimbra, 3000-456 Coimbra, Portugal
2
Department of Forestry Science, NEIKER-BRTA, 01080 Arkaute, Spain
3
Departament of Ciencias Agrarias y del Medio Natural, Instituto Universitario de Ciencias Ambientales, Universidad de Zaragoza, Escuela Politécnica Superior, 22071 Huesca, Spain
*
Authors to whom correspondence should be addressed.
Plants 2021, 10(11), 2333; https://doi.org/10.3390/plants10112333
Submission received: 3 September 2021 / Revised: 25 October 2021 / Accepted: 26 October 2021 / Published: 29 October 2021

Abstract

:
In the current context of climate change, plants need to develop different mechanisms of stress tolerance and adaptation to cope with changing environmental conditions. Temperature is one of the most important abiotic stresses that forest trees have to overcome. Recent research developed in our laboratory demonstrated that high temperatures during different stages of conifer somatic embryogenesis (SE) modify subsequent phases of the process and the behavior of the resulting ex vitro somatic plants. For this reason, Aleppo pine SE was induced under different heat stress treatments (40 °C for 4 h, 50 °C for 30 min, and 60 °C for 5 min) in order to analyze its effect on the global DNA methylation rates and the differential expression of four stress-related genes at different stages of the SE process. Results showed that a slight decrease of DNA methylation at proliferating embryonal masses (EMs) can correlate with the final efficiency of the process. Additionally, different expression patterns for stress-related genes were found in EMs and needles from the in vitro somatic plants obtained; the DEHYDRATION INDUCED PROTEIN 19 gene was up-regulated in response to heat at proliferating EMs, whereas HSP20 FAMILY PROTEIN and SUPEROXIDE DISMUTASE [Cu–Zn] were down-regulated in needles.

Graphical Abstract

1. Introduction

As long-lived sessile organisms with complex life cycles, plants need to develop different mechanisms of protection and adaptation for a broad range of biotic and abiotic stresses in order to maximize growth, reproduction, and survival [1,2]. As well as genetics, epigenetics has become an emerging and promising research field to understand tree phenotypic plasticity and adaptive responses [3,4,5,6,7].
Alterations in epigenetic marks are reversible enzyme-mediated modifications of DNA and/or associated histones that regulate transcriptional activity of genes as well as their sequences [8,9,10]. The most studied epigenetic mark is DNA methylation because of its stability, its incidence in both plants and mammals, and its influence on gene expression and genome structure regulation [11].
Most epigenetics marks are reverted when the environmental constraints that triggered them are no longer present. However, higher plants appear to be able to retain some “stress memory” or “stress imprinting”, since a first stress exposure often leads to an enhanced resistance to a later stress [12,13]. What is also known as “priming” or “hardening”, prior exposure to the eliciting factors leads to a faster and stronger induction of basal resistance mechanisms (or greater tolerance) against them. In some cases, crosstalk among different stimuli can happen, leading to multiple stress memory attainments (cross-priming) [14]. Heat is sometimes accompanied by other stresses, such as drought, and recent studies showed that both stresses have overlapping roles [15,16,17].
Variations in epigenetic marks were revealed to be involved in morphological and physiological changes in trees in a large number of processes, including embryogenesis, organ maturation, phase change, and bud set or burst [18,19,20,21]. Moreover, epigenetic regulation plays a critical role in modulation of multiple aspects of plant development through the adjustment of gene expression in response to environmental factors [8].
Somatic embryogenesis (SE) is a worldwide studied biotechnology tool that allows large-scale propagation for many conifers [22]. The first report of SE in Aleppo pine was carried out in our laboratory [23]. Later, a subsequently developed experiment showed that changes in temperature and water availability at the induction phase of SE affects the success of the process in this species [15]. Considering those results, we focused on the application of higher temperatures (40 °C (4 h), 50 °C (30 min), and 60 °C (5 min)) at the initial stage of SE to see if a “priming” effect could be obtained. It was already found that this stress application during SE induction can modulate the morphology and hormonal profiles of embryonal masses (EMs) as well as the efficiency of the process itself [24].
Notwithstanding that stress affects plants at different levels and several defense mechanisms are activated, it leads to the accumulation of reactive oxygen species (ROS) that can reach toxic levels and cause cell damage and death [25]. To avoid this, plants have evolved antioxidant machinery consisting of enzymatic components, such as superoxide dismutase [26]. Additionally, as previously mentioned, heat is sometimes accompanied with drought, and the dehydration-induced 19 family of proteins appears to be involved in the response to both stresses [27,28]. Finally, the induction of heat shock proteins (HSPs) seems to be essential to facilitate continued homeostasis and survival against heat stress [29,30].
In this sense, in the present study, the effects of different temperature treatments applied during the initial stage of Aleppo pine SE on the epigenetic patterns were evaluated. For this reason, levels of cytosine residues 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), both at proliferating EMs and at needles from the in vitro somatic plants produced, were measured. At the same time, the expression of stress-related genes involved in the defense mechanisms mentioned above was analyzed to assess if the initial heat stress triggered long-lasting modifications at the transcriptome along the different stages of the process.

2. Results

2.1. Global DNA Methylation/Hydroxymethylation Analysis

No statistically significant differences were found regarding global DNA methylation rates (%) between samples from the different induction temperature treatments applied at the initiation stage of P. halepensis SE (Table S1).
Regarding the results obtained for EMs, the treatment that presented the lowest levels of 5mC was 60 °C (5 min) (37.52%), followed by the control (23 °C) (38.01%), and the highest methylation rate was obtained at 50 °C (30 min) with a difference of 3.3% with respect to the lowest (40.82%) (Table S2). The values obtained in needles were very similar between treatments (Table S2). All samples presented high rates of global DNA methylation (between 37.52 and 41.56%) (Table S2).
The detection of 5hmC was not possible in proliferating EMs, and it was only achieved, at very low levels, in five out of twelve analyzed samples of needles from in vitro somatic plants. Thus, no further analysis was carried out concerning hydroxymethylation data. Even so, it should be noticed that three of the five detected values at needles corresponded to the samples from the control (23 °C).

2.2. Relative Expression of Stress-Related Genes

The different induction temperature treatments applied during the initial stage of P. halepensis SE led to changes in the relative expression patterns of the stress-related genes with respect to the control (23 °C).
The results found in proliferating EMs showed that statistically significant differences were found for the relative expression of DEHYDRATION INDUCED PROTEIN 19 (DI19) (Table S3); in both 50 °C (30 min) and 60 °C (5 min) treatments, despite the lower fold change obtained (Figure 1a), the relative expression of this gene was considerably higher than in the control treatment (23 °C). CHLOROPLAST SMALL HEAT PROTEIN (P439), HSP20 FAMILY PROTEIN (P444), and SUPEROXIDE DISMUTASE [Cu–Zn] (SOD) were slightly higher expressed in samples induced at higher temperatures when compared to the control, except for SOD at the 50 °C treatment (Figure 1a), but these differences were not significant.
In contrast, in needles from in vitro somatic plants, statistically significant differences were found for the relative expression of P444 and SOD genes. A gradual decrease of P444 relative expression was observed when priming temperature increased, being significantly repressed compared to control (23 °C) for the 60 °C treatment. Expression of the SOD gene was lower in relation to the control (23 °C), but significant differences were only detected in samples coming from the 50 °C treatment, which attained the highest fold change (Figure 1b). P439 and DI19 were slightly down-regulated at samples induced at higher temperatures when compared to the control (23 °C), but these differences were not significant. In conclusion, induction of SE at higher temperatures resulted in different relative expression patterns of these stress-related genes in proliferating EMs and in needles from in vitro somatic plants. At proliferating EMs, the stress-related genes studied were generally overexpressed in primed plant material, whilst in the in vitro somatic plants obtained, their expression was lower in primed than in control.

3. Discussion

Accumulating data have shown that epigenetics changes are involved in many physiological processes, and the understanding of these mechanisms is crucial for forest tree management and breeding in the context of climate change [31,32,33]. One of the most important epigenetic marks is DNA methylation. It occurs through the addition of a methyl group at position 5 of the pyrimidine ring of cytosine in the CG, CHG, and CHH (where H = A, T, or C) contexts [34,35], and understanding how its levels change under stressful conditions can lead to a better knowledge of the plant response to environmental changes [11,36]. Considering the use of epigenetic variations for breeding applications relies on their transmission features, and since 5mC patterns can be transmitted through mitosis as well as meiosis, this DNA methylation mark could be valuable in all crops regardless of their propagation method [37]. In this sense, the concentrations of 5mC and 5hmC were assessed in order to perform a global DNA methylation study in samples from proliferating EMs and needles from in vitro somatic plants of Aleppo pine induced under high temperatures.
No statistically significant differences in methylation status were found between the induction treatments applied. Nonetheless, at proliferating EMs, a difference of 3.3% was found between the 60 °C (5 min) treatment that presented the lowest concentrations of 5mC with respect to the highest methylation rate obtained at 50 °C (30 min). In Pinus nigra, specific DNA methylation levels were analyzed in tissues with different embryogenic potentials, and the lowest levels were found in those with higher embryogenic capacity [38]. In Coffea canephora indirect SE, 5mC rates were stable during cell differentiation, and a significant increase occurred during somatic embryos regeneration [39]. Several studies suggest that DNA methylation is critical for SE success, and it is common to find lower values in embryogenic tissues with respect to the non-embryogenic tissues [5,40,41]. Taking this into account and correlating the results presented here with the ones obtained in a previous study, where the effect of Pinus halepensis SE induction under high temperatures on the efficiency of the process was assessed, we can reaffirm this plausible hypothesis. The highest number of somatic embryos produced was obtained in EMs from the 60 °C (5 min) treatment (a mean of 317.4 somatic embryos per gram) followed by those from the control; the lowest number was gained when a priming treatment of 50 °C (30 min) was applied (a mean of 135.2 somatic embryos per gram) [24].
Considering the effect of heat stress on the variation of DNA methylation levels, different species and cell types display different responses. In a similar study, performed with P. radiata, the same decrease in samples induced at the highest temperature of 60 °C (5 min) compared to the control temperature and the 40 °C (4 h) treatment was found [42]. On the contrary, in Norway spruce, seedlings originating from a warm embryonic environment [43] and the exposure of Arabidopsis plants to heat [44] resulted in an increase of global DNA methylation. Interestingly, also in Arabidopsis, DNA methylation increased during heat stress, followed by the reduction of DNA methylation levels after transfer to the control conditions [36].
Needles presented slightly higher concentrations of 5mC relative to EMs, and the global DNA methylation values were stabilized between treatments at this stage. This is in accordance with the fact that a gradual increase in DNA methylation throughout ageing was previously reported for several forest trees [45,46]. In P. radiata needles from field trees, DNA methylation increased with both ageing and phase change [47]. In contrast, 5mC levels in proliferating EMs and needles from one-year-old somatic plants of P. radiata were similar [42].
The 5hmC was detected in needles from in vitro somatic plants of Aleppo pine. As reviewed in [48], numerous studies indicate that 5hmC acts not only as an intermediate during 5mC demethylation but also plays important roles during maintenance of pluripotency in animal embryonic stem cells; however, in plants the roles of 5hmC during development are still unknown. The first discovery of its presence in conifers was in Norway spruce [49], and it was found that its concentration can fluctuate not only at different tissues but also at different temperatures in P. radiata [42].
The differential expression of stress-related genes (P439, P444, SOD, and DI19) in primed plant material was also assessed in this study. To respond and adapt to different stresses, plants have developed a complex of molecular mechanisms by modulating the expression of a specific set of genes [11,33,50,51], and we found that differential expression patterns changed between the two stages of the propagation process analyzed.
Regarding the relative expression of the stress-related genes in proliferating EMs, statistically significant differences were only found for DI19. A gradual increase at the relative expression of this gene along the higher temperature treatments was found and in EMs initiated at 50 °C (30 min) and 60 °C (5 min) was considerably higher than in the control (23 °C). The DI19 family of proteins is a novel type of Cys2/His2 zinc-finger proteins involved in the response to several abiotic stresses, especially drought [52]. In Populus simonii, it was found that there is an overlapping heat–drought response [53], and our results suggest that DI19 overexpression is involved in the abiotic stress response in P. halepensis. In accordance, in rice, the overexpression of the DI19-4 resulted in significantly increased tolerance to drought stress [54], and in Arabidopsis, DI19-1 overexpressing lines presented higher tolerance to drought stress than the wild-type lines [55]. On the other hand, in P. radiata emerging EMs induced at higher temperatures, DI19 initially presented similar levels to control; however, after 4 weeks, actively proliferating EMs induced at 40 °C (4 h) were significantly underexpressed [42]. In another study in Arabidopsis, transgenic lines with DI19–3 overexpression were more sensitive to salinity and drought than wild-types [52]. It appears that different protein members of this family, despite being clearly involved in stress response mechanisms, may have different functions unknown yet.
No differences were obtained in proliferating EMs for P439, P444, and SOD. However, needles from in vitro grown somatic plants showed statistically significant differences in transcript abundance for both P444 and SOD genes. The relative expression of these genes was lower in plants coming from EMs initiated at high temperatures when compared to the control plants. Heat stress severely affects the stability of various cellular components, causing a state of metabolic imbalance and a cascade of cellular reactions. This disruption of the steady-state flux of cellular metabolites usually leads to the accumulation of toxic products, such as ROS [50,53]. Furthermore, HSPs are closely involved in cellular protection, its structures and responses to heat stress being highly conserved amongst several organisms [29,56]. According to this, in maritime pine, when heat priming was performed at immature megagametophytes through SE, similar levels of expression were observed in primed and control EMs for HSP70 and SOD, but the expression of HSP70 at the derived in vitro somatic plants was higher in control plants [57]. In P. radiata, the gene coding for a heat shock protein (HSP20) was down-regulated in proliferating EMs and somatic plants coming from EMs initiated at high temperatures [42]. In contrast, P. radiata one-year-old seedlings subjected to heat treatments showed significantly higher short-term expression of P439 and P444 [58]. Different sampling times, as well as different tissues, seem to lead to different relative expressions of stress-related genes. Taking this into account, we suggest that in P. halepensis stress-related gene overexpression may happen during and/or shortly after the heat stress occurs. With time and at different tissues, the priming effect leads to their stabilization and their lower expression relative to controls.
Finally, it is important to note that different patterns of stress-related genes relative expression have been observed. In this sense, DI19 presented a gradual increase along the higher temperatures at EMs and P444 a gradual decrease in the needles from in vitro somatic plants. In its turn, SOD was significantly repressed in needles sampled in plants derived from priming at 50 °C (30 min), while those from plants primed at 60 °C (5 min) showed expression rates similar to non-primed plants. When the effect of these temperature treatments in EMs endogenous cytokinin (CK) profiles was assessed in a previous work on this species, diverse patterns were also found for different CKs [24]. Those results were related to the different induction times between treatments and the concept that for some CKs, temperature treatments acted as short or mild stress, while for others they were sensed as a prolonged or more severe stress. In Arabidopsis, heat stress has been connected with fluctuations in the endogenous levels of CKs and ABA that seem to be involved in HSPs regulation [59]. It is possible that an interaction between CKs and the regulation of other stress-related genes, in response to heat stress, also occurs in P. halepensis.

4. Materials and Methods

4.1. SE Temperature Experiment and Plant Material Collection

Induction of Pinus halepensis EMs under different temperatures was performed as described in [24]. Briefly, one-year-old green female cones, enclosing immature seeds of P. halepensis from five open pollinated trees were used; storage and preparation of plant material was the same as described in [60]. Whole megagametophytes were placed horizontally on DCR initiation medium [61]. For temperature treatments, closed Petri dishes containing initiation medium were preheated for 30 min, and immature megagametophytes were cultured at 40, 50, and 60 °C for 4 h, 30 min, and 5 min, respectively. As control, 23 °C was used, and, after the application of the different treatments, all explants were kept at standard conditions in darkness.
After nine weeks on the initiation medium, proliferating EMs were detached from the megagametophyte and transferred to the proliferation medium. This medium had the same composition to that used in the initiation stage, but a higher gellan gum concentration (4.5 g L−1). EMs were subcultured every two weeks and kept in the dark.
Following 4 subcultures, fresh tissue from twenty proliferating embryogenic cell lines (ECLs) were immersed in liquid nitrogen and immediately stored at −80 °C until further analysis. Eight/nine ECLs per treatment were selected to maturation at DCR medium supplemented with 60 g L−1 sucrose, 75.0 µM abscisic acid, the EDM amino acid mixture [62], and 9 g L−1 Gelrite®.
For germination, somatic embryos were transferred to Petri dishes containing half-strength macronutrient LP medium [63,64] supplemented with 2 g L−1 of activated charcoal and 9.5 g L−1 Difco® granulated agar (Becton Dickinson, Franklin Lakes, NJ, USA). Cultures were cultured under dim light for 7 days and afterwards were kept under a 16:8 h photoperiod at 100 µmol m−2 s−1 provided by cool white fluorescent tubes (TFL 58 W/33, Philips, France). The obtained plantlets were transferred onto fresh medium of the same composition every 4 weeks. During the first 8 weeks, plantlets were cultured in Petri dishes and then transferred to glass culture vessels. After 6 months, needles from in vitro somatic plants from eleven ECLs were immersed in liquid nitrogen and immediately stored at −80 °C until further analysis.

4.2. Global DNA Methylation/Hydroxymethylation Analysis

Genomic DNA extraction and subsequent methylation analysis were performed both on samples from twenty proliferating EMs and needles from in vitro somatic plants from eleven ECLs (comprising five samples of EMs and three samples of needles, per treatment). Previously collected samples were lyophilized, and 15 mg of homogenized lyophilized tissue were used.
DNA extraction and its hydrolyzation were performed as described in [42].
Methylation and hydroxymethylation levels of cytosine were analyzed on a 1200 Series HPLC system coupled to a 6410 Triple Quad mass spectrometer (Agilent Technologies, Santa Clara, CA, USA). The chromatographic separation was performed on a Zorbax SB-C18 column (2.1 × 100 mm, 3.5 µm, Agilent Technologies). The mobile phase was 11% methanol and 0.1% formic acid in water, and 5 µL of samples were injected in the column at a flow rate of 0.1 mL min−1. The electrospray ionization source (ESI) was operated in the positive ion multiple reaction monitoring mode (MRM) set to an ion spray voltage of 3500 V, 40 psi for nebulizer, and source temperature at 350 °C. The intensities of specific MH+→ fragment ion transitions were recorded (5mC m/z 242→126, 5hmC m/z 258→142, and C m/z 228→112). Identification of cytosine, 5mC, and 5hmC was assessed by injection of commercial standards (5-Methylcytosine and 5-Hydroxymethylcytosine DNA Standard Set, Zymo Research, Irvine, CA, USA) under the same LC-ESI-MS/MS-MRM conditions. The percentage of 5mC and 5hmC at each sample was calculated from the MRM peak area divided by the combined peak areas for 5mC, 5hmC, and cytosine.

4.3. Relative Expression of Stress-Related Genes

RNA extraction and further analysis of expression patterns from stress-related genes were performed on samples from sixteen proliferating EMs and needles from in vitro somatic plants from eleven embryogenic lines (comprising four samples of EMs and three samples of needles, per treatment). As initial material, 10 mg of lyophilized tissue, previously grinded for homogenization in a TissueLyser II (Qiajen, Hilden, Germany), were used. The analysis was performed based on the protocol described in [42].
Total RNA extraction was carried out using a plant/fungi total RNA purification kit (Norgen Biotek Corp., Thorold, ON, Canada), and genomic DNA was degraded by using recombinant DNase I (RNase-free, Takara Bio Inc., Shiga, Japan), following manufacturer’s instructions. A NanodropTM 2000 was used for RNA quantification, and its integrity was assessed by agar gel electrophoresis.
cDNA was synthesized from 1000 ng of RNA using the PrimeScript RT Reagent Kit (Takara) and random hexamers as primers following the manufacturer’s instructions. Real time PCR amplifications were performed in StepOne Plus (Applied Biosystems, Carlsbad, CA, USA), using a final volume of 20 µL containing 0.8 µM of each primer and 10 µL of SYBR Green I Master mix (Takara Bio Inc., Shiga, Japan) in triplicate for each sample. The PCR conditions were an initial denaturation at 95 °C for 20 s, followed by 40 cycles of 95 °C for 3 s and 60 °C for 30 s.
Primers previously described [28,58] were used and their efficiencies estimated using the qPCR Efficiency Calculator available at Thermofisher.com, based on the standard curve previously developed with four dilution points for each primer. Analyzed genes as well as primers details are summarized in Table 1.
The relative transcript levels were normalized using ACTIN (ACT), and the relative expression of each gene (R) was calculated on the basis of ΔCt values using the following formula: R = 2−ΔCt [65]. Finally, the fold changes between expression values obtained at control treatment (23 °C) and different temperature treatments were calculated in logarithmic scale.

4.4. Statistical Analysis

A one-way analysis of variance, through the application of the non-parametric Kruskal–Wallis test, was carried out to assess the effect of different temperature treatments both for total DNA methylation rates (%) and for genes’ relative expression values. When significant differences were found (p < 0.05), Dunn’s multiple comparison test was carried out to find out which treatments were statistically different.

5. Conclusions

As far as it is known, this is the first report concerning the effect in DNA methylation and expression of stress-related genes in response to heat stress application during SE induction in Aleppo pine.
Regarding the DNA methylation results, the temperatures treatments applied were not enough to provoke significantly different levels of DNA methylations at the analyzed samples. Nonetheless, it is important to note that the treatment that presented the lowest methylation level was also the one that produced the highest number of somatic embryos (60 °C, 5 min) [24]. It appears that lower levels of DNA methylation/DNA hypomethylation are associated with higher embryogenic capacity and, therefore, with a higher number of somatic embryos produced.
Concerning the application of high temperatures at the early stage of P. halepensis SE at the expression of stress-related genes, long-term changes in the differential expression of stress-related genes, specifically DI19, P444, and SOD, at different stages of SE were found. Despite the fact that the pattern of overexpression of primed EMs and lower expression of primed needles were consistent, it appears that different induction times between treatments had an effect on the relative expression of the stress-related genes studied.
Further analysis concerning the effect of heat stress application during the induction phase of P. halepensis SE regarding its effect on metabolomics and protein profiles will be performed. This data can lead to a better understanding of all the mechanisms involved on the heat stress response in this species.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10112333/s1, Table S1: One-way analysis of variance for methylation rates (%) detected in P. halepensis embryonal masses (EMs) and needles from in vitro somatic plants induced under different temperature treatments (23 °C, 9 weeks; 40 °C, 4 h; 50 °C, 30 min; 60 °C, 5 min). Table S2: Total methylation rates (%) detected in P. halepensis proliferating embryonal masses (EMs) and needles from in vitro somatic plants induced under different temperature treatments (23 °C (control); 40 °C, 4 h; 50 °C, 30 min; 60 °C, 5 min). Table S3: One-way analysis of variance for expression of different genes detected in P. halepensis embryonal masses (EMs) and needles from in vitro somatic plants induced under different temperature treatments (23 °C, 9 weeks; 40 °C, 4 h; 50 °C, 30 min; 60 °C, 5 min).

Author Contributions

Conceptualization, P.M. and J.C.; methodology, P.M., J.C. and I.A.M.; formal analysis, C.P., A.C.-O. and E.S.; investigation, C.P., A.C.-O. and I.A.M.; data curation, C.P.; writing—original draft preparation, C.P.; writing—review and editing, C.P., J.C., P.M. and I.A.M.; supervision, P.M. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MINECO (Spanish Government) project (AGL2016-76143-C4-3R), BIOALI-CYTED (P117RT0522), DECO (Basque government, Ayudas de formación a jóvenes investigadores y tecnólogos), F4F-Forest for the future (CENTRO-08-5864-FSE-000031, Programa Operacional Regional do Centro, Fundo Social Europeu), the Portuguese Foundation for Science and Technology (FCT) (SFRH/BD/123702/2016), Fundo Social Europeu (FSE), Programa Operacional Regional do CENTRO-Centro 2020 (EU), and MULTIFOREVER (Project MULTIFOREVER is supported under the umbrella of ERA-NET cofund Forest Value by ANR (FR), FNR (DE), MINCyT (AR), MINECO-AEI (ES), MMM (FI), and VINNOVA (SE)). Forest Value has received funding from the European Union’s Horizon 2020 research and innovation programmed under agreement No 773324.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

It is all original.

Acknowledgments

We want to thank Mónica Meijón and Maria Jesús Cañal for their help and facilities availability for DNA extraction. This work was carried out at the R&D Unit Center for Functional Ecology-Science for People and the Planet (CFE), with reference UIDB/04004/2020, financed by FCT/MCTES through national funds (PIDDAC).

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. Baránek, M.; Křižan, B.; Ondrušíková, E.; Pidra, M. DNA-methylation changes in grapevine somaclones following in vitro culture and thermotherapy. Plant Cell Tissue Organ. Cult. 2010, 101, 11–22. [Google Scholar] [CrossRef]
  2. Edreva, A.; Velikova, V.; Tsonev, T.; Dagnon, S.; Gürel, A.; Aktaş, L.; Gesheva, E. Stress-protective role of secondary metabolites: Diversity of functions and mechanisms. Gen. Appl. Plant Physiol. 2008, 34, 67–78. [Google Scholar]
  3. Arnholdt-Schmitt, B. Stress-induced cell reprogramming. A role for global genome regulation? Plant Physiol. 2004, 136, 2579–2586. [Google Scholar] [CrossRef] [Green Version]
  4. Correia, B.; Valledor, L.; Meijón, M.; Rodriguez, J.L.; Dias, M.C.; Santos, C.; Cañal, M.J.; Rodriguez, R.; Pinto, G. Is the interplay between epigenetic markers related to the acclimation of Cork oak plants to high temperatures? PLoS ONE 2013, 8, e53543. [Google Scholar] [CrossRef] [Green Version]
  5. Lee, K.; Seo, P.J. Dynamic epigenetic changes during plant regeneration. Trends Plant Sci. 2018, 23, 235–247. [Google Scholar] [CrossRef] [PubMed]
  6. Ribeiro, T.; Viegas, W.; Morais-Cecílio, L. Epigenetic marks in the mature pollen of Quercus suber L. (Fagaceae). Sex. Plant Reprod. 2009, 22, 1–7. [Google Scholar] [CrossRef] [Green Version]
  7. Lira-Medeiros, C.F.; Parisod, C.; Fernandes, R.A.; Mata, C.S.; Cardoso, M.A.; Ferreira, P.C.G. Epigenetic variation in mangrove plants occurring in contrasting natural environment. PLoS ONE 2010, 5, e10326. [Google Scholar] [CrossRef]
  8. Hauser, M.T.; Aufsatz, W.; Jonak, C.; Luschnig, C. Transgenerational epigenetic inheritance in plants. Biochim. Biophys. Acta Gene Regul. Mech. 2011, 1809, 459–468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Smulders, M.J.M.; de Klerk, G.J. Epigenetics in plant tissue culture. Plant Growth Regul. 2010, 63, 137–146. [Google Scholar] [CrossRef] [Green Version]
  10. Us-Camas, R.; Rivera-Solís, G.; Duarte-Aké, F.; De-la-Peña, C. In vitro culture: An epigenetic challenge for plants. Plant Cell Tissue Organ. Cult. 2014, 118, 187–201. [Google Scholar] [CrossRef]
  11. Amaral, J.; Ribeyre, Z.; Vigneaud, J.; Dia Sow, R.; Fichot, R.; Messier, C.; Pinto, G.; Nolet, P.; Maury, S. Advances and promises of epigenetics for forest trees. Forests 2020, 11, 976. [Google Scholar] [CrossRef]
  12. Bruce, T.J.A.; Matthes, M.C.; Napier, J.A.; Pickett, J.A. Stressful “memories” of plants: Evidence and possible mechanisms. Plant Sci. 2007, 173, 603–608. [Google Scholar] [CrossRef]
  13. Conrath, U. Molecular aspects of defence priming. Trends Plant Sci. 2011, 16, 524–531. [Google Scholar] [CrossRef]
  14. Blödner, C.; Skroppa, T.; Johnsen, Ø.; Polle, A. Freezing tolerance in two Norway spruce (Picea abies [L.] Karst.) progenies is physiologically correlated with drought tolerance. J. Plant Physiol. 2005, 162, 549–558. [Google Scholar] [CrossRef] [PubMed]
  15. Pereira, C.; Montalbán, I.A.; García-Mendiguren, O.; Goicoa, T.; Ugarte, M.D.; Correia, S.; Canhoto, J.M.; Moncaleán, P. Pinus halepensis somatic embryogenesis is affected by the physical and chemical conditions at the initial stages of the process. J. For Res. 2016, 21, 143–150. [Google Scholar] [CrossRef]
  16. Moncaleán, P.; García-Mendiguren, O.; Novák, 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]
  17. 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 Struct. Funct. 2016, 30, 949–958. [Google Scholar] [CrossRef]
  18. Yakovlev, I.A.; Carneros, E.; Lee, Y.K.; Olsen, J.E.; Fossdal, C.G. Transcriptional profiling of epigenetic regulators in somatic embryos during temperature induced formation of an epigenetic memory in Norway spruce. Planta 2016, 243, 1237–1249. [Google Scholar] [CrossRef] [PubMed]
  19. Castander-Olarieta, A.; Montalbán, I.A.; De Medeiros Oliveira, E.; Dell’aversana, E.; D’amelia, L.; Carillo, P.; Steiner, N.; Fraga, H.P.D.F.; Guerra, M.P.; Goicoa, T.; et al. Effect of thermal stress on tissue ultrastructure and metabolite profiles during initiation of radiata pine somatic embryogenesis. Front Plant Sci. 2019, 9, 1–16. [Google Scholar] [CrossRef] [PubMed]
  20. Valledor, L.; Jorrín, J.V.; Rodríguez, J.L.; Lenz, C.; Meijón, M.; Rodríguez, R.; Cañal, M.J. Combined proteomic and transcriptomic analysis identifies differentially expressed pathways associated to Pinus radiata needle maturation. J. Proteome Res. 2010, 9, 3954–3979. [Google Scholar] [CrossRef] [PubMed]
  21. Johnsen, Ø.; Fossdal, C.G.; Nagy, N.; MØlmann, J.; Dæhlen, O.G.; SkrØppa, T. Climatic adaptation in Picea abies progenies is affected by the temperature during zygotic embryogenesis and seed maturation. Plant Cell Environ. 2005, 28, 1090–1102. [Google Scholar] [CrossRef]
  22. Von Arnold, S.; Clapham, D.; Abrahamsson, M. Embryology in Conifers. In Advances in Botanical Research, 1st ed.; Cánovas, F.N., Ed.; Elsevier: Cambridge, MA, USA, 2019; Volume 89, pp. 157–184. [Google Scholar]
  23. Montalbán, I.A.; Setién-Olarra, A.; Hargreaves, C.L.; Moncaleán, P. Somatic embryogenesis in Pinus halepensis Mill.: An important ecological species from the Mediterranean forest. Trees 2013, 27, 1339–1351. [Google Scholar] [CrossRef]
  24. Pereira, C.; Castander-Olarieta, A.; Montalbán, I.A.; Pěnčík, A.; Petřík, I.; Pavlović, I.; Oliveira, E.D.M.; Freitas Fraga, H.P.d.; Guerra, M.P.; Novák, O.; et al. Embryonal masses induced at high temperatures in Aleppo pine: Cytokinin profile and cytological characterization. Forests 2020, 11, 807. [Google Scholar] [CrossRef]
  25. Feher, A.; Ötvös, K.; Pasternak, T.P.; Pettkó-Szandtner, A. The involvement of reactive oxygen species (ROS) in the cell cycle activation (G 0 -to-G 1 transition) of plant cells. Plant Signal Behav. 2008, 3, 823–826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Almeselmani, M.; Deshmukh, P.S.; Sairam, R.K.; Kushwaha, S.R.; Singh, T.P. Protective role of antioxidant enzymes under high temperature stress. Plant Sci. 2006, 171, 382–388. [Google Scholar] [CrossRef]
  27. Liu, W.X.; Zhang, F.C.; Zhang, W.Z.; Song, L.F.; Wu, W.H.; Chen, Y.F. Arabidopsis Di19 functions as a transcription factor and modulates PR1, PR2, and PR5 expression in response to drought stress. Mol. Plant 2013, 6, 1487–1502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Alvarez, C.; Valledor, L.; Sáez, P.; Hasbún, R.; Sánchez-Olate, M.; Cañal, M.J.; Ríos, D. Changes in gene expression in needles and stems of Pinus radiata rootstock plants of different ontogenic age. Amer. J. Plant Sci. 2016, 07, 1205–1216. [Google Scholar] [CrossRef] [Green Version]
  29. Neilson, K.A.; Gayani Gammulla, C.; Mirzaei, M.; Imin, N.; Haynes, P.A. Proteomic analysis of temperature stress in plants. Wiley Online Libr. 2010, 10, 828–845. [Google Scholar] [CrossRef]
  30. Ling, Y.; Serrano, N.; Gao, G.; Atia, M.; Mokhtar, M.; Woo, Y.H.; Bazin, J.; Veluchamy, A.; Benhamed, M.; Crespi, M.; et al. Thermopriming triggers splicing memory in Arabidopsis. J. Exp. Bot. 2018, 69, 2659–2675. [Google Scholar] [CrossRef] [Green Version]
  31. Dia-Sow, M.; Allona, I.; Ambroise, C.; Conde, D.; Fichot, R.; Gribkova, S.; Jorge, V.; Le-Provost, G.; Pâques, L.; Plomion, C.; et al. Epigenetics in forest trees: State of the art and potential implications for breeding and management in a context of climate change. Adv. Bot. Res. 2018, 88, 387–453. [Google Scholar] [CrossRef]
  32. Boyko, A.; Kovalchuk, I. Epigenetic control of plant stress response. Environ. Mol. Mutag. 2008, 49, 61–72. [Google Scholar] [CrossRef]
  33. Lämke, J.; Bäurle, I. Epigenetic and chromatin-based mechanisms in environmental stress adaptation and stress memory in plants. Genome Biol. 2017, 18, 124. [Google Scholar] [CrossRef] [PubMed]
  34. Viejo, M.; Santamaría, M.E.; Rodríguez, J.L.; Valledor, L.; Meijón, M.; Pérez, M.; Pascual, J.; Hasbún, R.; Fraga, M.F.; Berdasco, M.; et al. Epigenetics, the role of DNA methylation in tree development. Methods Mol. Biol. 2012, 877, 277–301. [Google Scholar] [CrossRef] [PubMed]
  35. Kumar, S.; Kumari, R.; Sharma, V.; Sharma, V. Roles, and establishment, maintenance and erasing of the epigenetic cytosine methylation marks in plants. J. Genet. 2013, 92, 629–666. [Google Scholar] [CrossRef] [PubMed]
  36. Korotko, U.; Chwiałkowska, K.; Sańko-Sawczenko, I.; Kwasniewski, M. DNA demethylation in response to heat stress in Arabidopsis thaliana. Int. J. Mol. Sci. 2021, 22, 1555. [Google Scholar] [CrossRef]
  37. Gallusci, P.; Dai, Z.; Génard, M.; Gauffretau, A.; Leblanc-Fournier, N.; Richard-Molard, C.; Vile, D.; Brunel-Muguet, S. Epigenetics for plant improvement: Current knowledge and modeling avenues. Trends Plant Sci. 2017, 22, 610–623. [Google Scholar] [CrossRef]
  38. Noceda, C.; Salaj, T.; Pérez, M.; Viejo, M.; Cañal, M.J.; Salaj, J.; Rodriguez, R. DNA demethylation and decrease on free polyamines is associated with the embryogenic capacity of Pinus nigra Arn. cell culture. Trees 2009, 23, 1285–1293. [Google Scholar] [CrossRef]
  39. Amaral-Silva, P.; Clarindo, W.; Guilhen, J.; de Jesus Passos, A.; Sanglard, N.; Ferreira, A. Global 5-methylcytosine and physiological changes are triggers of indirect somatic embryogenesis in Coffea canephora. Protoplasma 2021, 258, 45–57. [Google Scholar] [CrossRef]
  40. De-la-Peña, C.; Nic-Can, G.I.; Galaz-Ávalos, R.M.; Avilez-Montalvo, R.; Loyola-Vargas, V.M. The role of chromatin modifications in somatic embryogenesis in plants. Front Plant Sci. 2015, 6, 635. [Google Scholar] [CrossRef] [Green Version]
  41. Bravo, S.; Bertín, A.; Turner, A.; Sepúlveda, F.; Jopia, P.; Parra, M.J.; Castillo, R.; Hasbún, R. Differences in DNA methylation, DNA structure and embryogenesis-related gene expression between embryogenic and non-embryogenic lines of Pinus radiata D. don. Plant Cell Tissue Organ. Cult. 2017, 130, 521–529. [Google Scholar] [CrossRef]
  42. 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]
  43. Johnsen, Ø.; Dæhlen, O.G.; Østreng, G.; Skrøppa, T. Daylength and temperature during seed production interactively affect adaptive performance of Picea abies progenies. New Phytol. 2005, 168, 589–596. [Google Scholar] [CrossRef] [PubMed]
  44. Boyko, A.; Blevins, T.; Yao, Y.; Golubov, A.; Bilichak, A.; Ilnytskyy, Y.; Hollander, J.; Meins, F., Jr.; Kovalchuk, I. Transgenerational adaptation of Arabidopsis to stress requires DNA methylation and the function of Dicer-Like proteins. PLoS ONE 2010, 5, e9514. [Google Scholar] [CrossRef]
  45. Bräutigam, K.; Vining, K.J.; Lafon-Placette, C.; Fossdal, C.G.; Mirouze, M.; Marcos, J.G.; Fluch, S.; Fraga, M.F.; Guevara, M.Á.; Abarca, D.; et al. Epigenetic regulation of adaptive responses of forest tree species to the environment. Ecol. Evol. 2013, 3, 399–415. [Google Scholar] [CrossRef] [PubMed]
  46. Valledor, L.; Hasbún, R.; Meijón, M.; Rodríguez, J.L.; Santamaría, E.; Viejo, M.; Berdasco, M.; Feito, I.; Fraga, M.F.; Cañal, M.J.; et al. Involvement of DNA methylation in tree development and micropropagation. Plant Cell Tissue Organ. Cult. 2007, 91, 75–86. [Google Scholar] [CrossRef] [Green Version]
  47. Fraga, M.F.; Cañal, M.; Rodríguez, R. Phase-change related epigenetic and physiological changes in Pinus radiata D. Don. Planta 2002, 215, 672–678. [Google Scholar] [CrossRef]
  48. Shi, D.; Ali, I.; Tang, J.; Yang, W. New insights into 5hmC DNA modification: Generation, distribution and function. Front. Genet. 2017, 8, 1–11. [Google Scholar] [CrossRef]
  49. Yakovlev, I.A.; Gackowski, D.; Abakir, A.; Viejo, M.; Ruzov, A.; Olinski, R.; Starczak, M.; Fossdal, C.G.; Krutovsky, K.V. Mass spectrometry reveals the presence of specific set of epigenetic DNA modifications in the Norway spruce genome. Sci. Rep. 2019, 9, 19314. [Google Scholar] [CrossRef]
  50. Mittler, R.; Finka, A.; Goloubinoff, P. How do plants feel the heat? Trends Biochem. Sci. 2012, 37, 118–125. [Google Scholar] [CrossRef]
  51. Escandón, M.; Cañal, M.J.; 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. 2015, 36, 63–77. [Google Scholar] [CrossRef]
  52. Milla, M.A.; Townsend, J.; Chang, I.F.; Cushman, J.C. The Arabidopsis AtDi19 gene family encodes a novel type of Cys2/His2 zinc-finger protein implicated in ABA-independent dehydration, high-salinity stress and light signaling pathways. Plant Mol. Biol. 2006, 61, 13–30. [Google Scholar] [CrossRef] [PubMed]
  53. Jia, J.; Zhou, J.; Shi, W.; Cao, X.; Luo, J.; Polle, A.; Luo, Z.B. Comparative transcriptomic analysis reveals the roles of overlapping heat-/drought-responsive genes in poplars exposed to high temperature and drought. Sci. Rep. 2017, 7, 43215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, L.; Yu, C.; Chen, C.; He, C.; Zhu, Y.; Huang, W. Identification of rice Di19 family reveals OsDi19—4 involved in drought resistance. Plant Cell Rep. 2014, 33, 2047–2062. [Google Scholar] [CrossRef]
  55. Qin, L.-X.; Li, Y.; Li, D.-D.; Xu, W.-L.; Zheng, Y.; Li, X.-B. Arabidopsis drought-induced protein Di19-3 participates in plant response to drought and high salinity stresses. Plant Mol. Biol. 2014, 86, 609–625. [Google Scholar] [CrossRef]
  56. Efeoǧlu, B. Heat shock proteins and heat shock response in plants. Gazi. Univ. J. Sci. 2009, 22, 67–75. [Google Scholar]
  57. 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]
  58. 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]
  59. Dobrá, J.; Černý, M.; Štorchová, H.; Dobrev, P.; Skalák, J.; Jedelský, P.L.; Lukšanová, H.; Gaudinová, A.; Pešek, B.; Malbecka, J.; et al. The impact of heat stress targeting on the hormonal and transcriptomic response in Arabidopsis. Plant Sci. 2015, 231, 52–61. [Google Scholar] [CrossRef]
  60. Montalbán, I.A.; García-Mendiguren, O.; Goicoa, T.; Ugarte, M.D.; Moncaleán, P. Cold storage of initial plant material affects positively somatic embryogenesis in Pinus radiata. New For. 2015, 46, 309–317. [Google Scholar] [CrossRef]
  61. Gupta, P.K.; Durzan, D.J. Plantlet regeneration via somatic embryogenesis from subcultured callus of mature embryos of Picea abies (Norway spruce). Vitr. Cell. Dev. Biol. 1986, 22, 685–688. [Google Scholar] [CrossRef]
  62. Walter, C.; Find, J.I.; Grace, L.J. Somatic embryogenesis and genetic transformation in Pinus radiata. In Protocol for Somatic Embryogenesis in Woody Plants; Jain, S.M., Gupta, P.K., Eds.; Springer: Dordrecht, The Netherlands, 2005; pp. 11–24. [Google Scholar] [CrossRef]
  63. Quoirin, M.; Lepoivre, P. Improved media for in vitro culture of Prunus sp. Acta Hortic. 1977, 78, 437–442. [Google Scholar] [CrossRef]
  64. Aitken-Christie, J.; Singh, A.P.; Davies, H. Multiplication of meristematic tissue: A new tissue culture system for radiata pine. In Genetic Manipulation of Woody Plants; Hanover, J.W., Keathley, D.E., Wilson, C.M., Kuny, G., Eds.; Springer: Boston, MA, USA, 1988; pp. 413–432. [Google Scholar] [CrossRef]
  65. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Fold-relative gene expression of four stress-related genes (P439, P444, DI19, and SOD) between P. halepensis samples induced under different temperature treatments (40 °C, 4 h; 50 °C, 30 min; 60 °C, 5 min) from (a) proliferating embryonal masses (EMs); (b) needles from in vitro somatic plants. Data are presented as mean values ± SE, and * represents statistically significant differences at p < 0.05 of different temperature treatments with respect to the control (23 °C).
Figure 1. Fold-relative gene expression of four stress-related genes (P439, P444, DI19, and SOD) between P. halepensis samples induced under different temperature treatments (40 °C, 4 h; 50 °C, 30 min; 60 °C, 5 min) from (a) proliferating embryonal masses (EMs); (b) needles from in vitro somatic plants. Data are presented as mean values ± SE, and * represents statistically significant differences at p < 0.05 of different temperature treatments with respect to the control (23 °C).
Plants 10 02333 g001
Table 1. List of primers used in quantitative real time PCR (qRT-PCR) for relative expression analysis. Names of the genes, forward and reverse primer sequences, and melting temperatures of primers are described.
Table 1. List of primers used in quantitative real time PCR (qRT-PCR) for relative expression analysis. Names of the genes, forward and reverse primer sequences, and melting temperatures of primers are described.
IDNameForward (5′ → 3′)Reverse (5′ → 3′)Tm (°C)
ACTACTINCACTGCACTTGCTCCCAGTAAACCTCCGATCCAAACACTG60
P439CHLOROPLAST SMALL HEAT PROTEINAAGTTGTCGGTTCGAACCCCCAGAACACCGTCCTCCACAG62
P444HSP20 FAMILY PROTEINTTTCCGACTTCTTCACGGGGTTTGACAGTCCCGGCATGTC62
DI19DEHIDRATION INDUCED PROTEIN 19ATAGATGCCCATGCTGTGTAGCTTCCCTCTGTTCCCACTTG54
SODSUPEROXIDE DISMUTASE [Cu–Zn]ACAAAACGGGTGCATGTCAACCCCATCCGCTCCTACAGTTAC66
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/plants10112333

AMA Style

Pereira C, Castander-Olarieta A, Sales E, Montalbán IA, 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(11):2333. https://doi.org/10.3390/plants10112333

Chicago/Turabian Style

Pereira, Cátia, Ander Castander-Olarieta, Ester Sales, Itziar A. Montalbán, Jorge Canhoto, and Paloma Moncaleán. 2021. "Heat Stress in Pinus halepensis Somatic Embryogenesis Induction: Effect in DNA Methylation and Differential Expression of Stress-Related Genes" Plants 10, no. 11: 2333. https://doi.org/10.3390/plants10112333

APA Style

Pereira, C., Castander-Olarieta, A., Sales, E., Montalbán, I. A., Canhoto, J., & Moncaleán, P. (2021). Heat Stress in Pinus halepensis Somatic Embryogenesis Induction: Effect in DNA Methylation and Differential Expression of Stress-Related Genes. Plants, 10(11), 2333. https://doi.org/10.3390/plants10112333

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