Next Article in Journal
Dopamine D3 Receptor Modulates Akt/mTOR and ERK1/2 Pathways Differently during the Reinstatement of Cocaine-Seeking Behavior Induced by Psychological versus Physiological Stress
Next Article in Special Issue
Abiotic Stress-Induced Leaf Senescence: Regulatory Mechanisms and Application
Previous Article in Journal
The Role of Autophagy and Apoptosis in Affected Skin and Lungs in Patients with Systemic Sclerosis
Previous Article in Special Issue
Transcriptomic, Physiological, and Metabolomic Response of an Alpine Plant, Rhododendron delavayi, to Waterlogging Stress and Post-Waterlogging Recovery
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 13
10.3390/ijms241311211
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of Diurnal and Nocturnal-Warmed Plants, the Molecular Mechanism Underlying Cold Deacclimation Response in Deschampsia antarctica

by
Dariel López
1,
Giovanni Larama
2,
Patricia L. Sáez
1 and
León A. Bravo
1,*
1
Departamento de Ciencias Agronómicas y Recursos Naturales, Facultad de Ciencias Agropecuarias y Medioambiente and Center of Plant, Soil Interactions and Natural Resources Biotechnology, Scientific and Technological Bioresource Nucleus, Universidad de La Frontera, Temuco 4811230, Chile
2
Biocontrol Research Laboratory and Scientific and Technological Bioresource Nucleus, Universidad de La Frontera, Temuco 4811230, Chile
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(13), 11211; https://doi.org/10.3390/ijms241311211
Submission received: 17 April 2023 / Revised: 5 May 2023 / Accepted: 7 May 2023 / Published: 7 July 2023
(This article belongs to the Special Issue Abiotic Stresses in Plants: From Molecules to Environment)
Browse Figures
Versions Notes

Abstract

