Next Article in Journal
Population Structure and Association Mapping for Agronomical and Biochemical Traits of a Large Spanish Apple Germplasm
Previous Article in Journal
Plant Species’ Capacity for Range Shifts at the Habitat and Geographic Scales: A Trade-Off-Based Framework
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 6
10.3390/plants12061250
plants-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

De Novo Transcriptome Assembly and Comparative Analysis of Differentially Expressed Genes Involved in Cold Acclimation and Freezing Tolerance of the Arctic Moss Aulacomnium turgidum (Wahlenb.) Schwaegr

by
Pilsung Kang
1,2,†,
Yo-Han Yoo
1,†,
Dong-Il Kim
2,
Joung Han Yim
1 and
Hyoungseok Lee
1,3,*
1
Division of Life Sciences, Korea Polar Research Institute, Incheon 21990, Republic of Korea
2
Department of Biological Engineering, Inha University, Incheon 22212, Republic of Korea
3
Polar Science, University of Science and Technology, Incheon 21990, Republic of Korea
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2023, 12(6), 1250; https://doi.org/10.3390/plants12061250
Submission received: 30 December 2022 / Revised: 6 March 2023 / Accepted: 7 March 2023 / Published: 9 March 2023
(This article belongs to the Topic Temperature Stress and Responses in Plants)
Browse Figures
Versions Notes

Abstract

:
Cold acclimation refers to a phenomenon in which plants become more tolerant to freezing after exposure to non-lethal low temperatures. Aulacomnium turgidum (Wahlenb.) Schwaegr is a moss found in the Arctic that can be used to study the freezing tolerance of bryophytes. To improve our understanding of the cold acclimation effect on the freezing tolerance of A. turgidum, we compared the electrolyte leakage of protonema grown at 25 °C (non-acclimation; NA) and at 4 °C (cold acclimation; CA). Freezing damage was significantly lower in CA plants frozen at −12 °C (CA-12) than in NA plants frozen at −12 °C (NA-12). During recovery at 25 °C, CA-12 demonstrated a more rapid and greater level of the maximum photochemical efficiency of photosystem II than NA-12, indicating a greater recovery capacity for CA-12 compared to NA-12. For the comparative analysis of the transcriptome between NA-12 and CA-12, six cDNA libraries were constructed in triplicate, and RNA-seq reads were assembled into 45,796 unigenes. The differential gene expression analysis showed that a significant number of AP2 transcription factor genes and pentatricopeptide repeat protein-coding genes related to abiotic stress and the sugar metabolism pathway were upregulated in CA-12. Furthermore, starch and maltose concentrations increased in CA-12, suggesting that cold acclimation increases freezing tolerance and protects photosynthetic efficiency through the accumulation of starch and maltose in A. turgidum. A de novo assembled transcriptome can be used to explore genetic sources in non-model organisms.

1. Introduction

Mosses are one of the oldest land plants and are spread all around the world [1]. They typically inhabit humid environments like rainforests, wetlands, and alpine ecosystems. However, they are also distributed in dry and cold polar regions [2]. Aulacomnium turgidum (Wahlenb.) Schwaegr is a moss widely found above the Arctic Circle, including in Svalbard and Alaska, with extreme climates [3]. It has been reported to be capable of tissue regeneration after 400 years of ice entombment [4], suggesting that it has a unique mechanism for freezing tolerance. In previous studies, its transcriptome information was revealed through the construction and analysis of a small-scale complementary DNA (cDNA) library [5], and a complete mitogenome sequence of A. turgidum has been reported [6]. However, even though A. turgidum is suitable for studying freezing tolerance, our understanding of A. turgidum remains insufficient.
Due to their sessile lifestyle, plants must cope with abiotic stresses such as cold, heat, salinity, and drought. Freezing stress is a significant challenge for plants to survive [7]. When exposed to freezing environments, organelles and cell membranes are physically damaged by ice crystals and osmotic stress occurring during the freezing process [8]. In addition, reactive oxygen species (ROS) such as singlet oxygen (1O2), superoxide (O2∙), hydroxyl radical (OH∙) and H2O2 are excessively generated and accumulated by electrons leaked from physically damaged mitochondria [9]. Enzymes, DNA, and unsaturated fatty acids are then denatured by accumulated ROS [10], which leads to increased cell damage and, ultimately, cell death [11].
Plant cells exposed to freezing stress can increase their tolerance using various cellular compounds and molecular events. These include anti-freeze compounds, osmoprotectants (proline, sucrose, fructan, mannitol, and glycine betaine), secondary messenger (calcium ion flux), and cold-regulated (COR) genes, including the LEA family [10,12]. For example, when plant cells are exposed to freezing-induced dehydration, COR15a protein increases the cryostability of the plasma membrane, which helps increase cell survival even at low freezing temperatures [13,14]. Fructan acts as an anti-freezing agent protecting cells by lowering the freezing point of water in the cells [15]. Proline, fructan, and glycine betaine play a role in regulating osmosis under osmotic stress due to dehydration, and they also have enzyme degradation inhibition and antioxidant effects [10,16]. The cytosolic calcium concentration of plant cells increases as plants are exposed to cold stress. This triggers phosphorylation cascades through calcium-binding proteins (CBPs) interacting with other proteins, ultimately leading to cold stress resistance. Interestingly, some species are able to increase freezing tolerance after exposure to low non-freezing temperatures, termed cold acclimation [17,18].
Genome-wide transcriptome analysis aids in understanding protective cellular responses. Quantification and comparison of transcript expression are widely used in molecular biology research. In particular, RNA-Seq technology and de novo transcriptome assembly are utilized to reconstruct and quantify whole transcriptomes in non-model organisms without a reference genome [19]. These technological advances can be employed to assess differential gene expression for non-model organisms and explore genetic diversity. Although mosses comprise a large phylum containing 13,000 species [20], only 21 genomes from eight genera have been reported, including model species, Antarctic moss, and peat mosses [21,22,23,24]. Additionally, RNA-seq-based de novo transcriptome studies have been reported in several non-seed plants, including Funaria hygrometrica (habitat: nutrient-rich soils including old bonfire sites) [25], Bryum argenteum (habitat: very dry and usually rich in nutrients such as nitrates) [26], Ceratodon purpureus (habitat: common and cosmopolitan in healthland and grassland) [27], Leptobryum pyriforme (habitat: easily found as a weed in glasshouses) [28], and Dicranum scoparium (habitat: frequent on the ground in woodland, and trees and logs) [29]. In 2019, the 1KP Plants project released transcriptomes of 1173 plant species, including 41 mosses [30]. However, there is still a need to accumulate more genomic resources through genome-wide transcriptome analysis of various moss species.
In the present study, we prepared two groups of axenic protonema of A. turgidum, with and without cold acclimation. After freezing stress treatment, we aimed to investigate transcriptome changes that mediate the cold acclimation process. To this end, we obtained Illumina technology-based RNA sequence data and created a transcriptome reference for A. turgidum through de novo assembly. Differentially expressed genes were selected from the cold acclimation treatment sample, and information about their biological functions was provided. Finally, we measured the concentration of starch and sugars to determine whether cold acclimation affects carbohydrate metabolism. The overall process for this study is summarized in Figure 1.

2. Results

2.1. Physiological Changes Due to Cold Acclimation Treatment under Freezing Stress

As electrolyte leakage (EL) is considered a good indicator of cell membrane integrity [31], we performed an electrolyte leakage assay to compare membrane integrity and freezing sensitivity of A. turgidum protonema with cold acclimation (CA) and without (non-acclimation; NA). The protonema of A. turgidum cultured at 25 °C for two weeks was transferred to a fresh BCDAT agar plate and acclimated at 4 °C (CA) and 25 °C (NA) under continuous light for 48 h. Then, changes in EL of protonema of A. turgidum at different sub-zero temperatures (0 to −20 °C) were examined. When the 50% electrolyte leakage (EL50) of the NA and CA were compared, the temperature of CA EL50 (−10.6 °C) was 3.3 °C lower than the temperature of NA EL50 (−7.3 °C). Furthermore, the greatest difference in EL of CA and NA was observed at −12 °C (Figure 2A).
Protonema of both CA and NA exposed to freezing at −12 °C were designated as CA-12 and NA-12, respectively. To assess the effect of freezing on the photosynthetic capacity of A. turgidum, they were transferred to the 25 °C culture room for recovery. Immediately after freezing exposure, FV/FM values of CA-12 and NA-12 were 0.28 ± 0.07 and 0.19 ± 0.02, respectively. On the second day of recovery, the value of CA-12 was higher than that of NA-12 and lower than that of the unfrozen control. By the fifth day of recovery, the value of CA-12 (0.66 ± 0.01) had recovered to the same level as the control (0.69 ± 0.01), whereas NA-12 (0.54 ± 0.04) had only recovered to 79.5% of the control (Figure 2B). This result suggested that cold acclimation resulted in physiological changes in the protonema of A. turgidum under freezing stress. Therefore, RNA from CA-12 and NA-12 samples was isolated for RNA sequencing.

2.2. De Novo Assembly and Annotation of Assembled A. turgidum Protonema Transcripts

