Next Article in Journal
Foraging Site Selection of the Middle Spotted Woodpecker (Leiopicus medius L.) in Primeval Oak-Lime-Hornbeam Forest of the Białowieża National Park: Comparison of Breeding and Non-Breeding Seasons
Next Article in Special Issue
Identification and Characterization of circRNAs under Drought Stress in Moso Bamboo (Phyllostachys edulis)
Previous Article in Journal
Revealing Changes in the Stem Form and Volume Allocation in Diverse Boreal Forests Using Two-Date Terrestrial Laser Scanning
Previous Article in Special Issue
Transcriptome Analysis Reveals Regulatory Framework for Salt and Drought Tolerance in Hibiscus hamabo Siebold & Zuccarini
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 12
Issue 7
10.3390/f12070836
forests-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

Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake

by
Hao Chen
1,2,3,
Jiang Li
1,2,3,
Bingfa Qiu
4,
Yingwei Zhao
4,
Zhao Liu
1,2,3,
Jun Yang
1,2,3,* and
Xiangyang Kang
1,2,3
1
Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing 100083, China
2
National Engineering Laboratory for Tree Breeding, Beijing Forestry University, Beijing 100083, China
3
Beijing Laboratory of Urban and Rural Ecological Environment, Beijing Forestry University, Beijing 100083, China
4
Guangxi Dongmen Forest Farm, Fusui 532199, China
*
Author to whom correspondence should be addressed.
Forests 2021, 12(7), 836; https://doi.org/10.3390/f12070836
Submission received: 11 May 2021 / Revised: 19 June 2021 / Accepted: 22 June 2021 / Published: 25 June 2021
(This article belongs to the Special Issue Resistance Genetics of Forest Trees)
Browse Figures
Versions Notes

Abstract

:
Long non-coding RNA (lncRNA) plays an important regulatory role in plant growth and development, but its systematic identification and analysis in Eucalyptus has not yet been reported. Cold stress has a huge impact on the survival and yield of Eucalyptus seedlings, but the regulatory mechanism of lncRNA in Eucalyptus in response to cold stress is still unclear. In this study, the transcriptomes of young leaves of Eucalyptus urophylla S.T.Blake under low-temperature treatment and restoration were analyzed by RNA-seq. A total of 11,394 lncRNAs and 46,276 mRNAs were identified, of which 300 were differentially expressed lncRNAs (DE_lncRNAs) and 5606 were differentially expressed target genes of lncRNAs under cold stress, with the total number of target genes of DE_lncRNAs being 1681. A total of 677 differentially expressed transcription factors (TFs) were also identified, mainly including ERF, MYB and the NAC transcription factor family. Gene ontology (GO) analysis of the differentially expressed genes (DEGs) and target genes of DE_lncRNAs was mostly related to the response to cold stress and external stimuli. Furthermore, lncRNA–miRNA–mRNA regulatory networks were constructed, and 22 DE_lncRNAs were predicted to be targets or targeting mimics of 20 miRNAs. A qRT-PCR was used to verify the relative expression of genes in the regulatory EuGBF3-EUC_00002677-MSTRG.7690 network, and it matched the transcriptome data, indicating that it may play an important role in the response to cold stress in E. urophylla. This study provides a new insight into lncRNA and its regulatory network under abiotic stress, especially cold stress in E. urophylla.

1. Introduction

Eucalyptus is the general term for Eucalyptus L’Herit in the Myrtaceae family. It is usually a perennial evergreen tree, mainly distributed naturally in Australia and its northern islands [1]. Eucalyptus is the most important plantation tree species in the world, providing a large amount of wood for industrial production [2,3] and sufficient raw material for pulping and papermaking [4,5,6]. Eucalyptus oil contained in leaves also has a variety of biological activities and is an important and effective ingredient in pesticides and herbicides [7,8]. As the most significant source of forest products in the world, Eucalyptus has been widely planted in more than 200 countries or regions, with a cumulative planting area of more than 20 million hectares [9,10]. Since it mostly originates in areas with hot climates, Eucalyptus does not tolerate low temperatures well. Cold stress limits the distribution and productivity of Eucalyptus plantations. When the temperature drops to 8 °C or below, Eucalyptus globulus exhibits various symptoms of chilling injury due to its inability to adapt to cold stress [11]. Especially after afforestation in mountainous regions in spring, occasional low temperatures may damage the buds and young leaves of seedlings, ultimately causing afforestation to fail. Therefore, research to reveal the cytological and molecular mechanisms of Eucalyptus in response to cold stress can provide important theoretical guidance for its future resistance breeding by molecular design or gene editing.
Cold stress is a kind of abiotic stress which is usually not conducive to plant growth and development [12]. In recent years, the signal pathways in plants that respond to external stress have been revealed. Many studies have shown that plants can respond to a short period of low but non-freezing temperatures through a complex transcriptional regulatory network which increases the plant’s cold tolerance. This process is called cold acclimation [13]. It changes the physiological state by altering the membrane lipid composition, photosynthetic pigments, affinity permeants (such as soluble sugars, betaine and proline) and antioxidant levels and by causing extensive rearrangement of gene expression [14,15,16,17,18,19]. The physiological and biochemical changes that occur during plant cold acclimation are mainly due to the fluctuations in cold-responsive gene expression [19]. For example, atypical R2R3 MYB transcription factors positively regulate cold hardiness and cold-responsive gene expression under cold stress and then increase the cold tolerance of apples through CBF-dependent and CBF-independent pathways [20]. After low-temperature treatment of two pepper varieties with different cold tolerances, ERF family genes had a specific regulatory effect on the cold response pathways through transcriptome sequencing [21]. miRNAs and lncRNAs are both involved in cold stress responses. For instance, overexpressing miR397a in Arabidopsis thaliana (L.) Heynh. increases the tolerance to cold stress [22], and a cold-responsive lncRNA named SVALKA plays a role in regulating CBF1 expression [23]. However, there is no detailed study on how lncRNA interacts with mRNA and miRNA to respond to cold stress.
In the past decade, many studies have shown that the genomes of plants are generally transcribed, while some transcripts do not encode proteins, producing a large number of long non-coding RNAs [24]. Long non-coding RNAs (lncRNAs) characterize a highly heterogeneous class of RNA molecules defined as transcripts, with over 200 nucleotides not translated into proteins [25]. lncRNAs contain natural anti-sense RNA (NATs) of protein-coding genes, small RNA precursors (microRNA, small nucleolar RNA/SnoRNA and other small regulatory RNA) and structural RNAs (tRNA and ribosomal RNA) [26]. They can originate from exonic, intronic, intragenic, intergenic, promoter regions, 3′- and 5′-UTR or enhancer sequences transcribed in either the sense or the antisense strand [27]. Based on genomic locations, lncRNAs can generally be divided into sense lncRNAs, antisense lncRNAs, bidirectional lncRNAs, intergenic lncRNAs and intronic lncRNAs [27]. Long non-coding RNAs act as significant managers in gene expression networks by governing structure and transcription in the nucleus and modulating mRNA steadiness, translation and post-translational modifications in the cytoplasm [28]. Because of the continuous developments in high-throughput sequencing technology, lncRNAs of various species have been identified, such as Arabidopsis [29], rice [30], maize [31], poplar [32] and Eucalyptus grandis W.Hill [33]. lncRNA in plants are involved in responses to salinity, osmotic stress, drought, heat stress and other stress. For example, lncRNAs are differentially and specifically expressed in roots and leaves in response to salinity in A. thaliana [34]. There are also studies on the response of lncRNA to low-temperature stress. For example, some cold-responsive lncRNAs were identified in Medicago truncatula Gaertn. seedlings, and potential regulatory networks between lncRNA MtCIR1 and MtCBFs were discovered [35]. However, no study has provided a comprehensive and detailed explanation of the role of lncRNAs in regulating the response to cold stress in plants.
Recently, many studies have investigated the physiological and gene expression changes in Eucalyptus in a low-temperature environment. Oberschelp et al. [36] conducted metabolomics and proteomics studies to compare cold acclimation and freezing tolerance in three Eucalyptus species and found that E. benthamii showed a higher cold tolerance than E. dunnii and E. grandis, possibly due to its higher accumulation of phenolics, anthocyanins and soluble sugars, as well as lower levels of photosynthetic pigments and related proteins. Gaete-Loyola et al. [37] used transcriptome analysis to reveal the molecular mechanism of cold acclimation and deacclimation in E. nitens, providing a comprehensive view of the total gene expression and revealing meaningful information about the dynamic and complex nature of gene expression.
Eucalyptus urophylla S.T.Blake is widely distributed in disjunct populations that occupy the islands of Indonesia and Timor [38]. Eucalyptus urophylla is fast-growing and easy to cross-breed. It is not only a vital cultivated tree species in China and other Asia-Pacific regions, but it is also an important parent of high-yielding hybrids such as Eucalyptus urophylla x E. grandis and Eucalyptus grandis x E. urophylla. This species grows well and is widely planted in China, especially southern China. It provides an indispensable source for the production of paper, pulp, cellulose products and lignocellulosic biofuels. The study revealed that the cold resistance mechanism of E. urophylla is important for the genetic improvement of Eucalyptus.
However, lncRNAs and mRNAs under cold stress in E. urophylla have not been identified and characterized, and the mechanisms by which lncRNA and mRNA participate in cold tolerance are still unclear. Therefore, this study subjects E. urophylla to low-temperature treatment to simulate damage of its young leaves by cold stress and uses next-generation DNA sequencing (NGS) to perform RNA-seq, exploring the response of lncRNA and mRNA to cold stress. Furthermore, an lncRNA–miRNA–mRNA regulatory network under cold stress in E. urophylla is constructed, and it provides a new perspective for Eucalyptus to better adapt to abiotic stress.

2. Materials and Methods

2.1. Plant Materials

