Genome-wide identification and characterization of novel lncRNAs in Populus under nitrogen deficiency

Genome-wide analysis of long non-coding RNAs responsive to multiple nutrient stresses in Arabidopsis thaliana

Nutrient stress is the most important environmental stress that limits plant growth and development. Although recent evidence highlights the vital functions of long non-coding RNAs (lncRNA) in response to single nutrient stress in some model plants, a comprehensive investigation of the effect of lncRNAs in response to nutrient stress has not been performed in Arabidopsis thaliana. Here, we presented the identification and characterization of lncRNAs under seven nutrient stress conditions. The expression pattern analysis revealed that aberrant expression of lncRNAs is a stress-specific manner under nutrient stress conditions and that lncRNAs are more sensitive to nutrient stress than protein-coding genes (PCGs). Moreover, competing endogenous RNA (ceRNA) network and lncRNA-mRNA co-expression network (CEN) were constructed to explore the potential function of these lncRNAs under nutrient stress conditions. We further combined different expressed lncRNAs with ceRNA network and CEN to select key lncRNAs in response to nutrient stress. Together, our study provides important information for further insights into the role of lncRNAs in response to stress in plants.

Long non-coding RNAs: emerging players regulating plant abiotic stress response and adaptation

BACKGROUND

The immobile nature of plants means that they can be frequently confronted by various biotic and abiotic stresses during their lifecycle. Among the various abiotic stresses, water stress, temperature extremities, salinity, and heavy metal toxicity are the major abiotic stresses challenging overall plant growth. Plants have evolved complex molecular mechanisms to adapt under the given abiotic stresses. Long non-coding RNAs (lncRNAs)-a diverse class of RNAs that contain > 200 nucleotides(nt)-play an essential role in plant adaptation to various abiotic stresses.

RESULTS

LncRNAs play a significant role as 'biological regulators' for various developmental processes and biotic and abiotic stress responses in animals and plants at the transcription, post-transcription, and epigenetic level, targeting various stress-responsive mRNAs, regulatory gene(s) encoding transcription factors, and numerous microRNAs (miRNAs) that regulate the expression of different genes. However, the mechanistic role of lncRNAs at the molecular level, and possible target gene(s) contributing to plant abiotic stress response and adaptation, remain largely unknown. Here, we review various types of lncRNAs found in different plant species, with a focus on understanding the complex molecular mechanisms that contribute to abiotic stress tolerance in plants. We start by discussing the biogenesis, type and function, phylogenetic relationships, and sequence conservation of lncRNAs. Next, we review the role of lncRNAs controlling various abiotic stresses, including drought, heat, cold, heavy metal toxicity, and nutrient deficiency, with relevant examples from various plant species. Lastly, we briefly discuss the various lncRNA databases and the role of bioinformatics for predicting the structural and functional annotation of novel lncRNAs.

CONCLUSIONS

Understanding the intricate molecular mechanisms of stress-responsive lncRNAs is in its infancy. The availability of a comprehensive atlas of lncRNAs across whole genomes in crop plants, coupled with a comprehensive understanding of the complex molecular mechanisms that regulate various abiotic stress responses, will enable us to use lncRNAs as potential biomarkers for tailoring abiotic stress-tolerant plants in the future.

Third-generation sequencing found LncRNA associated with heat shock protein response to heat stress in Populus qiongdaoensis seedlings

Background

As air temperatures increase globally, more and more plants are exposed to heat-stress conditions. Although many studies have explored regulation networks in plants with the aim of improving their heat-stress tolerance, only few have revealed them in trees. Here, individuals of Populus qiongdaoensis seedlings, which grows naturally in tropical areas, exposed to heat at 40 °C and the non-coding regulation networks were explored using the PacBio RSII and the Illumina sequencing platform.

Results

In total, we obtained 88,161 full-length transcripts representing 39,343 genes using 5,498,988 long reads and 350,026,252 clean reads, and also 216 microRNAs (miRNAs) via 95,794,107 reads. We then identified 928 putative long non-coding RNAs (lncRNAs), consisting of 828 sense lncRNAs (89.22%), 34 long intergenic non-coding RNAs (3.66%), 16 antisense (1.72%), and 50 sense intronic lncRNAs (5.39%). Under the dual criteria of |log₂fold-change| ≥ 1 and P-value < 0.05, 1690 genes, 25 lncRNAs, and 15 miRNAs were found differentially expressed under the heat stress treatment. Furthermore, 563 and 595 mRNAs were detected as target genes of 14 differently expressed miRNAs and 26 differentially expressed lncRNAs. Functional annotation analysis of these target genes demonstrated they were related to cell membrane stability, plant hormone signal transduction, antioxidation, and aldarate metabolism. Lastly, we uncovered a key interaction network of lncRNAs, miRNAs and mRNAs that consisted of miR1444d, miR482a.1, miR530a, lncHSP18.2, HSP18.1, and HSP18.2. Expression level analysis showed that miRNAs in the network were up-regulated, while mRNAs and lncRNA were down-regulated, and also found that lncHSP18.2 may cis-regulate HSP18.2.

Conclusions

Functional enrichment analysis of target genes of miRNAs and lncRNAs indicated that miRNAs and lncRNAs play an important role in the response to heat stress P. qiongdaoensis. Lastly, by investigating the miRNA–lncRNA–mRNA network of this species, we revealed that miRNAs may negatively regulate both lncRNAs and mRNAs in tree responses to heat stress, and found that lncHSP18.2 may cis-regulate HSP18.2.

A long non-coding apple RNA, MSTRG.85814.11, acts as a transcriptional enhancer of SAUR32 and contributes to the Fe-deficiency response

Iron (Fe) is an essential plant nutrient and its deficiency typically limits plant growth. Long non-coding (lnc) RNAs are involved in adaptive responses to nutrient stress; however, it is not known whether they function in the regulation of the canonical Fe-deficiency response. The expression of Malus domestica (apple) lncRNA MSTRG.85814 is induced by Fe deficiency, as identified by high-throughput strand-specific RNA-seq analysis of an apple homograft system. MSTRG.85814 has a complex structure, with 13 predicted RNA sequence variants, four of which are upregulated in the roots of plants experiencing Fe deficiency. We found that one MSTRG.85814 splice variant (MSTRG.85814.11) positively modulated its cis target mRNA derived from the small auxin upregulated gene SAUR32. This in turn promoted the expression of SAUR32 and caused an increase in the expression of a plasma membrane proton ATPase, AHA10. Using a pH imaging technique, a significant decrease in the apoplastic pH was observed to occur in the root tips of MSTRG.85814.11 or SAUR32-overexpressing apple plants. Thus MSTRG.85814.11 was shown to positively promote SAUR32 expression, which then activated proton extrusion involved in the Fe-deficiency response. These results reveal a mechanism by which lncRNA promotes environmental Fe-deficiency stress adaption.

Third-Generation Sequencing Reveals LncRNA-Regulated HSP Genes in the Populus x canadensis Moench Heat Stress Response