To comprehensively understand the transcriptome profile related to enhanced freezing tolerance due to cold acclimation, NA and CA plants were sampled in triplicate after treatment at −12 °C and then sequenced. The total number of raw reads from six libraries was 21,339,944, and 19,990,976 trimmed reads were generated using Trimmomatic (i.e., 93.65% of trimmed rate; Table 1). De novo assembly was carried out using Trinity, and 116,979 contigs and 45,796 unigenes with N50 (bp) values of 2481 and 1017 were obtained, respectively (Table 2).
Functional annotation of each unigene was obtained by BLASTx against the non-redundant protein database (Nr). The similarity analysis of the Nr database revealed that 56.12% of unigenes had a significant homology (E-value < 1 × 10−50) (Figure 3). The Nr annotation species distribution analysis revealed that a majority of A. turgidum unigenes had a significant match to Physcomitrella sp. (47.20%), followed by Selaginella sp. (10.74%), Picea sp. (3.38%), Amborella sp. (3.04%), Nelumbo sp. (2.18%), Vitis sp. (1.45%), Elaeis sp. (1.34%), Phoenix sp. (1.20%), Gossypium sp. (1.12%), Musa sp. (0.97%), and others (27.37%). Because of research into A. turgidum has been rarely carried out. Only three (0.0040%) unigenes were annotated in the Nr database against A. turgidum (Figure 3, Table S1).

2.3. Gene Ontology (GO) Distribution and KEGG Pathway Analysis Using Whole Unigenes

The GO distribution of the A. turgidum unigenes was analyzed. A total of 10,228 (22.33%) unigenes were assigned to at least one GO term, and 6992, 4965, and 6375 unigenes were classified into the groups of “biological processes”, “molecular functions”, and “cellular components”, respectively (Table S2). Within the “biological processes” category, “cellular processes”, “metabolic processes”, and “single-organism processes” were the most abundant sub-categories. In the “molecular functions” category, “binding” and “catalytic activity” sub-categories made up the main proportion of the unigenes. In the sub-category of “cellular components”, “cell parts”, “cells”, “membranes”, and “organelles” were abundant (Table S2). In addition, we performed KEGG pathway annotations with the acquired unigenes to gain insight into the molecular interactions, reactions, and relationship networks of A. turgidum. 2171 (4.74%) unigenes were assigned to at least one KEGG pathway annotation into 42 KEGG pathways. They were identified in five categories: “metabolism”, “genetic information processing”, “environmental information processing”, “cellular processes”, and “organismal systems”. The most abundant categories in the KEGG Orthology hierarchies were “ribosomes”, “protein processing in the endoplasmic reticulum”, and “spliceosome” (Table S3).

2.4. RNA-Seq Analysis Identified Candidate Unigenes Associated with Cold Acclimation and Freezing Stress

To identify the genes showing differential expression patterns by cold acclimation in A. turgidum, we compared CA-12 and NA-12 with three biological replicates. In CA-12, there were 2137 unigenes upregulated and 1907 unigenes downregulated when compared with NA-12 (p-values < 0.05 and |log2 fold-change| > 1) (Figure 4A, Table S4). To verify the DEGs, we selected four unigenes (unigene_02942, unigene_05085, unigene_06026, and unigene_13443) that were upregulated and four (unigene_21963, unigene_25485, unigene_29992, and unigene_38286) that were downregulated in CA-12. Then, we confirmed the expression patterns of the eight unigenes by real-time quantitative PCR (qPCR). The RNA-Seq data showed a significant positive correlation with the qPCR data (Figure 4B), verifying that DEGs selected through RNA-Seq analysis had high accuracy.

2.5. Transcription Factors Analysis of the Upregulated Group in the CA-12 Treatment

Transcription factors (TFs) are major regulators that control gene expression by binding to specific DNA sequences. We identified 25 upregulated unigenes encoding TFs out of 2137 unigenes upregulated by cold acclimation followed by freezing. These included 13 APETALA2/ethylene responsive factor (AP2/ERF) TFs (unigene_04492, unigene_08384, unigene_08902, unigene_09280, unigene_09431, unigene_20367, unigene_25566, unigene_26682, unigene_27557, unigene_27793, unigene_28723, unigene_29433, and unigene_34993), two basic helix-loop-helix (bHLH) TFs (unigene_02361 and unigene_39888), two MYB TFs (unigene_11015 and unigene_45663), two WRKY TFs (unigene_25745 and unigene_40956), and 6 other upregulated TF encoding unigenes under CA-12 conditions (Figure 5). These results indicate that AP2/ERF are likely pivotal regulators for the enhanced freezing tolerance due to cold acclimation in A. turgidum.

2.6. GO Distribution and KEGG Pathway Analysis of the Upregulated Group in CA-12

To identify the biological functions of the 2137 unigenes upregulated by CA-12 treatment, we performed a GO term analysis of these genes in the “biological processes” category. This revealed that 16% related to “cellular process”, 15% to “metabolic process”, 13% to “single-organism process”, 8% to “biological regulation”, and 7% to “regulation of biological processes”. In addition, the sub-category “response to stimulus” (6%) was significantly associated with cold acclimation (Figure 6, Table S6). Interestingly, the pentatricopeptide repeat protein was most commonly found in this category along with “ATP-binding cassette transporter”, “lipoxygenase”, and “cellulose synthase” (Table S6).
Subsequent KEGG pathway analysis (corrected p-Value < 0.01; Table 3) of these 2137 upregulated unigenes identified 22 pathways, 21 of which were included in the “metabolism” category. Interestingly, 10 of the 21 subcategories were associated with carbohydrate metabolisms, such as starch and sucrose metabolism (ppp00500; 32 unigenes), glycolysis/gluconeogenesis (ppp00010; 26 unigenes), amino sugar and nucleotide sugar metabolism (ppp00520; 23 unigenes), the pentose phosphate pathway (ppp00030; 16 unigenes), galactose metabolism (ppp00052; 13 unigenes), fructose and mannose metabolism (ppp00051; 12 unigenes), and others (Table 3 and Table S7). These results suggest a close relationship between cold acclimation and carbohydrate metabolism, which could potentially affect the freezing tolerance of A. turgidum.

2.7. Differences in Starch and Sugar Concentrations According to Cold Acclimation and Freezing Stress

To determine if cold acclimation affects carbohydrate metabolism, we measured the starch concentrations and detected a higher starch content in CA-12 than in NA-12 (Figure 7A). In addition, concentrations of maltose, sucrose, glucose, and fructose were measured using HPLC. The content of maltose and glucose in CA-12 was significantly higher than that of NA-12. However, the fructose content was lower in CA-12, and there was no significant difference in sucrose between NA-12 and CA-12 (Figure 7B–E). These results suggest that cold acclimation in A. turgidum increases tolerance to freezing stress through the accumulation of starch, maltose, and glucose.

3. Discussion

Although it is an extreme environment, many plant species live in the Arctic. Several ecological studies have reported that these plants have evolved in the direction of acquiring a freezing tolerance mechanism to adapt to extreme environments [2,32,33,34]. Still, explanations using physiological and molecular biological analysis analyses are limited.

3.1. Enhanced Freezing Tolerance by Cold Acclimation

Freezing can decrease photosynthetic efficiency by damaging chloroplasts, and the FV/FM ratio serves as an indicator to easily measure the degree of thylakoid membrane damage [35]. Researchers have found significant differences among winter wheat cultivars in terms of FV/FM recovery from freezing [36]. They observed that the degree of recovery of FV/FM in sweet cherry leaves depends on the number of freezing events experienced [37]. We assessed the level of freezing tolerance in A. turgidum protonema both before and after cold acclimation using a combination of electrolyte leakage and chlorophyll fluorescence measurement. The results indicate that CA-12 has an enhanced tolerance to freezing and a greater capacity for photosynthetic recovery than NA-12, showing more rapid and higher FV/FM recovering rates (Figure 2). Furthermore, we performed RNA-sequencing to generate transcriptome reference of A. turgidum and selected unigenes whose expression increased following cold acclimation under freezing stress. The transcriptome acquired in this study provides insight into the gene regulatory network of Arctic mosses and other non-seed plants for freezing stress tolerance mechanisms.

3.2. AP2/ERF Family Transcription Factors Are Predicted to Contribute to Freezing Stress Tolerance of A. turgidum

Transcription factors (TFs) are important regulators that activate or suppress gene expression to modulate signal transduction and play a pivotal role in plant development, cell signaling, and stress response [38]. TFs regulate downstream target gene expression by binding to cis-acting elements in the promoter region [39]. In this research, we identified 2137 unigenes that responded to cold acclimation under freezing stress, of which 25 were classified as TFs (Figure 5). Various major TF families, such as AP2/ERF, MYB, NAC, and WRKY, contribute to improved plant resistance in response to various stimuli [40]. Interestingly, 13 ou 25 TFs belonged to the AP2/ERF family (Figure 5). The AP2/ERF family proteins have an AP2/ERF DNA-binding domain that interacts directly with GCC box and/or dehydration-responsive element (DRE)/C-repeat element (CRT) cis-acting elements of the promoter [41]. This family can be further divided into five subfamilies, depending on the number and similarity of DNA-binding domains: AP2 (APETALA2), RAV (related to ABI3/VP1), DREB (dehydration-responsive element binding protein), ERF (ethylene-responsive factor), and others [42]. Among these subfamilies, the DREB-A1 subgroup, including several C-Repeat-Binding Factors (CBFs), is known to increase freezing stress tolerance. CBF activates DRE containing Cold Responsive Genes (CORs), along with the Inducer of CBF Expression (ICE) [43]. CORs encode the Late Embryogenesis Abundant (LEA) protein which enables resistance to cold stress in plants by modifying sugar metabolism and fatty acid desaturation [44]. Studies on the improvement of cold stress tolerance by DREBs have been conducted in various plants: TaDREB1 in wheat [45], BpERF13 in birch [46], and OsDREB1G in rice [47]. For example, expression of the DREB1A gene in Arabidopsis thaliana was induced by low-temperature stress, and transgenic plants overexpressing DREB1A showed enhanced tolerance to freezing and dehydration [48]. In rice, the OsDREB1A and OsDREB1B expression was induced by cold stress, and overexpression of OsDREB1A in transgenic Arabidopsis led to higher tolerance to drought, high salt, and freezing stresses [49]. Furthermore, AP2/ERF TFs have been known to protect plants from multiple stresses by the action of various plant hormones, protein interaction partners, and plant epigenetics, such as DNA methylation and histone modification [50,51]. Thus, the role of the AP2/ERF family TFs may be one of the reasons why A. turgidum can survive in the Arctic region. The 13 AP2/ERF unigenes presented in Figure 5 can be used to understand the biological processes that A. turgidum undergoes to acquire tolerance to freezing stress.

