Next Article in Journal
Beneficial Effect of Fenofibrate and Silymarin on Hepatic Steatosis and Gene Expression of Lipogenic and Cytochrome P450 Enzymes in Non-Obese Hereditary Hypertriglyceridemic Rats
Next Article in Special Issue
Novel Therapeutic Effects in Rat Spinal Cord Injuries: Recovery of the Definitive and Early Spinal Cord Injury by the Administration of Pentadecapeptide BPC 157 Therapy
Previous Article in Journal
The Treatment with Interleukin 17 Inhibitors and Immune-Mediated Inflammatory Diseases
Previous Article in Special Issue
N-Terminus-Mediated Solution Structure of Dimerization Domain of PRC1
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 44
Issue 5
10.3390/cimb44050128
cimb-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

Transcriptomic Analysis Reveals LncRNAs Associated with Flowering of Angelica sinensis during Vernalization

by
Xiaoxia Liu
1,
Mimi Luo
1,
Mengfei Li
1,* and
Jianhe Wei
2,*
1
State Key Laboratory of Aridland Crop Science, College of Life Science and Technology, Gansu Agricultural University, Lanzhou 730070, China
2
Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100193, China
*
Authors to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2022, 44(5), 1867-1888; https://doi.org/10.3390/cimb44050128
Submission received: 24 March 2022 / Revised: 19 April 2022 / Accepted: 23 April 2022 / Published: 26 April 2022
(This article belongs to the Collection Feature Papers in Current Issues in Molecular Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Angelica sinensis is a “low-temperature and long-day” perennial plant that produces bioactive compounds such as phthalides, organic acids, and polysaccharides for various types of clinical agents, including those with cardio-cerebrovascular, hepatoprotective, and immunomodulatory effects. To date, the regulatory mechanism of flowering under the photoperiod has been revealed, while the regulatory network of flowering genes during vernalization, especially in the role of lncRNAs, has yet to be identified. Here, lncRNAs associated with flowering were identified based on the full-length transcriptomic analysis of A. sinensis at vernalization and freezing temperatures, and the coexpressed mRNAs of lncRNAs were validated by qRT-PCR. We obtained a total of 2327 lncRNAs after assessing the protein-coding potential of coexpressed mRNAs, with 607 lncRNAs aligned against the TAIR database of model plant Arabidopsis, 345 lncRNAs identified, and 272 lncRNAs characterized on the SwissProt database. Based on the biological functions of coexpressed mRNAs, the 272 lncRNAs were divided into six categories: (1) chromatin, DNA/RNA and protein modification; (2) flowering; (3) stress response; (4) metabolism; (5) bio-signaling; and (6) energy and transport. The differential expression levels of representatively coexpressed mRNAs were almost consistent with the flowering of A. sinensis. It can be concluded that the flowering of A. sinensis is positively or negatively regulated by lncRNAs, which provides new insights into the regulation mechanism of the flowering of A. sinensis.

1. Introduction

Angelica sinensis (Oliv.) Diels is a “low-temperature and long-day” perennial plant that is native to Gansu Province, Northwest China [1,2]. The roots have been used as a traditional Chinese medicine for over 2000 years [3,4]. Currently, the roots are also applied in cardio-cerebrovascular, hepatoprotective, and immunomodulatory clinical agents, largely relying on bioactive compounds such as phthalides, organic acids, and polysaccharides [5,6,7].
Recently, the planting area of A. sinensis has exceeded 40,000 ha to satisfy the increasing demand; however, a higher rate (>40%) of early bolting and flowering (EBF) in commercial cultivation increases the lignified rate of roots and decreases the yield accordingly [8,9]. In order to inhibit the EBF, effective measures that have been taken include selecting excellent germplasm resources with a lower rate of EBF [8], controlling the seedling size (0.4 to 0.6 cm) to delay the transition from vegetative to reproductive growth [10], storing seedlings below freezing temperatures (−3 to −10 °C) to avoid vernalization (0 to 10 °C) [2,11,12], and shading the plants with sunshade nets (40% to 60%) to avoid the long-day conditions during the adult stages [13].
Regarding the regulatory mechanism of flowering in A. sinensis, key genes and the regulatory network during the photoperiod have been identified and mapped based on transcriptomic analysis. Specifically, 13 genes associated with the photoperiod, vernalization, sucrose, and GA pathways were identified from plants at the vegetative stage compared with the EBF stage [14]; 38 genes associated with the photoperiod, carbohydrates, hormone signaling, and floral development were identified from different development stages [15]; and 40 genes associated with the photoperiod, sucrose, GA, and floral development were identified from the EBF compared with Un-EBF [16]. In summary, key genes such as FLC, SOC1, FT, PHYA, AP1, and GA2OX1 were identified during the transition from the vegetative to the reproductive stage [14,15,16]. To date, physiological changes in the levels of carbohydrates, proteins, and hormones during vernalization have been investigated [17,18], while the regulatory network of flowering has yet to be identified.
Generally, ncRNAs include two categories, housekeeping and regulatory ncRNAs, and the latter can be further divided into sRNAs (i.e., miRNA, siRNA, and piRNA) and lncRNAs (>200 nucleotides) [19]. Both miRNAs and lncRNAs can influence plant developmental processes and stress responses [20], with the former being negative regulators functioning as specificity determinants, or guides, within complexes that promote the degradation of mRNA targets, and the latter acting either as precursors of miRNAs or endogenous target mimics (TMs), which mimic the real targets of miRNAs, thus rendering the corresponding miRNAs ineffective [21]. For example, previous studies on resistance against leaf rust in wheat found that 50 miRNAs and 1178 lncRNAs were identified and 49 lncRNAs were found to be the targets for miRNAs, with 1 lncRNA acting as a precursor of 2 miRNAs, and 3 lncRNAs acting as TMs [22].
Extensive investigations have demonstrated that lncRNAs regulate their downstream targets’ expression through the changing of epigenetic modification at the level of transcription and post-transcription by interacting with DNA, RNA, and proteins; thus, they are involved in various biological processes [23,24,25]. In this study, lncRNAs associated with flowering were identified based on the transcriptomic analysis of A. sinensis seedlings treated at vernalization and freezing temperatures (avoiding vernalization). We found that 272 lncRNAs directly or indirectly participate in regulating the flowering of A. sinensis and stress responses.

2. Materials and Methods

2.1. Plant Material

The seedlings (root–shoulder diameter 0.4–0.5 cm; Figure S1) of Angelica sinensis (cultivar Mingui 1) were stored at 0 (vernalization temperature) and −3 °C (freezing temperature), respectively. After storage at 0 °C for 14 (T1) and 60 days (T2), as well as at −3 °C for 125 days (T3), the shoot apical meristem (SAM) was cut from the root shoulder of the seedlings for transcriptomic analysis and qRT-PCR validation. Three biological replicates were performed for each treatment of T1, T2, and T3. Herein, the treatment of T1, T2, and T3 represents uncompleted, completed, and avoided vernalization, respectively, based on the EBF rate (Figure S2) when the stored seedlings were cultivated and grown in a long-day condition.

2.2. Full-Length Isoform Sequencing and Analysis

Total RNA of the SAM samples was extracted using Trizol reagent (Omega Bio-Tek, Norcross, GA, USA). The integrity of the RNA was determined using an Agilent 2100 Bioanalyzer (Agilent Technol., California, CA, USA) and agarose gel electrophoresis, and the purity and concentration of the RNA were determined using a microspectrophotometer (NanoDrop Technol., Wilmington, DE, USA). The high-quality RNAs were sequenced on a Pacific Biosciences Sequel platform (Gene Denovo Biotechnology Co., Ltd., Guangzhou, China). Raw reads of cDNA library were analyzed using a SMRT Link (V8.0.0) [26]. Briefly, high-quality CCS were extracted from the subreads BAM file; the integrity of transcripts (full-length sequences) was judged based on whether CCS reads contained primers (5′ and 3’) and polyAs; then, FLNC reads were generated by removing primers, barcodes, and polyAs; finally, FLNC reads were assembled to obtain the entire isoform [27].

2.3. Analysis of Long Noncoding RNAs (lncRNAs)

Isoforms that were not annotated against the four databases—NR, Swiss-Prot, Kyoto KEGG, and KOG—were used for the analysis of lncRNAs. The isoform that was assessed as a noncoding transcript by both CNCI and CPC software was finally confirmed as a lncRNA [28,29].

2.4. Characterization of LncRNAs

To date, the genome of A. sinensis has not been sequenced. Thus, the lncRNA analysis of A. sinensis was performed via a BLAST search with an E-value cut-off of ≤1 × 10−5 against the known lncRNAs from the TAIR database (https://www.arabidopsis.org accessed on 30 March 2022) [30]. The function of lncRNAs was annotated based on their coexpression mRNAs [31,32,33]. Herein, the biological functions of the coexpressed mRNAs were searched on the UniProt database (https://www.uniprot.org accessed on 30 March 2022).

2.5. qRT-PCR Validation

Based on the coding sequences (CDS) of coexpressed mRNA of lncRNA, 49 primer sequences of representatively coexpressed mRNAs (Table 1) were designed using the NCBI primer-blast tool. First-strand cDNA synthesis and qRT-PCR reaction were carried out using SuperRealPreMix Plus (FP205; Tiangen Biotech., Beijing, China) according to the manufacturer’s instructions; specifically, the cDNA was synthesized successively with one cycle (95 °C, 15 min) and 40 cycles (95 °C, 10 s; 60 °C, 20 s; and 72 °C, 30 s), and the qRT-PCR reaction was incubated successively at 95 °C for 15 s, 60 °C for 1 min and 95 °C for 1 s. The Actin (ACT) gene was used as a reference control gene with forward: TGGTATTGTGCTGGATTCTGGT and reverse: TGAGATCACCACCAGCAAGG (amplicon size 109 bp) [34]. Herein, the cycle threshold (Ct) values and standard curves of the ACT gene at different volumes (0.25, 0.5, 1.0, 1.5, 2.0, and 3.0 μL) was built to correct the gene expression level (Figure S3 and Figure S4), and the expression levels of the 49 candidate genes and their standard deviations for every variant were added to the Supplementary Materials (Table S1). The REL of coexpressed mRNAs was calculated using the 2−△△Ct method [35] according to the following formula.
△CtTest gene = CtTest gene − CtReference gene
△CtTarget gene = CtTarget gene − CtReference gene
−△△Ct(T2 vs. T1) = −(△CtT2 − △CtT1)
−△△Ct(T3 vs. T1) = −(△CtT3 − △CtT1)
Relative expression level (REL) = 2−△△Ct.

2.6. Statistical Analysis

In order to obtain the precise estimation of PCR efficiency, each experiment for qRT-PCR validation was performed with three biological replicates, along with three technical replicates [36]. A t-test in SPSS 22.0 was performed for independent experiments, with p < 0.05 as the basis for statistical differences.

3. Results

3.1. LncRNAs Analysis

In total, 2327 lncRNAs were obtained after assessing the protein-coding potential of coexpressed mRNAs based on the two software programs, CNCI and CPC (Figure 1A), with 607 genes aligned against the known lncRNAs from the TAIR database of model plant Arabidopsis (Figure 1B), and 345 lncRNAs with coexpressed mRNAs of A. sinensis identified (Figure 1C) based on the SwissProt database. Based on the biological functions, the 272 characterized lncRNAs (Figure 1D) were divided into six categories: chromatin, DNA/RNA and protein modification (29); flowering (36); stress response (24); metabolism (117); biosignaling (23), and energy and transport (43) (Figure 1E). The base sequences of the 272 lncRNAs are shown in Table S2.

3.2. LncRNAs Linked with Chromatin, DNA/RNA and Protein Modification, as well as Expression Levels of Their Coexpressed mRNAs

Based on the biological functions of coexpressed mRNAs, 29 lncRNAs were linked with chromatin (HMGB2 and HMGB3), DNA/RNA (At1g05910, RID2 and At4g26600) and protein modification (H2AV, At2g28720, HTR2, HOS15, REF6, SKP1A, ASK21, SRK2G, SRK2H, DET1, BOPAt4g295601, At1g45180, At3g50840, At2g36630, At3g47890, UBP7, DER2.1, GRP3, MDH9.13, At3g24715, At3g16560, At3g62260, ESMD1, and At1g27930) (Table 2). The expression levels of 14 representative coexpressed mRNAs were confirmed by qRT-PCR, with 3 mRNAs (HMGB2, HMGB3 and At1g05910) showing down-regulation at T2 versus (vs.) T1, and 11 mRNAs showing lower levels at T2 vs. T1 than T3 vs. T1, with the exception of 2 mRNAs (H2AV and ASK21) (Figure 2).

3.3. LncRNAs Linked with Flowering and Expression Levels of Their Coexpressed mRNAs

In total, 36 lncRNAs were linked with flowering based on the biological functions of their coexpressed mRNAs, with 12 lncRNAs directly associated with flowering, namely SRR1, PHL, PHYA, AGL62, AGL79, ATJ3, BBX29, CLE13, CLE44, MXC17.10, At1g06515, and BHLH30 (Table 3), and 24 lncRNAs indirectly associated with flowering, such as cell division, embryo development, and cell wall organization (Table S3). The expression levels of six representative coexpressed mRNAs (SRR1, PHL, PHYA, AGL62, AGL79, and ATJ3) were confirmed by qRT-PCR, with all six mRNAs showing down-regulation at T2 vs. T1 and T3 vs. T1, and five mRNAs showing lower levels at T2 vs. T1 than T3 vs. T1, with the exception of the gene AGL62 (Figure 3).

3.4. LncRNAs Linked with Stress Response and Expression Levels of Their Coexpressed mRNAs

In total, 24 lncRNAs were linked with the stress response based on the biological functions of their coexpressed mRNAs, with 14 lncRNAs associated with the temperature response, namely ACBP6, ENO2, ADH1, CSP2, RH20, RH52, RH53, RAB18, XERO1, MED14, HSP17.8, HSP70-3, HSP70-10, and HSP90-3 (Table 4), and 10 lncRNAs associated with other stresses responses, such as water, salt, and oxidative stress (Table S4). The expression levels of eight representative coexpressed mRNAs involved in the temperature response were confirmed by qRT-PCR, with four mRNAs (ACBP6, ENO2, CSP2, and HSP90-3) showing up-regulation at T3 vs. T1, and three mRNAs (ACBP6, ENO2, and CSP2) showing lower levels at T2 vs. T1 than T3 vs. T1 (Figure 4).

3.5. LncRNAs Linked with Metabolism and Expression Levels of Their Coexpressed mRNAs

In total, 117 lncRNAs were linked with metabolism based on the biological functions of their coexpressed mRNAs, with 19 lncRNAs associated with carbohydrate metabolism, namely CSY4, FBA3, GAPC1, GAPCP2, At3g52990, PGM2, USP, GLCNAC1PUT2, UXS2, XYLA, GALS1, OFUT31, At5g67460, RSW3, At1g59950, At5g25970, UGT76E7, At1g26850, and BAM1 (Table 5), and 98 lncRNAs associated with other types of metabolism, such as nucleotide, protein, and lipid metabolism (Table S5). The expression levels of five representative coexpressed mRNAs (CSY4, GAPCP2, At3g52990, PGM2, and BAM1) involved in carbohydrate metabolism were confirmed by qRT-PCR, with all five mRNAs showing down-regulation at T3 vs. T1, and three mRNAs showing higher levels at T2 vs. T1 than T3 vs. T1, with the exception of the two genes GAPCP2 and BAM1 (Figure 5).

3.6. LncRNAs Linked with Biosignaling and Expression Levels of Their Coexpressed mRNAs

In total, 23 lncRNAs were linked with metabolism based on the biological functions of their coexpressed mRNAs, with 13 lncRNAs associated with hormone signaling, namely ARF1, CUL1, T4L20.330, SOFL4, SOFL5, GRF2, GRF11, ERF3, CIPK20, SF1, AGD9, TIFY4B, and At2g34810 (Table 6), and 10 lncRNAs associated with other types of signaling, such as protein kinase, phosphatidylinositol-mediated, and cell surface receptor signaling (Table S6). The expression levels of eight representative coexpressed mRNAs associated with hormone signaling were confirmed by qRT-PCR, with three mRNAs (ARF1, CUL1, and GRF11) showing up-regulation at T2 vs. T1 and T3 vs. T1, and four mRNAs (ARF1, SOFL4, GRF2, and ERF3) showing lower levels at T2 vs. T1 than T3 vs. T1 (Figure 6).

3.7. LncRNAs Linked with Energy and Transport

In all, 43 lncRNAs were found to be linked with energy and transport based on the biological functions of their coexpressed mRNAs, with 5 lncRNAs associated with energy, such as PURU1, ndhB1, and NDB1, and 38 lncRNAs associated with transport, such as At5g11230, GPT1 and SMXL5 (Table 7). The expression levels of eight representative coexpressed mRNAs (PURU1, MES16, GPT1, ABCF5, SECA2, VPS26A, CML19, and NHX6) were confirmed by qRT-PCR, with three mRNAs (PURU, GPT1, and VPS26A) showing up-regulation at T2 vs. T1 and T3 vs. T1, and six mRNAs (PURU1, GPT1, ABCF5, SECA2, VPS26A, and CML19) showing higher levels at T2 vs. T1 than T3 vs. T1, with the exception of the two genes MES16 and NHX6 (Figure 7).

4. Discussion

Vernalization is a process considered to be an epigenetic switch whereby flowering is promoted by prolonged exposure to cold (0 to 10 °C); meanwhile, it can be lost at high temperatures or avoided below freezing temperatures [37,38]. Epigenetic regulation involves diverse molecular mechanisms including chromatin remodeling, DNA methylation, histone modification, histone variants, and ncRNAs [39]. Studies on Brassica rapa found that 127 differentially expressed lncRNAs were coexpressed with 128 differentially expressed genes, indicating that lncRNAs played an important role during vernalization [40]. In this study, 272 characterized lncRNAs were identified from A. sinensis and divided into six categories, namely (1) chromatin, DNA/RNA, and protein modification; (2) flowering; (3) stress response; (4) metabolism; (5) biosignaling; and (6) energy and transport, based on their coexpressed mRNAs.
FLC is a MADS-box transcriptional regulator that acts as a potent repressor of flowering [38]. In Arabidopsis, the epigenetic silencing of the floral repressor gene FLC is a well-characterized key event of vernalization [41]. In this study, 29 lncRNAs linked with chromatin, DNA/RNA, and protein modification were identified in A. sinensis during vernalization. For the chromatin modification, two coexpressed mRNAs, HMGB2 and HMGB3, are involved in binding preferentially double-stranded DNA and up-regulated in response to cold stress [42]. For the DNA/RNA modification, At1g05910 is involved in DNA demethylation and the negative regulation of chromatin silencing [43]; RID2 is involved in rRNA methyltransferase activity [44]; and At4g26600 is involved in RNA methylation [45]. For protein modification, 24 lncRNAs were involved and the roles of nine coexpressed mRNAs were represented as follows. H2AV plays a central role in regulating transcription, repairing DNA, replicating DNA, and stabilizing the nucleus chromosome [46]; At2g28720 and HTR2 are involved in compacting DNA into chromatin [47,48]; HOS15 promotes the deacetylation of histone H4 [49]; REF6 is involved in demethylating ‘Lys-27’ of histone H3, regulating flowering time by repressing FLC expression and interacting with the NF-Y complex to regulate SOC1 [50,51]; SKP1A and ASK21 are involved in ubiquitination and form an SCF E3 ubiquitin ligase complex together with CUL1, RBX1, and an F-box protein [52,53]; and SRK2G and SRK2H are involved in protein phosphorylation [54,55]. In previous studies, HMGB2 and HMGB3 were up-regulated in response to cold stress but down-regulated in response to drought and salt stresses [42]; H2AV, At2g28720 and HTR2 exhibited a high level in response to osmotic and drought stresses [56]; HOS15 was found to act as a repressor of cold stress-regulated gene expression and played a role in gene regulation for plant acclimation and tolerance to cold stress [49]; SKP1A and ASK21 were overexpressed in the host stress response [57]; and SRK2G and SRK2H were found to be positive regulators in stress responses such as drought, salt, and cold [54,55,58]. Currently, although most of the lncRNAs have been reported to be involved in the stress response in many plants, their roles in the regulation of flowering time have been studied in model plants [23]. Previous studies have found that FLC in Arabidopsis is epigenetically regulated by lncRNAs COOLAIR and COLDAIR [59,60,61]. In this study, a lower expression level was noted for almost all coexpressed mRNAs of lncRNAs involved in chromatin, DNA/RNA, and protein modification at vernalization (T2, 0 °C) compared with freezing temperature (T3, −3 °C), as well as down-regulation at T2, indicating that these lncRNAs participate in epigenetic silencing by transferring euchromatin to heterochromatin and confer early bolting and flowering of A. sinensis. Representatively, the down-regulation of mRNAs HMGB2 and HMGB3 during vernalization transfers the heterochromatin to the euchromatin [42], which generates the ability of flowering; in contrast, their up-regulation at freezing temperatures inhibits flowering by keeping the heterochromatin, which indicates that their coexpressed lncRNAs play positive roles in regulating the flowering of A. sinensis. The down-regulation of mRNA REF6 below 0 °C weakens the demethylation of histone H3 and delays flowering time by inducing FLC expression [50,51]; meanwhile, there is an increased expression level with decreased temperatures, which indicates that this coexpressed lncRNA plays a negative role in regulating the flowering of A. sinensis.
Flower formation occurs at the SAM and is a complex morphological event that is required not only for the circadian clock to measure the passage of time but also the regulation of meristem identity genes [42]. For the 6 representative coexpressed mRNAs of the 12 lncRNAs directly linked with flowering, SRR1 is involved in a circadian clock input pathway and regulation of the expression of clock-regulated genes such as CCA1 and TOC1 [62]; PHL is involved in the circadian rhythm and the regulation of the timing of transition [63]; PHYA is involved in the regulation of flowering time and expression of its own gene as negative feedback [64]; AGL62 is required for promoting the nuclear proliferation of early endosperm [65]; AGL79 is involved in positively regulating the transition of the meristem from the vegetative to reproductive phase [66]; and ATJ3 plays a continuous role in plant development, such as in photoperiodism, flowering, and positive regulation of flower development [67]. Previous studies found that mRNAs (e.g., miR156, miR169 and miR172) play a crucial role in developmental processes in rice, wheat, and maize, especially in the formation of the floral meristem, with miR172 controlling AP2-like genes [23,68,69]. Studies on the flowering of Chenopodium quinoa found that pivotal flowering homologs, including photoreceptor genes PHYA and CRY1, as well as genes associated with florigen-encoding genes (FT and TWIN SISTER of FT) and circadian clock-related genes (ELF3, LHY, and HY5), were specifically affected by night-break and competed with the positive- and negative-flowering lncRNAs [70]. In this study, down-regulation of all the coexpressed mRNAs involved in circadian clock and meristem identity genes was observed, which can be considered acceptable and reasonable, because these mRNAs are often highly expressed at the plant development stage (at photoperiod), while their expression levels were examined during vernalization. In addition, increased expression with decreased temperatures indicates that their coexpressed lncRNAs play negative roles in regulating the flowering of A. sinensis.
The expression of numerous lncRNAs has been demonstrated to be significantly affected by various stresses [23,30]. During vernalization, plants have to face and adapt to low temperature [38]. To date, extensive studies have reported that lncRNAs participate in defense responses associated with plant immunity and adaptation to the environment [22]. Heat-responsive lncRNAs have been found to be differentially expressed in Brassica juncea, and cold-responsive lncRNAs have been identified in grape and Arabidopsis [71,72,73]. lncRNAs could regulate HSP family genes (HSP82 and HSP83) in response to heat stress in Populus x canadensis Moench, and HSP18.1 in response to Cd stress Betula platyphylla [74,75]. For the eight representative coexpressed mRNAs of the 14 lncRNAs linked with the temperature response, ACBP6 confers resistance to cold and freezing [76]; ENO2 acts as a positive regulator in response to cold stress [77]; ADH1 is required for survival and acclimation in response to abiotic stresses (e.g., cold, salt, and drought) [78,79]; CSP2 contributes to the enhancement of cold and freezing tolerance [80]; and HSP17.8, HSP70-3, HSP70-10, and HSP90-3 play vital roles in adapting to biotic and abiotic stresses [81,82]. In this study, increased expression for cold-tolerated mRNAs (ACBP6, ENO2, ADH1, and CSP2), and decreased expression for heat-tolerated mRNAs (HSP17.8, HSP70-3, HSP70-10, and HSP90-3), were observed with decreased temperatures, which indicates that their coexpressed lncRNAs play positive roles in adapting to low temperatures.
For vernalization to occur, sources of energy (sugars) and carbohydrate metabolism are required [37]. In recent years, the roles of lncRNAs in regulating metabolism in cancer, insulin, and chicken have been reported [83,84,85], while, in plants, studies are still limited. For the five representative coexpressed mRNAs of the 19 lncRNAs linked with carbohydrate metabolism, CSY4 is involved in the synthesis of socitrate from oxaloacetate [86]; GAPCP2 plays a critical role in glycolysis and exhibits up-regulation under drought stress [87,88]; At3g52990 is involved in the synthesis of pyruvate from D-glyceraldehyde 3-phosphate [89]; PGM2 participates in the synthesis of glucose [90]; and BAM1 is required for starch breakdown [91]. In this study, the differential expression of these coexpressed mRNAs regulating sucrose and starch metabolism provided energy for the morphogenesis of seedlings and adaptation to low temperatures during vernalization. Representatively, the decreased expression of metabolite-synthesized mRNAs (CSY4, At3g52990, and PGM2), and increased expression of energy-produced mRNA GAPCP2 and metabolite-degraded mRNA BAM1, were observed with decreased temperatures, which indicates that their coexpressed lncRNAs play positive roles in carbohydrate metabolism.
Endogenous hormones such as gibberellin, auxin, cytokinin, brassinosteroid and abscisic acid can either inhibit or promote flowering [37]. In previous studies, a pre-miRNA of miR393 was identified in Brassica rapa during vernalization, and the overexpression of an miR393-resistant form of TIR1 (mTIR1) could enhance auxin sensitivity, thus leading to pleiotropic effects on plant development [92]. For the 8 representative coexpressed mRNAs of the 13 lncRNAs linked with hormone signaling, ARF1 is involved in the recruitment of COPI and GDAP1 to membranes and various auxin-dependent developmental processes [93]; CUL1 participates in forming a SCF complex together with SKP1, RBX1, and a F-box protein and is involved in floral organ development, the auxin signaling pathway and ethylene signaling [94]; SOFL4 and SOFL5 are involved in cytokinin-mediated development [95]; GRF2 and GRF11 participate in the brassinosteroid (BR)-mediated signaling pathway [96,97]; ERF3 is found to be differentially expressed in response to stresses and also regulates other ERFs [98]; and SF1 is required for development and is involved in the alternative splicing of FLM pre-mRNA [99]. In this study, the differential expression of these mRNAs involved in hormone signaling also played certain roles in regulating the flower-bud differentiation of seedlings and cold tolerance during vernalization. In previous studies, GRF2 and ERF3 were found to be down-regulated, while GRF11 was up-regulated in response to cold stress [100,101]; here, contrary findings for mRNAs GRF2 and GRF11 were observed with temperatures decreased, which indicate that their coexpressed lncRNAs may play negative roles in hormone signaling. The interaction between SF1 and FLM pre-mRNA controls flowering time in response to temperature [102]; here, decreased expression of mRNA SF1 was observed with decreased temperatures, which indicates that this coexpressed lncRNA may play a positive role in hormone signaling.
Energy generation and the transport of energy, nutrients and metabolites play essential roles in growth and development and stress tolerance [103]. For the 8 representative coexpressed mRNAs of the 43 lncRNAs linked with energy and transport, PURU1 is involved in photorespiration and participates in preventing the excessive accumulation of 5-formyl tetrahydrofolate [104]; MES16 is involved in chlorophyll breakdown by its demethylation [105]; GPT1 is involved in NADPH generation through a series of processes including Glc6P transport, starch biosynthesis, fatty acid biosynthesis, and oxidative pentose phosphate [106]; ABCF5 belongs to the ABC transporter superfamily and is involved in protein transport [107]; SECA2 is involved in protein export coupling ATP hydrolysis [108]; VPS26A plays a role in vesicular protein sorting and is essential in endosome-to-Golgi retrograde transport [109]; CML19 is a potential calcium sensor that binds calcium and is involved in the early response to stress [110]; and NHX6 is involved in trafficking to the vacuole and exchanging the low-affinity electroneutral Na(+) or K(+)/H(+) [111]. In this study, the differential expression of these coexpressed mRNAs associated with energy and transport provided energy, nutrients, and metabolites for A. sinensis seedlings to obtain the capacity for vernalization and, meanwhile, to adapt to low temperatures.
Based on the above functional analysis of lncRNAs identified in this study, flowering pathways proposed in the previous literature [14,15,16] and a general model of stress-responsive regulation by regulatory lncRNAs [23], a model of vernalization-induced bolting and flowering by regulatory lncRNAs in A. sinensis is proposed (Figure 8). Briefly, the vernalization of seedlings firstly triggers the differential expression of lncRNAs; then, the lncRNAs either act as a precursor of miRNAs or as a miRNA target mimic, which further binds their related targets; then, the binding of targets regulates the expression of their downstream mRNAs that are involved in various biological processes, including the temperature response, flowering pathways (i.e., epigenetic modification, flowering induction, carbohydrate metabolism, and hormone signaling), as well as energy and transport; finally, these biological processes promote the transition of the meristem from the vegetative phase to the bolting and flowering of A. sinensis.

5. Conclusions

From the above observations, we found that the lncRNAs positively or negatively regulated the expression of their downstream mRNAs through epigenetic changes at the level of transcription and post-transcription for the flowering of A. sinensis during vernalization. This coexpressed mRNA analysis of lncRNAs focused on five pathways, namely (1) chromatin, DNA/RNA, and protein modification; (2) floral development; (3) temperature response; (4) carbohydrate metabolism; and (5) hormone signaling. While several candidate lncRNAs were identified, their causative roles require further investigations.

Supplementary Materials

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

Author Contributions

X.L.: data curation, formal analysis, validation, and writing—original draft preparation; M.L. (Mimi Luo): formal analysis and validation; M.L. (Mengfei Li): conceptualization, project administration, supervision, and writing—review and editing; J.W.: funding acquisition and resources. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the State Key Laboratory of Aridland Crop Science/Gansu Agricultural University (GSCS-2021-Z03), National Natural Science Foundation of China (32160083), China Agriculture Research System of MOF and MARA (CARS-21), and Assurance Project of Ecological Planting and Quality of Daodi Herbs (202103003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets are available at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA789039 (release date on 13 January 2022).

Conflicts of Interest

All the authors declare no conflicts of interest.

Abbreviations

AP1APETALA 1
CCSCircular consensus sequence
CNCICoding–Noncoding Index
COLDAIRCOLD-ASSISTED INTRONIC NON-CODING RNA
COOLAIRCOLD-INDUCED LONG ANTISENSE INTRAGENIC RNAs
CPCCoding-Potential Calculator
FLCFLOWERING LOCUS C
FLMFLOWERING LOCUS M
FLNCFull-length nonchimeric
FTFLOWERING LOCUS T
GAGibberellin
GA2OX1Gibberellin 2-β-dioxygenase 1
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes
KOGeuKaryotic orthologous groups of proteins
lncRNAsLong noncoding RNAs
NCBINational Center for Biotechnology Information
ncRNAsNoncoding RNAs
NRNCBI nonredundant protein
PHYAPHYOCHROME A
RELRelative expression level
SOC1SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1
sRNAsSmall RNAs
TAIRThe Arabidopsis Information Resource

References

  1. Zhang, H.Y.; Bi, W.G.; Yu, Y.; Liao, W.B. Angelica sinensis (Oliv.) Diels in China: Distribution, cultivation, utilization and variation. Genet. Resour. Crop Evol. 2012, 59, 607–613. [Google Scholar] [CrossRef]
  2. Wang, W.J. Analysis and control of early bolting characteristic of Angelica sinensis. J. Northwest Univ. 1977, 7, 32–39. [Google Scholar]
  3. Committee for the Pharmacopoeia of PR China. Pharmacopoeia of the People’s Republic of China; Chinese Medical Science and Technology Press: Beijing, China, 2015; p. 133. [Google Scholar]
  4. Hook, I.L. Danggui to Angelica sinensis root: Are potential benefits to European women lost in translation? A review. J. Ethnopharmacol. 2014, 152, 1–13. [Google Scholar] [CrossRef]
  5. Upton, R. American Herbal Pharmacopoeia and Therapeutic Compendium: Dang Gui Root-Angelica sinensis (Oliv.); American Herbal Pharmacopoeia: Scotts Valley, CA, USA, 2003; pp. 1–41. [Google Scholar]
  6. Ma, J.P.; Guo, Z.-B.; Jin, L.; Li, Y.D. Phytochemical progress made in investigations of Angelica sinensis (Oliv.) Diels. Chin. J. Nat. Med. 2015, 13, 241–249. [Google Scholar] [CrossRef]
  7. Wei, W.L.; Zeng, R.; Gu, C.M.; Qu, Y.; Huang, L.F. Angelica sinensis in China-A review of botanical profile, ethnopharmacology, phytochemistry and chemical analysis. J. Ethnopharmacol. 2016, 190, 116–141. [Google Scholar] [CrossRef]
  8. Huang, L.Q.; Jin, L. Suitable Technology for Production and Processing of Angelica sinensis; China Pharmaceutical Science and Technology Press: Beijing, China, 2018; pp. 1–14. [Google Scholar]
  9. Li, M.L.; Cui, X.W.; Jin, L.; Li, M.F.; Wei, J.H. Bolting reduces ferulic acid and flavonoid biosynthesis and induces root lignification in Angelica sinensis. Plant Physiol. Biochem. 2022, 170, 171–179. [Google Scholar] [CrossRef] [PubMed]
  10. Lin, H.M.; Qiu, D.Y.; Chen, Y. Effect of root diameter on early bolting rate and yield in seedling of Angelica sinensis. Chin. Tradit. Herb. Drugs 2007, 38, 1386–1389. [Google Scholar]
  11. Wang, W.J. Technology and principle of seedling frozen storage of Angelica sinensis. J. Chin. Med. Mat. 1979, 3, 1–4. [Google Scholar]
  12. Jia, Z.; Di, S.Q.; Zhao, F.Y.; Li, S.X.; Wang, S.C. Effects of different low overwintering temperatures on Angelica vernalization and premature bolting. Agr. Sci. Tech. 2018, 19, 55–62. [Google Scholar]
  13. Yao, L. Effect of shading during the nursery of Angelica sinensis on bolting rate and economic characters. Gansu Agr. Sci. Tech. 2005, 10, 54–55. [Google Scholar]
  14. Yu, G.; Zhou, Y.; Yu, J.J.; Hu, X.Q.; Tang, Y.; Yan, H.; Duan, J.A. Transcriptome and digital gene expression analysis unravels the novel mechanism of early flowering in Angelica sinensis. Sci. Rep. 2019, 9, 10035. [Google Scholar] [CrossRef] [PubMed]
  15. Li, J.; Li, M.L.; Zhu, T.T.; Zhang, X.N.; Li, M.F.; Wei, J.H. Integrated transcriptomics and metabolites at different growth stagesreveals the regulation mechanism of bolting and flowering of Angelica sinensis. Plant Biol. 2021, 23, 574–582. [Google Scholar] [CrossRef] [PubMed]
  16. Li, M.F.; Li, J.; Wei, J.H.; Paré, P.W. Transcriptional controls for early bolting and flowering in Angelica sinensis. Plants 2021, 10, 1931. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, E.H.; Huang, P. Effects of vernalization treatment on physiological character of Angelica sinensis seedlings. J. Gansu Agric. Univ. 1998, 33, 240–243. [Google Scholar]
  18. Chen, H.G.; Du, T.; Zhu, T.T.; Gao, S.F.; Chai, L.; He, W.W. Study on physiological mechanisms in frozen storage to reduce early bolting of Angelica sinensis. Mod. Tradit. Chin. Med. Mater. Med. World Sci. Tech. 2014, 16, 203–206. [Google Scholar]
  19. 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]
  20. Axtell, M.J.; Westholm, J.O.; Lai, E.C. Vive la différence: Biogenesis and evolution of microRNAs in plants and animals. Genome Biol. 2011, 12, 221. [Google Scholar] [CrossRef] [Green Version]
  21. Carrington, J.C.; Ambros, V. Role of microRNAs in plant and animal development. Science 2003, 301, 336–338. [Google Scholar] [CrossRef] [Green Version]
  22. Jain, N.; Sinha, N.; Krishna, H.; Singh, P.K.; Gautam, T.; Prasad, P.; Balyan, H.S.; Gupta, P.K. A study of miRNAs and lncRNAs during Lr28-mediated resistance against leaf rust in wheat (Triticum aestivum L.). Physiol. Mol. Plant Pathol. 2020, 112, 101552. [Google Scholar] [CrossRef]
  23. Sharma, Y.; Sharma, A.; Madhu; Shumayla; Singh, K.; Upadhyay, S.K. Long non-coding RNAs as emerging regulators of pathogen response in plants. Non Coding RNA 2022, 8, 4. [Google Scholar] [CrossRef]
  24. Qi, L.; Li, X.; Zhang, S.; An, D. Genetic regulation by non-coding RNAs. Sci. China Ser. C Life Sci. 2006, 49, 201–217. [Google Scholar] [CrossRef]
  25. Mercer, T.R.; Dinger, M.E.; Mattick, J.S. Long non-coding RNAs: Insights into functions. Nat. Rev. Genet. 2009, 10, 155–159. [Google Scholar] [CrossRef] [PubMed]
  26. Gordon, S.P.; Tseng, E.; Salamov, A.; Zhang, J.; Meng, X.; Zhao, Z.; Kang, D.; Underwood, J.; Grigoriev, I.V.; Figueroa, M.; et al. Widespread polycistronic transcripts in fungi revealed by single-molecule mRNA sequencing. PLoS ONE 2015, 10, e0132628. [Google Scholar] [CrossRef] [Green Version]
  27. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.D.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Sun, L.; Luo, H.; Bu, D.; Zhao, G.; Yu, K.; Zhang, C.; Liu, Y.; Chen, R.; Zhao, Y. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. [Google Scholar] [CrossRef] [PubMed]
  29. Kong, L.; Zhang, Y.; Ye, Z.Q.; Liu, X.Q.; Zhao, S.Q.; Wei, L.; Gao, G. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35, W345–W349. [Google Scholar] [CrossRef] [PubMed]
  30. Poole, R.L. The TAIR database. Methods Mol. Biol. 2007, 406, 179–212. [Google Scholar] [CrossRef] [PubMed]
  31. Liao, Q.; Xiao, H.; Bu, D.C.; Xie, C.Y.; Miao, R.Y.; Luo, H.T.; Zhao, G.G.; Yu, K.T.; Zhao, H.T.; Skogerbø, G.; et al. ncFANs: A web server for functional annotation of long non-coding RNAs. Nucleic Acids Res. 2011, 39, w118–w124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Mattick, J.S.; Rinn, J.L. Discovery and annotation of long noncoding RNAs. Nat. Struct. Mol. Biol. 2015, 22, 5–7. [Google Scholar] [CrossRef]
  33. Shumayla; Sharma, S.; Taneja, M.; Tyagi, S.; Singh, K.; Upadhyay, S.K. Survey of high throughput RNA-Seq data reveals potential roles for lncRNAs during development and stress response in bread wheat. Front. Plant Sci. 2017, 8, 1019. [Google Scholar] [CrossRef] [Green Version]
  34. Xu, R.; Xu, J.; Li, Y.C.; Dai, Y.T.; Zhang, S.P.; Wang, G.; Liu, Z.G.; Dong, L.L.; Chen, S.L. Integrated chemical and transcriptomic analyses unveils synthetic characteristics of different medicinal root parts of Angelica sinensis. Chin. Herb. Med. 2020, 12, 19–28. [Google Scholar] [CrossRef]
  35. Willems, E.; Leyns, L.; Vandesompele, J. Standardization of real-time PCR gene expression data from independent biological replicates. Anal. Biochem. 2008, 379, 127–129. [Google Scholar] [CrossRef] [PubMed]
  36. Svec, D.; Tichopad, A.; Novosadova, V.; Pfaffl, M.W.; Kubista, M. How good is a PCR efficiency estimate: Recommendations for precise and robust qPCR efficiency assessments. Biomol. Detect. Quantif. 2015, 3, 9–16. [Google Scholar] [CrossRef] [Green Version]
  37. Taiz, L.; Zeiger, E. Plant physiology. In The Control of Flowering, 5th ed.; Fosket, D.E., Amasino, R., Eds.; Sinauer Associates, Inc.: Sunderland, MA, USA, 2010; pp. 559–590. [Google Scholar]
  38. Amasino, R.M. Vernalization and flowering time. Curr. Opin. Biotechnol. 2005, 16, 154–158. [Google Scholar] [CrossRef] [PubMed]
  39. He, K.; Cao, X.F.; Deng, X. Histone methylation in epigenetic regulation and temperature responses. Curr. Opin. Plant Biol. 2021, 61, 102001. [Google Scholar] [CrossRef] [PubMed]
  40. Liu, T.K.; Wu, P.; Wang, Q.; Wang, W.L.; Zhang, C.W.; Sun, F.F.; Liu, Z.K.; Li, Y.; Hou, X.L. Comparative transcriptome discovery and elucidation of the mechanism of long noncoding RNAs during vernalization in Brassica rapa. Plant Growth Regul. 2018, 85, 27–39. [Google Scholar] [CrossRef]
  41. Velanis, C.N.; Goodrich, J. Vernalization and epigenetic inheritance: A game of histones. Curr. Biol. 2017, 27, R1324–R1326. [Google Scholar] [CrossRef]
  42. Kwak, K.J.; Kim, J.Y.; Kim, Y.O.; Kang, H. Characterization of transgenic Arabidopsis plants overexpressing high mobility group b proteins under high salinity, drought or cold stress. Plant Cell Physiol. 2007, 48, 221–231. [Google Scholar] [CrossRef]
  43. Theologis, A.; Ecker, J.R.; Palm, C.J.; Federspiel, N.A.; Kaul, S.; White, O.; Alonso, J.; Altafi, H.; Araujo, R.; Bowman, C.L.; et al. Sequence and analysis of chromosome 1 of the plant Arabidopsis thaliana. Nature 2000, 408, 816–820. [Google Scholar] [CrossRef]
  44. Ohbayashi, I.; Konishi, M.; Ebine, K.; Sugiyama, M. Genetic identification of Arabidopsis RID2 as an essential factor involved in pre-rRNA processing. Plant J. 2011, 67, 49–60. [Google Scholar] [CrossRef]
  45. Cheng, C.Y.; Krishnakumar, V.; Chan, A.P.; Thibaud-Nissen, F.; Schobel, S.; Town, C.D. Araport11: A complete reannotation of the Arabidopsis thaliana reference genome. Plant J. 2017, 89, 789–804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Dai, X.Z.; Bai, Y.H.; Zhao, L.H.; Dou, X.Y.; Liu, Y.H.; Wang, L.L.; Li, Y.; Li, W.M.; Hui, Y.N.; Huang, X.Y.; et al. H2A.Z represses gene expression by modulating promoter nucleosome structure and enhancer histone modifications in Arabidopsis. Mol. Plant 2017, 10, 1274–1292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Jiang, D.; Borg, M.; Lorković, Z.J.; Montgomery, S.A.; Osakabe, A.; Yelagandula, R.; Axelsson, E.; Berger, F. The evolution and functional divergence of the histone H2B family in plants. PLoS Genet. 2020, 16, e1008964. [Google Scholar] [CrossRef] [PubMed]
  48. Sridhar, V.V.; Kapoor, A.; Zhang, K.L.; Zhu, J.J.; Zhou, T.; Hasegawa, P.M.; Bressan, R.A.; Zhu, J.K. Control of DNA methylation and heterochromatic silencing by histone H2B deubiquitination. Nature 2007, 447, 735–738. [Google Scholar] [CrossRef]
  49. Zhu, J.; Jeong, J.C.; Zhu, Y.; Sokolchik, I.; Miyazaki, S.; Zhu, J.K.; Hasegawa, P.M.; Bohnert, H.J.; Shi, H.; Yun, D.J.; et al. Involvement of Arabidopsis HOS15 in histone deacetylation and cold tolerance. Proc. Natl. Acad. Sci. USA 2008, 105, 4945–4950. [Google Scholar] [CrossRef] [Green Version]
  50. Lu, F.L.; Cui, X.; Zhang, S.B.; Jenuwein, T.; Cao, X.F. Arabidopsis REF6 is a histone H3 lysine 27 demethylase. Nat. Genet. 2011, 43, 715–719. [Google Scholar] [CrossRef]
  51. Hou, X.L.; Zhou, J.N.; Liu, C.; Liu, L.; Shen, L.S.; Yu, H. Nuclear factor Y-mediated H3K27me3 demethylation of the SOC1 locus orchestrates flowering responses of Arabidopsis. Nat. Commun. 2014, 5, 4601. [Google Scholar] [CrossRef] [Green Version]
  52. Zhao, D.Z.; Yu, Q.L.; Chen, M.; Ma, H. The ASK1 gene regulates B function gene expression in cooperation with UFO and LEAFY in Arabidopsis. Development 2001, 128, 2735–2746. [Google Scholar] [CrossRef]
  53. Liu, F.; Ni, W.; Griffith, M.E.; Huang, Z.; Chang, C.; Peng, W.; Ma, H.; Xie, D. The ASK1 and ASK2 genes are essential for Arabidopsis early development. Plant Cell 2004, 16, 5–20. [Google Scholar] [CrossRef] [Green Version]
  54. Yoshida, R.; Umezawa, T.; Mizoguchi, T.; Takahashi, S.; Takahashi, F.; Shinozaki, K. The regulatory domain of SRK2E/OST1/SnRK2.6 interacts with ABI1 and integrates abscisic acid (ABA) and osmotic stress signals controlling stomatal closure in Arabidopsis. J. Biol. Chem. 2006, 281, 5310–5318. [Google Scholar] [CrossRef] [Green Version]
  55. Yan, J.; Wang, P.C.; Wang, B.S.; Hsu, C.C.; Tang, K.; Zhang, H.R.; Hou, Y.J.; Zhao, Y.; Wang, Q.M.; Zhao, C.Z.; et al. The SnRK2 kinases modulate miRNA accumulation in Arabidopsis. PLoS Genet. 2017, 13, e1006753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Sura, W.; Kabza, M.; Karlowski, W.M.; Bieluszewski, T.; Kus-Slowinska, M.; Pawełoszek, Ł.; Sadowski, J.; Ziolkowski, P.A. Dual role of the histone variant H2A.Z in transcriptional regulation of stress-response genes. Plant Cell 2017, 29, 791–807. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hamorsky, K.T.; Kouokam, J.C.; Jurkiewicz, J.M.; Nelson, B.; Moore, L.J.; Husk, A.S.; Kajiura, H.; Fujiyama, K.; Matoba, N. N-Glycosylation of cholera toxin B subunit in Nicotiana benthamiana: Impacts on host stress response, production yield and vaccine potential. Sci. Rep. 2015, 5, 8003. [Google Scholar] [CrossRef] [PubMed]
  58. Fujii, H.; Verslues, P.E.; Zhu, J.K. Arabidopsis decuple mutant reveals the importance of SnRK2 kinases in osmotic stress re-sponses in vivo. Proc. Natl. Acad. Sci. USA 2011, 108, 1717–1722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Heo, J.B.; Sung, S. Vernalization-mediated epigenetic silencing by a long intronic noncoding RNA. Science 2011, 331, 76–79. [Google Scholar] [CrossRef] [Green Version]
  60. Liu, F.; Marquardt, S.; Lister, C.; Swiezewski, S.; Dean, C. Targeted 3′ processing of antisense transcripts triggers Arabidopsis FLC chromatin silencing. Science 2010, 327, 94–97. [Google Scholar] [CrossRef]
  61. Sun, Q.; Csorba, T.; Skourti-Stathaki, K.; Proudfoot, N.J.; Dean, C. R-loop stabilization represses antisense transcription at the Arabidopsis FLC locus. Science 2013, 340, 619–621. [Google Scholar] [CrossRef] [Green Version]
  62. Johansson, M.; Staiger, D. SRR1 is essential to repress flowering in non-inductive conditions in Arabidopsis thaliana. J. Exp. Bot. 2014, 65, 5811–5822. [Google Scholar] [CrossRef] [Green Version]
  63. Endo, M.; Tanigawa, Y.; Murakami, T.; Araki, T.; Nagatani, A. Phytochrome-dependent late-flowering accelerates flowering through physical interactions with phytochrome B and CONSTANS. Proc. Natl. Acad. Sci. USA 2013, 110, 18017–18022. [Google Scholar] [CrossRef] [Green Version]
  64. Yang, S.W.; Jang, I.C.; Henriques, R.; Chua, N.H. FAR-RED ELONGATED HYPOCOTYL1 and FHY1-LIKE associate with the Arabidopsis transcription factors LAF1 and HFR1 to transmit phytochrome a signals for inhibition of hypocotyl elongation. Plant Cell 2009, 21, 1341–1359. [Google Scholar] [CrossRef] [Green Version]
  65. Kang, I.H.; Steffen, J.G.; Portereiko, M.F.; Lloyd, A.; Drews, G.N. The AGL62 MADS domain protein regulates cellularization during endosperm development in Arabidopsis. Plant Cell 2008, 20, 635–647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Gao, R.M.; Wang, Y.; Gruber, M.Y.; Hannoufa, A. miR156/SPL10 modulates lateral root development, branching and leaf morphology in Arabidopsis by silencing AGAMOUS-LIKE 79. Front. Plant Sci. 2017, 8, 2226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Zhou, R.G.; Kroczynska, B.; Miernyk, J.A. The genomic sequence encoding the Arabidopsis thaliana molecular chaperone AtJ3. Plant Physiol. 2000, 122, 291. [Google Scholar] [CrossRef] [Green Version]
  68. Zhu, Q.H.; Helliwell, C.A. Regulation of flowering time and floral patterning by miR172. J. Exp. Bot. 2011, 62, 487–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Meng, F.; Liu, H.; Wang, K.; Liu, L.; Wang, S.; Zhao, Y.; Yin, J.; Li, Y. Development-associated microRNAs in grains of wheat (Triticum aestivum L.). BMC Plant Biol. 2013, 13, 140. [Google Scholar] [CrossRef] [Green Version]
  70. Wu, Q.; Luo, Y.; Wu, X.; Bai, X.; Ye, X.L.; Liu, C.Y.; Wan, Y.; Xiang, D.B.; Li, Q.; Zou, L.; et al. Identification of the specific long-noncoding RNAs involved in night-break mediated flowering retardation in Chenopodium quinoa. BMC Genom. 2021, 22, 284. [Google Scholar] [CrossRef] [PubMed]
  71. Calixto, C.P.G.; Tzioutziou, N.A.; James, A.B.; Hornyik, C.; Guo, W.B.; Zhang, R.X.; Nimmo, H.G.; Brown, J.W.S. Cold-dependent expression and alternative splicing of Arabidopsis long non-coding RNAs. Front. Plant Sci. 2019, 10, 235. [Google Scholar] [CrossRef] [Green Version]
  72. Wang, P.; Dai, L.; Ai, J.; Wang, Y.; Ren, F. Identification and functional prediction of cold-related long non-coding RNA (lncRNA) in grapevine. Sci. Rep. 2019, 9, 1–15. [Google Scholar] [CrossRef]
  73. Bhatia, G.; Singh, A.; Verma, D.; Sharma, S.; Singh, K. Genome-wide investigation of regulatory roles of lncRNAs in response to heat and drought stress in Brassica juncea (Indian mustard). Environ. Exp. Bot. 2020, 171, 103922. [Google Scholar] [CrossRef]
  74. Wen, X.J.; Ding, Y.; Tan, Z.L.; Wang, J.X.; Zhang, D.Y.; Wang, Y.C. Identification and characterization of cadmium stress-related LncRNAs from Betula platyphylla. Plant Sci. 2020, 299, 110601. [Google Scholar] [CrossRef]
  75. Xu, J.H.; Fang, M.; Li, Z.H.; Zhang, M.N.; Liu, X.Y.; Peng, Y.Y.; Wan, Y.L.; Chen, J.H. Third-generation sequencing reveals lncRNA-regulated HSP genes in the Populus x Canadensis moench heat stress response. Front. Genet. 2020, 11, 249. [Google Scholar] [CrossRef] [PubMed]
  76. Chen, Q.F.; Xiao, S.; Chye, M.L. Overexpression of the Arabidopsis 10-kilodalton acyl-coenzyme A-binding protein ACBP6 enhances freezing tolerance. Plant Physiol. 2008, 148, 304–315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Lee, H.; Guo, Y.; Ohta, M.; Xiong, L.; Stevenson, B.; Zhu, J.K. LOS2, a genetic locus required for cold-responsive gene transcription encodes a bifunctional enolase. EMBO J. 2002, 21, 2692–2702. [Google Scholar] [CrossRef] [PubMed]
  78. Ismond, K.P.; Dolferus, R.; de Pauw, M.; Dennis, E.S.; Good, A.G. Enhanced low oxygen survival in Arabidopsis through increased metabolic flux in the fermentative pathway. Plant Physiol. 2003, 132, 1292–1302. [Google Scholar] [CrossRef] [Green Version]
  79. Zhang, F.; Maeder, M.L.; Unger-Wallace, E.; Hoshaw, J.P.; Reyon, D.; Christian, M.; Li, X.; Pierick, C.J.; Dobbs, D.; Peterson, T.; et al. High frequency targeted mutagenesis in Arabidopsis thaliana using zinc finger nucleases. Proc. Natl. Acad. Sci. USA 2010, 107, 12028–12033. [Google Scholar] [CrossRef] [Green Version]
  80. Sasaki, K.; Kim, M.; Imai, R. Arabidopsis cold shock domain protein 2 is a RNA chaperone that is regulated by cold and developmental signals. Biochem. Biophys. Res. Commun. 2007, 364, 633–638. [Google Scholar] [CrossRef]
  81. Anaraki, Z.E.; Tafreshi, S.A.H.; Shariati, M. Transient silencing of heat shock proteins showed remarkable roles for HSP70 during adaptation to stress in plants. Environ. Exp. Bot. 2018, 155, 142–157. [Google Scholar] [CrossRef]
  82. Huang, S.; Monaghan, J.; Zhong, X.H.; Lin, L.; Sun, T.J.; Dong, O.X.; Li, X. HSP90s are required for NLR immune receptor accumulation in Arabidopsis. Plant J. 2014, 79, 427–439. [Google Scholar] [CrossRef]
  83. Lin, W.; Zhou, Q.; Wang, C.Q.; Zhu, L.; Bi, C.; Zhang, S.; Wang, X.; Jin, H. LncRNAs regulate metabolism in cancer. Int. J. Biol. Sci. 2020, 16, 1194–1206. [Google Scholar] [CrossRef]
  84. Karimi, P.; Bakhtiarizadeh, M.R.; Salehi, A.; Izadnia, H.R. Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken. Sci. Rep. 2022, 12, 2558. [Google Scholar] [CrossRef]
  85. Tello-Flores, V.A.; Beltrán-Anaya, F.O.; Ramírez-Vargas, M.A.; Esteban-Casales, B.E.; Navarro-Tito, N.; Alarcón-Romero, L.d.C.; Luciano-Villa, C.A.; Ramírez, M.; del Moral-Hernández, Ó.; Flores-Alfaro, E. Role of long non-coding RNAs and the molecular mechanisms involved in insulin resistance. Int. J. Mol. Sci. 2021, 22, 7256. [Google Scholar] [CrossRef] [PubMed]
  86. Slabas, A.R.; Ndimba, B.; Simon, W.J.; Chivasa, S. Proteomic analysis of the Arabidopsis cell wall reveals unexpected proteins with new cellular locations. Biochem. Soc. Trans. 2004, 32, 524–528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Muñoz-Bertomeu, J.; Cascales-Miñana, B.; Mulet, J.M.; Baroja-Fernández, E.; Pozueta-Romero, J.; Kuhn, J.M.; Segura, J.; Ros, R. Plastidial glyceraldehyde-3-phosphate dehydrogenase deficiency leads to altered root development and affects the sugar and amino acid balance in Arabidopsis. Plant Physiol. 2009, 151, 541–558. [Google Scholar] [CrossRef] [Green Version]
  88. Zhang, L.; Song, Z.; Li, F.; Li, X.; Ji, H.; Yang, S. The specific MYB binding sites bound by TaMYB in the GAPCp2/3 promoters are involved in the drought stress response in wheat. BMC Plant Biol. 2019, 19, 366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Rius, S.P.; Casati, P.; Iglesias, A.A.; Gomez-Casati, D.F. Characterization of Arabidopsis lines deficient in GAPC-1, a cytosolic NAD-dependent glyceraldehyde-3-phosphate dehydrogenase. Plant Physiol. 2008, 148, 1655–1667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Egli, B.; Kolling, K.; Kohler, C.; Zeeman, S.C.; Streb, S. Loss of cytosolic phosphoglucomutase compromises gametophyte development in Arabidopsis. Plant Physiol. 2010, 154, 1659–1671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Fulton, D.C.; Stettler, M.; Mettler, T.; Vaughan, C.K.; Li, J.; Francisco, P.; Gil, M.; Reinhold, H.; Eicke, S.; Messerli, G.; et al. Beta-AMYLASE4, a noncatalytic protein required for starch breakdown, acts upstream of three active beta-amylases in Arabidopsis chloroplasts. Plant Cell 2008, 20, 1040–1058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Oono, Y.; Seki, M.; Satou, M.; Iida, K.; Akiyama, K.; Sakurai, T.; Fujita, M.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Monitoring expression profiles of Arabidopsis genes during cold acclimation and deacclimation using DNA microarrays. Funct. Integr. Genom. 2006, 6, 212–234. [Google Scholar] [CrossRef]
  93. Tanaka, H.; Nodzylski, T.; Kitakura, S.; Feraru, M.I.; Sasabe, M.; Ishikawa, T.; Kleine-Vehn, J.; Kakimoto, T.; Friml, J. BEX1/ARF1A1C is required for BFA-sensitive recycling of PIN auxin transporters and auxin-mediated development in Ara-bidopsis. Plant Cell Physiol. 2014, 55, 737–749. [Google Scholar] [CrossRef] [Green Version]
  94. Feng, S.H.; Shen, Y.P.; Sullivan, J.A.; Rubio, V.; Xiong, Y.; Sun, T.P.; Deng, X.W. Arabidopsis CAND1, an unmodified CUL1-interacting protein, is involved in multiple developmental pathways controlled by ubiquitin/proteasome-mediated protein degradation. Plant Cell 2004, 16, 1870–1882. [Google Scholar] [CrossRef] [Green Version]
  95. Tayengwa, R.; Zhao, J.; Pierce, C.F.; Werner, B.E.; Neff, M.M. Synopsis of the SOFL plant-specific gene family. G3 Genes Genomes Genet. 2018, 8, 1281–1290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Gampala, S.S.; Kim, T.W.; He, J.X.; Tang, W.; Deng, Z.; Bai, M.Y.; Guan, S.; Lalonde, S.; Sun, Y.; Gendron, J.M.; et al. An essential role for 14-3-3 proteins in brassinosteroid signal transduction in Arabidopsis. Dev. Cell 2007, 13, 177–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Yang, J.L.; Chen, W.W.; Chen, L.Q.; Qin, C.; Jin, C.W.; Shi, Y.Z.; Zheng, S.J. The 14-3-3 protein general regulatory factor11 (GRF 11) acts downstream of nitric oxide to regulate iron acquisition in Arabidopsis thaliana. New Phytol. 2013, 197, 815–824. [Google Scholar] [CrossRef] [PubMed]
  98. Ohta, M.; Matsui, K.; Hiratsu, K.; Shinshi, H.; Ohme-Takagi, M. Repression domains of class II ERF transcriptional repressors share an essential motif for active repression. Plant Cell 2001, 13, 1959–1968. [Google Scholar] [CrossRef] [Green Version]
  99. Jang, Y.H.; Park, H.Y.; Lee, K.C.; Thu, M.P.; Kim, S.K.; Suh, M.C.; Kang, H.; Kim, J.K. A homolog of splicing factor SF1 is essential for development and is involved in the alternative splicing of pre-mRNA in Arabidopsis thaliana. Plant J. 2014, 78, 591–603. [Google Scholar] [CrossRef]
  100. Shang, S.; Wu, C.; Huang, C.; Tie, W.; Yan, Y.; Ding, Z.; Xia, Z.; Wang, W.; Peng, M.; Tian, L.; et al. Genome-wide analysis of the GRF family reveals their involvement in abiotic stress response in cassava. Genes 2018, 9, 110. [Google Scholar] [CrossRef] [Green Version]
  101. Fujimoto, S.Y.; Ohta, M.; Usui, A.; Shinshi, H.; Ohme-Takagi, M. Arabidopsis ethylene-responsive element binding factors act as transcriptional activators or repressors of GCC box-mediated gene expression. Plant Cell 2000, 12, 393–404. [Google Scholar] [CrossRef] [Green Version]
  102. Lee, K.C.; Chung, K.S.; Lee, H.T.; Park, J.H.; Lee, J.H.; Kim, J.K. Role of Arabidopsis splicing factor SF1 in tempera-ture-responsive alternative splicing of FLM pre-mRNA. Front. Plant Sci. 2020, 11, 596354. [Google Scholar] [CrossRef]
  103. Johnson, X.; Alric, J. Central carbon metabolism and electron transport in Chlamydomonas reinhardtii: Metabolic constraints for carbon partitioning between oil and starch. Eukaryot Cell 2013, 12, 776–793. [Google Scholar] [CrossRef] [Green Version]
  104. Collakova, E.; Goyer, A.; Naponelli, V.; Krassovskaya, I.; Gregory, J.F.; Hanson, A.D.; Shachar-Hill, Y. Arabidopsis 10-formyl tetrahydrofolate deformylases are essential for photorespiration. Plant Cell 2008, 20, 1818–1832. [Google Scholar] [CrossRef] [Green Version]
  105. Christ, B.; Schelbert, S.; Aubry, S.; Suessenbacher, I.; Mueller, T.; Kraeutler, B.; Hoertensteiner, S. MES16, a member of the methylesterase protein family, specifically demethylates fluorescent chlorophyll catabolites during chlorophyll breakdown in Arabidopsis. Plant Physiol. 2012, 158, 628–641. [Google Scholar] [CrossRef] [Green Version]
  106. Baune, M.C.; Lansing, H.; Fischer, K.; Meyer, T.; Charton, L.; Linka, N.; von Schaewen, A. The Arabidopsis plastidial glucose-6-phosphate transporter GPT1 is dually targeted to peroxisomes via the endoplasmic reticulum. Plant Cell 2020, 32, 1703–1726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Sánchez-Fernández, R.; Davies, T.G.E.; Coleman, J.O.; Rea, P.A. The Arabidopsis thaliana ABC protein superfamily, a complete inventory. J. Biol. Chem. 2001, 276, 30231–30244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Skalitzky, C.A.; Martin, J.R.; Harwood, J.H.; Beirne, J.J.; Adamczyk, B.J.; Heck, G.R.; Cline, K.; Fernandez, D.E. Plastids contain a second sec translocase system with essential functions. Plant Physiol. 2011, 155, 354–369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. Kikuchi, S.; Asakura, Y.; Imai, M.; Nakahira, Y.; Kotani, Y.; Hashiguchi, Y.; Nakai, Y.; Takafuji, K.; Bédard, J.; Hirabayashi-Ishioka, Y.; et al. A Ycf2-ftshi heteromeric AAA-ATPase complex is required for chloroplast protein import. Plant Cell 2018, 30, 2677–2703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Reddy, V.S.; Day, I.S.; Thomas, T.; Reddy, A.S.N. KIC, a novel Ca2+ binding protein with one EF-hand motif, interacts with a microtubule motor protein and regulates trichome morphogenesis. Plant Cell 2004, 16, 185–200. [Google Scholar] [CrossRef] [Green Version]
  111. Bassil, E.; Ohto, M.A.; Esumi, T.; Tajima, H.; Zhu, Z.; Cagnac, O.; Belmonte, M.; Peleg, Z.; Yamaguchi, T.; Blumwald, E. The Arabidopsis intracellular Na+/H+ antiporters NHX5 and NHX6 are endosome associated and necessary for plant growth and development. Plant Cell 2011, 23, 224–239. [Google Scholar] [CrossRef] [Green Version]
Cimb 44 00128 g001 550
Figure 1. Distribution and classification of lncRNAs in Angelica sinensis during vernalization, based on the biological functions of coexpressed mRNAs. Abbreviations: CNCI, coding–noncoding index; CPC, coding potential calculator. Images (A), (B), (C), (D) and (E) represent total, aligned, identified, characterized and classified lncRNAs, respectively.
Figure 1. Distribution and classification of lncRNAs in Angelica sinensis during vernalization, based on the biological functions of coexpressed mRNAs. Abbreviations: CNCI, coding–noncoding index; CPC, coding potential calculator. Images (A), (B), (C), (D) and (E) represent total, aligned, identified, characterized and classified lncRNAs, respectively.
Cimb 44 00128 g001
Cimb 44 00128 g002 550
Figure 2. The expression levels of coexpressed mRNAs of lncRNAs linked with chromatin, DNA/RNA and protein modification in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. T1, T2, and T3 represent uncompleted, completed, and avoided vernalization, respectively. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Figure 2. The expression levels of coexpressed mRNAs of lncRNAs linked with chromatin, DNA/RNA and protein modification in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. T1, T2, and T3 represent uncompleted, completed, and avoided vernalization, respectively. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Cimb 44 00128 g002
Cimb 44 00128 g003 550
Figure 3. The expression levels of coexpressed mRNAs of lncRNAs directly linked with flowering in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Figure 3. The expression levels of coexpressed mRNAs of lncRNAs directly linked with flowering in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Cimb 44 00128 g003
Cimb 44 00128 g004 550
Figure 4. The expression levels of coexpressed mRNAs of lncRNAs linked with temperature response in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Figure 4. The expression levels of coexpressed mRNAs of lncRNAs linked with temperature response in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Cimb 44 00128 g004
Cimb 44 00128 g005 550
Figure 5. The expression levels of coexpressed mRNAs of lncRNAs linked with carbohydrate metabolism in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Figure 5. The expression levels of coexpressed mRNAs of lncRNAs linked with carbohydrate metabolism in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Cimb 44 00128 g005
Cimb 44 00128 g006 550
Figure 6. The expression levels of coexpressed mRNAs of lncRNAs linked with biosignaling in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Figure 6. The expression levels of coexpressed mRNAs of lncRNAs linked with biosignaling in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Cimb 44 00128 g006
Cimb 44 00128 g007 550
Figure 7. The expression levels of coexpressed mRNAs of lncRNAs linked with energy and transport in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Figure 7. The expression levels of coexpressed mRNAs of lncRNAs linked with energy and transport in A. sinensis at T2 vs. T1 and T3 vs. T1, as determined by qRT-PCR. The “*” represents a significant difference at p < 0.05 level between T2 vs. T1 and T3 vs. T2 for the same gene.
Cimb 44 00128 g007
Cimb 44 00128 g008 550
Figure 8. A proposed model of vernalization-induced bolting and flowering by regulatory lncRNAs in A. sinensis.
Figure 8. A proposed model of vernalization-induced bolting and flowering by regulatory lncRNAs in A. sinensis.
Cimb 44 00128 g008
Table
Table 1. Primer sequences used in qRT-PCR validation.
Table 1. Primer sequences used in qRT-PCR validation.
lncRNA IDCoexpressed mRNAsmRNA IDPrimer Sequences (5’ to 3’)Amplicon Size (bp)
Isoform0062250HMGB2NM_001035997.1Forward: CAAAGCTGCTGCTAAGGAC155
Reverse: GGACTTCCACTTGTCTCCAGC
Isoform0061796HMGB3NM_001035998.1Forward: CCTTCCAGTGCCTTCTTCGT174
Reverse: CTCAACCTTGCGCTTGTCAG
Isoform0001498At1g05910NM_100472.2Forward: AGACCACTCTCTCCGGTTGT109
Reverse: TCGTCAACTCCGATGACGTG
Isoform0062769RID2NM_125110.6Forward: GCAGGGCTTAGGTCTTCGTT100
Reverse: ACGAGGTTCATGCGATGACT
Isoform0061049At4g26600NM_001341820.1Forward: TTCCGATTGGTGCAACTCCT118
Reverse: GCCATGTCCACAACTCGTTC
Isoform0034756H2AVNM_001339683.1Forward: CAGTTGGACGAATTCACAGGC176
Reverse: CAGATGCCTTGGCGTTATCC
Isoform0050517At2g28720NM_128433.4Forward: GCAAGAAGCTTCAAAATTAGC107
Reverse: TGCTTAGCAAGTTCACCAGG
Isoform0062503HTR2NM_113651.2Forward: CACCGGAGGAGTGAAGAAGC189
Reverse: TCTTGAAGAGCTGCGACTGC
Isoform0027210HOS15NM_126132.4Forward: TACAGGCGCAGAACCTATGG164
Reverse: CTGTTGCATCACCAGACCCT
Isoform0062048REF6NM_148863.4Forward: AGGGAACACAGCTTCTGGTG124
Reverse: TTCCCCAAGTGAACGGTCTG
Isoform0062818SKP1ANM_106245.5Forward: GTGCTGCTACCTCCGATGAC181
Reverse: GTGCGGATCTCTTCTGGAGT
Isoform0061474ASK21NM_001125404.1Forward: CCTGATGACCTTACTGAGGAG178
Reverse: CAGGTCATCCACTGAACGCT
Isoform0061497SRK2GNM_120946.5Forward: ACATCGAGAGAGGTCGCAAG110
Reverse: AGGTGTCAGGATCACCTCCTT
Isoform0062575SRK2HNM_125760.2Forward: TGGTCCTGTGGTGTGACTCT164
Reverse: GAGAGAAGGTGTCTGCACTCC
Isoform0062220SRR1NM_125348.4Forward: ATCGCATTGTTTGGGAACAGC117
Reverse: AGCAAACTCGCTTGTGACTCT
Isoform0061284PHLNM_001334526.1Forward: CAAAGTCCTCGTTTGTCGGC104
Reverse: GCAACTGCTCCATAGTGGGT
Isoform0057927PHYANM_001123784.1Forward: GTGCGATATGCTGATGCGTG149
Reverse: CCTGCAGGTGGAACTCACTT
Isoform0041956AGL62NM_125437.5Forward: CTCCTCACCAACACAACAAC197
Reverse: AACGCAAGTTCCTCAACGGG
Isoform0045502AGL79NM_113925.3Forward: AATCACCCCATGAGCTTCGC107
Reverse: TAGGGTTCCGGCAGCTACTT
Isoform0063170ATJ3NM_114279.4Forward: GAATACGCTCACGGAGTTGC135
Reverse: GCATCCCACTTGGCTCTCTC
Isoform0062470ACBP6NM_102916.4Forward: AATCACCCCATGAGCTTCGC107
Reverse: TAGGGTTCCGGCAGCTACTT
Isoform0061783ENO2NM_129209.4Forward: CACTGAGTGTGGAACCGAGG190
Reverse: GGTCATCACTCCCCAACCTG
Isoform0063049ADH1NM_106362.3Forward: TGTGACCGAGTGTGTGAACC123
Reverse: TGAATCATGGCCTGAACGCT
Isoform0062198CSP2NM_120029.3Forward: GATCTGGAGGTGGATACGGC115
Reverse: CAGTCTCTCGCCATGTGACC
Isoform0062617HSP17.8NM_100614.3Forward: AACATCGGCGATAACGAACG154
Reverse: CTCCACGTGTCTCTCTCCAC
Isoform0009507HSP70-3NM_001202918.1Forward: CGACTGCAGGAGACACTCAT144
Reverse: TCTCACAGGCGGTTCTCAAC
Isoform0034676HSP70-10NM_120996.4Forward: CGTGTCCCCAAGGTTCAGTC167
Reverse: CCGAGCGATAGAGGTGTGAC
Isoform0042993HSP90-3NM_124983.4Forward: AACAAGGAGGAGTACGCTGC193
Reverse: AGACACGACGGACATAGAGC
Isoform0061974CSY4NM_001337082.1Forward: GATGCAGAGCTCTACCGACC197
Reverse: CCTCTTCCGGGTCAAGCAAT
Isoform0062375GAPCP2NM_101496.3Forward: CATTTCTGCACCTTCAGCGG198
Reverse: TTCTGAGTAGCTGTGGTCGC
Isoform0062268At3g52990NM_115159.5Forward: GACAACTTGCGACCAACTCG101
Reverse: AATCCACGAAGGGTCTCAGC
Isoform0062585PGM2NM_001160993.2Forward: TGAACTGCGTACCCAAGGAG158
Reverse: TCGGTCTGCATCACCATCAG
Isoform0062370BAM1NM_113297.3Forward: AACTCTCTCGCTGTTCCTCG165
Reverse: GGAGAAGCCCGTCTCACAAT
Isoform0062586ARF1NM_001337250.1Forward: GTGACCGTGTTGTTGAAGCC148
Reverse: TGAAGCCCAAGCTTGTCAGT
Isoform0061377CUL1NM_001036498.3Forward: GTGCCGTGCATTGCTAAGAG153
Reverse: TCTTCGGCCTGTTGGACAAG
Isoform0057235SOFL4NM_123240.2Forward: AGGTCGTGGATGAGGACTAC144
Reverse: GAACCGCTGATAATTTGGCCC
Isoform0043114SOFL5NM_001342234.1Forward: TGCGAGTCAGGATGGACTCT193
Reverse: TCCTTGGACCAGAAGAAGCAT
Isoform0062152GRF2NM_106479.3Forward: AACTCTCCGGAATCTGCGAC192
Reverse: GAGCAGATTTGTAAGCGGCG
Isoform0062828GRF11NM_001084180.2Forward: GGTGCTAGGAGAGCATCGTG198
Reverse: GACGGTGGATTCTCCCGAAG
Isoform0061395ERF3NM_103946.3Forward: ATCGTTTAGCGGACCCAGAC101
Reverse: CGCAATCGCTGTGACAATCC
Isoform0022533SF1NM_001036978.2Forward: GGCTTAGGGTCAACTCCGAC164
Reverse: CCAGTCACACGGTCCTTGAT
Isoform0044730PURU1NM_124115.4Forward: ACGTCTTCTACTCTCGCAGC125
Reverse: AGGCACACGCACAACTGAAT
Isoform0047216MES16NM_117770.5Forward: CCATCCCTTCTCCGCATCTT195
Reverse: TCATAGGAGCAGGACGCAAC
Isoform0063248GPT1NM_124861.5Forward: CGCTGGTTCGTTGATGATGC193
Reverse: AAACGCAGGTTCACCACTCT
Isoform0062571ABCF5NM_125882.3Forward: TGCTGATAGGCTTGTGGCTT103
Reverse: CGGCTCATCAAGTAGCAGCA
Isoform0008194SECA2NM_001198130.1Forward: ACTGTGAGGCCCATTGTCTG117
Reverse: CTCTGCCACGAAGCTGGTTA
Isoform0062373VPS26ANM_124733.4Forward: TGTTCCGCTTCCTCCAATCAA196
Reverse: TGCTCCAGTTGATTCTCGCC
Isoform0062449CML19NM_119864.5Forward: CGAGCTCAACGTTGCTATGAG160
Reverse: GTCTATGGAGTCTCGTTCTCCG
Isoform0061307NHX6NM_106609.4Forward: GGCATTTGCTCTTGCTCTGC112
Reverse: TCCTCCAATCAGCAACACCG
Table
Table 2. Twenty-nine lncRNAs linked with chromatin, DNA/RNA, and protein modification.
Table 2. Twenty-nine lncRNAs linked with chromatin, DNA/RNA, and protein modification.
lncRNA IDCoexpressed mRNAsmRNA IDProteins Encoded by Coexpressed mRNAs
Chromatin modification (2)
Isoform0062250HMGB2AT1G20693.1High mobility group B protein 2
Isoform0061796HMGB3AT1G20696.1High mobility group B protein 3
DNA/RNA modification (3)
Isoform0001498At1g05910AT1G05910.1ATPase family AAA domain-containing protein At1g05910
Isoform0062769RID2AT5G57280.118S rRNA (guanine-N(7))-methyltransferase RID2 e
Isoform0061049At4g26600AT4G26600.8S-adenosyl-L-methionine-dependent methyltransferases superfamily protein
Protein modification (24)
Isoform0034756H2AVAT3G54560.1Histone H2A variant 1
Isoform0050517At2g28720AT2G28720.1Histone H2B.3
Isoform0062503HTR2AT1G09200.1Histone H3.2
Isoform0027210HOS15AT5G67320.1WD40 repeat-containing protein HOS15
Isoform0062048REF6AT3G48430.1Lysine-specific demethylase REF6
Isoform0062818SKP1AAT1G75950.1SKP1-like protein 1A
Isoform0061474ASK21AT3G61415.1SKP1-like protein 21
Isoform0061497SRK2GAT5G08590.1Serine/threonine-protein kinase SRK2G
Isoform0062575SRK2HAT5G63650.1Serine/threonine-protein kinase SRK2H
Isoform0053138DET1AT4G10180.1Light-mediated development protein DET1
Isoform0030135BOPAt4g295601AT3G57130.2Ankyrin repeat family protein/BTB/POZ domain-containing protein
Isoform0063715At1g45180AT1G45180.1F27F5.26
Isoform0044146At3g50840AT3G50840.3Phototropic-responsive NPH3 family protein
Isoform0034163At2g36630AT2G36630.1Sulfite exporter TauE/SafE family protein 4
Isoform0058941At3g47890AT3G47890.2Ubiquitin carboxyl-terminal hydrolase-related protein
Isoform0007659UBP7A0A1I9LL79Ubiquitin carboxyl-terminal hydrolase
Isoform0045576DER2.1AT3G21280.2Derlin-2.1
Isoform0062510GRP3AT2G05520.1Glycine-rich protein 3
Isoform0062575MDH9.13AT5G42440.1At5g42440
Isoform0051826At3g24715AT3G24715.1Kinase superfamily with octicosapeptide/Phox/Bem1p domain-containing protein
Isoform0061548At3g16560AT3G16560.4Probable protein phosphatase 2C 40
Isoform0049189At3g62260AT3G62260.1Probable protein phosphatase 2C 49
Isoform0060667ESMD1AT2G01480.1Protein ESMERALDA 1
Isoform0062775At1g27930AT1G27930.1Probable methyltransferase At1g27930
Table
Table 3. Twelve lncRNAs directly linked with flowering.
Table 3. Twelve lncRNAs directly linked with flowering.
lncRNA IDCoexpressed mRNAsmRNA IDProteins Encoded by Coexpressed mRNAs
Isoform0062220SRR1AT5G59560.2Protein SENSITIVITY TO RED LIGHT REDUCED 1
Isoform0061284PHLAT1G72390.1Protein PHYTOCHROME-DEPENDENT LATE-FLOWERING
Isoform0057927PHYAAT1G09570.6Phytochrome A
Isoform0041956AGL62AT5G60440.1Agamous-like MADS-box protein AGL62
Isoform0045502AGL79AT3G30260.1AGAMOUS-like 79
Isoform0063170ATJ3AT3G44110.1Chaperone protein dnaJ 3
Isoform0028325BBX29AT5G54470.1At5g54470
Isoform0054894CLE13AT1G73965.1CLAVATA3/ESR (CLE)-related protein 13
Isoform0061298CLE44AT4G13195.1CLAVATA3/ESR (CLE)-related protein 44
Isoform0056751MXC17.10AT5G24710.1Transducin/WD40 repeatlike superfamily protein
Isoform0062524At1g06515AT1G06515.2Transmembrane protein, putative (DUF3317)
Isoform0061573BHLH30AT1G68810.1Transcription factor bHLH30
Table
Table 4. Fourteen lncRNAs directly linked with temperature response.
Table 4. Fourteen lncRNAs directly linked with temperature response.
lncRNA IDCoexpressed mRNAsmRNA IDProteins Encoded by Coexpressed mRNAs
Isoform0062470ACBP6AT1G31812.1Acyl-CoA-binding domain-containing protein 6
Isoform0061783ENO2AT2G36530.1Bifunctional enolase 2/transcriptional activator
Isoform0063049ADH1AT1G77120.1Alcohol dehydrogenase class-P
Isoform0062198CSP2AT4G38680.1Cold shock protein 2
Isoform0035932RH20AT1G55150.2DEA(D/H)-box RNA helicase family protein
Isoform0019851RH52AT3G58570.1DEAD-box ATP-dependent RNA helicase 52
Isoform0062484RH53AT3G22330.1DEAD-box ATP-dependent RNA helicase 53, mitochondrial
Isoform0058095RAB18AT5G66400.1Dehydrin Rab18
Isoform0061968XERO1AT3G50980.1Dehydrin Xero 1
Isoform0020919MED14AT3G04740.1Mediator of RNA polymerase II transcription subunit 14
Isoform0062617HSP17.8AT1G07400.117.8 kDa class I heat shock protein
Isoform0009507HSP70-3AT3G09440.2Heat shock 70 kDa protein 3
Isoform0034676HSP70-10AT5G09590.1Heat shock 70 kDa protein 10, mitochondrial
Isoform0042993HSP90-3AT5G56010.1Heat shock protein 90-3
Table
Table 5. Nineteen lncRNAs directly linked with carbohydrate metabolism.
Table 5. Nineteen lncRNAs directly linked with carbohydrate metabolism.
lncRNA IDCoexpressed mRNAsmRNA IDProteins Encoded by Coexpressed mRNAs
Isoform0061974CSY4AT2G44350.2Citrate synthase 4, mitochondrial
Isoform0062251FBA3AT2G01140.1Fructose-bisphosphate aldolase 3, chloroplastic
Isoform0062131GAPC1AT3G04120.1Glyceraldehyde-3-phosphate dehydrogenase GAPC1, cytosolic
Isoform0062375GAPCP2AT1G16300.1Glyceraldehyde-3-phosphate dehydrogenase GAPCP2, chloroplastic
Isoform0062268At3g52990AT3G52990.1Pyruvate kinase
Isoform0062585PGM2AT1G70730.3Phosphoglucomutase (alpha-D-glucose-1,6-bisphosphate-dependent)
Isoform0063034USPAT5G52560.1UDP-sugar pyrophosphorylase
Isoform0042435GLCNAC1PUT2AT2G35020.1UDP-N-acetylglucosamine diphosphorylase 2
Isoform0062018UXS2AT3G62830.1UDP-glucuronic acid decarboxylase 2
Isoform0062750XYLAAT5G57655.2Xylose isomerase
Isoform0060700GALS1AT2G33570.1Galactan beta-1,4-galactosyltransferase GALS1
Isoform0062032OFUT31AT4G24530.1O-fucosyltransferase 31
Isoform0039872At5g67460AT5G67460.1Glucan endo-1,3-beta-D-glucosidase
Isoform0045893RSW3AT5G63840.2Glycosyl hydrolases family 31 protein
Isoform0036200At1g59950AT1G59950.1Aldo/keto reductase
Isoform0042053At5g25970AT5G25970.2Core-2/I-branching beta-1,6-N-acetylglucosaminyl-transferase family protein
Isoform0059181UGT76E7AT5G38040.1UDP-glycosyltransferase 76E7
Isoform0037698At1g26850AT1G26850.2Probable methyltransferase PMT2
Isoform0062370BAM1AT3G23920.1Beta-amylase 1, chloroplastic
Table
Table 6. Thirteen lncRNAs directly linked with hormone signaling.
Table 6. Thirteen lncRNAs directly linked with hormone signaling.
lncRNA IDCoexpressed mRNAsmRNA IDProteins Encoded by Coexpressed mRNAs
Isoform0062586ARF1AT2G47170.1ADP-ribosylation factor 1
Isoform0061377CUL1AT4G02570.1Cullin-1
Isoform0015752T4L20.330AT4G34750.1SAUR-like auxin-responsive protein family
Isoform0057235SOFL4AT5G38790.1Protein SOB FIVE-LIKE 4
Isoform0043114SOFL5AT4G33800.2Protein SOB FIVE-LIKE 5
Isoform0062152GRF2AT1G78300.114-3-3-like protein GF14 omega
Isoform0062828GRF11AT1G34760.114-3-3-like protein GF14 omicron
Isoform0061395ERF3AT1G50640.1Ethylene-responsive transcription factor 3
Isoform0007735CIPK20AT5G45820.1CBL-interacting serine/threonine-protein kinase 20
Isoform0022533SF1AT5G51300.2Splicing factorlike protein 1
Isoform0046311AGD9AT5G46750.1Probable ADP-ribosylation factor GTPase-activating protein AGD9
Isoform0057859TIFY4BAT4G14720.1Protein TIFY 4B
Isoform0062437At2g34810AT2G34810.1Berberine bridge enzyme-like 16
Table
Table 7. Forty-three lncRNAs linked with energy and transport.
Table 7. Forty-three lncRNAs linked with energy and transport.
lncRNA IDCoexpressed mRNAsmRNA IDProteins Encoded by Coexpressed mRNAs
Energy (5)
Isoform0044730PURU1AT5G47435.2Formyltetrahydrofolate deformylase 1, mitochondrial
Isoform0031698ndhB1ATCG01250.1NAD(P)H-quinone oxidoreductase subunit 2 A, chloroplastic
Isoform0029144NDB1AT4G28220.2NADH:ubiquinone reductase (nonelectrogenic)
Isoform0046995WNK9AT5G28080.3Nonspecific serine/threonine protein kinase
Isoform0047216MES16AT4G16690.1pFDCC methylesterase MES16
Transport (38)
Isoform0047603At5g11230AT5G11230.1Probable sugar phosphate/phosphate translocator At5g11230
Isoform0063248GPT1AT5G54800.1Glucose-6-phosphate/phosphate translocator 1, chloroplastic
Isoform0053163SMXL5AT5G57130.1Protein SMAX1-LIKE 5
Isoform0062571ABCF5AT5G64840.1ABC transporter F family member 5
Isoform0062182SFH8AT2G21520.2Phosphatidylinositol/phosphatidylcholine transfer protein SFH8
Isoform0022607BASS6AT4G22840.1Probable sodium/metabolite cotransporter BASS6, chloroplastic
Isoform0008194SECA2AT1G21650.3Protein translocase subunit SECA2, chloroplastic
Isoform0014041At4g14160AT4G14160.1Protein transport protein SEC23
Isoform0062660ycf2-AATCG01280.1Protein Ycf2
Isoform0036773At4g22990AT4G22990.1SPX domain-containing membrane protein At4g22990
Isoform0032676At3g49350AT3G49350.1At3g49350
Isoform0062443TMT2AT4G35300.5Tonoplast monosaccharide transporter2
Isoform0003274VPS24-1AT5G22950.1Vacuolar protein sorting-associated protein 24 homolog 1
Isoform0062373VPS26AAT5G53530.1Vacuolar protein sorting-associated protein 26A
Isoform0028575VPS52AT1G71270.1Vacuolar protein sorting-associated protein 52 A
Isoform0038654VPS60-2AT5G04850.1Vacuolar protein sorting-associated protein 60.2
Isoform0062996At5g19500AT5G19500.1At5g19500
Isoform0062800PIP1-5AT4G23400.1Probable aquaporin PIP1-5
Isoform0061992SULTR2;2AT1G77990.1Sulfate transporter 2.2
Isoform0062606MT2AAT3G09390.2Metallothionein-like protein 2A
Isoform0057136At3g52300AT3G52300.1ATP synthase subunit d, mitochondrial
Isoform0061450ABCB23AT4G28630.1ABC transporter B family member 23, mitochondrial
Isoform0062449CML19AT4G37010.1Calcium-binding protein CML19
Isoform0061763JJJ1AT1G74250.1DNAJ protein JJJ1 homolog
Isoform0061955PAP1AT1G13750.1Probable inactive purple acid phosphatase 1
Isoform0062074HIPP04AT1G29000.2Heavy metal-associated isoprenylated plant protein 4
Isoform0062972HIPP26AT4G38580.1Heavy metal-associated isoprenylated plant protein 26
Isoform0014712KINUAAT1G12430.1Kinesin-like protein KIN-UA
Isoform0063051IQD30AT1G18840.2Protein IQ-DOMAIN 30
Isoform0016777At3g18430AT3G18430.1Calcineurin b subunit (Protein phosphatase 2b regulatory subunit)-like protein
Isoform0063921CBL10AT4G33000.6Calcineurin B-like protein 10
Isoform0061307NHX6AT1G79610.1Sodium/hydrogen exchanger 6
Isoform0024682CNGC15AT2G28260.2Cyclic nucleotide-gated channel 15
Isoform00346382-OctAT1G79360.1Organic cation/carnitine transporter 2
Isoform0033810CNGC13AT4G01010.2Putative cyclic nucleotide-gated ion channel 13
Isoform0051075MTG13.4AT5G16680.1RING/FYVE/PHD zinc finger superfamily protein
Isoform0051563SUF4AT1G30970.3Zinc finger (C2H2 type) family protein
Isoform0037483AHA4AT3G47950.2ATPase 4, plasma membranetype
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, X.; Luo, M.; Li, M.; Wei, J. Transcriptomic Analysis Reveals LncRNAs Associated with Flowering of Angelica sinensis during Vernalization. Curr. Issues Mol. Biol. 2022, 44, 1867-1888. https://doi.org/10.3390/cimb44050128

AMA Style

Liu X, Luo M, Li M, Wei J. Transcriptomic Analysis Reveals LncRNAs Associated with Flowering of Angelica sinensis during Vernalization. Current Issues in Molecular Biology. 2022; 44(5):1867-1888. https://doi.org/10.3390/cimb44050128

Chicago/Turabian Style

Liu, Xiaoxia, Mimi Luo, Mengfei Li, and Jianhe Wei. 2022. "Transcriptomic Analysis Reveals LncRNAs Associated with Flowering of Angelica sinensis during Vernalization" Current Issues in Molecular Biology 44, no. 5: 1867-1888. https://doi.org/10.3390/cimb44050128

Article Metrics

Back to TopTop