Long non-coding RNAs (lncRNAs) regulate plant responses to abiotic stresses. However, the short reads produced by second-generation sequencing technology make it difficult to accurately explore full-length transcripts, limiting the study of lncRNAs. In this study, we used third-generation long-read sequencing technology with the PacBio Sequel and Illumina platform to explore the role of lncRNAs in the heat stress response of Populus x canadensis Moench trees. We using 382,034,416 short reads to correct 4,297,179 long reads by resulted in 66,657 full-length transcripts, representing 33,840 genes. Then, 753 putative lncRNAs were identified, including 658 sense lncRNAs (87.38%), 41 long intervening/intergenic non-coding RNAs (lincRNAs) (5.44%), 12 antisense lncRNAs (1.59%), and 42 sense intronic lncRNAs (5.58%). Using the criteria | log₂FC| ≥ 1 and q-value < 0.05, 3,493 genes and 78 lncRNAs were differentially expressed under the heat treatment. Furthermore, 923 genes were detected as targets of 43 differently expressed lncRNAs by cis regulation. Functional annotation demonstrated that these target genes were related to unfolded protein binding, response to stress, protein folding, and response to stimulus. Lastly, we identified a lncRNA–gene interaction network consisting of four lncRNAs and six genes [Heat Shock Protein 82 (HSP82), HSP83, Disease Resistance Protein 27 (DRL27), DnaJ family protein (DNJH), and two other predicted protein-coding genes], which showed that lncRNAs could regulate HSP family genes in response to heat stress in Populus. Therefore, our third-generation sequencing has improved the description of the P. canadensis transcriptome. The potential lncRNAs and HSP family genes identified here present a genetic resource to improve our understanding of the heat-adaptation mechanisms of trees.

A Comprehensive Transcriptomics Analysis Reveals Long Non-Coding RNA to be Involved in the Key Metabolic Pathway in Response to Waterlogging Stress in Maize

Waterlogging stress (WS) in a dynamic environment seriously limits plant growth, development, and yield. The regulatory mechanism underlying WS conditions at an early stage in maize seedlings is largely unknown. In the present study, the primary root tips of B73 seedlings were sampled before (0 h) and after (2 h, 4 h, 6 h, 8 h, 10 h, and 12 h) WS and then subjected to transcriptome sequencing, resulting in the identification of differentially expressed protein-coding genes (DEpcGs) and long non-coding RNAs (DElncRs) in response to WS. These DEpcGs were classified into nine clusters, which were significantly enriched in several metabolic pathways, such as glycolysis and methionine metabolism. Several transcription factor families, including AP2-EREBP, bZIP, NAC, bHLH, and MYB, were also significantly enriched. In total, 6099 lncRNAs were identified, of which 3190 were DElncRs. A co-expression analysis revealed lncRNAs to be involved in 11 transcription modules, 10 of which were significantly associated with WS. The DEpcGs in the four modules were enriched in the hypoxia response pathways, including phenylpropanoid biosynthesis, MAPK signaling, and carotenoid biosynthesis, in which 137 DElncRs were also co-expressed. Most of the co-expressed DElncRs were co-localized with previously identified quantitative trait loci associated with waterlogging tolerance. A quantitative reverse transcription-polymerase chain reaction analysis of DEpcG and DElncR expression among the 32 maize genotypes after 4 h of WS verified significant expression correlations between them as well as significant correlation with the phenotype of waterlogging tolerance. Moreover, the high proportion of hypoxia response elements in the promoter region increased the reliability of the DElncRs identified in this study. These results provide a comprehensive transcriptome in response to WS at an early stage of maize seedlings and expand our understanding of the regulatory network involved in hypoxia in plants.

Identification of lncRNAs and their regulatory relationships with target genes and corresponding miRNAs in melon response to powdery mildew fungi

Melon (Cucumis melo L.), an economically beneficial crop widely cultivated around the world, is vulnerable to powdery mildew (PM). However, the studies on molecular mechanism of melon response to PM fungi is still limited. Long non coding RNAs (lncRNAs) have emerged as new regulators in plants response to biotic stresses. We predicted and identified the intricate regulatory roles of lncRNAs in melon response to PM fungi. A total of 539 lncRNAs were identified from PM-resistant (MR-1) and susceptible melon (Top Mark), in which 254 were significantly altered after PM fungi infection. Multiple target genes of lncRNAs were found to be involved in the hydrolysis of chitin, callose deposition and cell wall thickening, plant-pathogen interaction and plant hormone signal transduction pathway. Additionally, a total of 42 lncRNAs possess the various functions with microRNAs (miRNAs), including lncRNAs that are targeted by miRNAs and function as miRNA precursors or miRNA sponges. These findings provide a comprehensive view of potentially functional lncRNAs, corresponding target genes and related lncRNA-miRNA pairs, which will greatly increase our knowledge of the mechanism underlying susceptibility and resistance to PM in melon.

Novel Aspects of Nitrate Regulation in Arabidopsis

Nitrogen (N) is one of the most essential macronutrients for plant growth and development. Nitrate (NO₃ ^-), the major form of N that plants uptake from the soil, acts as an important signaling molecule in addition to its nutritional function. Over the past decade, significant progress has been made in identifying new components involved in NO₃ ^- regulation and starting to unravel the NO₃ ^- regulatory network. Great reviews have been made recently by scientists on the key regulators in NO₃ ^- signaling, NO₃ ^- effects on plant development, and its crosstalk with phosphorus (P), potassium (K), hormones, and calcium signaling. However, several novel aspects of NO₃ ^- regulation have not been previously reviewed in detail. Here, we mainly focused on the recent advances of post-transcriptional regulation and non-coding RNA (ncRNAs) in NO₃ ^- signaling, and NO₃ ^- regulation on leaf senescence and the circadian clock. It will help us to extend the general picture of NO₃ ^- regulation and provide a basis for further exploration of NO₃ ^- regulatory network.

Transcriptional profiling of long noncoding RNAs associated with leaf-color mutation in Ginkgo biloba L

BACKGROUND

Long noncoding RNAs (lncRNAs) play an important role in diverse biological processes and have been widely studied in recent years. However, the roles of lncRNAs in leaf pigment formation in ginkgo (Ginkgo biloba L.) remain poorly understood.

RESULTS

In this study, lncRNA libraries for mutant yellow-leaf and normal green-leaf ginkgo trees were constructed via high-throughput sequencing. A total of 2044 lncRNAs were obtained with an average length of 702 nt and typically harbored 2 exons. We identified 238 differentially expressed lncRNAs (DELs), 32 DELs and 49 differentially expressed mRNAs (DEGs) that constituted coexpression networks. We also found that 48 cis-acting DELs regulated 72 target genes, and 31 trans-acting DELs regulated 31 different target genes, which provides a new perspective for the regulation of the leaf-color mutation. Due to the crucial regulatory roles of lncRNAs in a wide range of biological processes, we conducted in-depth studies on the DELs and their targets and found that the chloroplast thylakoid membrane subcategory and the photosynthesis pathways (ko00195) were most enriched, suggesting their potential roles in leaf coloration mechanisms. In addition, our correlation analysis indicates that eight DELs and 68 transcription factors (TFs) might be involved in interaction networks.

CONCLUSIONS

This study has enriched the knowledge concerning lncRNAs and provides new insights into the function of lncRNAs in leaf-color mutations, which will benefit future selective breeding of ginkgo.

Integrative Analysis of the lncRNA and mRNA Transcriptome Revealed Genes and Pathways Potentially Involved in the Anther Abortion of Cotton (Gossypium hirsutum L.)