3.3. AP2/ERF Family TFs Are Associated with Pentatricopeptide Repeat (PPR) Proteins and Carbohydrate Metabolism

Based on GO and KEGG analysis of upregulated unigenes in CA-12, we found that pentatricopeptide repeat (PPR) proteins and sugar metabolism were enriched by cold acclimation and freezing (Figure 6; Table 3). PPR proteins are RNA-binding proteins characterized by tandem arrays of a degenerate 35-amino-acid (PPR motifs) [52]. The PPR protein not only participates in post-transcriptional processes, including RNA editing, splicing, stability, cleavage, degradation, and translation but also appears to play an important role in response to abiotic stresses [53,54]. PPR proteins constitute a large family of land plants—with 450 members in Arabidopsis and 477 in rice, while the moss Physcomitrella patens have only 105 PPR genes [55,56]. Interestingly, PPR proteins are involved in the expression of AP2/ERF TFs. For example, GUN1, a DNA-binding chloroplast PPR protein, regulates the expression of ABI4 (ABSCISIC ACID INSENSITIVE-4), one of the AP2/ERF TFs [57]. In addition, the cytosol-nucleus dual-localized PPR protein SOAR1 is known to regulate the expression of cold-responsive genes in the C-repeat binding factor/DRE-binding factor (CBF/DREB) transcriptional regulatory cascade. SOAR1-overexpression lines were reported to increase the expression of CBF1/DREB1B, CBF2/DREB1C, and CBF3/DREB1A, as well as CBF downstream regulon genes COR15A, COR15B, COR414, and KIN1, during cold stress [58]. This implies that PPR might influence the acquisition of tolerance to freezing stress by regulating the expression of AP2/ERF TFs.
The AP2/ERF family also participates in carbohydrate metabolism. For instance, the galactinol synthase-encoding gene has CRT/DRE regulatory elements in its promoter region and is highly expressed in CBF-overexpressing Arabidopsis and rice [59,60]. In contrast, Arabidopsis GOLS3 in cbf123-1 triple mutants was reduced by about 55-fold, and the expression of 28 other genes involved in carbohydrate metabolisms, such as sucrose synthase and β-glucosidases, was also significantly decreased [44]. In addition, SUN6 (sucrose uncoupled-6)/ABI4 mutation has been reported to be insensitive to sugars, a substrate of hexokinase. These findings suggest that SUN6/ABI4 may play a role in the hexokinase-dependent sugar responses [61]. There have been multiple reports that freezing tolerance is improved by AP2/ERF [43,62] and soluble sugars [63,64]. We further verified that the concentrations of starch and maltose increased, and various unigenes related to sugar metabolism were upregulated by cold acclimation (Figure 7 and Figure S1). Taken together, the role of AP2/ERF family TFs is presumed to be one of the reasons for this result.

3.4. Starch and Soluble Sugars Play an Important Role in Freezing Stress Tolerance

Starch is the major storage unit of carbohydrates and the primary products of photosynthesis in plants. Recently, its importance in energy metabolism, developmental processes, and temperature acclimation has been recognized [65]. Interestingly, the role of starch in adaptation to low temperatures varies between different plant species [66]. For instance, Arabidopsis increases the content of glucose, fructose, sucrose, maltose, and starch under cold stress (4 °C) [67]. In this study, A. turgidum CA plants whose freezing tolerance was increased due to cold acclimation also had higher starch content than control NA plants. On the other hand, Physcomitrella patens was reported to have reduced starch content at low temperatures (−16 °C to 4 °C) [68]. These results confirm that the dynamics of synthesis and breakdown pathways play a more pivotal role than the absolute amount of starch under abiotic stresses [69].
Accumulation of soluble sugars correlates with the stabilization of biological components under cold stress [70]. Maltose, a disaccharide produced through the decomposition of starch by beta-amylase (βAM) [64], serves as an intermediate component when starch provides hexose to synthesize sucrose or when it is broken down into glucose [71,72]. Under cold stress, photosynthesis is maintained by maltose as it helps to protect chloroplasts from the osmotic stress caused by dehydration [69]. Freezing tolerance was enhanced by increasing maltose metabolism in MAL62 (a maltase encoding gene) overexpressing yeast [73]. Glucose is a well-known substance to assist in the regulation of intracellular osmotic pressure [74]. The higher glucose content of cold-tolerant Pinus halepensis is associated with increased cold tolerance compared to cold-sensitive P. halepensis [75]. This is in accordance with changes in the sugar contents of A. turgidum CA-12 plants observed in the current study. Starches, maltose and glucose concentrations may have increased in cold-acclimated A. turgidum plants due to the enhanced transcription of unigenes associated with the synthesis, decomposition, conversion, and transport of starch, maltose, sucrose, glucose, and fructose among the 2137 unigenes which were upregulated when the plants were exposed to freezing stress (Figure S1).
Plants are exposed to both biotic (insect herbivores and microbial pathogens) and abiotic (extreme temperature and inappropriate water supply) stresses throughout their lifetimes. Acclimation is a key strategy employed to survive these environmental stresses. Priming is, like acclimation, the strategy of remembering the first biotic attack to respond effectively to the second attack [76]. Research has shown that a “primed” state can increase defence responses and enhance resistance to stress [77]. An analysis of expression changes of unigenes linked to cold acclimation in A. turgidum is expected to provide biological information on molecular strategies of polar mosses to defend against various external stimuli.

4. Materials and Methods

4.1. Sample Preparation

A. turgidum specimens were collected in the vicinity of the Korean Dasan Arctic Station (78°54′ N/11°57′ W) near Ny-Ålesund, Svalbard, in August 2006. A small patch of A. turgidum was washed using 30 mL 1.2% NaOCl containing 1 drop Tween-80 at 180 rpm for 15 min, then washed twice with sterilized water. Axenic A. turgidum was cultured in BCDAT medium containing 250 mg/L KH2PO4, 250 mg/L MgSO4, 12.5 mg/L FeSO4·7H2O, 920.5 mg/L ammonium tartrate, 1010 mg/L KNO3, 147 mg/L CaCl2·2H2O, 0.614 mg/L H3BO3, 0.389 mg/L MnCl2·4H2O, 0.055 mg/L CuSO4·5H2O, 0.055 mg/L CoCl·6H2O, 0.055 mg/L ZnSO4·7H2O, 0.028 mg/L KI, 0.025 mg/L Na2MoO4·2H2O, and 0.8% phytoagar at pH 5.8.

4.2. Cold Acclimation and Freezing Stress Treatment

The two-week-old A. turgidum protonema grown on fresh BCDAT agar plates were transferred to a 4 °C cold room (cold acclimation; CA), and the others were left in the 25 °C culture room (non-acclimation; NA) under continuous light for 2 days. The refrigerated circulating water bath program and freezing stress treatment were carried out as described [78,79]. Briefly, 100 mg (dry weight) protonema from each of the CA and NA groups were transferred to tubes and rinsed carefully three times using deionized water (DW). Then, 100 μL DW was added to each tube. The tubes were placed in the refrigerated circulating bath (Gaon Science Instrument, South Korea) with the temperature preset at 0 °C for 30 min. Tiny pieces of ice were added to each tube for ice nucleation, and the circulating bath was programmed to cool at a rate of −2 °C per hour. The tubes were removed from the bath when the designated temperature was reached. Electrolyte leakage measurement and total RNA isolation were performed with samples exposed to freezing conditions (0, −4, −8, −12, −16, and −20 °C), frozen in liquid nitrogen, and then stored at −80 °C, as illustrated in Figure 1.

4.3. Electrolyte Leakage Measurement

The electrolyte leakage percentage of A. turgidum (NA and CA) exposed to different freezing conditions (0, −4, −8, −12, −16, and −20 °C) was measured as previously described, with minor modification [78]. Samples were removed when they reached their designated temperatures and then immediately placed on ice to allow gradual thawing. After complete thawing, 2.4 mL DW was added to the tubes, followed by overnight shaking at room temperature. The conductivity of the solution in each tube was measured with a conductivity meter B-173 (Horiba, Japan). The tubes were then autoclaved. After shaking at room temperature, the conductivity of the solution in each tube was measured again. Finally, the percentage of electrolyte leakage was calculated as the percentage of the conductivity before autoclaving divided by that after autoclaving [79].