:
Warming in the Antarctic Peninsula is one of the fastest on earth, and is predicted to become more asymmetric in the near future. Warming has already favored the growth and reproduction of Antarctic plant species, leading to a decrease in their freezing tolerance (deacclimation). Evidence regarding the effects of diurnal and nocturnal warming on freezing tolerance-related gene expression in D. antarctica is negligible. We hypothesized that freezing tolerance-related gene (such as CBF-regulon) expression is reduced mainly by nocturnal warming rather than diurnal temperature changes in D. antarctica. The present work aimed to determine the effects of diurnal and nocturnal warming on cold deacclimation and its associated gene expression in D. antarctica, under laboratory conditions. Fully cold-acclimated plants (8 °C/0 °C), with 16h/8h thermoperiod and photoperiod duration, were assigned to four treatments for 14 days: one control (8 °C/0 °C) and three with different warming conditions (diurnal (14 °C/0 °C), nocturnal (8 °C/6 °C), and diurnal-nocturnal (14 °C/6 °C). RNA-seq was performed and differential gene expression was analyzed. Nocturnal warming significantly down-regulated the CBF transcription factors expression and associated cold stress response genes and up-regulated photosynthetic and growth promotion genes. Consequently, nocturnal warming has a greater effect than diurnal warming on the cold deacclimation process in D. antarctica. The eco-physiological implications are discussed.

1. Introduction

The phenomenon known as climate change, implying shifts in temperature, precipitation, and CO2 levels among other factors, represents a major threat to agricultural and natural plant populations. Deforestation and fossil fuel burning have increased the concentration of greenhouse gases in the atmosphere, contributing to increases in the global average surface temperature of 1.0–3.7 °C during the past century [1,2]. It has been estimated that global warming will increase the frequency and intensity of stochastic extreme temperature events (freezing or heat waves), with changes in precipitation patterns during this century [2]. In addition, an asymmetric warming trend is expected, with a greater increase in minimum temperatures (nocturnal) with respect to maximum temperatures (diurnal) [2,3]. Furthermore, the warming process has been amplified in polar and high-mountain regions [4].
The climate of the Antarctic Peninsula and its associated islands has shown increases in the average annual temperature ca. 3.7 °C century−1 [5]. With longer and warmer growing seasons for the plant species that inhabit the Antarctic Peninsula, reproductive capacity and population sizes have increased [6,7,8,9]. Despite the warming trend, Antarctica is still among the coldest regions on the planet [10], and Antarctic plants must continually cope with freezing temperatures, even during the growing season [11,12]. This implies that Antarctic plants should maintain their freezing tolerance even during the growing season (summer).
Deschampsia antarctica Desv. (Poaceae), the only monocot that has naturally colonized Antarctica, is used to dealing with occasional daily warm temperatures since leaf temperatures can reach up to 20 °C on sunny days, even when the air temperature does not exceed 4 °C [13]. Although the rise of diurnal temperatures appears to have no effect on the freezing tolerance in this plant species [14], the foreseeable increase in the minimum nocturnal temperature (nocturnal warming) could lead to cold deacclimation [14].
Cold deacclimation entails the partial or total loss of previously acquired freezing tolerance traits [15], and is not just a passive reversal of cold acclimation, but rather a genetically and functionally distinct process [16,17]. Even though various studies have analyzed the cold deacclimation profile after a temperature increase in several species (Arabidopsis thaliana, Brassica napus, Hordeum vulgare, Rhododendron anthopogon, among others [17,18,19,20]), and identified up- and down-regulated genes during this process, they have not differentiated between the effects of diurnal or nocturnal temperature increments. In addition, the cold deacclimation dynamic could change among species [21], as well as the gene expression associated with this process.
The ‘C-repeat binding factor’ (CBF/DREB1) transcription factors are key in the regulation of a group of genes important in plants’ freezing tolerance [22,23,24]. As plants should not be expressing their cryoprotective mechanisms during cold deacclimation, it is expected that the CBF genes would be down-regulated. However, even though some of the physiological and molecular mechanisms of freezing tolerance seem to be similar in the main angiosperm clades, the genes and genetic pathways can vary widely among species [25]. For example, A. thaliana has 4 CBF genes, but only 3 respond to low temperatures [26], while in pooids, the number of CBF genes can exceed 20, with variable numbers in different species [27]. Thus, given the complexity of this pathway in cereals and grasses, more studies are needed to understand its functioning and interaction with variable environmental factors.
The paradox between improved growth vs. cold acclimation faces D. antarctica to a likely differential response against diurnal and nocturnal warming, this prompts the question about what happens with the expression of freezing tolerance-related genes. We propose that nocturnal warming down-regulates the expression of freezing tolerance-related genes in D. antarctica plants, while the decrease in temperatures close to 0 °C during the night is essential for transient CBF and associated gene expression, allowing D. antarctica plants to maintain freezing tolerance, even with a daytime temperature increase. Differentially expressed genes during the day and night experimental warming were analyzed to assess the differential effect of nocturnal warming in cold deacclimation regulation.

2. Results

2.1. RNA-Seq Analysis of Nocturnal Expressed Genes in D. antarctica

RNA-seq transcriptome analysis was carried out using RNA from D. antarctica leaves collected 2 h after the beginning of the nocturnal period. Three biological replicates for each treatment with and without warming were sent to the Illumina Novaseq 6000 platform for deep sequencing. The raw sequence data can be found at BioProject ID PRJNA941125.
A total of 275.3 million paired-end reads were obtained after removing the low-quality sequence and adaptor sequence. The percentage of Q30 bases was above 95.99%, while the mapping rates of all 12 libraries ranged from 86.2% to 87.6%, regarding high-quality reads (Supplementary Table S1).
After error correction and redundancy removal, de novo transcriptome assembly resulted in 193,333 transcripts corresponding to 84,162 genes, with an average GC percentage of 47.65%. The median contig length was 1127 bp, while the average contig length was 1533 bp (Supplementary Table S2). Based on benchmarking universal single-copy ortholog (BUSCO) analysis, 1549 (95.9%) of the 1614 expected embryophyte genes were identified as complete after removing redundancy, and only 4.1% were fragmented or missing (Supplementary Table S3).

2.2. Differentially Expressed Genes (DEGs) of Warmed vs. Cold-Acclimated Plants

To analyze the differences of the D. antarctica transcription levels during the nocturnal period, FPKM was used to calculate the amount of differentially expressed genes. The fold change was calculated based on a comparison of the FPKM between the cold-acclimated (CA) plants and warming treatments. The genes were considered differentially expressed when the fold change of the gene expression level was at least a two-fold change and Chi-square test (p < 0.05) FDR < 0.01.
From the total of genes found in the transcriptome, only 7.77% (6537 genes) were differentially expressed 2 h after the drop in temperatures during the night period, in warmed plants. A total of 871, 4692 and 2403 DEGs were identified in plants with nocturnal warming (NW+), diurnal warming (DW+), and diurnal-nocturnal warming (DNW+), respectively, at 2 h after the temperature drop in the nocturnal period (Figure 1). Among these DEGs, 263 were down-regulated in plants with nocturnal warming, of which 52.1% was common for NW+ and DNW+ treatments only (Figure 2A), while 190 were up-regulated in plants with nocturnal warming, and 91.1% were also common for NW+ and DNW+ treatments (Figure 2B).

2.3. Analysis of Co-Expressed DEGs in Response to Nocturnal Temperatures

The analysis of principal components (PCA) and self-organized maps (SOM), together with the grouping strategy used, allowed for the creation of eight super-nodes. Among them, two super-nodes show co-expressed genes profiles related to the presence or absence of nocturnal warming (Figure S1). The gene ontology enrichment analysis in these super-nodes showed that nocturnal-warmed plants presented a significant abundance of genes related to photosynthesis, chloroplast structures, and hormonal responses 2 h after the nocturnal period began (Figure 3A). Meanwhile, the plants without nocturnal warming presented a significant enrichment in the genes related to transporter activity and membrane structure, UDP- glycosyltransferase activity, and response to the cold (Figure 3B).
Likewise, the expression patterns of 1137 genes belonging to these 2 super-nodes show abiotic stress-related genes. Those related to cold tolerance are repressed in treatments with nocturnal warming. Among them, noteworthy genes include the ‘C-repeat binding factor’ (CBF or CBF/DREB1), sugar metabolism, dehydrins, and other protection-related proteins (Figure 4). Moreover, genes related to photosynthesis and vegetative growth, such as photosystems and antenna complex components, photosynthetic enzymes, and auxin regulation, among others, were induced in treatments with nocturnal warming (Figure 4).

2.4. Validation of RNA-Seq Gene Expression by qRT-PCR

The reliability of the transcriptome data was evaluated by quantitative real-time PCR (qRT-PCR) technology. The correlation between the log2 fold change of FPKM and the relative abundance determined by qRT-PCR was highly significant (R2 = 0.92; p < 0.05) for the 16 genes selected (Figure S2).

2.5. Nocturnal Kinetics of Interest DEGs

The relative gene expression during the nocturnal period was assessed with respect to the expression values average measured in CA plants, just before the nocturnal period started. The transcripts corresponding to ‘Dehydration-responsive element-binding protein 1A’ (DRE1A/CBF3), ‘Dehydration-responsive element-binding protein 1D’ (DRE1D/CBF4), ‘Dehydrin 1’ (DHN1), ‘E3 ubiquitin-protein ligase RGLG5’ (RGLG5), ‘Ornithine aminotransferase’ (OAT), ‘Sucrose phosphate synthase 4’ (SPS4F), ‘Sucrose synthase 4’ (SUS4), ‘UDP-glycosyltransferase 73B3’ (U73B3), and ‘Vacuolar cation/proton exchanger 1’ (CAX1) presented a significantly lower expression in the treatments with nocturnal warming, with respect to the CA plants, 2 h after beginning the night period (Figure 5A–I). Likewise, the aforementioned transcripts also presented significantly lower expressions in the treatments with nocturnal warming, compared to the plants with DW+, 2 h after beginning the nocturnal period, with the exception of SPS4F, (Figure 5A–I). However, in the case of SPS4F, both the CA and DW+ plants significantly increased their expression at 2 h with respect to their starting values (Figure 5I).
The expression of DRE1A/CBF3 was highly up-regulated during the night in CA and DW+ plants, but this was not modified significantly in the treatments with nocturnal warming (Figure 5A). On the other hand, the expression of DRE1D/CBF4 and RGLG5 remained stable overnight in CA plants, while their values increased in DW+ plants. Conversely, the above-mentioned genes in NW+ exhibited a reduced expression (only significantly for CBF4) to DNW+ similar values at 2 h, and remained low until the nocturnal period ended (Figure 5B,D). At the same time, DHN1 and SUS4 were the most down-regulated in the nocturnal warming treatments among the analyzed genes (Figure 5C,H).
In contrast, the transcripts corresponding to Fructose-bisphosphate aldolase 5 (FAB5), IAA-amino acid hydrolase ILR1-like 2 (ILL2), Ribulose bisphosphate carboxylase small subunit 1 (RBCS1), Oxygen-evolving enhancer protein 2 (PSBP2), Triose phosphate/phosphate translocator (TPT), and 1-deoxy-D-xylulose-5-phosphate synthase 1 (DXS1) significantly up-regulated their expression in treatments with nocturnal warming, at 2 h at the start of the nocturnal period, compared to CA and DW+ plants, (Figure 6A–F). Highlighting that, ILL2, DXS1, and FBA5 presented a significantly higher up-regulation in NW+ plants (Figure 6A,B,E). The transcripts corresponding to ILL2 and TPT1 maintained their significantly higher expression until the overnight period ended (Figure 6A,C). The ILL2 expression presented the higher increment during the night in plants with nocturnal warming (Figure 6A), with 9 time increments in NW+ plants at 2 h after the nocturnal period began, compared to its initial value.

3. Discussion

In the current global warming context, with more frequent drastic and stochastic temperature changes, such as heat waves and spring frosts, the early cold deacclimation of plants becomes more relevant, since leaf tissues or the whole plant may become vulnerable to subsequent freezing events, threatening their survival [21,28]. Our findings suggest that only nocturnal warming is capable of down-regulating freezing tolerant-related genes and up-regulating growth promotion and carbon assimilation-related genes, inducing cold deacclimation in D. antarctica.
Among the freezing tolerance-related genes that were only repressed by nocturnal warming were the CBF transcription factors, which regulate the expression of a group of genes essential in plant freezing tolerance, the CBF-regulon [22,23,24,29]. In particular, the ‘Dehydration-responsive element-binding protein 1A’ (DREB1A/CBF3) transcription factor was the most induced overnight, increasing its expression 17.6 times in CA plants 2 h into the nocturnal period, and up to 2457 times in plants with CD at the end of the night. This could suggest that DREB1A/CBF3 has a very important role in maintaining freezing tolerance in plants with nocturnal temperatures below the cold deacclimation threshold, even when CBF transcription factors show a high functional redundancy [23,24]. It also has been postulated that the expression of CBFs needs to be transient, and their transcripts return to basal levels a few hours after their induction at low temperatures [30]. In accordance with the above, the treatment with CD, which presented the highest levels of transcription for most of the CBF analyzed, also presented the highest levels of coding transcripts for the 14-3-3 protein (results not shown), which has been postulated to function as an intermediate in the CBFs’ ubiquitination modification for their subsequent degradation via the proteasome [30].
The CBF-regulon controls the expression of a diverse group of genes, and among them are other transcription factors, channel proteins and membrane transporters, enzymes related to sugar and proline metabolism, etc. [22,23,29]. Among these genes is the ‘Vacuolar cation/proton exchanger 1’ (CAX1), which has been reported to increase expression in cold-acclimated plants, such as the CA and DW+ plants in this study, and is key in the regulation of intracellular Ca2+ [31], one of the second messengers responsible for transmitting the cold signaling [32,33]. Some gene isoforms coding for sucrose synthase (SUS) enzymes [24,34] are also part of this regulon, and they degrade the sucrose molecule yielding UDP-glucose through a reversible reaction, allowing for its incorporation into more complex molecules synthesis, such as cellulose during the thickening of cell walls [35], or moving onto catabolic pathways as a result of an increased energy demand in the face of abiotic stress [36]. Sucrose synthase has also been found to have a predominant role in sucrose loading and unloading in phloem transport [37,38]. Both situations could be occurring in CA and DW+ plants, which would explain the induction of SUS4 in these treatments, since the cell walls’ rearrangement is an expected process in D. antarctica cold-acclimated plants [39], and cold acclimation is also an energetically demanding process, needing to mobilize sucrose to sinks organs.
Dehydrins gene expression are also regulated by CBF transcription factors [23,24]. This agrees with dehydrins transcript’s drastic reduction in plants with nocturnal warming, and with their protein levels reduction during the cold deacclimation of D. antarctica and other pooids (Poa annua, Agrostis stolonifera, and Cynodon spp.) [14,40,41,42,43]. Dehydrins are associated with freezing tolerance due to their ability to stabilize cell membranes, preventing the lipid bilayer transition to the hexagonal phase at freezing temperatures, as a result of the lipid interaction in a medium with less polarity, resulting leakage of electrolytes, and other essential cytoplasm components [44].
Proline accumulation has also been associated with cold acclimation and tolerance to low temperatures in Arabidopsis thaliana, Chrysanthemum dichrum, and Lolium perenne, among others [45,46,47,48]. In this sense, it has been described that the gene coding for ‘δ-1-pyrroline-5-carboxylate synthase 2’ (P5CS2), a key enzyme for proline synthesis, is regulated by CBFs [23]. However, the P5CS2 gene reduced its expression in all warming treatments (results not shown), while the ‘Ornithine aminotransferase’ (OAT) gene reduced its expression in plants that received nocturnal warming. The mitochondrial enzyme Ornithine aminotransferase is part of an alternative proline synthesis pathway usually activated by abiotic stress, dehydration, and oxidative stress mainly [49,50]. The OAT enzyme transforms ornithine into glutamyl-5-semi-aldehyde (GSA) for its subsequent conversion to pyrroline-5-carboxylate (P5C), which is exported to the cytoplasm and continues its conversion to proline [49]. Although the OAT genes’ expression is not significantly affected during Chrysanthemum dichrum cold acclimation [47], its expression was reduced in D. antarctica cold-deacclimated plants. This could be related to the possible existing interspecific differences, but also reaffirms the postulate that cold deacclimation is more than a passive regression from a previous cold-acclimated state [21].
Both tissue dehydration and oxidative stress can be the result of low temperature exposure, since the freezing of apoplastic fluid results in severe cell dehydration, the formation of ice crystals that cause cell membrane rupture [51], and the release of electrolytes and free radicals [52]. In this sense, the synthesis of flavonols, flavonoids, and phenylpropanoids has been proposed in D. antarctica as a way to increase the antioxidant capacity of cold-acclimated plants [39]. Our results confirm that ‘flavonol 4′-sulfotransferase’ and various ‘UDP-glycosyltransferase’ genes (UGT73B3, UGT73C1, UGT73C5, UGT91C1) that encode enzymes related to these antioxidant molecules synthesis pathways [53,54,55,56,57] were co-expressed together in freezing tolerant plants (CA and DW+), while their expression was down-regulated during cold deacclimation as result of nocturnal warming.
Another gene whose expression is modified by the cold is ‘E3 ubiquitin-protein ligase RGLG5’ (RGLG5), which together with RGLG1, is responsible for phosphatase 2CA (PP2CA) ubiquitination and its degradation [58], and PP2CA is a recognized repressor of the hormonal abscisic acid (ABA) synthesis pathway [58,59]. Our results show that RGLG5 transcripts were down-regulated in cold-deacclimated plants, while their expression was higher and stable in CA plants, and it gradually increased overnight in DW+ plants. This could imply a reduction in D. antarctica leaf tissue ABA concentration in cold-deacclimated plants by nocturnal warming, while the ABA concentration of the leaf tissue from DW+ plants could be similar to that of cold-acclimated plants at the end of the night. However, more studies are needed to fully validate this idea. In accordance with the previous approach, although the CBF-regulon expression corresponds to an ABA-independent pathway in response to the cold [60], its expression, as well the expression of other cold-regulated genes (COR), requires ABA presence to potentiate its induction by the cold [61]. This has been related to the presence in the promoter regions of the freezing tolerance genes of the CRT/DRE elements, and the sequences called ‘abscisic acid-responsive element’ (ABRE) [60].
The CBFs’ expression is not only related to the induction of freeze tolerance genes, but also to the repression of a group of genes related to plant growth [24]. Arabidopsis thaliana plants overexpressing one or more CBF genes exhibit constitutive freezing tolerance, slow growth, and a dwarf phenotype [62]. In addition, rice plants overexpressing the DaCBF4 gene presented a greater cold and drought tolerance than wild type [63], favoring the up-regulation of dehydrins and other LEA proteins, while down-regulating the expression of several proteins related to photosynthesis such as Chlorophyll A-B binding protein, 10 kDa chloroplast precursor polypeptide of photosystem II, 44 kDa reaction center protein of photosystem II, and the D2 protein of photosystem II [63]. Similar results were found in our study of plants that maintained freezing tolerance compared to those that were cold-deacclimated by nocturnal warming. This could suggest that DaCBF4 and other DaCBFs directly repress or promote the transcription of other repressors over a photosynthesis-related group of genes, given the functional redundancy found among CBFs genes [23,24].
Even though the transcription of genes related to photosynthetic activity does not necessarily imply a higher photosynthetic rate, it has been reported that night-warmed plants exhibited an increase in CO2 assimilation during the day [64,65]. Additionally, the increased expression of carbohydrate metabolism-related genes, such as the triose phosphate/phosphate translocator and Fructose-bisphosphate aldolase 5 in plants with nocturnal warming, suggests that a remobilization process of previously accumulated photo-assimilates occurs at night to obtain energy and/or to be used as intermediates in biosynthetic processes, as proposed by Turnbull et al. [64]. This nocturnal remobilization of photo-assimilates does not seem to be directed towards sucrose synthesis, since the expression of sucrose phosphate synthase transcripts did not vary throughout the night in plants with nocturnal warming. Previous results in D. antarctica plants under experimental laboratory conditions with nocturnal warming show a reduction in the ratio between respiration and carbon assimilation, mainly due to a greater increase in CO2 assimilation [65], but this increase in CO2 assimilation is not reflected in the accumulation of soluble carbohydrates, since there is also a reduction in the amount of sucrose accumulated [14]. All the above suggests that the accumulated carbohydrates and the recently assimilated CO2 were used in other biological processes, potentially including vegetative growth. This idea is consistent with the increased expression of genes related to the hormonal response in plants with nocturnal warming, such as ‘IAA-amino acid hydrolase ILR1-like 2’, which allows for the release of previously accumulated auxins in its active form, promoting plant growth [66,67].
The relationship between the consumption of accumulated soluble carbohydrates, with the activation of photosynthetic metabolism and vegetative growth, has been previously associated with the cold deacclimation process [16,21,68]. Usually, cold deacclimation occurs during the spring season, and for some species, it could be a passive process [17]. However, D. antarctica plants should maintain their freezing tolerance for the whole year, since they need to deal with freezing temperatures during the growing season (summer) [11,12]. Therefore, when cold deacclimation occurs in D. antarctica (at least 7 °C increase in LT50), it was promoted by nocturnal warming [14], and this should be an active process, as was also proposed by Wójcik-Jagła et al. [17] for winter barley submitted to warm periods during the winter season. This active cold deacclimation process involves a complex gene regulation that was activated by nocturnal warming only, as has been shown in this laboratory study. Even D. antarctica plants cold deacclimated by nocturnal warming maintain a certain freezing resistance that prevents serious damage by the Antarctic spring and summer freezing events [14]. However, it is not possible to predict the effect that the vegetative growth promotion by nocturnal warming should have on other physiological processes, such as reproduction and seed maturity. While most of the researchers found that nocturnal warming reduces the reproductive success of cereals [69], others report opposite results [70]. Moreover, D. antarctica plants in the field are exposed to other variables beside temperature that can influence their physiological response, such as wind intensity, variable frequency, nutrient availability, and the intensity of freezing events, among others. Therefore, this raises questions about the relative contribution of nocturnal warming to induce the cold deacclimation process of D. antarctica plants under field conditions.

4. Materials and Methods

4.1. Plant Material and Warming Treatments

Deschampsia antarctica Desv. plants were collected near the Henryk Arctowski Polish Scientific Station, King George Island, Maritime Antarctic (62°09′41″ S–58°28′10″ W) during the 2018 growing season. Three individual plants (biological replicates) were vegetatively propagated, as previously described by Bravo et al. [71], in 400 mL pots with soil mixture 3:2:1 (vegetal soil:vermiculite:peat), with 18 h/6 h of light/dark photoperiod, in a controlled climate room at 13 °C constant temperature. These plants (four clones for each biological replicate) were cold-acclimated during 28 days with a gradual temperature reduction in the 18 h/6 h day/night thermoperiod, as described by López et al. [14]. Then, groups of pots were transferred to four thermoperiods treatments that were established with the following temperature regimes: (CA) maintained at 8 °C/0 °C day/night, the control treatment where plants remain in the final cold acclimation condition; (NW+) only nocturnal warming at 8 °C/6 °C day/night; (DW+) only diurnal warming at 14 °C/0 °C day/night; and (DNW+), both diurnal and nocturnal warming at 14 °C/6 °C day/night. Each treatment contained clones of the 3 biological replicates. The maximum and minimum temperatures of the above-mentioned thermoperiods were selected based on microenvironmental data recorded in the field, near the Henryk Arctowski Polish Antarctic Station, King George Island. It has been reported that an increase of 6 °C in the minimum nocturnal temperature leads to the cold deacclimation of D. antarctica, unlike the same increase in the maximum daytime temperature [14]. The plants remained in these conditions for 14 days.
During the fourteenth night, the plants were kept at the maximum temperature of their corresponding thermoperiod until the end of the diurnal photoperiod. Samples were collected from 3 plants for each treatment during the night. First, just before the nocturnal period began, which was considered time 0 h, the initial state before the beginning of the night and the temperature dropped. Next, the plants were moved into dark culture chambers with their corresponding nocturnal thermoperiod of each treatment (0 °C CA and DW+; 6 °C NW+ and DNW+). After 2 h, when the expression of the CBF genes should be the maximum [26], the same plants were sampled again, and also at 6 h, just before the diurnal photoperiod began and the temperature rose, in order to analyze the effect of nocturnal warming on CBFs and freezing tolerant-related genes.

4.2. Total RNA Extraction and Quantification

Total RNA extraction and purification was performed using a RNeasy ®® Plant Mini Kit (QIAGEN), according to the manufacturer’s recommendations. Total RNA concentration was determined using an Infinite M200 NanoQuant spectrophotometer measured at 260 nm, and the A260/A280 absorbance ratio was used to estimate its purity.
Additionally, from extracted RNA, cDNA synthesis was performed with the commercial AffinityScript qRT-PCR kit (Stratagene, Cedar Creek, TX, USA), following the manufacturer’s instructions. Reverse transcription was performed from 1 μg of total RNA, using oligos-dT and random primers as primers in a final volume of 20 μL. At the end of the reaction, the cDNA in solution was diluted in a 1:10 ratio with nuclease-free water, before being aliquoted and stored at −80 °C until use.

4.3. Sequencing and De Novo Transcriptome Assembly

A total of twelve samples were sequenced in a Novaseq 6000 platform, for 150 cycles in the paired-end mode at Novogene facilities. The resulting fastq files were processed using Trimmomatic v0.38 [72] to remove residual adapters and low-quality sequences. These high-quality reads were used to construct a de novo assembly of a Deschampsia antarctica transcriptome using Trinity v2.12 [73].

4.4. Transcriptome Assessment and Abundance Estimation

The raw assembly was clustered at 95% similarity using CD-HIT [74], and the completeness was evaluated using BUSCO v5 (Simao et al., 2015), which queried a collection of proteins (Embryophyta odb10) to the assembled transcripts by similarity. The relative abundance estimation was calculated using RSEM v1.2.26 [75].

4.5. Differential Gene Expression and Co-Expression Analysis

Changes in relative RNA abundance were analyzed with Bioconductor package DESeq2 [76] in the R statistical environment. The significance of the changes in gene expression were judged using a False Discovery Rate (FDR) of less than 0.05 and a minimum fold change (FC) of 2 as the threshold. The expression patterns in the differentially expressed genes were clustered using dimensionality reduction with Principal Component Analysis (PCA), followed by a 8 × 8 Self-Organizing Map (SOM) using one hundred training iterations, as implemented in the Kohonen package in R [77]. The resulting SOM nodes were grouped using a hierarchical clustering approach with a branch length cutoff of 2.2, resulting in 9 super-nodes. The abundance patterns of each super-node were visualized in a boxplot to associate their accumulation with any treatment.

4.6. Transcriptome Annotation and Enrichment Test

The resulting genes were compared by homology into the UniProt/SwissProtKB database using BLAST+ with an e-value of 1 × 10−10 as the threshold. The ontology assignment was performed using the PANTHER (version 2020_4) classification system [78], from EMBL. Each super-node was analyzed for the enrichment of gene ontologies with a false discovery rate (FDR) of 0.05 using GOSeq [79].

4.7. Quantitative Real-Time (qRT-PCR) Analysis

The expression patterns of DEGs were analyzed by qRT-PCR, performed in a Stratagene AriaMX thermocycler (Agilent, Santa Clara, CA, USA) using the previously obtained cDNA and Brilliant II SYBR Green QPCR Master mix reagent (Stratagene, Cedar Creek, TX, USA), according to the manufacturer’s instructions. The gene-specific primers (Table S4) used for the qRT-PCR analysis were designed using the Primer-BLAST tool (http://www.ncbi.nlm.nih.gov/tools/primer-blast accessed on 12 October 2022). The genes Polyadenylate-binding protein2 (PAB2) and Eukaryotic translation initiation factor (IF4E4) were used as internal references to normalize the expression data. The relative expression levels were calculated according to the method proposed by Pfaffl [80], based on the gene primer efficiency of the e−ΔΔCT (cycle threshold). To validate the repeatability of the RNA-seq data, 16 DEGs were selected for verification by qRT-PCR, using the cDNA obtained from the same sample that was sent to the sequence. The relative expression was estimated with respect to their values in the control treatment (CA), 2 h into the nocturnal period. In contrast, the kinetic variation of genes expression during the night was evaluated with respect to their values in the control treatment (CA) at 0 h just before the nocturnal temperature began.

5. Conclusions

The evidence obtained in the present study confirms that nocturnal warming is key in the cold deacclimation process through the down-regulation of freezing tolerance-related genes, and the up-regulation of photosynthetic and growth-related genes in D. antarctica. Furthermore, the expression of the CBF transcription factor and its regulon during the night is necessary to maintain the previously acquired freezing tolerance. However, more studies are necessary to understand the consequences that regional warming will have on D. antarctica field populations. In this sense, it would be interesting to evaluate how gene expression is affected by nocturnal warming under field experimental conditions.

Supplementary Materials

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

Author Contributions

Conceptualization, D.L., P.L.S. and L.A.B.; methodology, D.L., G.L. and L.A.B.; investigation, D.L. and L.A.B.; formal analysis, D.L., G.L. and L.A.B.; resources, D.L. and L.A.B.; writing—original draft preparation and editing, D.L.; writing—review and editing, G.L., P.L.S. and L.A.B.; visualization, D.L. and G.L.; project administration, L.A.B.; funding acquisition, D.L., P.L.S. and L.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by ANID through Fondecyt 1151173 and Scholarship Program/Beca Doctorado Nacional/2017—21170470, NEXER-UFRO NXR17-0002, and INACH RT 18-18.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequencing data was deposited in the National Center for Biotechnology Information (NCBI) and can be accessed via BioProject ID PRJNA941125.

Acknowledgments

The authors thank the 42nd crew of the Henryk Arctowski Station for their tremendous help and logistical support, and the Instituto Antártico Chileno (INACH) for the support provided (including the permits to work and sample close to ASPA 128) during the 54th Antarctic Scientific Expedition (ECA). The authors acknowledge the supercomputing infrastructure of Soroban (SATREPS MACH—JPM/JSA1705) at Centro de Modelación y Computación Científica, Universidad de La Frontera. The authors thank Charles L. Guy for proofreading this manuscript, and Fernanda Meneses−Ponzini for the graphical abstract drawing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meehl, G.A.; Stocker, T.F.; Collins, W.D.; Friedlingstein, P.; Gaye, A.T.; Gregory, J.M.; Kitoh, A.; Knutti, R.; Murphy, J.M.; Noda, A.; et al. Global Climate Projections Coordinating. In Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K.B., Tignor, M., Miller, H.L., Eds.; Cambridge University Press: Cambridge, UK, 2007; pp. 747–845. ISBN 978-0521-70596-7. [Google Scholar]
  2. Pachauri, R.K.; Allen, M.R.; Barros, V.R.; Broome, J.; Cramer, W.; Christ, R.; Church, J.A.; Clarke, L.; Dahe, Q.; Dasgupta, P.; et al. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  3. Vose, R.S.; Easterling, D.R.; Gleason, B. Maximum and minimum temperature trends for the globe: An update through 2004. Geophys. Res. Lett. 2005, 32, L23822. [Google Scholar] [CrossRef] [Green Version]
  4. Xie, H.; Id, M.Z.; Yu, Y.; Zeng, X.; Tang, G.; Duan, Y.; Wang, J.; Yu, Y. Comparative transcriptome analysis of the cold resistance of the sterile rice line 33S. PLoS ONE 2022, 17, e0261822. [Google Scholar] [CrossRef]
  5. Vaughan, D.G.; Marshall, G.J.; Connolley, W.M.; Parkinson, C.; Mulvaney, R.; Hodgson, D.A.; King, J.C.; Pudsey, C.J.; Turner, J. Recent Rapid Regional Climate Warming on the Antarctic Peninsula. Clim. Change 2003, 60, 243–274. [Google Scholar] [CrossRef]
  6. Fowbert, J.A.; Smith, R.I.L. Rapid population increases in native vascular plants in the Argentine Islands, Antarctic Peninsula. Arct. Alp. Res. 1994, 26, 290–296. [Google Scholar] [CrossRef]
  7. Smith, R.I.L. The enigma of Colobanthus quitensis and Deschampsia antarctica in Antarctica. In Antarctic Biology in a Global Context; Huiskes, A.H.L., Gieskes, W.W.C., Rozema, J., Schorno, R.M.L., van der Vies, S.M., Eds.; Backhuys: Leiden, The Netherlands, 2003; pp. 234–239. [Google Scholar]
  8. Torres-Mellado, G.A.; Jaña, R.; Casanova-Katny, M.A. Antarctic hairgrass expansion in the South Shetland archipelago and Antarctic Peninsula revisited. Polar Biol. 2011, 34, 1679–1688. [Google Scholar] [CrossRef]
  9. Cannone, N.; Guglielmin, M.; Convey, P.; Worland, M.R.; Favero Longo, S.E. Vascular plant changes in extreme environments: Effects of multiple drivers. Clim. Change 2016, 134, 651–665. [Google Scholar] [CrossRef] [Green Version]
  10. Robinson, S.A.; Wasley, J.; Tobin, A.K. Living on the edge—Plants and global change in continental and maritime Antarctica. Glob. Chang. Biol. 2003, 9, 1681–1717. [Google Scholar] [CrossRef] [Green Version]
  11. Convey, P. The influence of environmental characteristics on life history attributes of Antarctic terrestrial biota. Biol. Rev. 1996, 71, 191–225. [Google Scholar] [CrossRef]
  12. Convey, P.; Chown, S.L.; Clarke, A.; Barnes, D.K.A.; Bokhorst, S.; Cummings, V.; Ducklow, H.W.; Frati, F.; Green, T.G.A.; Gordon, S.; et al. The spatial structure of antarctic biodiversity. Ecol. Monogr. 2014, 84, 203–244. [Google Scholar] [CrossRef] [Green Version]
  13. Day, T.A.; Ruhland, C.T.; Grobe, C.W.; Xiong, F. Growth and reproduction of Antarctic vascular plants in response to warming and UV radiation reductions in the field. Oecologia 1999, 119, 24–35. [Google Scholar] [CrossRef]
  14. López, D.; Sanhueza, C.; Salvo-Garrido, H.; Bascunan-Godoy, L.; Bravo, L.A. How Does Diurnal and Nocturnal Warming Affect the Freezing Resistance of Antarctic Vascular Plants? Plants 2023, 12, 806. [Google Scholar] [CrossRef] [PubMed]
  15. Vyse, K.; Pagter, M.; Zuther, E.; Hincha, D.K. Deacclimation after cold acclimation- a crucial, but widely neglected part of plant winter survival. J. Exp. Bot. 2019, 70, 4595–4604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zuther, E.; Juszczak, I.; Ping Lee, Y.; Baier, M.; Hincha, D.K. Time-dependent deacclimation after cold acclimation in Arabidopsis thaliana accessions. Sci. Rep. 2015, 5, 12199. [Google Scholar] [CrossRef] [Green Version]
  17. Wójcik-Jagła, M.; Daszkowska-Golec, A.; Fiust, A.; Kopeć, P.; Rapacz, M. Identification of the genetic basis of response to de-acclimation in winter barley. Int. J. Mol. Sci. 2021, 22, 1057. [Google Scholar] [CrossRef] [PubMed]
  18. Pagter, M.; Alpers, J.; Erban, A.; Kopka, J.; Zuther, E.; Hincha, D.K. Rapid transcriptional and metabolic regulation of the deacclimation process in cold acclimated Arabidopsis thaliana. BMC Genom. 2017, 18, 731. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Horvath, D.P.; Zhang, J.; Chao, W.S.; Mandal, A.; Rahman, M.; Anderson, J.V. Genome-wide association studies and transcriptome changes during acclimation and deacclimation in divergent Brassica napus varieties. Int. J. Mol. Sci. 2020, 21, 9148. [Google Scholar] [CrossRef]
  20. Rathore, N.; Kumar, P.; Mehta, N.; Swarnkar, M.K.; Shankar, R.; Chawla, A. Time-series RNA-Seq transcriptome profiling reveals novel insights about cold acclimation and de-acclimation processes in an evergreen shrub of high altitude. Sci. Rep. 2022, 12, 15553. [Google Scholar] [CrossRef]
  21. Pagter, M.; Arora, R. Winter survival and deacclimation of perennials under warming climate: Physiological perspectives. Physiol. Plant. 2012, 147, 75–87. [Google Scholar] [CrossRef]
  22. Fowler, S.; Thomashow, M.F. Arabidopsis Transcriptome Profiling Indicates That Multiple Regulatory Pathways Are Activated during Cold Acclimation in Addition to the CBF Cold Response Pathway. Plant Cell 2002, 14, 1675–1690. [Google Scholar] [CrossRef] [Green Version]
  23. Jia, Y.; Ding, Y.; Shi, Y.; Zhang, X.; Gong, Z.; Yang, S. The cbfs triple mutants reveal the essential functions of CBFs in cold acclimation and allow the definition of CBF regulons in Arabidopsis. New Phytol. 2016, 212, 345–353. [Google Scholar] [CrossRef] [Green Version]
  24. Zhao, C.; Zhang, Z.; Xie, S.; Si, T.; Li, Y.; Zhu, J.-K. Mutational Evidence for the Critical Role of CBF Transcription Factors in Cold Acclimation in Arabidopsis. Plant Physiol. 2016, 171, 2744–2759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Sandve, S.R.; Kosmala, A.; Rudi, H.; Fjellheim, S.; Rapacz, M.; Yamada, T.; Rognli, O.A. Molecular mechanisms underlying frost tolerance in perennial grasses adapted to cold climates. Plant Sci. 2011, 180, 69–77. [Google Scholar] [CrossRef] [PubMed]
  26. Thomashow, M.F. Molecular Basis of Plant Cold Acclimation: Insights Gained from Studying the CBF Cold Response Pathway. Plant Physiol. 2010, 154, 571–577. [Google Scholar] [CrossRef] [Green Version]
  27. Badawi, M.; Danyluk, J.; Boucho, B.; Houde, M.; Sarhan, F. The CBF gene family in hexaploid wheat and its relationship to the phylogenetic complexity of cereal CBFs. Mol. Genet. Genom. 2007, 277, 533–554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kalberer, S.R.; Wisniewski, M.; Arora, R. Deacclimation and reacclimation of cold-hardy plants: Current understanding and emerging concepts. Plant Sci. 2006, 171, 3–16. [Google Scholar] [CrossRef]
  29. Campoli, C.; Matus-Cádiz, M.A.; Pozniak, C.J.; Cattivelli, L.; Fowler, D.B. Comparative expression of Cbf genes in the Triticeae under different acclimation induction temperatures. Mol. Genet. Genom. 2009, 282, 141–152. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, Z.; Jia, Y.; Ding, Y.; Shi, Y.; Li, Z.; Guo, Y.; Gong, Z.; Yang, S. Plasma Membrane CRPK1-Mediated Phosphorylation of 14-3-3 Proteins Induces Their Nuclear Import to Fine-Tune CBF Signaling during Cold Response. Mol. Cell 2017, 66, 117–128.e5. [Google Scholar] [CrossRef] [Green Version]
  31. Catalá, R.; Santos, E.; Alonso, J.M.; Ecker, J.R.; Martínez-Zapater, J.M.; Salinas, J. Mutations in the Ca2+/H+ Transporter CAX1 Increase CBF/DREB1 Expression and the Cold-Acclimation Response in Arabidopsis. Plant Cell 2003, 15, 2940–2951. [Google Scholar] [CrossRef] [Green Version]
  32. Iqbal, Z.; Memon, A.G.; Ahmad, A.; Iqbal, M.S. Calcium Mediated Cold Acclimation in Plants: Underlying Signaling and Molecular Mechanisms. Front. Plant Sci. 2022, 13, 855559. [Google Scholar] [CrossRef]
  33. Tong, T.; Li, Q.; Jiang, W.; Chen, G.; Xue, D.; Deng, F.; Zeng, F.; Chen, Z.-H. Molecular Evolution of Calcium Signaling and Transport in Plant Adaptation to Abiotic Stress. Int. J. Mol. Sci. 2021, 22, 12308. [Google Scholar] [CrossRef]
  34. Zhang, W.; Wang, J.; Huang, Z.; Mi, L.; Xu, K.; Wu, J.; Fan, Y.; Ma, S.; Jiang, D. Effects of Low Temperature at Booting Stage on Sucrose Metabolism and Endogenous Hormone Contents in Winter Wheat Spikelet. Front. Plant Sci. 2019, 10, 498. [Google Scholar] [CrossRef] [PubMed]
  35. Stein, O.; Granot, D. An Overview of Sucrose Synthases in Plants. Front. Plant Sci. 2019, 10, 95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Bilska-Kos, A.; Mytych, J.; Suski, S.; Magoń, J.; Ochodzki, P.; Zebrowski, J. Sucrose phosphate synthase (SPS), sucrose synthase (SUS) and their products in the leaves of Miscanthus × giganteus and Zea mays at low temperature. Planta 2020, 252, 23. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, A.; Kao, M.-H.; Yang, W.; Sayion, Y.; Liu, L.-F.; Lee, P.-D.; Jong-Ching, S. Differentially and Developmentally Regulated Expression of Three Rice Sucrose Synthase Genes. Plant Cell Physiol. 1999, 40, 800–807. [Google Scholar] [CrossRef] [PubMed]
  38. Nolte, K.D.; Koch, K.E. Companion-Cell Specific Localization of Sucrose Synthase in Zones of Phloem Loading and Unloading. Plant Physiol. 1993, 101, 899–905. [Google Scholar] [CrossRef] [Green Version]
  39. Clemente-Moreno, M.J.; Omranian, N.; Sáez, P.L.; Figueroa, C.M.; Del-Saz, N.; Elso, M.; Poblete, L.; Orf, I.; Cuadros-Inostroza, A.; Cavieres, L.A.; et al. Low-temperature tolerance of the Antarctic species Deschampsia antarctica: A complex metabolic response associated with nutrient remobilization. Plant Cell Environ. 2020, 43, 1376–1393. [Google Scholar] [CrossRef]
  40. Borovik, O.A.; Pomortsev, A.V.; Korsukova, A.V.; Polyakova, E.A.; Fomina, E.A.; Zabanova, N.S.; Grabelnych, O.I. Effect of Cold Acclimation and Deacclimation on the Content of Soluble Carbohydrates and Dehydrins in the Leaves of Winter Wheat. J. Stress Physiol. Biochem. 2019, 15, 62–67. [Google Scholar]
  41. Hoffman, L.; DaCosta, M.; Scott Ebdon, J. Examination of cold deacclimation sensitivity of annual bluegrass and creeping bentgrass. Crop Sci. 2014, 54, 413–420. [Google Scholar] [CrossRef] [Green Version]
  42. Vítámvás, P.; Prášil, I.T. WCS120 protein family and frost tolerance during cold acclimation, deacclimation and reacclimation of winter wheat. Plant Physiol. Biochem. 2008, 46, 970–976. [Google Scholar] [CrossRef]
  43. Zhang, X.; Wang, K.; Ervin, E.H.; Waltz, C.; Murphy, T. Metabolic changes during cold acclimation and deacclimation in five bermudagrass varieties. I. Proline, total amino acid, protein, and dehydrin expression. Crop Sci. 2011, 51, 838–846. [Google Scholar] [CrossRef]
  44. Banerjee, A.; Roychoudhury, A. Group II late embryogenesis abundant (LEA) proteins: Structural and functional aspects in plant abiotic stress. Plant Growth Regul. 2016, 79, 1–17. [Google Scholar] [CrossRef]
  45. Kaplan, F.; Kopka, J.; Sung, D.Y.; Zhao, W.; Popp, M.; Porat, R.; Guy, C.L. Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content. Plant J. 2007, 50, 967–981. [Google Scholar] [CrossRef] [PubMed]
  46. Hoermiller, I.I.; Funck, D.; Schönewolf, L.; May, H.; Heyer, A.G. Cytosolic proline is required for basal freezing tolerance in Arabidopsis. Plant Cell Environ. 2022, 45, 147–155. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, Y.; Jiang, J.; Chang, Q.; Gu, C.; Song, A.; Chen, S.; Dong, B.; Chen, F. Cold acclimation induces freezing tolerance via antioxidative enzymes, proline metabolism and gene expression changes in two chrysanthemum species. Mol. Biol. Rep. 2014, 41, 815–822. [Google Scholar] [CrossRef] [PubMed]
  48. Hoffman, L.; DaCosta, M.; Ebdon, J.S.; Watkins, E. Physiological changes during cold acclimation of perennial ryegrass accessions differing in freeze tolerance. Crop Sci. 2010, 50, 1037–1047. [Google Scholar] [CrossRef]
  49. Anwar, A.; She, M.; Wang, K.; Riaz, B.; Ye, X. Biological Roles of Ornithine Aminotransferase (OAT) in Plant Stress Tolerance: Present Progress and Future Perspectives. Int. J. Mol. Sci. 2018, 19, 3681. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. You, J.; Hu, H.; Xiong, L. Plant Science An ornithine δ-aminotransferase gene OsOAT confers drought and oxidative stress tolerance in rice. Plant Sci. 2012, 197, 59–69. [Google Scholar] [CrossRef]
  51. Steponkus, P.L. Role of the plasma membrane in freezing injury and cold acclimation. Annu. Rev. Plant Physiol. 1984, 35, 543–584. [Google Scholar] [CrossRef]
  52. Pearce, R.S. Plant freezing and damage. Ann. Bot. 2001, 87, 417–424. [Google Scholar] [CrossRef]
  53. Langlois-Meurinne, M.; Gachon, C.M.M.; Saindrenan, P. Pathogen-Responsive Expression of Glycosyltransferase Genes UGT73B3 and UGT73B5 Is Necessary for Resistance to Pseudomonas syringae pv tomato. Plant Physiol. 2005, 139, 1890–1901. [Google Scholar] [CrossRef] [Green Version]
  54. Mostek, A.; Börner, A.; Weidner, S. Plant Physiology and Biochemistry Comparative proteomic analysis of b-aminobutyric acid-mediated alleviation of salt stress in barley. Plant Physiol. Biochem. 2016, 99, 150–161. [Google Scholar] [CrossRef]
  55. Pan, L.; Gao, H.; Xia, W.; Zhang, T.; Dong, L. Establishing a herbicide-metabolizing enzyme library in Beckmannia syzigachne to identify genes associated with metabolic resistance. J. Exp. Bot. 2016, 67, 1745–1757. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, J.; Jia, J.; Sun, J.; Pang, X.; Li, B.; Yuan, J.; Chen, E.; Li, X. Trypsin preservation: CsUGT91C1 regulates Trilobatin Biosynthesis in Cucumis sativus during Storage. Plant Growth Regul. 2023. [Google Scholar] [CrossRef]
  57. Poppenberger, B.; Fujioka, S.; Soeno, K.; George, G.L.; Seto, H.; Takatsuto, S.; Adam, G.; Yoshida, S.; Bowles, D. The UGT73C5 of Arabidopsis thaliana glucosylates brassinosteroids. Proc. Natl. Acad. Sci. USA 2005, 102, 15253–15258. [Google Scholar] [CrossRef] [PubMed]
  58. Wu, Q.; Zhang, X.; Peirats-Llobet, M.; Belda-Palazon, B.; Wang, X.; Cui, S.; Yu, X.; Rodriguez, P.L.; Chengcai, A. Ubiquitin Ligases RGLG1 and RGLG5 Regulate Abscisic Acid Signaling by Controlling the Turnover of Phosphatase PP2CA. Plant Cell 2016, 28, 2178–2196. [Google Scholar] [CrossRef] [Green Version]
  59. Tähtiharju, S.; Palva, T. Antisense inhibition of protein phosphatase 2C accelerates cold acclimation in Arabidopsis thaliana. Plant J. 2001, 26, 461–470. [Google Scholar] [CrossRef]
  60. Thomashow, M.F. PLANT COLD ACCLIMATION: Freezing Tolerance Genes and Regulatory Mechanisms. Annu. Rev. Plant Physiol. Plant Mol. Biol. 1999, 50, 571–599. [Google Scholar] [CrossRef] [Green Version]
  61. Xiong, L.; Ishitani, M.; Lee, H.; Zhu, J.-K. The Arabidopsis LOS5/ABA3 locus encodes a molybdenum cofactor sulfurase and modulates cold stress- and osmotic stress- responsive gene expression. Plant Cell 2001, 13, 2063–2083. [Google Scholar]
  62. Gilmour, S.J.; Fowler, S.G.; Thomashow, M.F. Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol. Biol. 2004, 54, 767–781. [Google Scholar] [CrossRef]
  63. Byun, M.Y.; Cui, L.H.; Lee, J.; Park, H.; Lee, A.; Kim, W.T.; Lee, H. Identification of Rice Genes Associated With Enhanced Cold Tolerance by Comparative Transcriptome Analysis With Two Transgenic Rice Plants Isolated From Antarctic Flowering Plant Deschampsia antarctica. Front. Plant Sci. 2018, 9, 601. [Google Scholar] [CrossRef] [Green Version]
  64. Turnbull, M.H.; Murthy, R.; Griffin, K.L. The relative impacts of daytime and night-time warming on photosynthetic capacity in Populus deltoides. Plant Cell Environ. 2002, 25, 1729–1737. [Google Scholar] [CrossRef] [Green Version]
  65. Sanhueza, C.; Fuentes, F.; Cortes, D.; Bascuñan, L.; Saez, P.; Bravo, L.A.; Cavieres, L. Contrasting thermal acclimation of leaf dark respiration and photosynthesis of Antarctic vascular plant species exposed to nocturnal warming. Physiol. Plant. 2019, 167, 205–216. [Google Scholar] [CrossRef]
  66. Korasick, D.A.; Enders, T.A.; Strader, L.C. Auxin biosynthesis and storage forms. J. Exp. Bot. 2013, 64, 2541–2555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. LeClere, S.; Tellez, R.; Rampey, R.A.; Matsuda, S.P.T.; Bartel, B. Characterization of a Family of IAA-Amino Acid Conjugate Hydrolases from Arabidopsis. J. Biol. Chem. 2002, 277, 20446–20452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Wan, S.; Xia, J.; Liu, W.; Niu, S. Photosynthetic overcompensation under nocturnal warming enhances grassland carbon sequestration. Ecology 2009, 90, 2700–2710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Sadok, W.; Jagadish, S.V.K. The Hidden Costs of Nighttime Warming on Yields. Trends Plant Sci. 2020, 25, 644–651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Fan, Y.; Lv, Z.; Qin, B.; Yang, J.; Ren, K.; Liu, Q.; Jiang, F.; Zhang, W.; Ma, S.; Ma, C.; et al. Night warming at the vegetative stage improves pre-anthesis photosynthesis and plant productivity involved in grain yield of winter wheat. Plant Physiol. Biochem. 2022, 186, 19–30. [Google Scholar] [CrossRef]
  71. Bravo, L.A.; Ulloa, N.; Zuñiga, G.E.; Casanova, A.; Corcuera, L.J.; Alberdi, M. Cold resistance in antarctic angiosperms. Physiol. Plant. 2001, 111, 55–65. [Google Scholar] [CrossRef]
  72. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  73. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Chen, Z.; et al. Trinity: Reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 2013, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  74. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef] [PubMed]
  75. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Wehrens, R.; Buydens, L.M.C. Self- and super-organizing maps in R: The kohonen package. J. Stat. Softw. 2007, 21, 1–19. [Google Scholar] [CrossRef] [Green Version]
  78. Mi, H.; Muruganujan, A.; Huang, X.; Ebert, D.; Mills, C.; Guo, X.; Thomas, P.D. Protocol Update for Large-scale genome and gene function analysis with PANTHER Classification System (v.14.0). Nat. Protoc. 2019, 14, 703–721. [Google Scholar] [CrossRef]
  79. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [Green Version]
  80. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001, 29, 16–21. [Google Scholar] [CrossRef]
Ijms 24 11211 g001 550
Figure 1. Number of differentially expressed genes (DEGs). The amounts of differentially expressed genes in warming treatments, up-regulated (light yellow) and down-regulated (orange) in relation to the control condition of cold-acclimated plants, 8 °C/0 °C (CA). Warming treatments received the following thermoperiods: 8 °C/6 °C, nocturnal warming (NW+); 14 °C/0 °C, diurnal warming (DW+); and 14 °C/6 °C, diurnal-nocturnal warming (DNW+). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Figure 1. Number of differentially expressed genes (DEGs). The amounts of differentially expressed genes in warming treatments, up-regulated (light yellow) and down-regulated (orange) in relation to the control condition of cold-acclimated plants, 8 °C/0 °C (CA). Warming treatments received the following thermoperiods: 8 °C/6 °C, nocturnal warming (NW+); 14 °C/0 °C, diurnal warming (DW+); and 14 °C/6 °C, diurnal-nocturnal warming (DNW+). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Ijms 24 11211 g001
Ijms 24 11211 g002 550
Figure 2. Venn diagram of the number of differentially expressed genes (DEGs). The number and percentage of down-regulated genes (A) and up-regulated genes (B) in warming treatments in relation to the control condition of cold-acclimated plants, 8 °C/0 °C (CA). Warming treatments received the following thermoperiods: 8 °C/6 °C, nocturnal warming (NW+); 14 °C/0 °C, diurnal warming (DW+); and 14 °C/6 °C, diurnal-nocturnal warming (DNW+). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Figure 2. Venn diagram of the number of differentially expressed genes (DEGs). The number and percentage of down-regulated genes (A) and up-regulated genes (B) in warming treatments in relation to the control condition of cold-acclimated plants, 8 °C/0 °C (CA). Warming treatments received the following thermoperiods: 8 °C/6 °C, nocturnal warming (NW+); 14 °C/0 °C, diurnal warming (DW+); and 14 °C/6 °C, diurnal-nocturnal warming (DNW+). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Ijms 24 11211 g002
Ijms 24 11211 g003 550
Figure 3. Clusters obtained from principal component analysis and self-organizing maps (PCA-SOM). The nodes group transcripts up-regulated (A) and down-regulated (B) by nocturnal warming, and the enriched gene ontology categories with their corresponding p-values are shown. The boxes represent the interquartile range (middle 50% of the data), its midline indicates the median value, and the whiskers indicate the ranges of the upper and lower 25%, excluding outliers. Corresponding to treatments with the following thermoperiods (day/night): 8 °C/0 °C, plants acclimated to cold (AF); 8 °C/6 °C, night heating (CN); 14 °C/0 °C, daytime warming (CD); and 14 °C/6 °C, daytime-night warming (CDN). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Figure 3. Clusters obtained from principal component analysis and self-organizing maps (PCA-SOM). The nodes group transcripts up-regulated (A) and down-regulated (B) by nocturnal warming, and the enriched gene ontology categories with their corresponding p-values are shown. The boxes represent the interquartile range (middle 50% of the data), its midline indicates the median value, and the whiskers indicate the ranges of the upper and lower 25%, excluding outliers. Corresponding to treatments with the following thermoperiods (day/night): 8 °C/0 °C, plants acclimated to cold (AF); 8 °C/6 °C, night heating (CN); 14 °C/0 °C, daytime warming (CD); and 14 °C/6 °C, daytime-night warming (CDN). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Ijms 24 11211 g003
Ijms 24 11211 g004 550
Figure 4. Heat-map of co-expressed genes’ response to nocturnal warming. The genes’ relative abundance is presented by the z-score of the log10 FPKM, after 2 h of temperature drop during the nocturnal period, of plants with a cold acclimation of 8 °C/0 °C (CA); nocturnal warming of 8 °C/6 °C (NW+); diurnal warming of 14 °C/0 °C (DW+); and diurnal-nocturnal warming of 14 °C/6 °C (DNW+). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Figure 4. Heat-map of co-expressed genes’ response to nocturnal warming. The genes’ relative abundance is presented by the z-score of the log10 FPKM, after 2 h of temperature drop during the nocturnal period, of plants with a cold acclimation of 8 °C/0 °C (CA); nocturnal warming of 8 °C/6 °C (NW+); diurnal warming of 14 °C/0 °C (DW+); and diurnal-nocturnal warming of 14 °C/6 °C (DNW+). The plants were cultivated in chambers under controlled conditions at a photoperiod and thermoperiod of 18 h/6 h duration.
Ijms 24 11211 g004
Ijms 24 11211 g005 550
Figure 5. Relative expression of genes related to freezing tolerance. The variation on gene relative expression was evaluated by qRT-PCR. Values represent the fold change mean ± standard error (n = 3), in relation to cold-acclimated (CA) plants 8 °C/0 °C, at the beginning of the nocturnal period (0 h). Transcripts representing variation included: (A) ‘Dehydration-responsive element-binding protein 1A’ (DRE1A/CBF3), (B) ‘Dehydration-responsive element-binding protein 1D’ (DRE1D/CBF4), (C) ‘Dehydrin 1’ (DHN1), (D) ‘E3 ubiquitin-protein ligase RGLG5’ (RGLG5), (E) ‘UDP-glycosyltransferase 73B3’ (U73B3), (F) ‘Vacuolar cation/proton exchanger 1’ (CAX1), (G) ‘Ornithine aminotransferase’ (OAT), (H) ‘Sucrose synthase 4’ (SUS4), and (I) ‘Sucrose phosphate synthase 4’ (SPS4F). Warming treatments correspond to: nocturnal warming, 8 °C/6 °C (NW+); diurnal warming, 14 °C/0 °C (DW+); and diurnal-nocturnal warming, 14 °C/6 °C (DNW+), with 18 h/6 h photoperiod and thermoperiods duration. Significant differences among factors are shown as different lower cases (p < 0.05).
Figure 5. Relative expression of genes related to freezing tolerance. The variation on gene relative expression was evaluated by qRT-PCR. Values represent the fold change mean ± standard error (n = 3), in relation to cold-acclimated (CA) plants 8 °C/0 °C, at the beginning of the nocturnal period (0 h). Transcripts representing variation included: (A) ‘Dehydration-responsive element-binding protein 1A’ (DRE1A/CBF3), (B) ‘Dehydration-responsive element-binding protein 1D’ (DRE1D/CBF4), (C) ‘Dehydrin 1’ (DHN1), (D) ‘E3 ubiquitin-protein ligase RGLG5’ (RGLG5), (E) ‘UDP-glycosyltransferase 73B3’ (U73B3), (F) ‘Vacuolar cation/proton exchanger 1’ (CAX1), (G) ‘Ornithine aminotransferase’ (OAT), (H) ‘Sucrose synthase 4’ (SUS4), and (I) ‘Sucrose phosphate synthase 4’ (SPS4F). Warming treatments correspond to: nocturnal warming, 8 °C/6 °C (NW+); diurnal warming, 14 °C/0 °C (DW+); and diurnal-nocturnal warming, 14 °C/6 °C (DNW+), with 18 h/6 h photoperiod and thermoperiods duration. Significant differences among factors are shown as different lower cases (p < 0.05).
Ijms 24 11211 g005
Ijms 24 11211 g006 550
Figure 6. Relative expression of up-regulated genes by nocturnal warming. The variation in relative gene expression was evaluated by qRT-PCR. Values represent the fold change mean ± standard error (n = 3), in relation to cold-acclimated (CA) plants 8 °C/0 °C, at the beginning of the nocturnal period (0 h). Transcript variation shown includes: (A) ‘IAA-amino acid hydrolase ILR1-like 2’ (ILL2), (B) ‘1-deoxy-D-xylulose-5-phosphate synthase 1’ (DXS1), (C) ‘Triose phosphate/phosphate translocator’ (TPT), (D) ‘Oxygen-evolving enhancer protein 2’ (PSBP2), (E) ‘Fructose-bisphosphate aldolase 5’ (FAB5), and (F) ‘Ribulose bisphosphate carboxylase small subunit 1’ (RBCS1). Warming treatments correspond to: nocturnal warming, 8 °C/6 °C (NW+); diurnal warming, 14 °C/0 °C (DW+); and diurnal-nocturnal warming, 14 °C/6 °C (DNW+), with 18 h/6 h photoperiod and thermoperiods duration. Significant differences among factors are shown as different lower cases (p < 0.05).
Figure 6. Relative expression of up-regulated genes by nocturnal warming. The variation in relative gene expression was evaluated by qRT-PCR. Values represent the fold change mean ± standard error (n = 3), in relation to cold-acclimated (CA) plants 8 °C/0 °C, at the beginning of the nocturnal period (0 h). Transcript variation shown includes: (A) ‘IAA-amino acid hydrolase ILR1-like 2’ (ILL2), (B) ‘1-deoxy-D-xylulose-5-phosphate synthase 1’ (DXS1), (C) ‘Triose phosphate/phosphate translocator’ (TPT), (D) ‘Oxygen-evolving enhancer protein 2’ (PSBP2), (E) ‘Fructose-bisphosphate aldolase 5’ (FAB5), and (F) ‘Ribulose bisphosphate carboxylase small subunit 1’ (RBCS1). Warming treatments correspond to: nocturnal warming, 8 °C/6 °C (NW+); diurnal warming, 14 °C/0 °C (DW+); and diurnal-nocturnal warming, 14 °C/6 °C (DNW+), with 18 h/6 h photoperiod and thermoperiods duration. Significant differences among factors are shown as different lower cases (p < 0.05).
Ijms 24 11211 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

López, D.; Larama, G.; Sáez, P.L.; Bravo, L.A. Transcriptome Analysis of Diurnal and Nocturnal-Warmed Plants, the Molecular Mechanism Underlying Cold Deacclimation Response in Deschampsia antarctica. Int. J. Mol. Sci. 2023, 24, 11211. https://doi.org/10.3390/ijms241311211

AMA Style

López D, Larama G, Sáez PL, Bravo LA. Transcriptome Analysis of Diurnal and Nocturnal-Warmed Plants, the Molecular Mechanism Underlying Cold Deacclimation Response in Deschampsia antarctica. International Journal of Molecular Sciences. 2023; 24(13):11211. https://doi.org/10.3390/ijms241311211

Chicago/Turabian Style

López, Dariel, Giovanni Larama, Patricia L. Sáez, and León A. Bravo. 2023. "Transcriptome Analysis of Diurnal and Nocturnal-Warmed Plants, the Molecular Mechanism Underlying Cold Deacclimation Response in Deschampsia antarctica" International Journal of Molecular Sciences 24, no. 13: 11211. https://doi.org/10.3390/ijms241311211

APA Style

López, D., Larama, G., Sáez, P. L., & Bravo, L. A. (2023). Transcriptome Analysis of Diurnal and Nocturnal-Warmed Plants, the Molecular Mechanism Underlying Cold Deacclimation Response in Deschampsia antarctica. International Journal of Molecular Sciences, 24(13), 11211. https://doi.org/10.3390/ijms241311211

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