Cotton plays an important role in the economy of many countries. Many studies have revealed that numerous genes and various metabolic pathways are involved in anther development. In this research, we studied the differently expressed mRNA and lncRNA during the anther development of cotton between the cytoplasmic male sterility (CMS) line, C2P5A, and the maintainer line, C2P5B, using RNA-seq analysis. We identified 17,897 known differentially expressed (DE) mRNAs, and 865 DE long noncoding RNAs (lncRNAs) that corresponded to 1172 cis-target genes at three stages of anther development using gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of DE mRNAs; and cis-target genes of DE lncRNAs probably involved in the degradation of tapetum cells, microspore development, pollen development, and in the differentiation, proliferation, and apoptosis of the anther cell wall in cotton. Of these DE genes, LTCONS_00105434, LTCONS_00004262, LTCONS_00126105, LTCONS_00085561, and LTCONS_00085561, correspond to cis-target genes Ghir_A09G011050.1, Ghir_A01G005150.1, Ghir_D05G003710.2, Ghir_A03G016640.1, and Ghir_A12G005100.1, respectively. They participate in oxidative phosphorylation, flavonoid biosynthesis, pentose and glucuronate interconversions, fatty acid biosynthesis, and MAPK signaling pathway in plants, respectively. In summary, the transcriptomic data indicated that DE lncRNAs and DE mRNAs were related to the anther development of cotton at the pollen mother cell stage, tetrad stage, and microspore stage, and abnormal expression could lead to anther abortion, resulting in male sterility of cotton.

Genome-Wide Identification of Long Non-coding RNA in Trifoliate Orange (Poncirus trifoliata (L.) Raf) Leaves in Response to Boron Deficiency

Long non-coding RNAs (lncRNAs) play important roles in plant growth and stress responses. As a dominant abiotic stress factor in soil, boron (B) deficiency stress has impacted the growth and development of citrus in the red soil region of southern China. In the present work, we performed a genome-wide identification and characterization of lncRNAs in response to B deficiency stress in the leaves of trifoliate orange (Poncirus trifoliata), an important rootstock of citrus. A total of 2101 unique lncRNAs and 24,534 mRNAs were predicted. Quantitative real-time polymerase chain reaction (qRT-PCR) experiments were performed for a total of 16 random mRNAs and lncRNAs to validate their existence and expression patterns. Expression profiling of the leaves of trifoliate orange under B deficiency stress identified 729 up-regulated and 721 down-regulated lncRNAs, and 8419 up-regulated and 8395 down-regulated mRNAs. Further analysis showed that a total of 84 differentially expressed lncRNAs (DELs) were up-regulated and 31 were down-regulated, where the number of up-regulated DELs was 2.71-fold that of down-regulated. A similar trend was also observed in differentially expressed mRNAs (DEMs, 4.21-fold). Functional annotation of these DEMs was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, and the results demonstrated an enrichment of the categories of the biosynthesis of secondary metabolites (including phenylpropanoid biosynthesis/lignin biosynthesis), plant hormone signal transduction and the calcium signaling pathway. LncRNA target gene enrichment identified several target genes that were involved in plant hormones, and the expression of lncRNAs and their target genes was significantly influenced. Therefore, our results suggest that lncRNAs can regulate the metabolism and signal transduction of plant hormones, which play an important role in the responses of citrus plants to B deficiency stress. Co-expression network analysis indicated that 468 significantly differentially expressed genes may be potential targets of 90 lncRNAs, and a total of 838 matched lncRNA-mRNA pairs were identified. In summary, our data provides a rich resource of candidate lncRNAs and mRNAs, as well as their related pathways, thereby improving our understanding of the role of lncRNAs in response to B deficiency stress, and in symptom formation caused by B deficiency in the leaves of trifoliate orange.

Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression of Heteromorphic Leavesin Response to Environmental Heterogeneityin Populus euphratica

Heterophylly is the phenomenon of leaf forms varying along the longitudinal axis within a single plant. Populus euphratica, a heterophyllous woody plant, develops lanceolate leaves and dentate broad-ovate leaves on the bottom and top of the canopy, respectively, which are faced with different intensities of ambient solar radiation. However, the mechanism of the heteromorphic leaf response to the microenvironment in P. euphratica remains elusive. Here, we show that the dentate broad-ovate leaves have advantages in tolerating high light intensity, while lanceolate leaves are excellent at capturing light. Compared with lanceolate leaves, more trichomes, higher stomatal density, thicker lamina, and higher specific leaf weight were observed in dentate broad-ovate leaves. Furthermore, high-throughput RNA sequencing analysis revealed that the expression patterns of genes and long noncoding RNAs (lncRNAs) are different between the two heteromorphic leaves. A total of 36,492 genes and 1725 lncRNAs were detected, among which 586 genes and 54 lncRNAs were differentially expressed. Based on targets prediction, lncRNAs and target genes involved in light adaption, protein repair, stress response, and growth and development pathways were differentially expressed in heteromorphic leaves, 10 pairs of which were confirmed by quantitative real-time PCR. Additionally, the analysis of interactions indicated that lncRNA-mRNA interactions were involved in the response to the microenvironment of heteromorphic leaves. Taken together, these results suggest that the morphological features and joint regulation of lncRNA-mRNA in heteromorphic leaves may serve as survival strategies for P. euphratica, which could lead to optimal utilization of environmental factors.

Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots

BACKGROUND

Environmental stimuli can activate a series of physiological and biochemical responses in plants accompanied by extensive transcriptional reprogramming. Long non-coding RNAs (lncRNAs), as versatile regulators, control gene expression in multiple ways and participate in the adaptation to biotic and abiotic stresses.

RESULTS

In this study, soybean seedlings were continuously cultured for 15 days with high salinity solutions started from seed germination. Strand-specific whole transcriptome sequencing and stringent bioinformatic analysis led to the identification of 3030 long intergenic non-coding RNAs (lincRNAs) and 275 natural antisense transcripts (lncNATs) in soybean roots. In contrast to mRNAs, newly identified lncRNAs exhibited less exons, similar AU content to UTRs, even distribution across the genome and low evolutionary conservation. Remarkably, more than 75% of discovered lncRNAs that were activated or up-regulated by continuous salt stress mainly targeted proteins with binding and catalytic activities. Furthermore, two DNA methylation maps with single-base resolution were generated by using reduced representation bisulfite sequencing, offering a genome-wide perspective and important clues for epigenetic regulation of stress-associated lncRNAs and protein-coding genes.

CONCLUSIONS

Taken together, our findings systematically demonstrated the characteristics of continuous salt stress-induced lncRNAs and extended the knowledge of corresponding methylation profiling, providing valuable evidence for a better understanding of how plants cope with long-term salt stress circumstances.

The long noncoding RNA T5120 regulates nitrate response and assimilation in Arabidopsis

Long noncoding RNAs (lncRNAs) are crucial regulators in many plant biological processes. However, it remains unknown whether lncRNAs can respond to nitrate or function in nitrate regulation. We detected 695 lncRNAs, 480 known and 215 novel, in Arabidopsis seedling roots; six showed altered expression in response to nitrate treatment, among which T5120 showed the highest induction. Overexpression of T5120 in Arabidopsis promoted the response to nitrate, enhanced nitrate assimilation and improved biomass and root development. Biochemical and molecular analyses revealed that NLP7, a master nitrate regulatory transcription factor, directly bound to the nitrate-responsive cis-element (NRE)-like motif of the T5120 promoter and activated T5120 transcription. In addition, T5120 partially restored the nitrate signalling and assimilation phenotypes of nlp7 mutant, suggesting that T5120 is involved in NLP7-mediated nitrate regulation. Interestingly, the expression of T5120 was regulated by the nitrate sensor NRT1.1. Therefore, T5120 is modulated by NLP7 and NRT1.1 to regulate nitrate signalling. Our work reveals a new regulatory mechanism in which lncRNA T5120 functions in nitrate regulation, providing new insights into the nitrate signalling network. Importantly, lncRNA T5120 can promote nitrate assimilation and plant growth to improve nitrogen use efficiency.