Eucalyptus urophylla plants preserved in the Guangxi Dongmen Forest Farm were used. In vitro plantlets produced by tissue culture were transplanted to pots and cultivated for about 2 months in a greenhouse under natural light conditions (1250 µmol m−2s−1 of photosynthetically active radiation) at 25 ± 1 °C and 50 ± 1% relative humidity, with a 16/8 h light/dark cycle. For cold treatment, the seedlings were transferred to a chamber at 4 °C for different time periods: LT6, LT12, LT24 and LT48, referring to low-temperature treatment for 6, 12, 24 and 48 h, respectively, and RE24, referring to returning to a normal temperature of 25 °C for 24 h after 48 h of low-temperature treatment. Some seedlings were moved out of the chamber and used as the non-treatment control group (CK). After each treatment, young leaves were collected, immediately frozen in liquid nitrogen and stored at −80 °C until use. Three biological replicates were set for each treatment and sampling.

2.2. Measurement of Fv/Fm and Relative Electrolyte Leakage

After 20 min of dark adaptation, the first through third fully expanded young leaves were used to measure the maximal quantum efficiency of photosystem II (PSII) (Fv/Fm) [39]. These values were automatically monitored and calculated using a chlorophyll fluorometer (MultispeQ, East Lansing, MI, USA). Subsequently, 0.1 g of a young leaf cut into small pieces of 0.5 cm2 was used to measure the relative electrolyte leakage (REL) using a conductivity meter (Leici-DDS-307A, Shanghai, China). After 6 h of shock at 160 rpm, leaf samples packed in the centrifuge tube containing 10 mL of deionized water were used to determine the initial conductivity (REL1). Then, these leaf samples were boiled for about 20 min, and the final conductivity (REL2) was measured after cooling to a normal temperature. The REL value was calculated as (REL1/REL2) × 100% [40].

2.3. RNA Extraction and Sequencing

The total RNA was extracted and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA), following the manufacturer’s instructions. The total RNA quantity and purity were analyzed using Bioanalyzer 2100 and the RNA 1000 Nano LabChip Kit (Agilent, CA, USA) with an RIN >7.0. About 5 µg of the total RNA was used to deplete the ribosomal RNA according to the manual of the Ribo-Zero™ rRNA Removal Kit (Illumina, San Diego, USA). The enriched circRNAs were broken into small pieces with divalent cations in a high-temperature environment. After removing ribosomal RNAs, the remaining RNAs were reverse-transcribed to create cDNA, which was used to synthesize U-labeled, second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base was then added to the blunt ends of each strand, preparing them for ligation to the indexed adapters. Each adapter contained a T-base overhang for ligating the adapter to the A-tailed fragmented DNA. Single- or dual-index adapters were ligated to the fragments, and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled, second-stranded DNAs, the ligated products were amplified by PCR under the following conditions: initial denaturation at 95 °C for 3 min, 8 cycles of denaturation at 98 °C for 15 s, annealing at 60 °C for 15 s, extension at 72 °C for 30 s and final extension at 72 °C for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). Finally, paired-end sequencing was performed on an Illumina Novaseq™ 6000, following the vendor’s recommended protocol.

2.4. Transcript Assembly

Cutadapt [41] was used to remove the reads that contained adapter contamination, low-quality bases and undetermined bases. The clean reads obtained were aligned to the Eucalyptus genome using Hisat2 [42]. Transcriptome assembly and expression quantification were performed using StringTie [43].

2.5. Identification of lncRNA

Assembled transcripts with overlapping of known mRNAs or that were shorter than 200 bp were discarded. CPC and CNCI were used to predict the coding potential of the remaining transcripts with default parameters. The transcripts with a CPC score ≥−1 and a CNCI score ≥0 were considered lncRNAs.

2.6. Different Expression Analysis of mRNAs and lncRNAs

Edger R was used to identify differentially expressed mRNA or lncRNA genes. Genes with |log2 (fold change)| ≥ 1 and p < 0.05 were considered differentially expressed genes.

2.7. Target Gene Prediction of lncRNAs with Protein-Coding Genes

In this study, coding genes in 100,000 upstream and downstream components were selected by a Python script.

2.8. Functional Analysis of DEGs and DE_lncRNAs:

Functional annotation of differentially expressed lncRNA loci was based on their target genes. The target genes and differentially expressed genes were annotated based on the GO database (http://www.geneontology.org; accessed on 12 July 2020).

2.9. Identification of Transcription Factors

The Plant Transcription Factor Database was used to identify and classify TFs from differentially expressed genes.

2.10. Prediction of lncRNAs as miRNA Target Mimics

Eucalyptus grandis miRNA sequences were downloaded from a previous study [44]. The target mimics were predicted by submitting the lncRNAs and miRNAs to psRNATarget (http://plantgrn.noble.org/psRNATarget; accessed on 30 July 2020). Less than three mismatches and G/U pairs were allowed within the lncRNA and miRNA pairing regions [45]. Cytoscape was used to establish the lncRNA–miRNA–mRNA interaction network.

2.11. Quantitative Real-Time PCR

The total RNA was isolated from the young leaves of E. urophylla using an RNAprep Pure Plant Plus Kit (Tiangen, China, cat DP441). About 1.5 μg of RNA was reverse-transcribed into first-stand cDNA using a cDNA synthesis kit (Tiangen, China, cat KR106). qPCR was subsequently performed using a TransStart® Tip Green qPCR SuperMix (TransGen Biotech, Beijing, China) with EuGADPH as a reference gene. qPCR on the Applied Biosystems 7500 FAST real-time PCR system according to the manufacturer’s instructions. The primers used for qRT-PCR analysis are listed in Table S1.

2.12. Statistical Analysis

The Statistical Product and Service Solutions (SPSS) 20 program was used to conduct all statistical analyses of data, and data were presented as the mean ± SD (standard deviation). The significant differences between means (p < 0.05) were examined by a t-test and analysis of variance (ANOVA). The Waller-Ducan test was used for post hoc comparison after ANOVA.

3. Results

3.1. The Effect of Cold Stress on the Phenotype of Eucalyptus urophylla

To analyze the effect of cold stress on E. urophylla, the leaf morphology, relative electrolyte leakage (REL) and maximum photochemical efficiency of PSII (Fv/Fm) under different durations of low-temperature treatment (4 °C) were carefully observed. Compared with the control group (CK), after 6, 12, 24 and 48 h of low-temperature treatment (LT6, LT12, LT24 and LT48, respectively), the young leaves had a different extent of wilting (Figure 1a). In the first 6 h of low-temperature treatment (LT6), the morphological changes in the young leaves were not obvious. As the duration of low-temperature treatment increased from 12 h to 48 h (LT12, LT24 and LT48), the young leaves showed different levels of curling and wilting. After 48 h of treatment, the plants were placed in a normal-temperature environment for 24 h of recovery (RE24), but the leaves did not return to normal. The REL value is a key indicator of the severity of cell penetration. Compared with low temperature treatment for 6 h, the REL value when treated for 12 h did not change significantly (Figure 1c). After that, they gradually increased from 12.67% to 57.79%, especially in the 48 h and later treatment groups (LT48 and RE24), and they were significantly higher than in other treatments. Chlorophyll fluorescence is a key parameter for evaluating photosynthesis. The Fv/Fm value is used to represent the maximum photochemical efficiency of PSII, which is often used to determine the response degree of plants to various stress treatments. Compared with the control group (CK), Fv/Fm decreased from 0.68 to 0.56 in the initial 6 h of low-temperature treatment (LT6) (Figure 1b) and then rose slightly after adaptation at 12 h of treatment (LT12). With the increase in the cold stress duration, Fv/Fm further gradually decreased to 0.53 at 48 h of treatment (LT48). After that, with 24 h of recovery (RE24), Fv/Fm directly decreased to 0.09, which was significantly lower than in other treatments. The phenotype of E. urophylla underwent major changes under different durations of low-temperature treatment.

3.2. Genome-Wide Identification and Characterization of lncRNAs and mRNAs of E. urophylla under Cold Stress

Young leaves of E. urophylla from six groups were collected, including one control group (CK), four low-temperature treatment groups chilled at 4 °C for 6, 12, 24 and 48 h (LT6, LT12, LT24 and LT48, respectively) and one recovery group (recovered to a normal-temperature environment for 24 h after 48 h of a low-temperature treatment (RE24)). RNA sequencing from 18 cDNA libraries with three biological replicates from these six groups was conducted to systematically identify the lncRNA response to cold stress. Low-quality reads were abandoned from the RNA-seq data, and adapter sequences from residual reads were trimmed to obtain clean reads (Table 1). We obtained 46,276 mRNAs, and the cDNA of each sample from the six groups produced an average of 53,175,489, 52,968,464, 46,365,606, 47,024,170, 50,875,522 and 39,763,655 raw reads, respectively, and 52,294,501, 52,130,830, 45,584,262, 46,220,078, 49,941,482 and 39,114,580 clean reads, respectively. The average mapping rate to the E. grandis genome was more than 89%, the Q30 value was more than 98%, and the GC content was above 49%, which indicated that the transcriptome data were reliable.
A series of screening processes was followed to identify the lncRNAs from CK, LT6–LT48 and RE24 groups (Figure 2). Numerous candidate lncRNAs were predicted using the Coding Potential Calculator (CPC) and Non-Coding Index (CNCI) databases, including 433 intronic lncRNAs, 188 bidirectional lncRNAs, 1429 sense lncRNAs, 7407 intergenic lncRNAs and 1937 antisense lncRNAs (Figure 3a). The mRNAs and lncRNAs involved in cold stress were characterized according to five aspects: the exon number, ORF length, isoform number per gene, transcript length and differential expression level. Most of the mRNAs and lncRNAs encompassed one or two exons, accounting for 34% and 91%, respectively. The maximum number of exons of lncRNAs was five, while the number of exons in mRNAs was in a wider distribution range (Figure 3b). The whole transcript length and ORF length of the lncRNAs were apparently much shorter than those of the mRNAs. The number of lncRNAs with ORF lengths less than 100 nt was 2081, occupying more than 83%, while the number of mRNAs with ORF lengths less than 100 nt was 2902. The ORF lengths of the mRNAs were primarily distributed between 100 and 500 nt. The longest ORF length of the mRNAs was more than 1200 nt, while the longest ORF length of the lncRNAs was shorter than 600 nt (Figure 3c). In terms of the transcript length, most of the mRNA transcripts were longer than 1000 nt, accounting for about 69%, while lncRNA transcripts not exceeding 300 nt accounted for a maximum ratio of 46% (Figure 3e). The expression levels of the mRNAs and lncRNAs were also compared using the fragments per kilobase of exon per million fragments mapped (FPKM). The lncRNAs and mRNAs were expressed at similar levels in all six groups (Figure 3f). However, the overall expression levels of the mRNAs were generally higher than those of lncRNAs. Around 8% of the lncRNAs were alternatively spliced, while about 23% of the mRNAs were spliced (Figure 3d).
Three lncRNAs and three mRNAs from the treatment and recovery groups were randomly selected and analyzed by qRT-PCR in order to explore their expression patterns. The results showed that the expression patterns of lncRNAs and mRNAs involved in cold stress were relatively consistent and comparable (Figure 4), suggesting that their expression based on RNA-seq data was reliable.

3.3. Identification of DEGs and DE_lncRNAs

To understand gene expression under cold stress, the transcript levels of the lncRNAs and mRNAs in terms of fragments per kilobase of exon per million fragments mapped (FPKM) based on RNA-Seq data were measured (Figure S1, Tables S2 and S3). The conditions to select differentially expressed mRNAs (DEGs) and differentially expressed lncRNAs (DE_lncRNAs) were set as |log2 (fold change)| ≥ 1 and p < 0.05. The heat maps of the mRNAs and lncRNAs demonstrated gene expression after normalizing the data of all six groups. The results showed significant differences in gene expression in E. urophylla under cold stress, especially comparing the CK group with the LT6, LT12, LT24, LT48 and RE24 groups, indicating a rapid response to low-temperature stimulus of gene regulation. The number of DEGs (DE_lncRNAs) of five comparison groups (LT6 vs. CK, LT12 vs. CK, LT24 vs. CK, LT48 vs. CK and RE24 vs. CK) were 1893 (57), 2108 (85), 3076 (159), 4707 (232) and 3123 (45), respectively. The expression patterns of the mRNAs and lncRNAs in each treatment group were generally similar, except in the LT6 and LT12 groups (Figure 5c,d). Compared with the control group, with the increase in cold stress time, the number of DEGs and DE_lncRNAs gradually increased. After 24 h of recovery, the number of DEGs and DE_lncRNAs decreased to a certain extent. Pairwise comparisons were performed among the treatment groups and the recovery group, and a large number of DEGs and DE_lncRNAs were obtained. Compared with the LT48 group, the RE24 group had 4359 DEGs and 74 DE_lncRNAs, suggesting the recovery process had a vital influence on the biological process in E. urophylla after low-temperature treatment. Detailed analysis of differentially expressed genes showed that the number downregulated DEGs was more than the upregulated DEGs with the increase in treatment time (Figure 5a), while the upregulated DE_lncRNA were mostly greater than the downregulated ones (Figure 5b). These results show that cold stress caused different changes in the mRNAs and lncRNAs in E. urophylla.

3.4. Prediction of Target Genes of lncRNAs

Since lncRNA plays an important role in regulating gene function, identification and analysis of their target genes could be helpful to explore the function of lncRNAs. In this study, the number of target genes was 19,418, of which 5606 were differentially expressed target genes, while 300 were DE_lncRNAs and 579 were differentially expressed target genes of DE_lncRNA (Table 2).

3.5. Identification of Transcription Factors in E. urophylla

A total of 677 differentially expressed TFs under cold stress were identified in each group (Table S4). These TFs belonged to 48 families, of which the top 12 TF families were ERF, MYB, NAC, bHLH, C2H2, WRKY, GRAS, HD-ZIP, G2-like, bZIP, Dof and MYB_related, consisting of 77, 66, 64, 53, 42, 41, 26, 23, 22, 21, 21 and 21 genes, respectively. Most of these gene families are related to the response to cold stress [46] (Figure 6). Among them, the TF family with the maximum differentially expressed genes was ERF, another subfamily of the AP2/ERF TF family that plays a vital role in the biological regulation of the response to stress [47,48]. Among these differentially expressed TFs, 14 were the target genes of DE_lncRNAs, belonging to the ERF, MYB_related, NAC, C3H, CO_like and WRKY TF families (Table S4). Of these 14, 4 differentially expressed TFs were the target genes of MSTRG.6543, namely Eucgr.B02398, Eucgr.B02399, Eucgr.B02396 and Eucgr.B02393, which all belonged to the C2H2 TF family.

3.6. Functional Analysis of DEGs and Target Genes of DE_lncRNAs

Gene ontology (GO) was used to evaluate the functions of DEGs and was divided into three parts: biological process, cell component and molecular function. To further analyze the effect of cold stress on the gene expression of E. urophylla, GO enrichment analysis of the DEGs and target genes of DE_lncRNAs was performed (Figure 7). The DEGs were mainly enriched in four GO terms in the cell component process under low-temperature stress in E. urophylla: anchored component of membrane, kinesin complex, nucleosome and cell surface. In general, the DEGs were mainly enriched in parts of the biological processes and molecular functions. Most of the enrichment pathways of the biological processes were related to the cold stress response, such as the response to cold being enriched within 134 genes. The pathway with most DEGs enriched was the response to salt stress, which was enriched within 164 genes. The results also show that the pathways including the response to absicic acid, response tosalicylic acid, enthylene-activated signaling pathway, response to ethylene and response to jasmonic acid were enriched within 160, 100, 99, 62 and 61 genes, respectively. In the GO terms of molecular functions, the most enriched pathways were reansferase activity, transferring glycosyl groups, protein kinase activity, quercetin 3-O-glucosyltransferase activity and quercetin 7-O-glucosyltransferase activity.
The GO enrichment of the target genes of DE_lncRNAs was similar to that of DEGs, with only one pathway, the cortical actin cytoskeleton, enriched in the GO terms of cellular components, which were also mainly enriched in biological processes and molecular functions. Most of the biological process enrichment pathways were related to cold stress and external stimuli, including the response to cold, cellular response to extracellular stimulus, cellular response to external stimulus, response to extracellular stimulus, cellular response to abiotic stimulus, cellular response to environmental stimulus and cellular response to cold. The cold-related and drought- and water-deficiency pathways were also enriched, such as the response to water deprivation, response to water and cellular response to water deprivation. The target genes of DE_lncRNAs were also enriched in the area of molecular functions, and the number of pathways enriched within more than 20 genes was 6. The most enriched pathway was protein homodimerization activity, containing 68 genes, while the second- and third-most enriched pathways were transition metal ion binding and lipase activity, including 41 and 38 genes, respectively.

3.7. DEG, DE_lncRNA, Differentially Expressed TF and miRNA Interaction Network Construction

A large number of lncRNAs that could be used as targets or target mimics of miRNA were recognized. Among these 300 cold stress responsive DE_lncRNAs, 22 were identified as targets or target mimics of 20 miRNAs, belonging to the miR166, miR167 and miR482 families containing 27 targeting relationships (Table S5). In the prediction of the miRNA targeting relationship, the results also show the interaction between the target genes of the DE_lncRNAs and miRNAs. A total of 11 target genes of the DE_lncRNAs were identified as target mimics of 14 miRNAs.
To analyze the interaction between lncRNA and miRNA, transcriptional regulatory networks among the DEGs, DE_lncRNAs, differentially expressed TFs (DE_TFs) and miRNAs under low-temperature stress were constructed in E. urophylla (Figure 8). One of the regulatory networks showed that EUC_00002677 from the miR167d family had four target genes (Figure 8I), one of which was a DEG named Eucgr.J01413, while the others were all TF, including EuGNC (Eucgr.C03895), EuFRSII (Eucgr.K03522) and EuGBF3 (Eucgr.E000078), which belong to the GATA, FAR1 and bZIP TF families, respectively. EUC_00002677 can also form a targeting relationship with two lncRNAs, MSTRG.48389.1 and MSTRG.7690.1, which respond to cold stress. Another regulatory network showed that (Figure 8VI) miRNA EUC_00008890, belonging to the miR167d family, can target the Dof zinc finger gene EuDAG2 (Eucgr.B0245) and lncRNA MSTRG.46425, which respond to cold stress. The third regulatory network showed that MSTRG.26750 has two target DEGs, and it forms a targeting relationship with two miRNAs, EUC_00005453 (miR482e) and EUC_00006406 (miR447B-3P). Furthermore, EUC_00005453 (miR482e) could also bind to four target mRNAs, all of which were TFs, including three members from the NAC family and one from the SRS family (Figure 8II).
In addition, qRT-PCR was used to determine the relative gene expression in the mRNA–miRNA–lncRNA regulatory network of E. urophylla under low-temperature treatment to verify their interaction. EUC_00002677 from the miR167d family could target an lncRNA named MSTRG7690 and a transcription factor named EuGBF3 (Eucgr.E00078) from the bZIP family. The relative expression of the above three genes was measured to explore the relevance of their changing trends. The expression of miRNA EUC_00002677 was downregulated under low-temperature treatment, which promoted the upregulation of lncRNAs MSTRG.7690 and EuGBF3. The change was particularly obvious after 6 h of low-temperature treatment (LT6). The upregulated expression of MSTRG.7690 and EuGBF3 was nearly 9 and 1.5 times that of the control group, respectively, which indicates that lncRNAs may be more sensitive to cold stress. The expression of lncRNAs MSTRG.7690 and EuGBF3 was downregulated after 12 h of treatment (LT12), but overall, it was higher in the LT6, LT12, LT24 and LT48 groups than in the control group, while the expression of miRNA EUC_00002677 was lower than in the control group. In addition, at 24 h of recovery after the low-temperature treatment, EUC_00002677 was upregulated, and its expression was approximately twice that of the control group, while the expression of MSTRG.7690 and EuGBF3 was downregulated and was lower than that of the control group. In general, under cold stress, the expression of EUC_00002677 in E. urophylla was downregulated, which promoted the upregulation of MSTRG.7690 and EuGBF3, showing a negative correlation, which was in line with the ceRNA theory.

4. Discussion

Many studies have shown that lncRNAs can interact with DNA, RNA or TFs to alter the expression of target genes in order to regulate different biological processes [49]. With the continuous development of high-throughput sequencing technology, lncRNAs of many species of plants have been identified, such as Arabidopsis thaliana [29], rice [30], corn [31] and poplar [32]. Although some lncRNAs have been reported in Eucalyptus [33], the lncRNA response to cold stress and its mRNA–miRNA–lncRNA regulatory network have not been identified in E. urophylla. In this study, E. urophylla seedlings were subjected to low-temperature treatments (4 °C) for different durations. Through transcriptome sequencing, 11,394 lncRNAs and 46,276 mRNAs were identified, of which 300 were DE_lncRNAs responding to low temperatures, and 5606 were garget genes of lncRNAs. The DE_lncRNAs and DEGs were mainly enriched in the abiotic stress response, hormone signal response and carbohydrate biosynthetic processes, among others. In the newly constructed mRNA–miRNA–lncRNA regulatory network, 22 DE_lncRNAs were identified as miRNA targets. The miRNAs in E. urophylla negatively regulate the expression of lncRNA and mRNA to respond to a low-temperature environment. In general, the regulatory network in E. urophylla under low-temperature stress was explored, which could provide a molecular basis for the research on resistance to cold stress and improve its adaptation to an external environment.
Plant responses to cold stress involve multifaceted processes, such as changes in physiology, metabolism and macromolecular processes [50]. Cold stress can cause obvious damage to plants, including wilting, drying of leaf edges, discoloration, incomplete ripening, accelerated aging and even death [51]. The phenotype changes significantly after cold treatment. Young leaves are more sensitive to cold stress than mature leaves [52]. Phenotypic changes in young leaves began to appear in the LT12 group, and the leaves began to wither in the LT24 group, indicating that low temperatures damage the cell membrane, especially in young leaves, and further cause water shortage [53]. In addition, the REL value increased after the cold treatment, which confirmed that low temperatures destroy the structures of the cell membranes and further cause a large amount of extravasation of soluble substances. The Fv/Fm value significantly reduced after cold treatment, indicating that low temperature can destroy photosynthetic components, thereby affecting the growth and development of plants. Both REL and Fv/Fm showed obvious changes after 6 h of cold treatment and a certain degree of recovery after 12 h of cold treatment, which may be related to the balance of oxidative stress and the antioxidant system in the plant [53,54]. Plants can induce reactive oxygen species (ROS) in a low-temperature environment [55,56,57,58,59,60]. To cope with the massive production of ROS, plants can run an efficient purification system, including increasing the activity of non-enzymatic antioxidants and antioxidant enzymes [61,62]. When the temperature was restored to room temperature for 24 h after low-temperature treatment, neither the leaf state nor the REL or Fv/Fm recovered, which indicated that low temperatures cause irreversible damage to plant leaves, especially young leaves.
To understand the effect of cold stress on the regulation of gene expression in E. urophylla, GO enrichment analysis of DEGs and the target genes of DE_lncRNAs was performed. In the cell components, the GO term with the maximum DEG enrichment was “anchored component of membrane”, suggesting that the main effect of cold stress on the cells in E. urophylla is a change in the composition of the cell membrane, especially the saturation of constituent lipids [63], eventually causing young leaves to lose water and wither. The response of DE_lncRNAs to cold stress at the cellular level is mainly through the influence of the cortical actin cytoskeleton. In the biological process, in addition to the “response to cold” term, the DEGs and the target genes of DE_lncRNAs were mainly enriched in response to salt, drought, osmotic stress and plant hormones, indicating crosstalk in the response to different abiotic stresses in E. urophylla under cold stress. This is because when plants are under abiotic stresses, such as cold stress, signal transduction, such as the calcium ion signal transduction pathway, could play a corresponding role to create various responses to external stimuli [50,64]. Plant hormones also play an important role in the response to stimuli and regulate signal transduction pathways [65], such as jasmonic acid [66], ethylene [67] and abscisic acid [68]. Some DEGs were also enriched in the flavonoid biosynthetic process, because cold stress can regulate flavonoid synthesis pathways through transcriptome modification, including anthocyanin synthesis and metabolism pathways [69]. In addition, some DE_lncRNAs were enriched in the carbohydrate synthesis pathway. Notably, DE_lncRNAs related to the carbohydrate biosynthetic process under cold stress are likely beneficial to homeostasis between cellular membranes and ROS [70,71], which could partly explain the increase in Fv/Fm and the decrease in REL in the LT12 group compared with the LT6 group. Regarding molecular function, most of the DEGs were enriched in pathways related to enzyme activity, such as transferase activity, and glycosyl group transfer, which may be because E. urophylla accumulates carbohydrates through the activation of specific enzymes under cold stress [72,73], and these carbohydrates further act as osmotic regulators and signaling molecules [74]. The DE_lncRNAs were mostly enriched in protein homodimerization activity, which shows that cold stress had a significant impact on protein structure formation. These analyses indicate that large amounts of DEGs and DE_lncRNAs could respond to cold stress to regulate biological activities, which ultimately led to changes in the plant phenotype and physiological indicators.
TFs are proteins that affect many biological processes, including growth, development, cell division and responses to environment stimuli such as cold stress [75]. Transcriptome data showed that 677 transcription factors were differentially expressed under cold stress in E. urophylla. Of these, 14 TFs were identified as the target genes of DE_lncRNAs. Further investigation of the DE_TFs showed that 77 DEGs belonged to the ERF transcription factor family, which was the most enriched TF family, followed by the MYB and NAC TF families. APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERFs) are key regulators of several abiotic stresses [76,77], and CBF/dehydration-responsive element-binding factor (DREB1) from the AP2/ERF family (CBF1/DREB1B, CBF2/DREB1C and CBF3/DREB1A) in Arabidopsis plays a prominent role in cold acclimation [78,79]. Analyzing changes in the expression of EuCBF1/EuCBF2/EuCBF3 after low-temperature treatment showed that, compared with the CK group, the expression of EuCBF1-3 in E. urophylla after low-temperature treatment significantly increased in the treatment groups, while its expression in the recovery group reduced compared with the treatment groups (Table S6). This indicates that the CBF transcription factor in E. urophylla could respond to cold stress. Low temperatures induce its expression, while normal conditions inhibit its expression.
In addition to the ERF family of TFs, MYB transcription factors also responded significantly to cold stress. We found 66 DEGs belonging to the MYB transcription factor family. MYB proteins are one of the largest TF families and are considered unique to plants. They are classified into four major subfamilies based on their structural features: 1R-MYB proteins, R2R3-MYB proteins, 3R-MYB proteins and 4R-MYB proteins [80]. They can promote or reduce the cold tolerance of plants by regulating the expression of downstream genes [81,82]. The NAC transcription factor family also showed significant differences in expression under cold stress in E. urophylla. All proteins in this family have a conserved N-terminal DNA-binding domain called the NAC structure domain [83]. Although most of the NAC genes are involved in regulating the salt and drought stress tolerance in plants, some are also involved in regulating temperature tolerance [46]. For example, overexpression of GmNAC20 and GmNAC11 in transgenic Arabidopsis enhances its salt and cold tolerance, respectively, due to the interactive effects of these genes with the DREB/CBF pathway and other stress-related genes [83]. Transcription factors themselves can be regulated at the transcription level by upstream components and also be subjected to various tiers of modifications at the post-transcription level, such as ubiquitination and sumoylation, thus forming a complex regulatory network to modulate the expression of stress-responsive genes [84,85]. Therefore, analyzing the roles of these TFs in various signal transduction pathways could help us better understand the mechanism of response to cold stress in E. urophylla.
Studies have shown that miRNA recognizes target mRNA through base pairing and participates in the degradation or inhibition of the silent complex that regulates the translation of target mRNA [86]. The competitive endogenous RNA (ceRNA) hypothesis proposed by Salmena et al. [87] is considered a new regulatory mechanism between non-coding RNA and mRNA. Based on the previous research on genome-wide miRNA identification, we analyzed the targeting relationship among miRNAs, DE_lncRNAs and mRNAs under cold stress to construct an mRNA–miRNA–lncRNA regulatory network in E. urophylla. In addition to the DE_lncRNAs identified in this study, miRNAs in E. urophylla also responded to cold stress. In the regulatory network constructed, miRNAs and lncRNAs constitute a targeting relationship, and 22 DE_lncRNAs are predicted to be targets or targeting mimics of 20 miRNAs. In this study, miRNA also showed differential expression under cold stress, which affected the expression of mRNA and lncRNA, indicating that cold stress can change the expression of miRNA and further regulate mRNA and lncRNA to respond to cold temperatures in E. urophylla (Figure 9). Cold stress can cause significant changes in the expression of these miRNAs. For example, miR167 in tomatoes [88] and miR193, miR169 and miR168 in Arabidopsis thaliana [89,90,91] are upregulated under low-temperature treatment, and the expression of miR166 in cabbage [92], miR482 in Hevea brasiliensis and miR858 in Citrullus lanatus change under cold stress [93,94].
miRNAs do not directly play a role in plant growth and development or in the response to abiotic stress but do participate in a plant’s response to abiotic stress by regulating key components of a complex regulatory network. Both lncRNAs and mRNAs can act as target genes of miRNAs. miR167 is a highly conserved family acting in a wide range of biological processes in plants [95]. ARF6 and ARF8 are the targets of miR167 that promote adventitious rooting. Cold stress induces the significant differential expression of miR167 in many plant species [96]. In this study, we found new regulatory networks that involved miR167 in E. urophylla under cold stress. For example, miRNAs EUC_00002677 and EUC_00008890, both belonging to the miR167 family, responded to cold stress by interacting with the target genes in E. urophylla. EUC_00002677 targeted novel lncRNAs MSTRG.7690 and EuGBF3, and EUC_00008890 targeted novel lncRNAs MSTRG.46425 and EuDAG2. EuGBF3 and lncRNA MSTRG.7690 were both involved in the response to cold stress regulated by miRNA EUC_00002677, inducing the gene expression difference in E. urophylla (Figure 9). DOF AFFECTING GERMINATION 2 (DAG2), which encodes Dof zinc finger transcription factors, is a positive regulator of light-mediated seed germination and is repressed by DAG1 [97]. This interaction network suggested that EuDAG2 and novel lncRNA MSTRG.46425 were both responsive to a low-temperature environment, which may be due to the regulation of miRNA EUC_00008890. In general, the mRNA–miRNA–lncRNA regulatory network we constructed revealed a new molecular mechanism of the response to cold stress in E. urophylla and may play a significant role in understanding the genetic improvement of cold tolerance in this tree species.

5. Conclusions

This study identified the lncRNAs’ response to cold stress in E. urophylla and found that lncRNAs exhibit different functions. We also constructed mRNA–miRNA–lncRNA regulatory networks in E. urophylla under cold stress, which may provide new insights into enhancing the cold tolerance in this tree species.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/f12070836/s1. Figure S1: Number of DE_lncRNAs and DEGs in the control group, four cold treatment groups and one recovery group; Table S1: List of primers for genes used for qPCR; Table S2: List of DEGs in Eucalyptus urophylla; Table S3: List of DE_lncRNAs and target genes in Eucalyptus urophylla; Table S4: Annotation of DE_TFs; Table S5: List of putative targets and target mimics of DE_lncRNAs for miRNA; Table S6: EuCBF gene expression under cold stress.

Author Contributions

Conceptualization, H.C., J.L. and J.Y.; methodology, H.C., J.L. and X.K.; formal analysis, H.C., J.L. and J.Y.; investigation, J.L., B.Q., Y.Z. and Z.L.; writing—original draft preparation, H.C. and J.L.; writing—review and editing, J.Y. and X.K.; funding acquisition, J.Y. and X.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (31901337).

Data Availability Statement

The data sets generated or analyzed during this study are included in this published article and its Supplementary Materials.

Acknowledgments

The authors are grateful to Jianzhong Wang from Guangxi Dongmen Forest Farm for collecting materials and additional help.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

lncRNALong non-coding RNA
DEGsDifferentially expressed genes
DE_lncRNADifferentially expressed long non-coding RNA
DE_TFDifferentially expressed transcription factor
GOGene ontology
FPKMFragments per kilo-base per million reads
qRT-PCRQuantitative real-time polymerase chain reaction
ROSReactive oxygen species
NGSNext-generation DNA sequencing
ORFOpen reading frame
Fv/FmThe maximal quantum efficiency of photosystem II
RELThe relative electrolyte leakage
ceRNACompetitive endogenous RNA

References

  1. Williams, J.; Woinarski, J. Eucalypt Ecology: Individuals to Ecosystems; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  2. Pérez, S.; Renedo, C.; Ortiz, A.; Mañana, M.; Silió, D. Energy evaluation of the Eucalyptus globulus and the Eucalyptus nitens in the north of Spain (Cantabria). Thermochim. Acta 2006, 451, 57–64. [Google Scholar] [CrossRef]
  3. Cherubini, F.; Strømman, A.H. Life cycle assessment of bioenergy systems: State of the art and future challenges. Bioresour. Technol. 2011, 102, 437–451. [Google Scholar] [CrossRef] [PubMed]
  4. Rockwood, D.L.; Rudie, A.W.; Ralph, S.A.; Zhu, J.Y.; Winandy, J.E. Energy product options for Eucalyptus species grown as short rotation woody crops. Int. J. Mol. Sci. 2008, 9, 1361–1378. [Google Scholar] [CrossRef] [PubMed]
  5. Carroll, A.; Somerville, C. Cellulosic biofuels. Annu. Rev. Plant Biol. 2009, 60, 165–182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hung, T.D.; Brawner, J.T.; Meder, R.; Lee, D.J.; Southerton, S.; Thinh, H.H.; Dieters, M.J. Estimates of genetic parameters for growth and wood properties in Eucalyptus pellita F. Muell. to support tree breeding in Vietnam. Ann. For. Sci. 2014, 72, 205–217. [Google Scholar] [CrossRef]
  7. Carrillo, I.; Vidal, C.; Elissetche, J.P.; Mendonça, R.T. Wood anatomical and chemical properties related to the pulpability of Eucalyptus globulus: A review. South. For. J. For. Sci. 2018, 80, 1–8. [Google Scholar] [CrossRef]
  8. Dhakad, A.K.; Pandey, V.V.; Beg, S.; Rawat, J.M.; Singh, A. Biological, medicinal and toxicological significance of Eucalyptus leaf essential oil: A review. J. Sci. Food Agric. 2018, 98, 833–848. [Google Scholar] [CrossRef] [PubMed]
  9. Forrester, D.I. Growth responses to thinning, pruning and fertiliser application in Eucalyptus plantations: A review of their production ecology and interactions. For. Ecol. Manag. 2013, 310, 336–347. [Google Scholar] [CrossRef]
  10. Bayle, G. Ecological and social impacts of Eucalyptus tree plantation on the environment. J. Biodivers. Conserv. Bioresour. Manag. 2019, 5, 93–104. [Google Scholar] [CrossRef]
  11. Sands, P.J.; Landsberg, J.J. Parameterisation of 3-PG for plantation grown Eucalyptus globulus. For. Ecol. Manag. 2002, 163, 273–292. [Google Scholar] [CrossRef]
  12. Chang, Y.N.; Zhu, C.; Jiang, J.; Zhang, H.; Zhu, J.K.; Duan, C.G. Epigenetic regulation in plant abiotic stress responses. J. Integr. Plant Biol. 2020, 62, 563–580. [Google Scholar] [CrossRef]
  13. Keller, G.; Cao, P.B.; Clemente, H.S.; El Kayal, W.; Marque, C.; Teulieres, C. Transcript profiling combined with functional annotation of 2,662 ESTs provides a molecular picture of Eucalyptus gunnii cold acclimation. Trees 2013, 27, 1713–1735. [Google Scholar] [CrossRef]
  14. Thomashow, M.F. Plant cold acclimation: Freezing tolerance genes and regulatory mechanisms. Annu. Rev. Plant Biol. 1999, 50, 571–599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tarkowski, Ł.P.; Ende, W.V.D. Cold tolerance triggered by soluble sugars: A multifaceted countermeasure. Front. Plant Sci. 2015, 6, 203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Krasensky, J.; Jonak, C. Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J. Exp. Bot. 2012, 63, 1593–1608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Liu, X.-G.; Xu, H.; Zhang, J.-Y.; Liang, G.-W.; Liu, Y.-T.; Guo, A.-G. Effect of low temperature on chlorophyll biosynthesis in albinism line of wheat (Triticum aestivum) FA85. Physiol. Plant. 2012, 145, 384–394. [Google Scholar] [CrossRef]
  18. Zhou, C.D.; Yang, C.D.; Zhan, L. Effects of low temperature stress on physiological and biochemical characteristics of Podocarpus nagi. Agric. Sci. Technol. 2012, 13, 533–536. [Google Scholar] [CrossRef]
  19. Liu, Y.; Jiang, Y.; Lan, J.; Zou, Y.; Gao, J. Comparative transcriptomic analysis of the response to cold acclimation in Eucalyptus dunnii. PLoS ONE 2014, 9, e113091. [Google Scholar] [CrossRef] [Green Version]
  20. Xie, Y.; Chen, P.; Yan, Y.; Bao, C.; Li, X.; Wang, L.; Shen, X.; Li, H.; Liu, X.; Niu, C.; et al. An atypical R2R3 MYB transcription factor increases cold hardiness by CBF-dependent and CBF-independent pathways in apple. New Phytol. 2018, 218, 201–218. [Google Scholar] [CrossRef]
  21. Lee, J.G.; Yi, G.; Seo, J.; Kang, B.-C.; Choi, J.H.; Lee, E.J. Jasmonic acid and ERF family genes are involved in chilling sensitivity and seed browning of pepper fruit after harvest. Sci. Rep. 2020, 10, 17949. [Google Scholar] [CrossRef]
  22. Dong, C.-H.; Pei, H. Over-expression of miR397 improves plant tolerance to cold stress in Arabidopsis thaliana. J. Plant Biol. 2014, 57, 209–217. [Google Scholar] [CrossRef]
  23. Kindgren, P.; Ard, R.; Ivanov, M.; Marquardt, S. Transcriptional read-through of the long non-coding RNA SVALKA governs plant cold acclimation. Nat. Commun. 2018, 9, 4561. [Google Scholar] [CrossRef] [Green Version]
  24. Kopp, F.; Mendell, J.T. Functional classification and experimental dissection of long noncoding RNAs. Cell 2018, 172, 393–407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Gil, N.; Ulitsky, I. Regulation of gene expression by cis-acting long non-coding RNAs. Nat. Rev. Genet. 2020, 21, 102–117. [Google Scholar] [CrossRef] [PubMed]
  26. Muthusamy, M.; Uma, S.; Backiyarani, S.; Saraswathi, M.S. Genome-wide screening for novel, drought stress-responsive long non-coding RNAs in drought-stressed leaf transcriptome of drought-tolerant and -susceptible banana (Musa spp.) cultivars using Illumina high-throughput sequencing. Plant Biotechnol. Rep. 2015, 9, 279–286. [Google Scholar] [CrossRef]
  27. Ponting, C.P.; Oliver, P.L.; Reik, W. Evolution and functions of long noncoding RNAs. Cell 2009, 136, 629–641. [Google Scholar] [CrossRef] [Green Version]
  28. Yao, R.; Wang, Y.; Chen, L.-L. Cellular functions of long noncoding RNAs. Nat. Cell Biol. 2019, 21, 542–551. [Google Scholar] [CrossRef]
  29. Liu, J.; Jung, C.; Xu, J.; Wang, H.; Deng, S.; Bernad, L.; Arenas-Huertero, C.; Chua, N.-H. Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis. Plant Cell 2012, 24, 4333–4345. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, Y.-C.; Liao, J.-Y.; Li, Z.-Y.; Yu, Y.; Zhang, J.-P.; Li, Q.-F.; Qu, L.-H.; Shu, W.-S.; Chen, Y.-Q. Genome-wide screening and functional analysis identify a large number of long noncoding RNAs involved in the sexual reproduction of rice. Genome Biol. 2014, 15, 512. [Google Scholar] [CrossRef] [Green Version]
  31. Li, L.; Eichten, S.R.; Shimizu, R.; Petsch, K.; Yeh, C.-T.; Wu, W.; Chettoor, A.M.; A Givan, S.; A Cole, R.; E Fowler, J.; et al. Genome-wide discovery and characterization of maize long non-coding RNAs. Genome Biol. 2014, 15, R40. [Google Scholar] [CrossRef] [Green Version]
  32. Shuai, P.; Liang, D.; Tang, S.; Zhang, Z.; Ye, C.-Y.; Su, Y.; Xia, X.; Yin, W. Genome-wide identification and functional prediction of novel and drought-responsive lincRNAs in Populus trichocarpa. J. Exp. Bot. 2014, 65, 4975–4983. [Google Scholar] [CrossRef] [PubMed]
  33. Lin, Z.; Long, J.; Yin, Q.; Wang, B.; Li, H.; Luo, J.; Wang, H.C.; Wu, A.-M. Identification of novel lncRNAs in Eucalyptus grandis. Ind. Crop. Prod. 2019, 129, 309–317. [Google Scholar] [CrossRef]
  34. Ben Amor, B.; Wirth, S.; Merchan, F.; Laporte, P.; D’Aubenton-Carafa, Y.; Hirsch, J.; Maizel, A.; Mallory, A.; Lucas, A.; Deragon, J.M.; et al. Novel long non-protein coding RNAs involved in Arabidopsis differentiation and stress responses. Genome Res. 2008, 19, 57–69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhao, M.; Wang, T.; Sun, T.; Yu, X.; Tian, R.; Zhang, W.-H. Identification of tissue-specific and cold-responsive lncRNAs in Medicago truncatula by high-throughput RNA sequencing. BMC Plant Biol. 2020, 20, 1–16. [Google Scholar] [CrossRef] [Green Version]
  36. Oberschelp, G.P.J.; Guarnaschelli, A.B.; Teson, N.; Harrand, L.; Podestá, F.E.; Margarit, E. Cold acclimation and freezing tolerance in three Eucalyptus species: A metabolomic and proteomic approach. Plant Physiol. Biochem. 2020, 154, 316–327. [Google Scholar] [CrossRef]
  37. Gaete-Loyola, J.; Lagos, C.; Beltrán, M.F.; Valenzuela, S.; Emhart, V.; Fernández, M. Transcriptome profiling of Eucalyptus nitens reveals deeper insight into the molecular mechanism of cold acclimation and deacclimation process. Tree Genet. Genomes 2017, 13, 37. [Google Scholar] [CrossRef]
  38. Pryor, L.D.; Williams, E.R.; Gunn, B.V. A morphometric analysis of Eucalyptus urophylla and related taxa with descriptions of two new species. Aust. Syst. Bot. 1995, 8, 57–70. [Google Scholar] [CrossRef]
  39. He, F.; Wang, H.-L.; Li, H.-G.; Su, Y.; Li, S.; Yang, Y.; Feng, C.-H.; Yin, W.; Xia, X. PeCHYR1, a ubiquitin E3 ligase from Populus euphratica, enhances drought tolerance via ABA-induced stomatal closure by ROS production in Populus. Plant Biotechnol. J. 2018, 16, 1514–1528. [Google Scholar] [CrossRef] [Green Version]
  40. Ban, Q.; Wang, X.; Pan, C.; Wang, Y.; Kong, L.; Jiang, H.; Xu, Y.; Wang, W.; Pan, Y.; Li, Y.; et al. Comparative analysis of the response and gene regulation in cold resistant and susceptible tea plants. PLoS ONE 2017, 12, e0188514. [Google Scholar] [CrossRef] [Green Version]
  41. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  42. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [Green Version]
  43. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [Green Version]
  44. Pappas, M.D.C.R.; Pappasjr, G.J.; Grattapaglia, D. Genome-wide discovery and validation of Eucalyptus small RNAs reveals variable patterns of conservation and diversity across species of Myrtaceae. BMC Genom. 2015, 16, 1113. [Google Scholar] [CrossRef] [Green Version]
  45. Dai, X.; Zhuang, Z.; Zhao, P.X. psRNATarget: A plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018, 46, W49–W54. [Google Scholar] [CrossRef] [Green Version]
  46. Erpen, L.; Devi, H.S.; Grosser, J.W.; Dutt, M. Potential use of the DREB/ERF, MYB, NAC and WRKY transcription factors to improve abiotic and biotic stress in transgenic plants. Plant Cell Tissue Organ Cult. 2018, 132, 1–25. [Google Scholar] [CrossRef]
  47. Büttner, M.; Singh, K.B. Arabidopsis thaliana ethylene-responsive element binding protein (AtEBP), an ethylene-inducible, GCC box DNA-binding protein interacts with an ocs element binding protein. Proc. Natl. Acad. Sci. USA 1997, 94, 5961–5966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Zarei, A.; Körbes, A.P.; Younessi, P.; Montiel, G.; Champion, A.; Memelink, J. Two GCC boxes and AP2/ERF-domain transcription factor ORA59 in jasmonate/ethylene-mediated activation of the PDF1.2 promoter in Arabidopsis. Plant Mol. Biol. 2011, 75, 321–331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Chen, L.-L.; Carmichael, G.G. Decoding the function of nuclear long non-coding RNAs. Curr. Opin. Cell Biol. 2010, 22, 357–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Janská, A.; Maršík, P.; Zelenková, S.; Ovesná, J. Cold stress and acclimation—What is important for metabolic adjustment? Plant Biol. 2010, 12, 395–405. [Google Scholar] [CrossRef]
  51. Guo, X.; Liu, D.; Chong, K. Cold signaling in plants: Insights into mechanisms and regulation. J. Integr. Plant Biol. 2018, 60, 745–756. [Google Scholar] [CrossRef] [Green Version]
  52. Wong, C.E.; Li, Y.; Labbe, A.; Guevara, D.; Nuin, P.; Whitty, B.; Díaz, C.; Golding, G.B.; Gray, G.R.; Weretilnyk, E.A.; et al. Transcriptional profiling implicates novel interactions between abiotic stress and hormonal responses in Thellungiella, a close relative of Arabidopsis. Plant Physiol. 2006, 140, 1437–1450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. de Araújo, N.O.; Santos, M.N.D.S.; de Araujo, F.F.; Véras, M.L.M.; Tello, J.P.D.J.; Arruda, R.D.S.; Fugate, K.K.; Finger, F.L. Balance between oxidative stress and the antioxidant system is associated with the level of cold tolerance in sweet potato roots. Postharvest Biol. Technol. 2021, 172, 111359. [Google Scholar] [CrossRef]
  54. Lee, H.J.; Lee, J.H.; Wi, S.; Jang, Y.; An, S.; Choi, C.K.; Jang, S. Exogenously applied glutamic acid confers improved yield through increased photosynthesis efficiency and antioxidant defense system under chilling stress condition in Solanum lycopersicum L. cv. Dotaerang Dia. Sci. Hortic. 2021, 277, 109817. [Google Scholar] [CrossRef]
  55. Zhou, J.; Wang, J.; Shi, K.; Xia, X.J.; Zhou, Y.H.; Yu, J.Q. Hydrogen peroxide is involved in the cold acclimation-induced chilling tolerance of tomato plants. Plant Physiol. Biochem. 2012, 60, 141–149. [Google Scholar] [CrossRef]
  56. Hetherington, S.E.; Oquist, G. Monitoring chilling injury: A comparison of chlorophyll fluorescence measurements, post-chilling growth and visible symptoms of injury in Zea Mays. Physiol. Plant. 1988, 72, 241–247. [Google Scholar] [CrossRef]
  57. Farooq, M.; Aziz, T.; Wahid, A.; Lee, D.-J.; Siddique, K.H.M. Chilling tolerance in maize: Agronomic and physiological approaches. Crop. Pasture Sci. 2009, 60, 501–516. [Google Scholar] [CrossRef]
  58. Baker, N.R.; East, T.M.; Long, S.P. Chilling damage to photosynthesis in young Zea Mays. J. Exp. Bot. 1983, 34, 189–197. [Google Scholar] [CrossRef]
  59. Apel, K.; Hirt, H. Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Annu. Rev. Plant Biol. 2004, 55, 373–399. [Google Scholar] [CrossRef] [Green Version]
  60. Allen, D.J.; Ort, D.R. Impacts of chilling temperatures on photosynthesis in warm-climate plants. Trends Plant Sci. 2001, 6, 36–42. [Google Scholar] [CrossRef]
  61. Suzuki, N.; Mittler, R. Reactive oxygen species and temperature stresses: A delicate balance between signaling and destruction. Physiol. Plant. 2006, 126, 45–51. [Google Scholar] [CrossRef]
  62. Li, H.; Jiang, X.; Lv, X.; Ahammed, G.J.; Guo, Z.; Qi, Z.; Yu, J.; Zhou, Y.H. Tomato GLR3.3 and GLR3.5 mediate cold acclimation-induced chilling tolerance by regulating apoplastic H2O2 production and redox homeostasis. Plant Cell Environ. 2019, 42, 3326–3339. [Google Scholar] [CrossRef]
  63. Zhang, C.; Tian, S. Crucial contribution of membrane lipids’ unsaturation to acquisition of chilling-tolerance in peach fruit stored at 0 °C. Food Chem. 2009, 115, 405–411. [Google Scholar] [CrossRef]
  64. Kudla, J.; Becker, D.; Grill, E.; Hedrich, R.; Hippler, M.; Kummer, U.; Parniske, M.; Romeis, T.; Schumacher, K. Advances and current challenges in calcium signaling. New Phytol. 2018, 218, 414–431. [Google Scholar] [CrossRef]
  65. Peleg, Z.; Blumwald, E. Hormone balance and abiotic stress tolerance in crop plants. Curr. Opin. Plant Biol. 2011, 14, 290–295. [Google Scholar] [CrossRef] [PubMed]
  66. Yuan, H.-M.; Sheng, Y.; Chen, W.-J.; Lu, Y.-Q.; Tang, X.; Ou-Yang, M.; Huang, X. Overexpression of Hevea brasiliensis HbICE1 enhances cold tolerance in Arabidopsis. Front. Plant Sci. 2017, 8, 1462. [Google Scholar] [CrossRef] [Green Version]
  67. Hu, Z.R.; Liu, A.; Bi, A.Y.; Amombo, E.; Gitau, M.M.; Huang, X.B.; Chen, L.; Fu, J.M. Identification of differentially expressed proteins in bermudagrass response to cold stress in the presence of ethylene. Environ. Exp. Bot. 2017, 139, 67–78. [Google Scholar] [CrossRef]
  68. Huang, G.-T.; Ma, S.-L.; Bai, L.-P.; Zhang, L.; Ma, H.; Jia, P.; Liu, J.; Zhong, M.; Guo, Z.-F. Signal transduction during cold, salt, and drought stresses in plants. Mol. Biol. Rep. 2012, 39, 969–987. [Google Scholar] [CrossRef] [PubMed]
  69. Crifò, T.; Puglisi, I.; Petrone, G.; Recupero, G.R.; Piero, A.R.L. Expression analysis in response to low temperature stress in blood oranges: Implication of the flavonoid biosynthetic pathway. Gene 2011, 476, 1–9. [Google Scholar] [CrossRef]
  70. Campbell, K.; Vowinckel, J.; Keller, M.A.; Ralser, M. Methionine metabolism alters oxidative stress resistance via the pentose phosphate pathway. Antioxid. Redox Signal. 2016, 24, 543–547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Moghaddam, M.R.B.; Le Roy, K.; Xiang, L.; Rolland, F.; Ende, W.V.D. Sugar signalling and antioxidant network connections in plant cells. FEBS J. 2010, 277, 2022–2037. [Google Scholar] [CrossRef]
  72. Couée, I.; Sulmon, C.; Gouesbet, G.; El Amrani, A. Involvement of soluble sugars in reactive oxygen species balance and responses to oxidative stress in plants. J. Exp. Bot. 2006, 57, 449–459. [Google Scholar] [CrossRef]
  73. Bohnert, H.J.; Sheveleva, E. Plant stress adaptations—Making metabolism move. Curr. Opin. Plant Biol. 1998, 1, 267–274. [Google Scholar] [CrossRef]
  74. Welling, A.; Palva, E.T. Molecular control of cold acclimation in trees. Physiol. Plant. 2006, 127, 167–181. [Google Scholar] [CrossRef]
  75. Carretero-Paulet, L.; Galstyan, A.; Villanova, I.R.; Martinez-Garcia, J.; Bilbao-Castro, J.R.; Robertson, D.L. Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis, poplar, rice, moss, and algae. Plant Physiol. 2010, 153, 1398–1412. [Google Scholar] [CrossRef] [Green Version]
  76. Mizoi, J.; Shinozaki, K.; Yamaguchi-Shinozaki, K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 2012, 1819, 86–96. [Google Scholar] [CrossRef] [PubMed]
  77. Dietz, K.-J.; Vogel, M.O.; Viehhauser, A. AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma 2010, 245, 3–14. [Google Scholar] [CrossRef] [PubMed]
  78. Novillo, F.; Medina, J.; Salinas, J. Arabidopsis CBF1 and CBF3 have a different function than CBF2 in cold acclimation and define different gene classes in the CBF regulon. Proc. Natl. Acad. Sci. USA 2007, 104, 21002–21007. [Google Scholar] [CrossRef] [Green Version]
  79. Cook, D.; Fowler, S.; Fiehn, O.; Thomashow, M.F. A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis. Proc. Natl. Acad. Sci. USA 2004, 101, 15243–15248. [Google Scholar] [CrossRef] [Green Version]
  80. Xing, C.; Liu, Y.; Zhao, L.; Zhang, S.; Huang, X. A novel MYB transcription factor regulates ascorbic acid synthesis and affects cold tolerance. Plant Cell Environ. 2018, 42, 832–845. [Google Scholar] [CrossRef]
  81. Kim, S.H.; Kim, H.S.; Bahk, S.; An, J.; Yoo, Y.; Kim, J.-Y.; Chung, W.S. Phosphorylation of the transcriptional repressor MYB15 by mitogen-activated protein kinase 6 is required for freezing tolerance in Arabidopsis. Nucleic Acids Res. 2017, 45, 6613–6627. [Google Scholar] [CrossRef] [Green Version]
  82. Aida, M.; Ishida, T.; Fukaki, H.; Fujisawa, H.; Tasaka, M. Genes involved in organ separation in Arabidopsis: An analysis of the cup-shaped cotyledon mutant. Plant Cell 1997, 9, 841–857. [Google Scholar] [CrossRef] [Green Version]
  83. Hao, Y.-J.; Wei, W.; Song, Q.-X.; Chen, H.-W.; Zhang, Y.-Q.; Wang, F.; Zou, H.-F.; Lei, G.; Tian, A.; Zhang, W.-K.; et al. Soybean NAC transcription factors promote abiotic stress tolerance and lateral root formation in transgenic plants. Plant J. 2011, 68, 302–313. [Google Scholar] [CrossRef]
  84. Samad, A.F.A.; Sajad, M.; Nazaruddin, N.; Fauzi, I.A.; Murad, A.M.A.; Zainal, Z.; Ismail, I. MicroRNA and transcription factor: Key players in plant regulatory network. Front. Plant Sci. 2017, 8, 565. [Google Scholar] [CrossRef] [Green Version]
  85. Moore, J.W.; Loake, G.; Spoel, S.H. Transcription dynamics in plant immunity. Plant Cell 2011, 23, 2809–2820. [Google Scholar] [CrossRef] [Green Version]
  86. Yang, J.; Gong, Y.; Cai, J.; Liu, Q.; Zhang, Z. lnc-3215 suppression leads to calcium overload in selenium deficiency-induced chicken heart lesion via the lnc-3215-miR-1594-TNN2 pathway. Mol. Ther. Nucleic Acids 2019, 18, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Jodder, J.; Das, R.; Sarkar, D.; Bhattacharjee, P.; Kundu, P. Distinct transcriptional and processing regulations control miR167a level in tomato during stress. RNA Biol. 2017, 15, 130–143. [Google Scholar] [CrossRef] [Green Version]
  89. Sunkar, R.; Zhu, J.-K. Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell 2004, 16, 2001–2019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Singh, A.; Gandhi, N.; Mishra, V.; Yadav, S.; Rai, V.; Sarkar, A.K. Role of abiotic stress responsive miRNAs in Arabidopsis root development. J. Plant Biochem. Biotechnol. 2020, 29, 733–742. [Google Scholar] [CrossRef]
  91. Liu, H.-H.; Tian, X.; Li, Y.-J.; Wu, C.-A.; Zheng, C.-C. Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA 2008, 14, 836–843. [Google Scholar] [CrossRef] [Green Version]
  92. Zeng, X.; Xu, Y.; Jiang, J.; Zhang, F.; Ma, L.; Wu, D.; Wang, Y.; Sun, W. Identification of cold stress responsive microRNAs in two winter turnip rape (Brassica rapa L.) by high throughput sequencing. BMC Plant Biol. 2018, 18, 52. [Google Scholar] [CrossRef]
  93. Kuruvilla, L.; Sathik, M.M.; Thomas, M.; Luke, L.; Kv, S. Identification and validation of cold responsive microRNAs of Hevea brasiliensis using high throughput sequencing. J. Crop. Sci. Biotechnol. 2017, 20, 369–377. [Google Scholar] [CrossRef]
  94. Li, H.; Dong, Y.; Chang, J.; Hejie, C.; Chen, H.; Liu, Q.; Wei, C.; Ma, J.; Zhang, Y.; Yang, J.; et al. High-throughput microRNA and mRNA sequencing reveals that microRNAs may be involved in melatonin-mediated cold tolerance in Citrullus lanatus L. Front. Plant Sci. 2016, 7, 1231. [Google Scholar] [CrossRef] [Green Version]
  95. Zhao, Z.; Feng, Q.; Cao, X.; Zhu, Y.; Wang, H.; Chandran, V.; Fan, J.; Zhao, J.; Pu, M.; Li, Y.; et al. Osa-miR167d facilitates infection of Magnaporthe oryzae in rice. J. Integr. Plant Biol. 2019, 62, 702–715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Zhang, B. MicroRNA: A new target for improving plant tolerance to abiotic stress. J. Exp. Bot. 2015, 66, 1749–1761. [Google Scholar] [CrossRef] [PubMed]
  97. Santopolo, S.; Boccaccini, A.; Lorrai, R.; Ruta, V.; Capauto, D.; Minutello, E.; Serino, G.; Costantino, P.; Vittorioso, P. DOF AFFECTING GERMINATION 2 is a positive regulator of light-mediated seed germination and is repressed by DOF AFFECTING GERMINATION 1. BMC Plant Biol. 2015, 15, 72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Forests 12 00836 g001 550
Figure 1. Cold responses of two-month-old Eucalyptus urophylla S.T.Blake subjected to cold stress (4 °C). (a) Plant morphology in CK, LT6, LT12, LT24, LT48 and RE24. The white box is a close-up of the terminal buds and young leaves. (b) Fv/Fm value in each group. (c) REL value in each group. Values shown in (b,c) are means ± SD with three biological replicates. The letters indicate significant differences (one way ANOVA followed by post hoc Waller-Duncan test, p < 0.05).
Figure 1. Cold responses of two-month-old Eucalyptus urophylla S.T.Blake subjected to cold stress (4 °C). (a) Plant morphology in CK, LT6, LT12, LT24, LT48 and RE24. The white box is a close-up of the terminal buds and young leaves. (b) Fv/Fm value in each group. (c) REL value in each group. Values shown in (b,c) are means ± SD with three biological replicates. The letters indicate significant differences (one way ANOVA followed by post hoc Waller-Duncan test, p < 0.05).
Forests 12 00836 g001
Forests 12 00836 g002 550
Figure 2. Informatics pipeline for the identification of lncRNAs in Eucalyptus urophylla under cold stress. i: intronic lncRNA; j: bidirectional lncRNA; o: sense lncRNA; u: intergenic lncRNA; x: antisense lncRNA.
Figure 2. Informatics pipeline for the identification of lncRNAs in Eucalyptus urophylla under cold stress. i: intronic lncRNA; j: bidirectional lncRNA; o: sense lncRNA; u: intergenic lncRNA; x: antisense lncRNA.
Forests 12 00836 g002
Forests 12 00836 g003 550
Figure 3. Characteristics of mRNAs and lncRNAs. (a) Compositions of different kinds of lncRNAs. (b) Number of exons per transcription for mRNAs and lncRNAs. (c) Open reading frame (ORF) lengths of all mRNAs and lncRNAs. (d) Isoform number per gene of lncRNAs and mRNAs. (e) Transcript size distributions for all mRNAs and lncRNAs. (f) Comparison of mRNA and lncRNA expression levels. The crosses represent the maximum and minimum values, circles represent the medium values and the horizontal lines represent the first quartile and third quartile values.
Figure 3. Characteristics of mRNAs and lncRNAs. (a) Compositions of different kinds of lncRNAs. (b) Number of exons per transcription for mRNAs and lncRNAs. (c) Open reading frame (ORF) lengths of all mRNAs and lncRNAs. (d) Isoform number per gene of lncRNAs and mRNAs. (e) Transcript size distributions for all mRNAs and lncRNAs. (f) Comparison of mRNA and lncRNA expression levels. The crosses represent the maximum and minimum values, circles represent the medium values and the horizontal lines represent the first quartile and third quartile values.
Forests 12 00836 g003
Forests 12 00836 g004 550
Figure 4. Confirmation of the expression patterns of three lncRNAs and three mRNAs of Eucalyptus urophylla under cold stress using qRT-PCR. The first layer is three randomly selected mRNAs, and the second layer is three lncRNAs. The values shown are the means ± SD of three replicates. EuGADPH was used as the reference gene.
Figure 4. Confirmation of the expression patterns of three lncRNAs and three mRNAs of Eucalyptus urophylla under cold stress using qRT-PCR. The first layer is three randomly selected mRNAs, and the second layer is three lncRNAs. The values shown are the means ± SD of three replicates. EuGADPH was used as the reference gene.
Forests 12 00836 g004
Forests 12 00836 g005 550
Figure 5. Summary of DEGs and DE_lncRNAs in Eucalyptus urophylla under cold stress. (a) Number of DEGs upregulated or downregulated by cold stress and recovery treatment. (b) Number of DE_lncRNAs upregulated or downregulated by cold stress and recovery treatment. (c) Heat map of DE_lncRNAs under cold stress. (d) Heat map of DEGs under cold stress. All expression levels were normalized by FPKMs.
Figure 5. Summary of DEGs and DE_lncRNAs in Eucalyptus urophylla under cold stress. (a) Number of DEGs upregulated or downregulated by cold stress and recovery treatment. (b) Number of DE_lncRNAs upregulated or downregulated by cold stress and recovery treatment. (c) Heat map of DE_lncRNAs under cold stress. (d) Heat map of DEGs under cold stress. All expression levels were normalized by FPKMs.
Forests 12 00836 g005
Forests 12 00836 g006 550
Figure 6. Distribution of transcription factors responsive to cold stress in Eucalyptus urophylla. Data are sorted by the number of DEGs, and the top 18 families are shown.
Figure 6. Distribution of transcription factors responsive to cold stress in Eucalyptus urophylla. Data are sorted by the number of DEGs, and the top 18 families are shown.
Forests 12 00836 g006
Forests 12 00836 g007 550
Figure 7. GO enrichment analysis of mRNAs and lncRNAs, showing differentially expressed genes in Eucalyptus urophylla under cold stress. (a) GO enrichment of DEGs. (b) GO enrichment of the target genes of DE_lncRNAs.
Figure 7. GO enrichment analysis of mRNAs and lncRNAs, showing differentially expressed genes in Eucalyptus urophylla under cold stress. (a) GO enrichment of DEGs. (b) GO enrichment of the target genes of DE_lncRNAs.
Forests 12 00836 g007
Forests 12 00836 g008 550
Figure 8. Representatives of the predicted regulatory networks among DEGs, DE_lncRNAs, DE_TFs and miRNAs under cold stress in Eucalyptus urophylla. The square, circular, octagonal and triangular nodes represent the DEGs, DE_lncRNAs, DE_TFs, and miRNAs, respectively.
Figure 8. Representatives of the predicted regulatory networks among DEGs, DE_lncRNAs, DE_TFs and miRNAs under cold stress in Eucalyptus urophylla. The square, circular, octagonal and triangular nodes represent the DEGs, DE_lncRNAs, DE_TFs, and miRNAs, respectively.
Forests 12 00836 g008
Forests 12 00836 g009 550
Figure 9. Relative expression of mRNA, miRNA and lncRNA in Eucalyptus urophylla under different cold treatments. Purple bar: Eucgr.E00078; blue bar: EUC_00002677; green bar: MSTRG.7690. EuGADPH was used as the reference gene. The values shown are the means ± SD of three replicates. Asterisks indicate statistically significant differences compared with the control group by a t-test (* p < 0.05; ** p < 0.01; *** p < 0.001).
Figure 9. Relative expression of mRNA, miRNA and lncRNA in Eucalyptus urophylla under different cold treatments. Purple bar: Eucgr.E00078; blue bar: EUC_00002677; green bar: MSTRG.7690. EuGADPH was used as the reference gene. The values shown are the means ± SD of three replicates. Asterisks indicate statistically significant differences compared with the control group by a t-test (* p < 0.05; ** p < 0.01; *** p < 0.001).
Forests 12 00836 g009
Table
Table 1. RNA-seq data for 18 samples.
Table 1. RNA-seq data for 18 samples.
SamplesRaw ReadsClean ReadsValid Ratio (%)Q30 (%)GC Content (%)Mapping Rate (%)
CK_154,799,92253,974,78098.4998.304989.17
CK_251,332,04850,536,53498.4598.3549.5089.51
CK_353,394,49852,372,19098.0998.4249.5088.59
LT6_155,730,66054,848,75498.4298.385089.61
LT6_251,042,17650,186,76498.3298.544989.61
LT6_352,132,55851,356,97298.5198.335090.25
LT12_144,118,45443,365,94098.2998.5549.5089.27
LT12_249,217,46648,395,61698.3398.695089.75
LT12_345,760,90044,991,23098.3298.625089.70
LT24_152,758,63451,927,87898.4398.6049.5089.61
LT24_246,546,27645,664,63898.1198.595089.02
LT24_341,767,60641,067,72098.3298.534989.41
LT48_153,199,23852,266,37698.2598.5149.5089.84
LT48_248,131,67647,281,36098.2398.5949.5089.87
LT48_351,295,65450,276,71298.0198.4949.5089.72
RE24_139,175,45238,513,80098.3198.505190.70
RE24_235,969,90635,382,10298.3798.2950.5090.51
RE24_344,145,60843,447,84098.4298.4350.5090.58
Table
Table 2. Statistical data of annotation lncRNAs and target genes.
Table 2. Statistical data of annotation lncRNAs and target genes.
Number of lncRNAsTarget GenesDifferentially Expressed Target GenesDE_lncRNAsTarget Genes of DE_lncRNAsDifferentially Expressed Target Genes of DE_lncRNAs
834719,41856063001681579
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, H.; Li, J.; Qiu, B.; Zhao, Y.; Liu, Z.; Yang, J.; Kang, X. Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake. Forests 2021, 12, 836. https://doi.org/10.3390/f12070836

AMA Style

Chen H, Li J, Qiu B, Zhao Y, Liu Z, Yang J, Kang X. Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake. Forests. 2021; 12(7):836. https://doi.org/10.3390/f12070836

Chicago/Turabian Style

Chen, Hao, Jiang Li, Bingfa Qiu, Yingwei Zhao, Zhao Liu, Jun Yang, and Xiangyang Kang. 2021. "Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake" Forests 12, no. 7: 836. https://doi.org/10.3390/f12070836

APA Style

Chen, H., Li, J., Qiu, B., Zhao, Y., Liu, Z., Yang, J., & Kang, X. (2021). Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake. Forests, 12(7), 836. https://doi.org/10.3390/f12070836

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