4.4. Chlorophyll Fluorescence Measurement

Immediately after freezing, A. turgidum plants were transferred to a BCDAT medium in a 25 °C culture room. Following 2 and 5 days of recovery, the maximum photochemical efficiency of photosystem II (calculated as FV/FM, where FM is the maximum fluorescence of the dark-adapted plants under the light-saturating flash, and FV is the maximum variable fluorescence, FM − F0) was measured with samples dark-adapted for 30 min, using IMAGING-PAM (Walz, Effeltrich, Germany).

4.5. RNA Extraction and Library Preparation and Sequencing

To determine the transcriptome of A. turgidum protonema depending on whether or not cold acclimation treatment was undertaken, the total RNA of each sample (NA treated with −12 °C freezing stress; NA-12 vs. CA treated with −12 °C freezing stress; CA-12) was isolated using the RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions and then RNA was treated with RNase-Free DNase I (Takara, Japan) to remove any possible DNA. The integrity and concentration of RNA were determined using a Bioanalyzer (RIN > 6) (Agilent Technologies, Santa Clara, CA, USA) and a Qubit RNA Broad-Range Assay Kit (Life Technologies, Carlsbad, CA, USA), respectively. To construct the sequencing library, 1.5 μg of total RNA from each sample was used as the input for the TruSeq RNA sample prep kit v2 (Illumina, San Diego, CA, USA) following the manufacturer’s recommended method. In total, six libraries composed of three biological replicates of NA-12 and CA-12 samples were validated and then quantified using the Bioanalyzer and the library qPCR quantification method. These libraries were paired-end sequenced on a MiSeq using a MiSeq Reagent Kit v3 (2 × 300 bp) (Illumina). A total of 5 Gb of sequence data were generated (Q30 > 93%). The RNA-Seq data were deposited into the Sequence Read Archive under accession numbers SRS16292302 and SRS16292303.

4.6. De Novo Transcriptome Sequence Assembly and Functional Annotation