Metabolic and physiological analyses reveal that Populus cathayana males adopt an energy-saving strategy to cope with phosphorus deficiency

Dioecious trees have evolved sex-specific adaptation strategies to cope with inorganic phosphorus (Pi) limitation. Yet, little is known about the effects of Pi limitation on plant metabolism, particularly in dioecious woody plants. To identify potential gender-specific metabolites appearing in response to Pi limitation in poplars, we studied the metabolic and ionomic responses in the roots and leaves of Populus cathayana Rehd males and females exposed to a 60-day period of Pi deficiency. Besides significant decreases in phosphorus contents in both Pi-deficient roots and leaves, the calcium level decreased significantly and the sulfur content increased significantly in Pi-deficient male roots, while the zinc and ferrum contents increased significantly in Pi-deficient female roots. Inorganic P deficiency caused a smaller change in the abscisic acid content, but a significant increase in the jasmonic acid content was detected in both leaves and roots. Salicylic acid significantly decreased under Pi deficiency in male leaves and female roots. Changes were found in phospholipids and phosphorylated metabolites (e.g., fructose-6-phosphate, glycerol-3-phosphate, glucose-6-phosphate, phosphoric acid and inositol-1-phosphate) in roots and leaves. Both P. cathayana males and females relied on inorganic pyrophosphate-dependent but not on Pi-dependent glycolysis under Pi-deficient conditions. Sex-specific metabolites in leaves were primarily in the category of primary metabolites (e.g., amino acids), while in roots primarily in the category of secondary metabolites (e.g., organic acids) and sugars. The metabolome analysis revealed that sexually different pathways occurred mainly in amino acid metabolism, and the tissue-related differences were in the shikimate pathway and glycolysis. We observed changes in carbon flow, reduced root biomass and increased amino acid contents in P. cathayana males but not in females, which indicated that males have adopted an energy-saving strategy to adapt to Pi deficiency. Thus, this study provides new insights into sex-specific metabolic responses to Pi deficiency.

Long noncoding RNAs in Brassica rapa L. following vernalization

Brassica rapa L. is an important agricultural crop that requires a period of prolonged cold for flowering. This process is known as vernalization. Studies have shown that long noncoding RNAs (lncRNAs) play important roles in abiotic stress responses and several cold-responsive noncoding RNAs have been suggested to be involved in vernalization. We examined the transcriptome of the Chinese cabbage inbred line (B. rapa L. var. pekinensis) RJKB-T24, and identified 1,444 long intergenic noncoding RNAs (lincRNAs), 551 natural antisense transcripts (NATs), and 93 intronic noncoding RNAs (incRNAs); 549 of the 2,088 lncRNAs significantly altered their expression in response to four weeks of cold treatment. Most differentially expressed lncRNAs did not lead to a change of expression levels in mRNAs covering or near lncRNAs, suggesting that the transcriptional responses to four weeks of cold treatment in lncRNA and mRNA are independent. However, some differentially expressed mRNAs had NATs with expression altered in the same direction. These genes were categorized as having an abiotic stress response, suggesting that the paired-expression may play a role in the transcriptional response to vernalization or cold treatment. We also identified short-term cold treatment induced NATs in BrFLC and BrMAF genes, which are involved in vernalization. The lncRNAs we identified differed from those reported in Arabidopsis thaliana, suggesting the role of lncRNAs in vernalization differ between these two species.

LncRNA33732-respiratory burst oxidase module associated with WRKY1 in tomato- Phytophthora infestans interactions

Our previous studies indicated that tomato WRKY1 transcription factor acts as a positive regulator during tomato resistance to Phytophthora infestans. However, the molecular mechanism of WRKY1-mediated resistance regulation remains unclear. Here, we used a comparative transcriptome analysis between wild-type and WRKY1-overexpressing tomato plants to identify differentially expressed genes (DEGs) and long non-coding RNAs (DELs), and we examined long non-coding RNA (lncRNA)-gene networks. The promoter sequences of the upregulated DEGs and DELs were analyzed. Among 1073 DEGs and 199 DELs, 1 kb 5'-upstream regions of 59 DEGs and 22 DELs contain the W-box, the target sequence of the WRKY1. The results of promoter-β-glucuronidase (GUS) fusion and yeast one-hybrid assay showed that lncRNA33732 was activated by WRKY1 through sequence-specific interactions with the W-box element in its promoter. The overexpression and silencing analysis of lncRNA33732 in tomato showed that lncRNA33732 acts as a positive regulator and enhanced tomato resistance to P. infestans by induction of the expression of respiratory burst oxidase (RBOH) and increase in the accumulation of H₂ O₂ . When the expression of RBOH gene was inhibited in tomato plants, H₂ O₂ accumulation decreased and resistance were impaired. These findings suggest that lncRNA33732 activated by WRKY1 induces RBOH expression to increase H₂ O₂ accumulation in early defense reaction of tomato to P. infestans attack. Our results provide insights into the WRKY1-lncRNA33732-RBOH module involved in the regulation of H₂ O₂ accumulation and resistance to P. infestans, as well as provide candidates to enhance broad-spectrum resistance to pathogens in tomato.

Transcriptomic and functional analyses unveil the role of long non-coding RNAs in anthocyanin biosynthesis during sea buckthorn fruit ripening

Fruit ripening is a developmental process regulated by a complex network of endogenous and exogenous cues. Sea buckthorn is an excellent material for fruit ripening studies due to its dramatic ripening process and high contents of nutritional and anti-oxidant compounds in berries. Here, the whole transcriptome of sea buckthorn fruit at three development stages were analysed using multiple high-throughput sequencings. We assembled and annotated 9,008 long non-coding RNAs (lncRNAs) in sea buckthorn fruits, and identified 118 differentially expressed lncRNAs (DE-lncRNAs) and 32 differentially expressed microRNAs in fruit developmental process. In addition, we predicted 1,061 cis-regulated and 782 trans-regulated targets of DE-lncRNAs, and these DE-lncRNAs are specifically enriched in the biosynthesis of ascorbic acid, carotenoids and flavonoids. Moreover, the silencing of two lncRNAs (LNC1 and LNC2) in vivo and expression analysis revealed that LNC1 and LNC2 can act as endogenous target mimics of miR156a and miR828a to reduce SPL9 and induce MYB114 expression, respectively, which lead to increased and decreased anthocyanin content as revealed by high-performance liquid chromatography analysis. Our results present the first global functional analysis of lncRNA in sea buckthorn and provide two essential regulators of anthocyanin biosynthesis, which provides new insights into the regulation of fruit quality.

Cold-Dependent Expression and Alternative Splicing of Arabidopsis Long Non-coding RNAs

Plants re-program their gene expression when responding to changing environmental conditions. Besides differential gene expression, extensive alternative splicing (AS) of pre-mRNAs and changes in expression of long non-coding RNAs (lncRNAs) are associated with stress responses. RNA-sequencing of a diel time-series of the initial response of Arabidopsis thaliana rosettes to low temperature showed massive and rapid waves of both transcriptional and AS activity in protein-coding genes. We exploited the high diversity of transcript isoforms in AtRTD2 to examine regulation and post-transcriptional regulation of lncRNA gene expression in response to cold stress. We identified 135 lncRNA genes with cold-dependent differential expression (DE) and/or differential alternative splicing (DAS) of lncRNAs including natural antisense RNAs, sORF lncRNAs, and precursors of microRNAs (miRNAs) and trans-acting small-interfering RNAs (tasiRNAs). The high resolution (HR) of the time-series allowed the dynamics of changes in transcription and AS to be determined and identified early and adaptive transcriptional and AS changes in the cold response. Some lncRNA genes were regulated only at the level of AS and using plants grown at different temperatures and a HR time-course of the first 3 h of temperature reduction, we demonstrated that the AS of some lncRNAs is highly sensitive to small temperature changes suggesting tight regulation of expression. In particular, a splicing event in TAS1a which removed an intron that contained the miR173 processing and phased siRNAs generation sites was differentially alternatively spliced in response to cold. The cold-induced reduction of the spliced form of TAS1a and of the tasiRNAs suggests that splicing may enhance production of the siRNAs. Our results identify candidate lncRNAs that may contribute to the regulation of expression that determines the physiological processes essential for acclimation and freezing tolerance.

Copaifera langsdorffii Novel Putative Long Non-Coding RNAs: Interspecies Conservation Analysis in Adaptive Response to Different Biomes

Long non-coding RNAs (lncRNAs) are involved in multiple regulatory pathways and its versatile form of action has disclosed a new layer in gene regulation. LncRNAs have their expression levels modulated during plant development, and in response to stresses with tissue-specific functions. In this study, we analyzed lncRNA from leaf samples collected from the legume Copaifera langsdorffii Desf. (copaíba) present in two divergent ecosystems: Cerrado (CER; Ecological Station of Botanical Garden in Brasília, Brazil) and Atlantic Rain Forest (ARF; Rio de Janeiro, Brazil). We identified 8020 novel lncRNAs, and they were compared to seven Fabaceae genomes and transcriptomes, to which 1747 and 2194 copaíba lncRNAs were mapped, respectively, to at least one species. The secondary structures of the lncRNAs that were conserved and differentially expressed between the populations were predicted using in silico methods. A few selected lncRNA were confirmed by RT-qPCR in the samples from both biomes; Additionally, the analysis of the lncRNA sequences predicted that some might act as microRNA (miRNA) targets or decoys. The emerging studies involving lncRNAs function and conservation have shown their involvement in several types of biotic and abiotic stresses. Thus, the conservation of lncRNAs among Fabaceae species considering their rapid turnover, suggests they are likely to have been under functional conservation pressure. Our results indicate the potential involvement of lncRNAs in the adaptation of C. langsdorffii in two different biomes.

Genome-wide analysis of long non-coding RNAs in Catalpa bungei and their potential function in floral transition using high-throughput sequencing

BACKGROUND

Long non-coding RNAs (lncRNAs) have crucial roles in various biological regulatory processes. However, the study of lncRNAs is limited in woody plants. Catalpa bungei is a valuable ornamental tree with a long cultivation history in China, and a deeper understanding of the floral transition mechanism in C. bungei would be interesting from both economic and scientific perspectives.

RESULTS

In this study, we categorized C. bungei buds from early flowering (EF) and normal flowering (NF) varieties into three consecutive developmental stages. These buds were used to systematically study lncRNAs during floral transition using high-throughput sequencing to identify molecular regulatory networks. Quantitative real-time PCR was performed to study RNA expression changes in different stages. In total, 12,532 lncRNAs and 26,936 messenger RNAs (mRNAs) were detected. Moreover, 680 differentially expressed genes and 817 differentially expressed lncRNAs were detected during the initiation of floral transition. The results highlight the mRNAs and lncRNAs that may be involved in floral transition, as well as the many lncRNAs serving as microRNA precursors. We predicted the functions of lncRNAs by analysing the relationships between lncRNAs and mRNAs. Seven lncRNA-mRNA interaction pairs may participate in floral transition.

CONCLUSIONS

This study is the first to identify lncRNAs and their potential functions in floral transition, providing a starting point for detailed determination of the functions of lncRNAs in C. bungei.

Long Non-Coding RNAs as Endogenous Target Mimics and Exploration of Their Role in Low Nutrient Stress Tolerance in Plants

Long non-coding RNA (lncRNA) research in plants has recently gained momentum taking cues from studies in animals systems. The availability of next-generation sequencing has enabled genome-wide identification of lncRNA in several plant species. Some lncRNAs are inhibitors of microRNA expression and have a function known as target mimicry with the sequestered transcript known as an endogenous target mimic (eTM). The lncRNAs identified to date show diverse mechanisms of gene regulation, most of which remain poorly understood. In this review, we discuss the role of identified putative lncRNAs that may act as eTMs for nutrient-responsive microRNAs (miRNAs) in plants. If functionally validated, these putative lncRNAs would enhance current understanding of the role of lncRNAs in nutrient homeostasis in plants.

Comparative Transcriptome Profiling Analysis of Red- and White-Fleshed Strawberry (Fragaria�ananassa) Provides New Insight into the Regulation of the Anthocyanin Pathway

Anthocyanins are water-soluble pigments in plants. They confer both economic and healthy profits for humans. To gain a deeper insight into the regulation of anthocyanin biosynthesis in octoploid strawberry (Fragaria�ananassa; Fa), a widely consumed economically important fruit, we performed comparative transcriptomic analysis of red- and white-fleshed strawberry cultivars in two ripening stages. In total, 365,455 non-redundant transcripts were assembled from the RNA sequencing (RNAseq) data. Of this collection, 377 were annotated as putative anthocyanin-related transcripts. Differential expression analysis revealed that 57 anthocyanin biosynthesis transcripts were down-regulated, and 89 transcription factors (TFs) were either down- or up-regulated under anthocyanin deficiency. Additionally, amongst the 50,601 putative long non-coding RNAs (lncRNAs) identified here, 68 lncRNAs were differentially expressed and co-expressed with differentially expressed anthocyanin-related mRNAs; 2,070 co-expressing lncRNA-mRNA pairs were generated. Expression profile analysis revealed that it was the limited expression of FaF3'H (flavonoid 3'-hydroxylase) that blocked the cyanidin 3-glucoside accumulation in the two investigated strawberry cultivars. This was further supported by a transient overexpression experiment with FaMYB10. The down-regulated lncRNAs might participate in anthocyanin regulation by acting as targets for microRNAs (miRNAs). The level of competitive intensity in miRNA and lncRNA for the same mRNA targets was probably lower in the white-fleshed strawberries, which can release the repression effect of the mRNAs in red-fleshed strawberry as a result. This study for the first time presents lncRNAs related to anthocyanins in strawberries, provides new insights into the anthocyanin regulatory network and also lays the foundation for identifying new anthocyanin regulators in strawberry.

ceRNA Cross-Talk in Paulownia Witches' Broom Disease