To improve the accuracy of the analysis, the quality check of the raw sequence data obtained after sequencing was performed using the fastQC (v0.11.1; http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and Trimmomatic (v0.33; http://www.usadellab.org/cms/?page=trimmomatic), with adaptor sequences and low-quality, reads removed. Filtered reads were merged to perform de novo assembly using Trinity (v2.1.1; http://trinityrnaseq.github.io/) with default parameters. Using TrinityStats, Transdecoder.LongOrfs (https://transdecoder.github.io/), and the CD-HIT-EST v4.6 (https://github.com/weizhongli/cdhit/blob/master/doc/cdhit-user-guide.wiki), assembled reads were made into high-quality, unique transcripts. The resulting unigene sequences were compared against the plant non-redundant (nr) protein database at the National Centre for Biotechnology Information (NCBI) using BLASTP with an E-value parameter not greater than 1 × 10−4 for identification of the best significant match. The BLASTP results were then imported into Blast2GO v3.1 [80] to retrieve the Gene Ontology (GO) terms of the assembled unigenes, and the annotation was further continued with unique enzyme codes (EC) and Kyoto Encyclopedia of Genes and Genomes (KEGG) maps. Moreover, the KEGG Automated Annotation Server (KAAS) was used for pathway mapping in addition to Blast2GO. GO terms are precisely defined as controlled vocabulary which can be used to describe the functions of genes or gene products. The assembled transcripts based on the retrieved GO terms were classified into three categories: biological processes, molecular functions, and cellular components. Pathway maps were determined from the KEGG database with an E-value of 1 × 10−5.

4.7. RNA-Seq Analysis

The high-quality trimmed raw reads were mapped to an assembled reference transcriptome by the HISAT v2.0.0 (Bowite2 v2.2.3; https://ccb.jhu.edu/software/hisat2/index.shtml) with default parameters, and the number of mapped reads were calculated by samtools v1.2 (http://www.htslib.org/). Differentially Expressed Genes (DEGs) were evaluated using Cuffdiff to compare treatment conditions. Unigenes with p-values < 0.05 and log2 fold-changes > 1 (i.e., fold-change > 2) were considered to be differentially expressed. A further screening of initial DEGs was based on fragments per kilobase per million fragments mapped (FPKM) values. The selected DEGs were used to create a heatmap using the Multi Experiment Viewer (MeV v.4.9.0) software tool [81].

4.8. Starch and Sugars Concentration Analysis

The starch extraction was performed as previously described [82]. Briefly, each sample was treated with 0.7 M perchloric acid. After centrifugation at 16,000× g for 5 min at room temperature, the supernatant was removed; the pellet was washed twice with 80% (v/v) ethanol and dried. DMSO was added, and the sample was incubated at 100 °C for 20 min to disperse the starch. The starch content was measured using Lugol’s solution at 620 nm. The absorbances of the sample were converted to concentration (ug/mL) using a standard curve of soluble starch (48425-1501, Junsei, Tokyo, Japan). To measure the soluble sugar content, the method described in previous research [83] was partially modified and used. One milliliter of 80% (v/v) methanol (MeOH) was added to 200 mg of the homogenized sample and mixed vigorously. Then, the mixed solution was incubated at 90 °C for 10 min and centrifuged at 16,000× g for 5 min. The supernatant was transferred to a new test tube, and the extraction process was repeated. The collected extract was lyophilized, dissolved in DW, and analyzed using HPLC with a Quaternay pump (flow rate: 0.5 mL/min), RI detector (Ultimate 3000; Dionex, Sunnyvale, CA, USA), and Aminex 87P ion exclusion column (10 µm, 7.8 × 300 mm). Distilled water was used as a mobile phase. The content of each detected sugar was calculated using the peak area of the sugar standards (maltose, sucrose, glucose, fructose) analyzed by HPLC (Figure S2).

4.9. Statistical Analysis

Data were presented as the mean ± standard deviation. Statistical analyses were performed using Student’s t-test in Microsoft excel (* p < 0.05; ** p < 0.01; and *** p < 0.001).

5. Conclusions

Svalbard, a Norwegian archipelago, has a typical Arctic environment. Summer in the archipelago lasts 3.5 months, from mid-May to September, with an average temperature of 5.8 °C. Winter lasts 8.5 months, from September to mid-May of the following year, and has an average temperature of −12 °C. The overall Arctic atmosphere is dry [84], though atmospheric humidity is relatively higher in summer than in winter. Here, we used Aulacomnium turgidum (Wahlenb.) Schwaegr, a moss found in this region, to investigate cold acclimation and freezing tolerance through de novo transcriptome assembly and RNA-seq analysis. AP2/ERF family transcription factors and PPR protein genes involved in carbohydrate metabolism were identified as DEGs induced by cold acclimation under freezing stress. In addition, it was found that the concentration of starch and maltose in A. turgidum increased with cold acclimation. The unigenes and DEGs associated with cold acclimation identified in this study provide a foundation for further studies on the mechanisms behind freezing stress regulation in A. turgidum and other non-seed plant species. Furthermore, lineage-specific stress-regulated genes not found in flowering plants could be important clues for protecting crops in a climate change caused by global warming.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12061250/s1. Table S1. Detailed information on non-redundant database blast results. Table S2. Gene ontology (GO) assignments in the categories of biological processes, molecular functions, and cellular components for whole unigenes. Table S3. KEGG pathway analysis in the cellular processes, environmental information processing, genetic information processing, metabolism, and organismal systems for whole unigenes. Table S4. List of differentially expressed genes between NA-12 and CA-12. Table S5. Primer sequences used for real-time RT-PCR analyses. Table S6. Detailed information for Gene Ontology (GO) assigned to unigenes upregulated in CA-12. Table S7. Detailed information of KEGG pathway analysis assigned to unigenes upregulated in CA-12. Figure S1. Lists of unigenes related to the synthesis, decomposition, conversion, and transport of starch, maltose, glucose, fructose, and sucrose. Figure S2. Chromatograms of the sugar standards, NA-12 and CA-12 using HPLC-RI.

Author Contributions

Conceptualization, D.-I.K., J.H.Y. and H.L.; methodology, formal analysis investigation and data curation, P.K. and Y.-H.Y.; software, Y.-H.Y.; resources, P.K., D.-I.K. and J.H.Y.; writing—original draft preparation, P.K., Y.-H.Y. and H.L.; writing—review and editing, P.K., Y.-H.Y., D.-I.K., J.H.Y. and H.L.; supervision, D.-I.K. and J.H.Y.; funding acquisition, J.H.Y. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Korea Polar Research Institute (KOPRI) grants funded by the Ministry of Oceans and Fisheries (Project no. KOPRI PE23150 and PE23160).

Data Availability Statement

The data that support the findings of this study are openly available in GenBank of NCBI at https://www.ncbi.nlm.nih.gov. The BioProject, SRA, and Bio-Sample numbers are PRJNA916379, SRS16292302, SRS16292303, SAMN32423464, and SAMN32423465, respectively.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Morris, J.L.; Puttick, M.N.; Clark, J.W.; Edwards, D.; Kenrick, P.; Pressel, S.; Wellman, C.H.; Yang, Z.; Schneider, H.; Donoghue, P.C.J. The timescale of early land plant evolution. Proc. Natl. Acad. Sci. USA 2018, 115, E2274–E2283. [Google Scholar] [CrossRef] [Green Version]
  2. Gornall, J.L.; Jónsdóttir, I.S.; Woodin, S.J.; Van der Wal, R. Arctic mosses govern below-ground environment and ecosystem processes. Oecologia 2007, 153, 931–941. [Google Scholar] [CrossRef]
  3. GBIF Secretariat. GBIF Backbone Taxonomy. Available online: https://www.gbif.org/species/2675975 (accessed on 13 December 2022).
  4. La Farge, C.; Williams, K.H.; England, J.H. Regeneration of Little Ice Age bryophytes emerging from a polar glacier with implications of totipotency in extreme environments. Proc. Natl. Acad. Sci. USA 2013, 110, 9839–9844. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, S.; Lee, H.; Kang, P.S.; Huang, X.; Yim, J.H.; Lee, H.K.; Kim, I.C. Complementary DNA library construction and expressed sequence tag analysis of an Arctic moss, Aulacomnium turgidum . Polar Biol. 2010, 33, 617–626. [Google Scholar] [CrossRef]
  6. Kang, P.; Cho, S.M.; Lee, J.; Yim, J.H.; Lee, J.H.; Lee, H. The complete mitogenome of the Arctic moss Aulacomnium turgidum (Wahlenb.) Schwaegr. Mitochondrial DNA Part B 2019, 4, 3446–3447. [Google Scholar] [CrossRef]
  7. Levitt, J. Chilling, Freezing, and High Temperature Stresses. In Responses of Plants to Environmental Stress, 2nd ed.; Academic Press: New York, NY, USA, 1980; Volume 1, p. 497. [Google Scholar]
  8. Jang, T.H.; Park, S.C.; Yang, J.H.; Kim, J.Y.; Seok, J.H.; Park, U.S.; Choi, C.W.; Lee, S.R.; Han, J. Cryopreservation and its clinical applications. Integr Med. Res. 2017, 6, 12–18. [Google Scholar] [CrossRef]
  9. Park, J.; Lee, J.; Choi, C. Mitochondrial network determines intracellular ROS dynamics and sensitivity to oxidative stress through switching inter-mitochondrial messengers. PLoS ONE 2011, 6, e23211. [Google Scholar] [CrossRef]
  10. Parvanova, D.; Ivanov, S.; Konstantinova, T.; Karanov, E.; Atanassov, A.; Tsvetkov, T.; Alexieva, V.; Djilianov, D. Transgenic tobacco plants accumulating osmolytes show reduced oxidative damage under freezing stress. Plant Physiol. Biochem. 2004, 42, 57–63. [Google Scholar] [CrossRef]
  11. Orrenius, S.; Gogvadze, V.; Zhivotovsky, B. Mitochondrial oxidative stress: Implications for cell death. Annu. Rev. Pharmacol. Toxicol. 2007, 47, 143–183. [Google Scholar] [CrossRef]
  12. Hashempour, A.; Ghasemnezhad, M.; Sohani, M.M.; Ghazvini, R.F.; Abedi, A. Effects of freezing stress on the expression of fatty acid desaturase (FAD2, FAD6 and FAD7) and beta-glucosidase (BGLC) genes in tolerant and sensitive olive cultivars. Russ. J. Plant Physiol. 2019, 66, 214–222. [Google Scholar] [CrossRef]
  13. Artus, N.N.; Uemura, M.; Steponkus, P.L.; Gilmour, S.J.; Lin, C.; Thomashow, M.F. Constitutive expression of the cold-regulated Arabidopsis thaliana COR15a gene affects both chloroplast and protoplast freezing tolerance. Proc. Natl. Acad. Sci. USA 1996, 93, 13404–13409. [Google Scholar] [CrossRef] [Green Version]
  14. Sowemimo, O.; Borcherds, W.; Knox-Brown, P.; Rindfleisch, T.; Thalhammer, A.; Daughdrill, G. Evolution of transient helicity and disorder in late embryogenesis abundant protein COR15A. Biophys. J. 2019, 116, 473a. [Google Scholar] [CrossRef] [Green Version]
  15. Livingston, D.P.; Hincha, D.K.; Heyer, A.G. Fructan and its relationship to abiotic stress tolerance in plants. Cell. Mol. Life Sci. 2009, 66, 2007–2023. [Google Scholar] [CrossRef] [Green Version]
  16. Konstantinova, T.; Parvanova, D.; Atanassov, A.; Djilianov, D. Freezing tolerant tobacco, transformed to accumulate osmoprotectants. Plant Sci. 2002, 163, 157–164. [Google Scholar] [CrossRef]
  17. Hughes, M.A.; Dunn, M.A. The molecular biology of plant acclimation to low temperature. J. Exp. Bot. 1996, 47, 291–305. [Google Scholar] [CrossRef] [Green Version]
  18. Thomashow, M.F. Plant cold acclimation: Freezing tolerance genes and regulatory mechanisms. Annu. Rev. Plant Biol. 1999, 50, 571–599. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  20. Goffinet, B.; Buck, W.R.; Shaw, A.J. Morphology, anatomy, and classification of the Bryophyta. In Bryophyte Biology; Shaw, A.J., Ed.; Cambridge University Press: Cambridge, UK, 2008; pp. 55–138. [Google Scholar]
  21. Fernandez-Pozo, N.; Haas, F.B.; Gould, S.B.; Rensing, S.A. An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes. J. Exp. Bot. 2022, 73, 4291–4305. [Google Scholar] [CrossRef]
  22. Liu, S.; Fang, S.; Cong, B.; Li, T.; Yi, D.; Zhang, Z.; Zhao, L.; Zhang, P. The Antarctic moss Pohlia nutans genome provides insights into the evolution of bryophytes and the adaptation to extreme terrestrial habitats. Front. Plant Sci. 2022, 13, 920138. [Google Scholar] [CrossRef]
  23. Meleshko, O.; Martin, M.D.; Korneliussen, T.S.; Schröck, C.; Lamkowski, P.; Schmutz, J.; Healey, A.; Piatkowski, B.T.; Shaw, A.J.; Weston, D.J.; et al. Extensive genome-wide phylogenetic discordance is due to incomplete lineage sorting and not ongoing introgression in a rapidly radiated bryophyte genus. Mol. Biol. Evol. 2021, 38, 2750–2766. [Google Scholar] [CrossRef]
  24. Rensing, S.A.; Lang, D.; Zimmer, A.D.; Terry, A.; Salamov, A.; Shapiro, H.; Nishiyama, T.; Perroud, P.F.; Lindquist, E.A.; Kamisugi, Y.; et al. The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science 2008, 319, 64–69. [Google Scholar] [CrossRef] [Green Version]
  25. Szövényi, P.; Rensing, S.A.; Lang, D.; Wray, G.A.; Shaw, A.J. Generation-biased gene expression in a bryophyte model system. Mol. Biol. Evol. 2011, 28, 803–812. [Google Scholar] [CrossRef]
  26. Gao, B.; Zhang, D.; Li, X.; Yang, H.; Zhang, Y.; Wood, A.J. De novo transcriptome characterization and gene expression profiling of the desiccation tolerant moss Bryum argenteum following rehydration. BMC Genom. 2015, 16, 416. [Google Scholar] [CrossRef] [Green Version]
  27. Szövényi, P.; Perroud, P.F.; Symeonidi, A.; Stevenson, S.; Quatrano, R.S.; Rensing, S.A.; Cuming, A.C.; McDaniel, S.F. De novo assembly and comparative analysis of the Ceratodon purpureus transcriptome. Mol. Ecol. Resour. 2015, 15, 203–215. [Google Scholar] [CrossRef]
  28. Liu, S.; Fang, S.; Liu, C.; Zhao, L.; Cong, B.; Zhang, Z. Transcriptomics integrated with metabolomics reveal the effects of ultraviolet-B radiation on flavonoid biosynthesis in Antarctic moss. Front. Plant Sci. 2021, 12, 788377. [Google Scholar] [CrossRef]
  29. Chen, K.H.; Liao, H.L.; Arnold, A.E.; Korotkin, H.B.; Wu, S.H.; Matheny, P.B.; Lutzoni, F. Comparative transcriptomics of fungal endophytes in co-culture with their moss host Dicranum scoparium reveals fungal trophic lability and moss unchanged to slightly increased growth rates. New Phytol. 2022, 234, 1832–1847. [Google Scholar] [CrossRef]
  30. Carpenter, E.J.; Matasci, N.; Ayyampalayam, S.; Wu, S.; Sun, J.; Yu, J.; Vieira, F.R.J.; Bowler, C.; Dorrell, R.G.; Gitzendanner, M.A.; et al. Access to RNA-Sequencing Data from 1173 Plant Species: The 1000 Plant Transcriptomes Initiative (1KP). GigaScience 2019, 8, giz126. [Google Scholar] [CrossRef] [Green Version]
  31. Lee, B.H.; Zhu, J.K. Response to the Abiotic Environment. In Arabidopsis: A Laboratory Manual; Weigel, D., Glazebrook, J., Eds.; Cold Spring Harbor Laboratory Press: Laurel Hollow, NY, USA, 2002; pp. 71–75. [Google Scholar]
  32. Kennedy, A.D. Photosynthetic response of the Antarctic moss Polytrichum alpestre Hoppe to low temperatures and freeze-thaw stress. Polar Biol. 1993, 13, 271–279. [Google Scholar] [CrossRef]
  33. Alberdi, M.; Bravo, L.A.; Gutiérrez, A.; Gidekel, M.; Corcuera, L.J. Ecophysiology of Antarctic vascular plants. Physiol. Plant. 2002, 115, 479–486. [Google Scholar] [CrossRef]
  34. Bokhorst, S.; Huiskes, A.; Convey, P.; Aerts, R. The effect of environmental change on vascular plant and cryptogam communities from the Falkland Islands and the Maritime Antarctic. BMC Ecol. 2007, 7, 15. [Google Scholar] [CrossRef] [Green Version]
  35. Fracheboud, Y.; Leipner, J. The application of chlorophyll fluorescence to study light, temperature, and drought stress. In Practical Applications of Chlorophyll Fluorescence in Plant Biology; De Ell, J.R., Toivonen, P.M.A., Eds.; Springer: Boston, MA, USA, 2013; pp. 125–150. [Google Scholar]
  36. Rapacz, M. Chlorophyll a fluorescence transient during freezing and recovery in winter wheat. Photosynthetica 2007, 45, 409–418. [Google Scholar] [CrossRef]
  37. Vosnjak, M.; Sircelj, H.; Vodnik, D.; Usenik, V. Physio-biochemical responses of sweet cherry leaf to natural cold conditions. Plants 2022, 11, 3507. [Google Scholar] [CrossRef]
  38. Gonzalez, D.H. Introduction to transcription factor structure and function. In Plant Transcription Factors; Gonzalez, D.H., Ed.; Academic Press: Cambridge, MA, USA, 2016; pp. 3–11. [Google Scholar]
  39. Inukai, S.; Kock, K.H.; Bulyk, M.L. Transcription factor–DNA binding: Beyond binding site motifs. Curr. Opin. Genet. Dev. 2017, 43, 110–119. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, H.; Wang, H.; Shao, H.; Tang, X. Recent advances in utilizing transcription factors to improve plant abiotic stress tolerance by transgenic technology. Front. Plant Sci. 2016, 7, 67. [Google Scholar] [CrossRef] [Green Version]
  41. Nakano, T.; Suzuki, K.; Fujimura, T.; Shinshi, H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006, 140, 411–432. [Google Scholar] [CrossRef] [Green Version]
  42. Sakuma, Y.; Liu, Q.; Dubouzet, J.G.; Abe, H.; Shinozaki, K.; Yamaguchi-Shinozaki, K. DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration-and cold-inducible gene expression. Biochem. Biophys. Res. Commun. 2002, 290, 998–1009. [Google Scholar] [CrossRef]
  43. Liu, J.; Shi, Y.; Yang, S. Insights into the regulation of C-repeat binding factors in plant cold signaling. J. Integr. Plant Biol. 2018, 60, 780–795. [Google Scholar] [CrossRef] [Green Version]
  44. 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] [Green Version]
  45. Shen, Y.G.; Zhang, W.K.; He, S.J.; Zhang, J.S.; Liu, Q.; Chen, S.Y. An EREBP/AP2-type protein in Triticum aestivum was a DRE-binding transcription factor induced by cold, dehydration and ABA stress. Theor. Appl. Genet. 2003, 106, 923–930. [Google Scholar] [CrossRef]
  46. Lv, K.; Li, J.; Zhao, K.; Chen, S.; Nie, J.; Zhang, W.; Liu, G.; Wei, H. Overexpression of an AP2/ERF family gene, BpERF13, in birch enhances cold tolerance through upregulating CBF genes and mitigating reactive oxygen species. Plant Sci. 2020, 292, 110375. [Google Scholar] [CrossRef]
  47. Moon, S.J.; Min, M.K.; Kim, J.; Kim, D.Y.; Yoon, I.S.; Kwon, T.R.; Byun, M.O.; Kim, B.G. Ectopic expression of OsDREB1G, a member of the OsDREB1 subfamily, confers cold stress tolerance in rice. Front. Plant Sci. 2019, 10, 297. [Google Scholar] [CrossRef]
  48. Liu, Q.; Kasuga, M.; Sakuma, Y.; Abe, H.; Miura, S.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought-and low-temperature-responsive gene expression, respectively, in Arabidopsis . Plant Cell 1998, 10, 1391–1406. [Google Scholar] [CrossRef] [Green Version]
  49. Dubouzet, J.G.; Sakuma, Y.; Ito, Y.; Kasuga, M.; Dubouzet, E.G.; Miura, S.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt-and cold-responsive gene expression. Plant J. 2003, 33, 751–763. [Google Scholar] [CrossRef]
  50. Licausi, F.; Ohme-Takagi, M.; Perata, P. APETALA 2/Ethylene Responsive Factor (AP 2/ERF) transcription factors: Mediators of stress responses and developmental programs. New Phytol. 2013, 199, 639–649. [Google Scholar] [CrossRef]
  51. Xie, Z.; Nolan, T.M.; Jiang, H.; Yin, Y. AP2/ERF transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis . Front. Plant Sci. 2019, 10, 228. [Google Scholar] [CrossRef] [Green Version]
  52. Barkan, A.; Small, I. Pentatricopeptide repeat proteins in plants. Annu. Rev. Plant Biol. 2014, 65, 415–442. [Google Scholar] [CrossRef]
  53. Lurin, C.; Andrés, C.; Aubourg, S.; Bellaoui, M.; Bitton, F.; Bruyère, C.; Caboche, M.; Debast, C.; Gualberto, J.; Hoffmann, B.; et al. Genome-wide analysis of Arabidopsis pentatricopeptide repeat proteins reveals their essential role in organelle biogenesis. Plant Cell 2004, 16, 2089–2103. [Google Scholar] [CrossRef] [Green Version]
  54. Okuda, K.; Myouga, F.; Motohashi, R.; Shinozaki, K.; Shikanai, T. Conserved domain structure of pentatricopeptide repeat proteins involved in chloroplast RNA editing. Proc. Natl. Acad. Sci. USA 2007, 104, 8178–8183. [Google Scholar] [CrossRef] [Green Version]
  55. O’Toole, N.; Hattori, M.; Andres, C.; Iida, K.; Lurin, C.; Schmitz-Linneweber, C.; Sugita, M.; Small, I. On the expansion of the pentatricopeptide repeat gene family in plants. Mol. Biol. Evol. 2008, 25, 1120–1128. [Google Scholar] [CrossRef]
  56. Sugita, M.; Ichinose, M.; Ide, M.; Sugita, C. Architecture of the PPR gene family in the moss Physcomitrella patens . RNA Biol. 2013, 10, 1439–1445. [Google Scholar] [CrossRef] [Green Version]
  57. Koussevitzky, S.; Nott, A.; Mockler, T.C.; Hong, F.; Sachetto-Martins, G.; Surpin, M.; Lim, J.; Mittler, R.; Chory, J. Signals from chloroplasts converge to regulate nuclear gene expression. Science 2007, 316, 715–719. [Google Scholar] [CrossRef]
  58. Jiang, S.C.; Mei, C.; Liang, S.; Yu, Y.T.; Lu, K.; Wu, Z.; Wang, X.F.; Zhang, D.P. Crucial roles of the pentatricopeptide repeat protein SOAR1 in Arabidopsis response to drought, salt and cold stresses. Plant Mol. Biol. 2015, 88, 369–385. [Google Scholar] [CrossRef] [Green Version]
  59. Cui, L.H.; Byun, M.Y.; Oh, H.G.; Kim, S.J.; Lee, J.; Park, H.; Lee, H.; Kim, W.T. Poaceae Type II galactinol synthase 2 from Antarctic flowering plant Deschampsia antarctica and rice improves cold and drought tolerance by accumulation of raffinose family oligosaccharides in transgenic rice plants. Plant Cell Physiol. 2020, 61, 88–104. [Google Scholar] [CrossRef]
  60. Park, S.; Lee, C.M.; Doherty, C.J.; Gilmour, S.J.; Kim, Y.; Thomashow, M.F. Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network. Plant J. 2015, 82, 193–207. [Google Scholar] [CrossRef] [Green Version]
  61. Huijser, C.; Kortstee, A.; Pego, J.; Weisbeek, P.; Wisman, E.; Smeekens, S. The Arabidopsis SUCROSE UNCOUPLED-6 gene is identical to ABSCISIC ACID INSENSITIVE-4: Involvement of abscisic acid in sugar responses. Plant J. 2000, 23, 577–585. [Google Scholar] [CrossRef]
  62. Najafi, S.; Sorkheh, K.; Nasernakhaei, F. Characterization of the APETALA2/Ethylene-responsive factor (AP2/ERF) transcription factor family in sunflower. Sci. Rep. 2018, 8, 11576. [Google Scholar] [CrossRef] [Green Version]
  63. Kaplan, F.; Guy, C.L. β-Amylase induction and the protective role of maltose during temperature shock. Plant Physiol. 2004, 135, 1674–1684. [Google Scholar] [CrossRef] [Green Version]
  64. Tarkowski, Ł.P.; Van den Ende, W. Cold tolerance triggered by soluble sugars: A multifaceted countermeasure. Front. Plant Sci. 2015, 6, 203. [Google Scholar] [CrossRef] [Green Version]
  65. Zeeman, S.C.; Kossmann, J.; Smith, A.M. Starch: Its metabolism, evolution, and biotechnological modification in plants. Annu. Rev. Plant Biol. 2010, 61, 209–234. [Google Scholar] [CrossRef] [Green Version]
  66. Thalmann, M.; Santelia, D. Starch as a determinant of plant fitness under abiotic stress. New Phytol. 2017, 214, 943–951. [Google Scholar] [CrossRef] [Green Version]
  67. Kaplan, F.; Guy, C.L. RNA interference of Arabidopsis beta-amylase8 prevents maltose accumulation upon cold shock and increases sensitivity of PSII photochemical efficiency to freezing stress. Plant J. 2005, 44, 730–743. [Google Scholar] [CrossRef]
  68. Nagao, M.; Minami, A.; Arakawa, K.; Fujikawa, S.; Takezawa, D. Rapid degradation of starch in chloroplasts and concomitant accumulation of soluble sugars associated with ABA-induced freezing tolerance in the moss Physcomitrella patens . J. Plant Physiol. 2005, 162, 169–180. [Google Scholar] [CrossRef]
  69. Fürtauer, L.; Weiszmann, J.; Weckwerth, W.; Nägele, T. Dynamics of plant metabolism during cold acclimation. Int. J. Mol. Sci. 2019, 20, 5411. [Google Scholar] [CrossRef] [Green Version]
  70. Santarius, K.A. The protective effect of sugars on chloroplast membranes during temperature and water stress and its relationship to frost, desiccation and heat resistance. Planta 1973, 113, 105–114. [Google Scholar] [CrossRef]
  71. Sauter, J.J.; Van Cleve, B. Occurrence of a maltose pool and of maltase in poplar wood (Populus x canadensis robusta) during fall. J. Plant Physiol. 1993, 141, 248–250. [Google Scholar] [CrossRef]
  72. Sharkey, T.D.; Laporte, M.; Lu, Y.; Weise, S.; Weber, A.P.M. Engineering plants for elevated CO2: A relationship between starch degradation and sugar sensing. Plant Biol. 2004, 6, 280–288. [Google Scholar] [CrossRef]
  73. Sun, X.; Zhang, J.; Fan, Z.H.; Xiao, P.; Liu, S.N.; Li, R.P.; Zhu, W.B.; Huang, L. MAL62 overexpression enhances freezing tolerance of baker’s yeast in lean dough by enhancing Tps1 activity and maltose metabolism. J. Agric. Food Chem. 2019, 67, 8986–8993. [Google Scholar] [CrossRef]
  74. Adam, S.; Murthy, S.D.S. Effect of cold stress on photosynthesis of plants and possible protection mechanisms. In Approaches to Plant Stress and Their Management; Gaur, R., Sharma, P., Eds.; Springer: New Delhi, India, 2014; pp. 219–226. [Google Scholar]
  75. Taïbi, K.; Del Campo, A.D.; Vilagrosa, A.; Bellés, J.M.; López-Gresa, M.P.; López-Nicolás, J.M.; Mulet, J.M. Distinctive physiological and molecular responses to cold stress among cold-tolerant and cold-sensitive Pinus halepensis seed sources. BMC Plant Biol. 2018, 18, 236. [Google Scholar] [CrossRef] [Green Version]
  76. Aranega-Bou, P.; de la O Leyva, M.; Finiti, I.; García-Agustín, P.; González-Bosch, C. Priming of plant resistance by natural compounds. Hexanoic acid as a model. Front. Plant Sci. 2014, 5, 488. [Google Scholar] [CrossRef] [Green Version]
  77. Conrath, U. Priming of induced plant defense responses. Adv. Bot. Res. 2009, 51, 361–395. [Google Scholar]
  78. NDong, C.; Danyluk, J.; Wilson, K.E.; Pocock, T.; Huner, N.P.; Sarhan, F. Cold-regulated cereal chloroplast late embryogenesis abundant-like proteins. Molecular characterization and functional analyses. Plant Physiol. 2002, 129, 1368–1381. [Google Scholar] [CrossRef] [Green Version]
  79. Minami, A.; Nagao, M.; Arakawa, K.; Fujikawa, S.; Takezawa, D. Abscisic acid-induced freezing tolerance in the moss Physcomitrella patens is accompanied by increased expression of stress-related genes. J. Plant Physiol. 2003, 160, 475–483. [Google Scholar] [CrossRef]
  80. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef]
  81. Howe, E.; Holton, K.; Nair, S.; Schlauch, D.; Sinha, R.; Quackenbush, J. Mev: Multiexperiment viewer. In Biomedical Informatics for Cancer Research; Ochs, M.F., Casagrande, J.T., Davuluri, R.V., Eds.; Springer: Boston, MA, USA, 2010; pp. 267–277. [Google Scholar]
  82. Kobylińska, A.; Borek, S.; Posmyk, M.M. Melatonin redirects carbohydrates metabolism during sugar starvation in plant cells. J. Pineal Res. 2018, 64, e12466. [Google Scholar] [CrossRef]
  83. Menezes, F.O.; Carvalho, V.; Moreira, V.A.; Rigui, A.P.; Gaspar, M.; Nievola, C.C. Juvenile plants of an epiphytic bromeliad exhibit rapid changes in carbohydrate content and antioxidant capacity to cope with water withholding. Theor. Exp. Plant Physiol. 2020, 32, 89–98. [Google Scholar] [CrossRef]
  84. Maturilli, M.; Herber, A.; König-Langlo, G. Climatology and time series of surface meteorology in Ny-Ålesund, Svalbard. Earth Syst. Sci. Data 2013, 5, 155–163. [Google Scholar] [CrossRef] [Green Version]
Plants 12 01250 g001 550
Figure 1. Workflow diagram summarizing the analysis process used in this study.
Figure 1. Workflow diagram summarizing the analysis process used in this study.
Plants 12 01250 g001
Plants 12 01250 g002 550
Figure 2. Physiological changes of A. turgidum protonema induced by cold acclimation and freezing stress. Freezing stress was applied to A. turgidum protonema at a rate of decreasing the temperature of 2 °C per hour. Electrolyte leakage change under freezing stress for CA and NA (A). The maximum photochemical efficiency of photosystem II (FV/FM) values of CA-12 and NA-12 during recovery from freezing (B). Data are means ± SD of three repeats (***, p < 0.001; **, p < 0.01, Student t-test). CA, cold-acclimated A. turgidum at 4 °C for 2 days; NA, non-acclimated A. turgidum at 25 °C for 2 days; CA-12, CA exposed to freezing until −12 °C; NA-12, NA exposed to freezing until −12 °C.
Figure 2. Physiological changes of A. turgidum protonema induced by cold acclimation and freezing stress. Freezing stress was applied to A. turgidum protonema at a rate of decreasing the temperature of 2 °C per hour. Electrolyte leakage change under freezing stress for CA and NA (A). The maximum photochemical efficiency of photosystem II (FV/FM) values of CA-12 and NA-12 during recovery from freezing (B). Data are means ± SD of three repeats (***, p < 0.001; **, p < 0.01, Student t-test). CA, cold-acclimated A. turgidum at 4 °C for 2 days; NA, non-acclimated A. turgidum at 25 °C for 2 days; CA-12, CA exposed to freezing until −12 °C; NA-12, NA exposed to freezing until −12 °C.
Plants 12 01250 g002
Plants 12 01250 g003 550
Figure 3. E-value distribution of the top five BLASTx hits for each unigene (less than 1.0 × 10−4) (A) and Nr annotation species distribution analysis (B).
Figure 3. E-value distribution of the top five BLASTx hits for each unigene (less than 1.0 × 10−4) (A) and Nr annotation species distribution analysis (B).
Plants 12 01250 g003
Plants 12 01250 g004 550
Figure 4. Heatmap of differentially expressed unigenes of CA-12 compared to NA-12. Using RNA-seq data analysis under the criteria of p-value < 0.05 and a log2 ratio of <−1 or >1 of NA-12 versus CA-12, we identified 4044 differentially expressed unigenes (A). In the left panel, red indicates upregulation in NA-12; green indicates upregulation in CA-12. The right panel shows average normalized FPKM values from RNA-seq experiments; blue indicates the lowest expression level, and yellow is the highest. Detailed data from the RNA-seq analysis are presented in Table S4. The association with the heatmap was further confirmed by monitoring eight genes using qPCR (B). The y-axis indicates the expression level relative to EF1/unigene_08701 (internal control), and the x axis indicates samples used for qPCR. Data are means ± SD of three repeats (***, p < 0.001; **, p < 0.01, Student t-test). The primer sequences used for qPCR analysis are shown in Table S5. CA-12, A. turgidum protonema were cold-acclimated at 4 °C for 2 days and exposed to freezing until −12 °C; NA-12, A. turgidum protonema were non-acclimated at 25 °C for 2 days and exposed to freezing until −12 °C.
Figure 4. Heatmap of differentially expressed unigenes of CA-12 compared to NA-12. Using RNA-seq data analysis under the criteria of p-value < 0.05 and a log2 ratio of <−1 or >1 of NA-12 versus CA-12, we identified 4044 differentially expressed unigenes (A). In the left panel, red indicates upregulation in NA-12; green indicates upregulation in CA-12. The right panel shows average normalized FPKM values from RNA-seq experiments; blue indicates the lowest expression level, and yellow is the highest. Detailed data from the RNA-seq analysis are presented in Table S4. The association with the heatmap was further confirmed by monitoring eight genes using qPCR (B). The y-axis indicates the expression level relative to EF1/unigene_08701 (internal control), and the x axis indicates samples used for qPCR. Data are means ± SD of three repeats (***, p < 0.001; **, p < 0.01, Student t-test). The primer sequences used for qPCR analysis are shown in Table S5. CA-12, A. turgidum protonema were cold-acclimated at 4 °C for 2 days and exposed to freezing until −12 °C; NA-12, A. turgidum protonema were non-acclimated at 25 °C for 2 days and exposed to freezing until −12 °C.
Plants 12 01250 g004
Plants 12 01250 g005 550
Figure 5. List of up-regulated unigenes encoding transcription factors in the CA-12 treatment. The average FPKM normalization and log2 fold change values for the identified 25 unigenes are described. In the log2 fold change column, red indicates the highest value, and white indicates the lowest value. CA-12, A. turgidum protonema were cold-acclimated at 4 °C for 2 days and exposed to freezing until −12 °C; NA-12, A. turgidum protonema were non-acclimated at 25 °C for 2 days and exposed to freezing until −12 °C.
Figure 5. List of up-regulated unigenes encoding transcription factors in the CA-12 treatment. The average FPKM normalization and log2 fold change values for the identified 25 unigenes are described. In the log2 fold change column, red indicates the highest value, and white indicates the lowest value. CA-12, A. turgidum protonema were cold-acclimated at 4 °C for 2 days and exposed to freezing until −12 °C; NA-12, A. turgidum protonema were non-acclimated at 25 °C for 2 days and exposed to freezing until −12 °C.
Plants 12 01250 g005
Plants 12 01250 g006 550
Figure 6. Gene ontology (GO) analysis in the “biological processes” category for unigenes upregulated in response to CA-12 treatment. Thirteen GO terms were over-represented, with the rest collectively indicated as “others”. In the circle graph, the percentage occupied by each GO term is expressed as “%”. Fifteen unigenes encoding pentatricopeptide repeat (PPR) protein were identified in “response to stimulus” (6%), which is a distinguished feature for A. turgidum. Details of GO information for groups upregulated in CA-12 treatment are presented in Table S6.
Figure 6. Gene ontology (GO) analysis in the “biological processes” category for unigenes upregulated in response to CA-12 treatment. Thirteen GO terms were over-represented, with the rest collectively indicated as “others”. In the circle graph, the percentage occupied by each GO term is expressed as “%”. Fifteen unigenes encoding pentatricopeptide repeat (PPR) protein were identified in “response to stimulus” (6%), which is a distinguished feature for A. turgidum. Details of GO information for groups upregulated in CA-12 treatment are presented in Table S6.
Plants 12 01250 g006
Plants 12 01250 g007 550
Figure 7. Concentration changes of starch and sugars in A. turgidum protonema of NA-12 and CA-12. Starch (A), Maltose (B), Sucrose (C), Glucose (D) and Fructose (E). The y-axis scales of each graph are different. CA-12, A. turgidum protonema were cold-acclimated at 4 °C for 2 days and exposed to freezing until −12 °C; NA-12, A. turgidum protonema were non-acclimated at 25 °C for 2 days and exposed to freezing until −12 °C. Data are means ± SD of three repeats (***, p < 0.001; **, p < 0.01; *, p < 0.05, Student t-test).
Figure 7. Concentration changes of starch and sugars in A. turgidum protonema of NA-12 and CA-12. Starch (A), Maltose (B), Sucrose (C), Glucose (D) and Fructose (E). The y-axis scales of each graph are different. CA-12, A. turgidum protonema were cold-acclimated at 4 °C for 2 days and exposed to freezing until −12 °C; NA-12, A. turgidum protonema were non-acclimated at 25 °C for 2 days and exposed to freezing until −12 °C. Data are means ± SD of three repeats (***, p < 0.001; **, p < 0.01; *, p < 0.05, Student t-test).
Plants 12 01250 g007
Table
Table 1. Sequencing output statistics of A. turgidum from this study.
Table 1. Sequencing output statistics of A. turgidum from this study.
AccessionNo. Raw ReadsNo. Trimmed Reads% of Reads
after Trimming		Size of Trimmed Reads (bp)GC (%)Q30 (%)
NA-12_13,529,1803,305,32093.66814,851,89449.1393.77
NA-12_23,759,0023,536,72694.09879,302,30749.6393.84
NA-12_33,446,1143,202,40492.93789,298,95149.593.61
CA-12_13,636,5523,421,01094.07848,381,53849.2493.8
CA-12_23,241,5443,006,90692.76732,304,41349.2293.52
CA-12_33,727,5523,518,61094.39874,412,69049.6693.86
Table
Table 2. Summary of contigs and unigenes.
Table 2. Summary of contigs and unigenes.
ContigUnigene
Total Length (bp)157,581,30324,753,852
No. Sequence116,97945,796
Max. Length (bp)19,62417,139
Mean Length (bp)1347541
N50 (bp)24811017
N90 (bp)507207
N90 Sequence No.73,45731,863
GC%46.549.2
Table
Table 3. KEGG analysis of DEGs upregulated in CA-12.
Table 3. KEGG analysis of DEGs upregulated in CA-12.
ClassPathwayIDDEGAll Genesp-ValueCorrected
p-Value
Category/Metabolism
Carbohydrate metabolismStarch and sucrose metabolismppp00500321981.46 × 10−114.57 × 10−10
Glycolysis/gluconeogenesisppp00010261724.06 × 10−97.63 × 10−8
Amino sugar and nucleotide sugar metabolismppp00520231375.49 × 10−98.61 × 10−8
Galactose metabolismppp0005213522.32 × 10−73.12 × 10−6
Pentose phosphate pathwayppp0003016863.36 × 10−73.91 × 10−6
Butanoate metabolismppp006509241.06 × 10−69.93 × 10−6
Fructose and mannose metabolismppp0005112881.46 × 10−47.64 × 10−4
Ascorbate and aldarate metabolismppp000538529.09 × 10−43.88 × 10−3
Pyruvate metabolismppp00620121119.86 × 10−44.03 × 10−3
Glyoxylate and dicarboxylate metabolismppp00630111072.27 × 10−38.54 × 10−3
Amino acid metabolismPhenylalanine metabolismppp0036012453.74 × 10−73.91 × 10−6
Alanine, aspartate, and glutamate metabolismppp0025012594.38 × 10−63.74 × 10−5
Tyrosine metabolismppp003508432.96 × 10−41.39 × 10−3
Lipid metabolismFatty acid biosynthesisppp0006111602.56 × 10−52.01 × 10−4
Glycerophospholipid metabolismppp0056413873.37 × 10−52.43 × 10−4
Cutin, suberine, and wax biosynthesisppp000735111.41 × 10−47.64 × 10−4
Biosynthesis of other secondary metabolitesPhenylpropanoid biosynthesisppp00940151174.34 × 10−52.84 × 10−4
Metabolism of other amino acidsTaurine and hypotaurine metabolismppp00430584.54 × 10−52.84 × 10−4
beta-Alanine metabolismppp004106362.76 × 10−39.97 × 10−3
Energy metabolismCarbon fixation in photosynthetic organismsppp00710151258.56 × 10−55.03 × 10−4
Metabolism of cofactors and vitaminsUbiquinone and other terpenoid-quinone biosynthesisppp001308401.92 × 10−49.48 × 10−4
Category/Cellular Processes
Transport and catabolismPeroxisomeppp0414611961.03 × 10−34.02 × 10−3
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

Kang, P.; Yoo, Y.-H.; Kim, D.-I.; Yim, J.H.; Lee, H. De Novo Transcriptome Assembly and Comparative Analysis of Differentially Expressed Genes Involved in Cold Acclimation and Freezing Tolerance of the Arctic Moss Aulacomnium turgidum (Wahlenb.) Schwaegr. Plants 2023, 12, 1250. https://doi.org/10.3390/plants12061250

AMA Style

Kang P, Yoo Y-H, Kim D-I, Yim JH, Lee H. De Novo Transcriptome Assembly and Comparative Analysis of Differentially Expressed Genes Involved in Cold Acclimation and Freezing Tolerance of the Arctic Moss Aulacomnium turgidum (Wahlenb.) Schwaegr. Plants. 2023; 12(6):1250. https://doi.org/10.3390/plants12061250

Chicago/Turabian Style

Kang, Pilsung, Yo-Han Yoo, Dong-Il Kim, Joung Han Yim, and Hyoungseok Lee. 2023. "De Novo Transcriptome Assembly and Comparative Analysis of Differentially Expressed Genes Involved in Cold Acclimation and Freezing Tolerance of the Arctic Moss Aulacomnium turgidum (Wahlenb.) Schwaegr" Plants 12, no. 6: 1250. https://doi.org/10.3390/plants12061250

APA Style

Kang, P., Yoo, Y. -H., Kim, D. -I., Yim, J. H., & Lee, H. (2023). De Novo Transcriptome Assembly and Comparative Analysis of Differentially Expressed Genes Involved in Cold Acclimation and Freezing Tolerance of the Arctic Moss Aulacomnium turgidum (Wahlenb.) Schwaegr. Plants, 12(6), 1250. https://doi.org/10.3390/plants12061250

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