Long noncoding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA) are important in the regulation of life activities. However, their function is unclear in Paulownia fortunei. To identify lncRNAs, circRNAs, and miRNA, and investigate their roles in the infection progress of Paulownia witches’ broom (PaWB) disease, we performed RNA sequencing of healthy and infected P. fortunei. A total of 3126 lncRNAs, 1634 circRNAs, and 550 miRNAs were identified. Among them, 229 lncRNAs, 65 circRNAs, and 65 miRNAs were differentially expressed in a significant manner. We constructed a competing endogenous RNA (ceRNA) network, which contains 5 miRNAs, 4 circRNAs, 5 lncRNAs, and 15 mRNAs, all of which were differentially expressed between healthy and infected P. fortunei. This study provides the first catalog of candidate ceRNAs in Paulownia and gives a revealing insight into the molecular mechanism responsible for PaWB.

Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies

BACKGROUND

Nitrogen (N) is a key macronutrient essential for plant growth, and its availability has a strong influence on crop development. The application of synthetic N fertilizers on crops has increased substantially in recent decades; however, the applied N is not fully utilized due to the low N use efficiency of crops. To overcome this limitation, it is important to understand the genome-wide responses and functions of key genes and potential regulatory factors in N metabolism.

RESULTS

Here, we characterized changes in the rice (Oryza sativa) transcriptome, including genes, newly identified putative long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) and their target mRNAs in response to N starvation using four different transcriptome approaches. Analysis of rice genes involved in N metabolism and/or transport using strand-specific RNA-Seq identified 2588 novel putative lncRNA encoding loci. Analysis of previously published RNA-Seq datasets revealed a group of N starvation-responsive lncRNAs showing differential expression under other abiotic stress conditions. Poly A-primed sequencing (2P-Seq) revealed alternatively polyadenylated isoforms of N starvation-responsive lncRNAs and provided precise 3' end information on the transcript models of these lncRNAs. Analysis of small RNA-Seq data identified N starvation-responsive miRNAs and down-regulation of miR169 family members, causing de-repression of NF-YA, as confirmed by strand-specific RNA-Seq and qRT-PCR. Moreover, we profiled the N starvation-responsive down-regulation of root-specific miRNA, osa-miR444a.4-3p, and Degradome sequencing confirmed MADS25 as a novel target gene.

CONCLUSIONS

In this study, we used a combination of multiple RNA-Seq analyses to extensively profile the expression of genes, newly identified lncRNAs, and microRNAs in N-starved rice roots and shoots. Data generated in this study provide an in-depth understanding of the regulatory pathways modulated by N starvation-responsive miRNAs. The results of comprehensive, large-scale data analysis provide valuable information on multiple aspects of the rice transcriptome, which may be useful in understanding the responses of rice plants to changes in the N supply status of soil.

Analysis of long-non-coding RNAs associated with ethylene in tomato

Long-Non-Coding RNAs (LncRNAs) are a class of non-coding endogenous RNAs contributing to numerous biological processes. LeERF1 is a tomato ethylene response factor (ERF) near the end of the ethylene signal transduction pathway. To identify lncRNAs in tomato and elucidate their roles in ethylene signaling, deep sequencing was deployed in over-expression and repression LeERF1 transgenic and control tomato fruits. A total of 397 lncRNAs were identified, including 169 tomato lncRNAs that had not previously been identified. Among these, 12 were differentially expressed between the transgenic and control tomato fruits. Numerous lncRNA target genes were identified including many associated with ethylene signaling including auxin response factors and auxin-induced proteins, F-box proteins, ERFs and MADS-box proteins. In addition, two lncRNAs were found to be the precursor of three miRNAs and four lncRNAs could be targeted by five miRNAs. We propose a regulatory model highlighting the relationships between lncRNAs and their targets involved in ethylene signal transduction which establishes a foundation for addressing the role of LncRNAs in ethylene response.

Genome-Wide Identification and Functional Prediction of Novel Drought-Responsive lncRNAs in Pyrus betulifolia

Increasing evidence shows that long noncoding RNAs (lncRNAs) play important roles in developmental regulation and many other biological processes in plants. However, identification of lncRNAs in Pyrus betulifolia is limited compared with studies of functional gene expression. Using high-throughput sequencing technology, the transcriptome of P. betulifolia under drought stress was analyzed to identify lncRNAs. A total of 14,478 lncRNAs were identified, of which 251 were found to be drought-responsive. The putative target genes of these differentially expressed lncRNAs were significantly enriched in metabolic processes, organic substance metabolic processes, macromolecule metabolic processes, and heterocyclic compound binding. Real-time quantitative polymerase chain reaction validation suggested that the results of the RNA sequencing data analysis were reliable. This study will provide genetic resources for pear breeding and provide reference to other pomological studies.

Widespread antisense transcription of Populus genome under drought

Antisense transcription is widespread in many genomes and plays important regulatory roles in gene expression. The objective of our study was to investigate the extent and functional relevance of antisense transcription in forest trees. We employed Populus, a model tree species, to probe the antisense transcriptional response of tree genome under drought, through stranded RNA-seq analysis. We detected nearly 48% of annotated Populus gene loci with antisense transcripts and 44% of them with co-transcription from both DNA strands. Global distribution of reads pattern across annotated gene regions uncovered that antisense transcription was enriched in untranslated regions while sense reads were predominantly mapped in coding exons. We further detected 1185 drought-responsive sense and antisense gene loci and identified a strong positive correlation between the expression of antisense and sense transcripts. Additionally, we assessed the antisense expression in introns and found a strong correlation between intronic expression and exonic expression, confirming antisense transcription of introns contributes to transcriptional activity of Populus genome under drought. Finally, we functionally characterized drought-responsive sense-antisense transcript pairs through gene ontology analysis and discovered that functional groups including transcription factors and histones were concordantly regulated at both sense and antisense transcriptional level. Overall, our study demonstrated the extensive occurrence of antisense transcripts of Populus genes under drought and provided insights into genome structure, regulation pattern and functional significance of drought-responsive antisense genes in forest trees. Datasets generated in this study serve as a foundation for future genetic analysis to improve our understanding of gene regulation by antisense transcription.

Re-analysis of long non-coding RNAs and prediction of circRNAs reveal their novel roles in susceptible tomato following TYLCV infection

BACKGROUND

Long Noncoding-RNAs (LncRNAs) are known to be involved in some biological processes, but their roles in plant-virus interactions remain largely unexplored. While circular RNAs (circRNAs) have been studied in animals, there has yet to be extensive research on them in a plant system, especially in tomato-tomato yellow leaf curl virus (TYLCV) interaction.

RESULTS

In this study, RNA transcripts from the susceptible tomato line JS-CT-9210 either infected with TYLCV or untreated, were sequenced in a pair-end strand-specific manner using ribo-zero rRNA removal library method. A total of 2056 lncRNAs including 1767 long intergenic non-coding RNA (lincRNAs) and 289 long non-coding natural antisense transcripts (lncNATs) were obtained. The expression patterns in lncRNAs were similar in susceptible tomato plants between control check (CK) and TYLCV infected samples. Our analysis suggested that lncRNAs likely played a role in a variety of functions, including plant hormone signaling, protein processing in the endoplasmic reticulum, RNA transport, ribosome function, photosynthesis, glulathione metabolism, and plant-pathogen interactions. Using virus-induced gene silencing (VIGS) analysis, we found that reduced expression of the lncRNA S-slylnc0957 resulted in enhanced resistance to TYLCV in susceptible tomato plants. Moreover, we identified 184 circRNAs candidates using the CircRNA Identifier (CIRI) software, of which 32 circRNAs were specifically expressed in untreated samples and 83 circRNAs in TYLCV samples. Approximately 62% of these circRNAs were derived from exons. We validated the circRNAs by both PCR and Sanger sequencing using divergent primers, and found that most of circRNAs were derived from the exons of protein coding genes. The silencing of these circRNAs parent genes resulted in decreased TYLCV virus accumulation.

CONCLUSION

In this study, we identified novel lncRNAs and circRNAs using bioinformatic approaches and showed that these RNAs function as negative regulators of TYLCV infection. Moreover, the expression patterns of lncRNAs in susceptible tomato plants were different from that of resistant tomato plants, while exonic circRNAs expression positively associated with their respective protein coding genes. This work provides a foundation for elaborating the novel roles of lncRNAs and circRNAs in susceptible tomatoes following TYLCV infection.

Integrative analysis of genome-wide lncRNA and mRNA expression in newly synthesized Brassica hexaploids

Polyploidization, as a significant evolution force, has been considered to facilitate plant diversity. The expression levels of lncRNAs and how they control the expression of protein-coding genes in allopolyploids remain largely unknown. In this study, lncRNA expression profiles were compared between Brassica hexaploid and its parents using a high-throughput sequencing approach. A total of 2,725, 1,672, and 2,810 lncRNAs were discovered in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. It was also discovered that 725 lncRNAs were differentially expressed between Brassica hexaploid and its parents, and 379 lncRNAs were nonadditively expressed in this hexaploid. LncRNAs have multiple expression patterns between Brassica hexaploid and its parents and show paternal parent-biased expression. These lncRNAs were found to implement regulatory functions directly in the long-chain form, and acted as precursors or targets of miRNAs. According to the prediction of the targets of differentially expressed lncRNAs, 109 lncRNAs were annotated, and their target genes were involved in the metabolic process, pigmentation, reproduction, exposure to stimulus, biological regulation, and so on. Compared with the paternal parent, differentially expressed lncRNAs between Brassica hexaploid and its maternal parent participated in more regulation pathways. Additionally, 61 lncRNAs were identified as putative targets of known miRNAs, and 15 other lncRNAs worked as precursors of miRNAs. Some conservative motifs of lncRNAs from different groups were detected, which indicated that these motifs could be responsible for their regulatory roles. Our findings may provide a reference for the further study of the function and action mechanisms of lncRNAs during plant evolution.

Long Non-Coding RNAs Responsive to Salt and Boron Stress in the Hyper-Arid Lluteño Maize from Atacama Desert

Long non-coding RNAs (lncRNAs) have been defined as transcripts longer than 200 nucleotides, which lack significant protein coding potential and possess critical roles in diverse cellular processes. Long non-coding RNAs have recently been functionally characterized in plant stress-response mechanisms. In the present study, we perform a comprehensive identification of lncRNAs in response to combined stress induced by salinity and excess of boron in the Lluteño maize, a tolerant maize landrace from Atacama Desert, Chile. We use deep RNA sequencing to identify a set of 48,345 different lncRNAs, of which 28,012 (58.1%) are conserved with other maize (B73, Mo17 or Palomero), with the remaining 41.9% belonging to potentially Lluteño exclusive lncRNA transcripts. According to B73 maize reference genome sequence, most Lluteño lncRNAs correspond to intergenic transcripts. Interestingly, Lluteño lncRNAs presents an unusual overall higher expression compared to protein coding genes under exposure to stressed conditions. In total, we identified 1710 putatively responsive to the combined stressed conditions of salt and boron exposure. We also identified a set of 848 stress responsive potential trans natural antisense transcripts (trans-NAT) lncRNAs, which seems to be regulating genes associated with regulation of transcription, response to stress, response to abiotic stimulus and participating of the nicotianamine metabolic process. Reverse transcription-quantitative PCR (RT-qPCR) experiments were performed in a subset of lncRNAs, validating their existence and expression patterns. Our results suggest that a diverse set of maize lncRNAs from leaves and roots is responsive to combined salt and boron stress, being the first effort to identify lncRNAs from a maize landrace adapted to extreme conditions such as the Atacama Desert. The information generated is a starting point to understand the genomic adaptabilities suffered by this maize to surpass this extremely stressed environment.

Concurrent isotope-assisted metabolic flux analysis and transcriptome profiling reveal responses of poplar cells to altered nitrogen and carbon supply

Reduced nitrogen is indispensable to plants. However, its limited availability in soil combined with the energetic and environmental impacts of nitrogen fertilizers motivates research into molecular mechanisms toward improving plant nitrogen use efficiency (NUE). We performed a systems-level investigation of this problem by employing multiple 'omics methodologies on cell suspensions of hybrid poplar (Populus tremula × Populus alba). Acclimation and growth of the cell suspensions in four nutrient regimes ranging from abundant to deficient supplies of carbon and nitrogen revealed that cell growth under low-nitrogen levels was associated with substantially higher NUE. To investigate the underlying metabolic and molecular mechanisms, we concurrently performed steady-state ¹³ C metabolic flux analysis with multiple isotope labels and transcriptomic profiling with cDNA microarrays. The ¹³ C flux analysis revealed that the absolute flux through the oxidative pentose phosphate pathway (oxPPP) was substantially lower (~threefold) under low-nitrogen conditions. Additionally, the flux partitioning ratio between the tricarboxylic acid cycle and anaplerotic pathways varied from 84%:16% under abundant carbon and nitrogen to 55%:45% under deficient carbon and nitrogen. Gene expression data, together with the flux results, suggested a plastidic localization of the oxPPP as well as transcriptional regulation of certain metabolic branchpoints, including those between glycolysis and the oxPPP. The transcriptome data also indicated that NUE-improving mechanisms may involve a redirection of excess carbon to aromatic metabolic pathways and extensive downregulation of potentially redundant genes (in these heterotrophic cells) that encode photosynthetic and light-harvesting proteins, suggesting the recruitment of these proteins as nitrogen sinks in nitrogen-abundant conditions.

A Novel LncRNA, MuLnc1, Associated With Environmental Stress in Mulberry (Morus multicaulis)

Environmental stresses are major constraints that limit the leaf productivity and quality of mulberry. LncRNAs have emerged as important regulators in response to biotic and abiotic stresses in plants. However, the functions and mechanisms of most lncRNAs remain largely unknown. A novel lncRNA designated as MuLnc1 was found to be cleaved by mul-miR3954 and produce secondary siRNAs in a 21 nt phase in mulberry. It was demonstrated that one of the siRNAs produced, si161579, can silence the expression of the calmodulin-like protein gene CML27 of mulberry (MuCML27). When MuCML27 was heterologously expressed in Arabidopsis, the transgenic plants exhibited enhanced resistance to Botrytis cinerea and Pseudomonas syringae pv tomato DC3000. In addition, the transgenic MuCML27-overexpressing Arabidopsis plants are more tolerant to salt and drought stresses. Furthermore, the network of mul-miR3954-MuLnc1-siRNAs-mRNAs was modeled to elucidate the interaction between lncRNAs and sRNAs with mRNAs. All of these, taken together, suggest that MuLnc1 was associated with environmental stress in mulberry and may be considered as a potential genetic improvement target gene of mulberry. The information provided may shed light on the complicated gene expression regulatory mechanisms in mulberry stress responses.

Non-coding RNAs and Their Roles in Stress Response in Plants

Eukaryotic genomes encode thousands of non-coding RNAs (ncRNAs), which play crucial roles in transcriptional and post-transcriptional regulation of gene expression. Accumulating evidence indicates that ncRNAs, especially microRNAs (miRNAs) and long ncRNAs (lncRNAs), have emerged as key regulatory molecules in plant stress responses. In this review, we have summarized the current progress on the understanding of plant miRNA and lncRNA identification, characteristics, bioinformatics tools, and resources, and provided examples of mechanisms of miRNA- and lncRNA-mediated plant stress tolerance.

Genomewide differential expression profiling of long non-coding RNAs in androgenetic alopecia in a Chinese male population

BACKGROUND

Androgenetic alopecia (AGA), or male pattern baldness (MPB), is the most common form of hair loss in males. A combination of genetic and androgen causes have been suggested as factors that contribute to the development of AGA. However, the specific molecular mechanisms that underly AGA remain largely unknown. Long non-coding RNAs (lncRNAs), a new class of regulatory non-coding RNAs that are longer than 200 nucleotides, have been shown to play important roles in a number of cellular processes, including transcription, chromosome remodelling and post-transcriptional processing. The dysregulation of lncRNAs is associated with many forms of diseases, but it remains unknown whether lncRNAs are associated with AGA.

OBJECTIVE

The aim of this study was to identify AGA-associated lncRNAs and predict the potential roles of these lncRNAs in AGA.

METHODS

A genomewide microarray was used to identify lncRNAs that are differentially expressed between AGA and adjacent normal tissues. Real-time qRT-PCR was used to validate the microarray data.

RESULTS

A large number of lncRNAs were differentially expressed (fold change >2.4) between AGA and adjacent normal tissues. Of these, 770 were upregulated and 1373 were downregulated. Moreover, pathway analysis revealed that 53 functional pathways were associated with the upregulated transcripts, while 11 pathways were associated with the downregulated transcripts.

CONCLUSION

To our knowledge, this is the first study to investigate AGA-associated lncRNAs. lncRNA profiles are altered in AGA, and these lncRNAs and their target genes may serve as novel candidates for preventing and treating AGA.

Present Scenario of Long Non-Coding RNAs in Plants

Small non-coding RNAs have been extensively studied in plants over the last decade. In contrast, genome-wide identification of plant long non-coding RNAs (lncRNAs) has recently gained momentum. LncRNAs are now being recognized as important players in gene regulation, and their potent regulatory roles are being studied comprehensively in eukaryotes. LncRNAs were first reported in humans in 1992. Since then, research in animals, particularly in humans, has rapidly progressed, and a vast amount of data has been generated, collected, and organized using computational approaches. Additionally, numerous studies have been conducted to understand the roles of these long RNA species in several diseases. However, the status of lncRNA investigation in plants lags behind that in animals (especially humans). Efforts are being made in this direction using computational tools and high-throughput sequencing technologies, such as the lncRNA microarray technique, RNA-sequencing (RNA-seq), RNA capture sequencing, (RNA CaptureSeq), etc. Given the current scenario, significant amounts of data have been produced regarding plant lncRNAs, and this amount is likely to increase in the subsequent years. In this review we have documented brief information about lncRNAs and their status of research in plants, along with the plant-specific resources/databases for information retrieval on lncRNAs.

Systematic comparison of lncRNAs with protein coding mRNAs in population expression and their response to environmental change

BACKGROUND

Long non-coding RNA (lncRNA) is a class of non-coding RNA with important regulatory roles in biological process of organisms. The systematic comparison of lncRNAs with protein coding mRNAs in population expression and their response to environmental change are still poorly understood. Here we identified 17,610 lncRNAs and calculated their expression levels based on RNA-seq of 80 individuals of Miscanthus lutarioriparius from two environments, the nearly native habitats and transplanted field, respectively.

RESULTS

LncRNAs had significantly higher expression diversity and lower expression frequency in population than protein coding mRNAs in both environments, which suggested that lncRNAs may experience more relaxed selection or divergent evolution in population compared with protein coding RNAs. In addition, the increase of expression diversity for lncRNAs was always significantly higher and the magnitude of fold change of expression in new stress environment was significantly larger than protein-coding mRNAs. These results suggested that lncRNAs may be more sensitive to environmental change than protein-coding mRNAs. Analysis of environment-robust and environment-specific lncRNA-mRNA co-expression network between two environments revealed the characterization of lncRNAs in response to environmental change. Furthermore, candidate lncRNAs contributing to water use efficiency (WUE) identified based on the WUE-lncRNA-mRNA co-expression network suggested the roles of lncRNAs in response to environmental change.

CONCLUSION

Our study provided a comprehensive understanding of expression characterization of lncRNAs in population for M. lutarioriparius under field condition, which would be useful to explore the roles of lncRNAs and could accelerate the process of adaptation in new environment for many plants.

Comparative transcriptome analysis between resistant and susceptible tomato allows the identification of lncRNA16397 conferring resistance to Phytophthora infestans by co-expressing glutaredoxin

The rapid development of omics sequencing technology has facilitated the identification of thousands of long non-coding (lnc)RNAs in plant species, but the role of lncRNAs in plant-pathogen interactions remains largely unexplored. We used comparative transcriptome analysis of Phytophthora infestans-resistant and -susceptible tomatoes to identify differentially expressed genes (DEGs) and lncRNAs (DELs), and examine lncRNA-mRNA networks. A total of 1037 DEGs and 688 DELs were identified between P. infestans-resistant and -susceptible tomatoes. The co-localization networks, including 128 DEGs and 127 DELs, were performed. We found that lncRNA16397 acted as an antisense transcript of SlGRX22 to regulate its expression, and also induced SlGRX21 expression when lncRNA16397 was overexpressed. In addition, disease symptoms and reactive oxygen species (ROS) accumulation in tomatoes overexpressing lncRNA16397 and SpGRX were fewer and lower than those in wild-type after P. infestans infection. This result suggests that tomato lncRNA16397 induces SlGRX expression to reduce ROS accumulation and alleviate cell membrane injury, resulting in enhanced resistance to P. infestans. Our results provide insight into lncRNAs involved in the response of tomato to P. infestans infection, demonstrate that the lncRNA16397-GRXs network is an important component of the P. infestans network in tomato, and provide candidates for breeding to enhance biotic stress-resistance in tomato.

Long noncoding RNA (lincRNA), a new paradigm in gene expression control

Long intergenic non-coding RNAs (lincRNAs) are defined as RNA transcripts that are longer than 200 nucleotides. By definition, these RNAs must not have open reading frames that encode proteins. Many of these transcripts are encoded by RNA polymerase II, are spliced, and are poly-adenylated. This final fact indicates that there is a trove of information about lincRNAs in databases such as the Gene Expression Omnibus (GEO), which is a repository for RNAseq and microarray data. Recent experiments indicate that there are upwards of 15,000 lincRNAs encoded by the human genome. The term "intergenic" refers to the identification of these transcripts from regions of the genome that do not contain protein-encoding genes. These regions coincide with what was once labeled as the "junk DNA" portions of our genomes, which, upon careful examination by whole genome RNA sequencing experiments, clearly encode RNA transcripts. LincRNAs also contain promoter- or enhancer-associated RNAs that are gene proximal and can be either in the sense or antisense orientation, relative to the protein-coding gene with which they are associated. In this review, we describe the functions of lincRNAs playing roles in biological processes such as gene expression control, scaffold formation, and epigenetic control.