Arabidopsis noncoding RNA mediates control of photomorphogenesis by red light

First Come, First Served: Sui Generis Features of the First Intron

Most of the transcribed genes in eukaryotic cells are interrupted by intervening sequences called introns that are co-transcriptionally removed from nascent messenger RNA through the process of splicing. In Arabidopsis, 79% of genes contain introns and more than 60% of intron-containing genes undergo alternative splicing (AS), which ostensibly is considered to increase protein diversity as one of the intrinsic mechanisms for fitness to the varying environment or the internal developmental program. In addition, recent findings have prevailed in terms of overlooked intron functions. Here, we review recent progress in the underlying mechanisms of intron function, in particular by focusing on unique features of the first intron that is located in close proximity to the transcription start site. The distinct deposition of epigenetic marks and nucleosome density on the first intronic DNA sequence, the impact of the first intron on determining the transcription start site and elongation of its own expression (called intron-mediated enhancement, IME), translation control in 5′-UTR, and the new mechanism of the trans-acting function of the first intron in regulating gene expression at the post-transcriptional level are summarized.

Splicing conservation signals in plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) have recently emerged as prominent regulators of gene expression in eukaryotes. LncRNAs often drive the modification and maintenance of gene activation or gene silencing states via chromatin conformation rearrangements. In plants, lncRNAs have been shown to participate in gene regulation, and are essential to processes such as vernalization and photomorphogenesis. Despite their prominent functions, only over a dozen lncRNAs have been experimentally and functionally characterized. Similar to its animal counterparts, the rates of sequence divergence are much higher in plant lncRNAs than in protein coding mRNAs, making it difficult to identify lncRNA conservation using traditional sequence comparison methods. Beyond this, little is known about the evolutionary patterns of lncRNAs in plants. Here, we characterized the splicing conservation of lncRNAs in Brassicaceae. We generated a whole-genome alignment of 16 Brassica species and used it to identify synthenic lncRNA orthologs. Using a scoring system trained on transcriptomes from A. thaliana and B. oleracea, we identified splice sites across the whole alignment and measured their conservation. Our analysis revealed that 17.9% (112/627) of all intergenic lncRNAs display splicing conservation in at least one exon, an estimate that is substantially higher than previous estimates of lncRNA conservation in this group. Our findings agree with similar studies in vertebrates, demonstrating that splicing conservation can be evidence of stabilizing selection. We provide conclusive evidence for the existence of evolutionary deeply conserved lncRNAs in plants and describe a generally applicable computational workflow to identify functional lncRNAs in plants.

Unique and contrasting effects of light and temperature cues on plant transcriptional programs

Plants have adapted to tolerate and survive constantly changing environmental conditions by reprogramming gene expression in response to stress or to drive developmental transitions. Among the many signals that plants perceive, light and temperature are of particular interest due to their intensely fluctuating nature which is combined with a long-term seasonal trend. Whereas specific receptors are key in the light-sensing mechanism, the identity of plant thermosensors for high and low temperatures remains far from fully addressed. This review aims at discussing common as well as divergent characteristics of gene expression regulation in plants, controlled by light and temperature. Light and temperature signaling control the abundance of specific transcription factors, as well as the dynamics of co-transcriptional processes such as RNA polymerase elongation rate and alternative splicing patterns. Additionally, sensing both types of cues modulates gene expression by altering the chromatin landscape and through the induction of long non-coding RNAs (lncRNAs). However, while light sensing is channeled through dedicated receptors, temperature can broadly affect chemical reactions inside plant cells. Thus, direct thermal modifications of the transcriptional machinery add another level of complexity to plant transcriptional regulation. Besides the rapid transcriptome changes that follow perception of environmental signals, plant developmental transitions and acquisition of stress tolerance depend on long-term maintenance of transcriptional states (active or silenced genes). Thus, the rapid transcriptional response to the signal (Phase I) can be distinguished from the long-term memory of the acquired transcriptional state (Phase II - remembering the signal). In this review we discuss recent advances in light and temperature signal perception, integration and memory in Arabidopsis thaliana, focusing on transcriptional regulation and highlighting the contrasting and unique features of each type of cue in the process.

Long non-coding RNAs in the alkaline stress response in sugar beet (Beta vulgaris L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) play crucial roles in regulating numerous biological processes in which complicated mechanisms are involved. Nonetheless, little is known about the number, features, sequences, and possible effects of lncRNAs on plant responses to alkaline stress.

RESULTS

Leaf samples collected based on the control Beta vulgaris L., as well as those under short-term and long-term alkaline treatments, were subjected to high-throughput RNA sequencing, through which a total of 8535 lncRNAs with reliable expression were detected. Of these lncRNAs, 102 and 49 lncRNA expression profiles were altered after short- and long-term alkaline stress, respectively. Moreover, 7 lncRNAs were recognized as precursors to 17 previously identified miRNAs. Four lncRNAs responsive to alkaline stress were estimated as targets for 8 miRNAs. Moreover, computational analysis predicted 4318 potential target genes as lncRNAs responsive to alkaline stress. Analysis of functional annotations showed that the abovementioned possible target genes were involved in various bioprocesses, such as kinase activity, structural constituents of ribosomes, the ribonucleoprotein complex and protein metabolic processes. Association analysis provided convincing proof of the interplay of specific candidate target genes with lncRNAs.

CONCLUSION

LncRNAs likely exert vital roles during the regulation of the alkaline stress response and adaptation in plants through interaction with protein-coding genes. The findings of this study contribute to comprehensively examining lncRNAs in Beta vulgaris L. and shed more light on the possible roles and modulating interplays of lncRNAs responsive to alkaline stress, thereby laying a certain basis for functional analyses of these types of Beta vulgaris L. lncRNAs in the future.

PlantRegMap: charting functional regulatory maps in plants

With the goal of charting plant transcriptional regulatory maps (i.e. transcription factors (TFs), cis-elements and interactions between them), we have upgraded the TF-centred database PlantTFDB (http://planttfdb.cbi.pku.edu.cn/) to a plant regulatory data and analysis platform PlantRegMap (http://plantregmap.cbi.pku.edu.cn/) over the past three years. In this version, we updated the annotations for the previously collected TFs and set up a new section, ‘extended TF repertoires’ (TFext), to allow users prompt access to the TF repertoires of newly sequenced species. In addition to our regular TF updates, we are dedicated to updating the data on cis-elements and functional interactions between TFs and cis-elements. We established genome-wide conservation landscapes for 63 representative plants and then developed an algorithm, FunTFBS, to screen for functional regulatory elements and interactions by coupling the base-varied binding affinities of TFs with the evolutionary footprints on their binding sites. Using the FunTFBS algorithm and the conservation landscapes, we further identified over 20 million functional TF binding sites (TFBSs) and two million functional interactions for 21 346 TFs, charting the functional regulatory maps of these 63 plants. These resources are publicly available at PlantRegMap (http://plantregmap.cbi.pku.edu.cn/) and a cloud-based mirror (http://plantregmap.gao-lab.org/), providing the plant research community with valuable resources for decoding plant transcriptional regulatory systems.

Long noncoding RNAs shape transcription in plants

The advent of novel high-throughput sequencing techniques has revealed that eukaryotic genomes are massively transcribed although only a small fraction of RNAs exhibits protein-coding capacity. In the last years, long noncoding RNAs (lncRNAs) have emerged as regulators of eukaryotic gene expression in a wide range of molecular mechanisms. Plant lncRNAs can be transcribed by alternative RNA polymerases, acting directly as long transcripts or can be processed into active small RNAs. Several lncRNAs have been recently shown to interact with chromatin, DNA or nuclear proteins to condition the epigenetic environment of target genes or modulate the activity of transcriptional complexes. In this review, we will summarize the recent discoveries about the actions of plant lncRNAs in the regulation of gene expression at the transcriptional level.

The Emerging Role of Long Non-Coding RNAs in Plant Defense Against Fungal Stress

Growing interest and recent evidence have identified long non-coding RNA (lncRNA) as the potential regulatory elements for eukaryotes. LncRNAs can activate various transcriptional and post-transcriptional events that impact cellular functions though multiple regulatory functions. Recently, a large number of lncRNAs have also been identified in higher plants, and an understanding of their functional role in plant resistance to infection is just emerging. Here, we focus on their identification in crop plant, and discuss their potential regulatory functions and lncRNA-miRNA-mRNA network in plant pathogen stress responses, referring to possible examples in a model plant. The knowledge gained from a deeper understanding of this colossal special group of plant lncRNAs will help in the biotechnological improvement of crops.

Genome-wide identification of Arabidopsis long noncoding RNAs in response to the blue light

Long non-coding RNAs (lncRNAs) have been shown in animals to play roles in a wide range of biological processes. In plant, light modulates the growth and development as a key external signal. However, little is known about the role of plant lncRNA in response to light. In this study, we sequenced the messenger RNAs (mRNAs), lncRNAs and microRNAs (miRNAs) in Arabidopsis seedlings under blue light for 2 h and 8 h. Compared to dark, we identified 4197 mRNAs, 375 miRNAs and 481 lncRNAs, or 5207 mRNAs, 286 miRNAs and 545 lncRNAs of differential expressions under blue light treatments for 2 h or 8 h respectively. Subsequently, a total of 407 competing endogenous RNA (ceRNA) pairs (lncRNA-mRNA-miRNA) were constructed. We identified a blue light-induced lncRNA which plays roles in blue light-directed plant photomorphogenesis and response to mannitol stress by serving as a ceRNA to sequester miR167 in a type of target mimicry. These results revealed previously unknown roles of the lncRNA in blue light signaling and mannitol stress, and provided useful resources of lncRNAs associated with miRNAs in response to blue light.

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein-encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue-specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis-related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker-assisted breeding efforts to genetically improve woody plants.

Genome-wide identification of lncRNAs during hickory (Carya cathayensis) flowering

Non-coding RNAs with lengths greater than 200 bp are known as long non-coding RNAs (lncRNAs), and these RNAs play important role in gene regulation and plant development. However, to date, little is known regarding the role played by lncRNAs during flowering in hickory (Carya cathayensis). Here, we performed whole transcriptome RNA-sequencing of samples from hickory female and male floral buds, in which three samples (H0311PF, H0318PF, and H0402PF) represent pre-flowering, flowering, and post-flowering, respectively, while eight male samples collected from May 8th to June 13th as this time course are the key stage for male floral bud differentiation. We identified 2163 lncRNAs in hickory during flowering, containing 213 intronic, 1488 intergenic, and 462 antisense lncRNAs. We noticed that 510 and 648 lncRNAs were differentially expressed corresponding to female and male floral buds, respectively. And some of the lncRNAs were in a tightly tissue-specific or stage-specific manner. To further understand the roles of the lncRNAs, we predicted the function of the lncRNAs in cis- and trans-acting modes. The results showed that 924 lncRNAs were cis-correlated with 1536 protein-coding genes, while 1207 lncRNAs co-expressed (trans-acting) with 7432 protein-coding genes (R > 0.95 or R < - 0.95). These lncRNAs were all enriched in flower development-associated biological processes, i.e., circadian rhythm, vernalization response, response to gibberellin, inflorescence development, floral organ development, etc. To further understand the relationships between lncRNAs and floral-core genes, we build a co-expressing lncRNA-mRNA flowering network. We classified these floral genes into different pathway (photoperiod, vernalization, gibberellin, autonomous, and sucrose pathway) according to their particular functions. We found a set of lncRNAs that preferentially expressed in these pathways. The network showed that some lncRNAs (i.e., XLOC_038669, XLOC_017938) functioned in a particular flowering time pathway, while others (i.e., XLOC_011251, XLOC_04110) were involved in multiple pathway. Furthermore, some lncRNAs (i.e., XLOC_038669, XLOC_009597, and XLOC_049539) played roles in single or multiple pathways via interaction with each other. This study provides a genome-wide survey of hickory flower-related lncRNAs and will contribute to further understanding of the molecular mechanism underpinning flowering in hickory.

Research Progress on Plant Long Non-Coding RNA

Non-coding RNAs (ncRNAs) that were once considered “dark matter” or “transcriptional noise” in genomes are research hotspots in the field of epigenetics. The most well-known microRNAs (miRNAs) are a class of short non-coding, small molecular weight RNAs with lengths of 20–24 nucleotides that are highly conserved throughout evolution. Through complementary pairing with the bases of target sites, target gene transcripts are cleaved and degraded, or translation is inhibited, thus regulating the growth and development of organisms. Unlike miRNAs, which have been studied thoroughly, long non-coding RNAs (lncRNAs) are a group of poorly conserved RNA molecules with a sequence length of more than 200 nucleotides and no protein encoding capability; they interact with large molecules, such as DNA, RNA, and proteins, and regulate protein modification, chromatin remodeling, protein functional activity, and RNA metabolism in vivo through cis- or trans-activation at the transcriptional, post-transcriptional, and epigenetic levels. Research on plant lncRNAs is just beginning and has gradually emerged in the field of plant molecular biology. Currently, some studies have revealed that lncRNAs are extensively involved in plant growth and development and stress response processes by mediating the transmission and expression of genetic information. This paper systematically introduces lncRNA and its regulatory mechanisms, reviews the current status and progress of lncRNA research in plants, summarizes the main techniques and strategies of lncRNA research in recent years, and discusses existing problems and prospects, in order to provide ideas for further exploration and verification of the specific evolution of plant lncRNAs and their biological functions.

Integrative analyses of noncoding RNAs reveal the potential mechanisms augmenting tumor malignancy in lung adenocarcinoma

Precise noncoding RNA (ncRNA)-based network prediction is necessary to reveal ncRNA functions and pathological mechanisms. Here, we established a systemic pipeline to identify prognostic ncRNAs, predict their functions and explore their pathological mechanisms in lung adenocarcinoma (LUAD). After in silico and experimental validation based on evaluations of prognostic value in multiple LUAD cohorts, we selected the PTTG3P pseudogene from among other prognostic ncRNAs (MIR497HG, HSP078, TBX5-AS1, LOC100506990 and C14orf64) for mechanistic studies. PTTG3P upregulation in LUAD cells shortens the metaphase to anaphase transition in mitosis, increases cell viability after cisplatin or paclitaxel treatment, facilitates tumor growth that leads to poor survival in orthotopic lung models, and is associated with a poor survival rate in LUAD patients in the TCGA cohort who received chemotherapy. Mechanistically, PTTG3P acts as an ncRNA that interacts with the transcription factor FOXM1 to regulate the transcriptional activation of the mitotic checkpoint kinase BUB1B, which augments tumor growth and chemoresistance and leads to poor outcomes for LUAD patients. Overall, we established a systematic strategy to uncover prognostic ncRNAs with functional prediction methods suitable for pan-cancer studies. Moreover, we revealed that PTTG3P, due to its upregulation of the PTTG3P/FOXM1/BUB1B axis, could be a therapeutic target for LUAD patients.

Identification and integrated analysis of glyphosate stress-responsive microRNAs, lncRNAs, and mRNAs in rice using genome-wide high-throughput sequencing

BACKGROUND

Glyphosate has become the most widely used herbicide in the world. Therefore, the development of new varieties of glyphosate-tolerant crops is a research focus of seed companies and researchers. The glyphosate stress-responsive genes were used for the development of genetically modified crops, while only the EPSPS gene has been used currently in the study on glyphosate-tolerance in rice. Therefore, it is essential and crucial to intensify the exploration of glyphosate stress-responsive genes, to not only acquire other glyphosate stress-responsive genes with clean intellectual property rights but also obtain non-transgenic glyphosate-tolerant rice varieties. This study is expected to elucidate the responses of miRNAs, lncRNAs, and mRNAs to glyphosate applications and the potential regulatory mechanisms in response to glyphosate stress in rice.

RESULTS

Leaves of the non-transgenic glyphosate-tolerant germplasm CA21 sprayed with 2 mg·ml- 1 glyphosate (GLY) and CA21 plants with no spray (CK) were collected for high-throughput sequencing analysis. A total of 1197 DEGs, 131 DELs, and 52 DEMs were identified in the GLY samples in relation to CK samples. Genes were significantly enriched for various biological processes involved in detoxification of plant response to stress. A total of 385 known miRNAs from 59 miRNA families and 94 novel miRNAs were identified. Degradome analysis led to the identification of 32 target genes, of which, the squamosa promoter-binding-like protein 12 (SPL12) was identified as a target of osa-miR156a_L + 1. The lncRNA-miRNA-mRNA regulatory network consisted of osa-miR156a_L + 1, two transcripts of SPL12 (LOC_Os06g49010.3 and LOC_Os06g49010.5), and 13 lncRNAs (e.g., MSTRG.244.1 and MSTRG.16577.1).

CONCLUSION

Large-scale expression changes in coding and noncoding RNA were observed in rice mainly due to its response to glyphosate. SPL12, osa-miR156, and lncRNAs (e.g., MSTRG.244.1 and MSTRG.16577.1) could be a novel ceRNA mechanism in response to glyphosate in rice by regulating transcription and metal ions binding. These findings provide a theoretical basis for breeding glyphosate-tolerant rice varieties and for further research on the biogenesis of glyphosate- tolerance in rice.

Cell lineage-specific transcriptome analysis for interpreting cell fate specification of proembryos

In Arabidopsis, a zygote undergoes asymmetrical cell division that establishes the first two distinct cell types of early proembryos, apical and basal cells. However, the genome-wide transcriptional activities that guide divergence of apical and basal cell development remain unknown. Here, we present a comprehensive transcriptome analysis of apical and basal cell lineages, uncovering distinct molecular pathways during cell lineage specification. Selective deletion of inherited transcripts and specific de novo transcription contribute to the establishment of cell lineage-specific pathways for cell fate specification. Embryo-related pathways have been specifically activated in apical cell lineage since 1-cell embryo stage, but quick transcriptome remodeling toward suspensor-specific pathways are found in basal cell lineage. Furthermore, long noncoding RNAs and alternative splicing isoforms may be involved in cell lineage specification. This work also provides a valuable lineage-specific transcriptome resource to elucidate the molecular pathways for divergence of apical and basal cell lineages at genome-wide scale.

Asymmetric division of the Arabidopsis zygote produces apical and basal cells that mainly develop into embryo and suspensor, respectively. Here, Zhou et al. show that de novo transcription and selective RNA turnover establish distinct apical and basal transcriptomes as early as the 1-cell stage of embryo development.

Translational Landscape of Protein-Coding and Non-Protein-Coding RNAs upon Light Exposure in Arabidopsis

Light is one of the most essential environmental clues for plant growth and morphogenesis. Exposure to blue monochromatic light from darkness is a turning point for plant biological activity, and as a result dramatic changes in gene expression occur. To understand the translational impacts of blue light, we have performed ribosome profiling analysis and called translated open reading frames (ORFs) de novo within not only mRNAs but also non-coding RNAs (ncRNAs). Translation efficiency of 3,823 protein-coding ORFs, such as nuclear chloroplast-related genes, was up-regulated by blue light exposure. Moreover, the translational activation of the microRNA biogenesis-related genes, DCL1 and HYL1, was induced by blue light. Considering the 3-nucleotide codon periodicity of ribosome footprints, a few hundred short ORFs lying on ncRNAs and upstream ORFs (uORFs) on mRNAs were found that had differential translation status between blue light and dark. uORFs are known to have a negative effect on the expression of the main ORFs (mORFs) on the same mRNAs. Our analysis suggests that the translation of uORFs is likely to be more stimulated than that of the corresponding mORFs, and uORF-mediated translational repression of the mORFs in five genes was alleviated by blue light exposure. With data-based annotation of the ORFs, our analysis provides insights into the translatome in response to environmental changes, such as those involving light.

Exogenously overexpressed intronic long noncoding RNAs activate host gene expression by affecting histone modification in Arabidopsis

Involvement of long non-coding RNAs (lncRNAs) in the regulation of gene expression in cis has been well studied in eukaryotes but relatively little is known whether and how lncRNAs affect gene expression in tans. In Arabidopsis thaliana, COLDAIR, a previously reported lncRNA, is produced from the first intron of FLOWERING LOCUS C (FLC), which encodes a repressor of flowering time. Our results indicated that the exogenously overexpressed COLDAIR enhances the expression of FLC in trans, resulting in a late-flowering phenotype. In 35S-COLDAIR lines, the enhanced expression of FLC is correlated with the down-regulation of the repressive histone mark H3K27me3 and with the up-regulation of the active histone mark H3K4me3 at the FLC chromatin. Furthermore, we demonstrated that overexpression of intronic lncRNAs from several other H3K27me3-enriched MADS-box genes also activates the expression of their host genes. This study suggests that the involvement of overexpressed intronic lncRNAs in gene activation may be conserved in H3K27me3-enriched genes in eukaryotes.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering

BACKGROUND

Drought stress is a major abiotic factor that affects rapeseed (Brassica napus L.) productivity. Though previous studies indicated that long non-coding RNAs (lncRNAs) play a key role in response to drought stress, a scheme for genome-wide identification and characterization of lncRNAs' response to drought stress is still lacking, especially in the case of B. napus. In order to further understand the molecular mechanism of the response of B. napus to drought stress, we compared changes in the transcriptome between Q2 (a drought-tolerant genotype) and Qinyou8 (a drought-sensitive genotype) responding drought stress and rehydration treatment at the seedling stage.

RESULTS

A total of 5546 down-regulated and 6997 up-regulated mRNAs were detected in Q2 compared with 7824 and 10,251 in Qinyou8, respectively; 369 down-regulated and 108 up- regulated lncRNAs were detected in Q2 compared with 449 and 257 in Qinyou8, respectively. LncRNA-mRNA interaction network analysis indicated that the co-expression network of Q2 was composed of 145 network nodes and 5175 connections, while the co-expression network of Qinyou8 was composed of 305 network nodes and 22,327 connections. We further identified 34 transcription factors (TFs) corresponding to 126 differentially expressed lncRNAs in Q2, and 45 TFs corresponding to 359 differentially expressed lncRNAs in Qinyou8. Differential expression analysis of lncRNAs indicated that up- and down-regulated mRNAs co-expressed with lncRNAs participated in different metabolic pathways and were involved in different regulatory mechanisms in the two genotypes. Notably, some lncRNAs were co-expressed with BnaC07g44670D, which are associated with plant hormone signal transduction. Additionally, some mRNAs co-located with XLOC_052298, XLOC_094954 and XLOC_012868 were mainly categorized as signal transport and defense/stress response.

CONCLUSIONS

The results of this study increased our understanding of expression characterization of rapeseed lncRNAs in response to drought stress and re-watering, which would be useful to provide a reference for the further study of the function and action mechanisms of lncRNAs under drought stress and re-watering.

PncStress: a manually curated database of experimentally validated stress-responsive non-coding RNAs in plants

Non-coding RNAs (ncRNAs) are recognized as key regulatory molecules in many biological processes. Accumulating evidence indicates that ncRNA-related mechanisms play important roles in plant stress responses. Although abundant plant stress-responsive ncRNAs have been identified, these experimentally validated results have not been gathered into a single public domain archive. Therefore, we established PncStress by curating experimentally validated stress-responsive ncRNAs in plants, including microRNAs, long non-coding RNAs and circular RNAs. The current version of PncStress contains 4227 entries from 114 plants covering 48 biotic and 91 abiotic stresses. For each entry, PncStress has biological information and network visualization. Serving as a manually curated database, PncStress will become a valuable resource in support of plant stress response research.

Long Non-coding RNA in Plants in the Era of Reference Sequences

The discovery of non-coding RNAs (ncRNAs), and the subsequent elucidation of their functional roles, was largely delayed due to the misidentification of non-protein-coding parts of DNA as "junk DNA," which forced ncRNAs into the shadows of their protein-coding counterparts. However, over the past decade, insight into the important regulatory roles of ncRNAs has led to rapid progress in their identification and characterization. Of the different types of ncRNAs, long non-coding RNAs (lncRNAs), has attracted considerable attention due to their mRNA-like structures and gene regulatory functions in plant stress responses. While RNA sequencing has been commonly used for mining lncRNAs, a lack of widespread conservation at the sequence level in addition to relatively low and highly tissue-specific expression patterns challenges high-throughput in silico identification approaches. The complex folding characteristics of lncRNA molecules also complicate target predictions, as the knowledge about the interaction interfaces between lncRNAs and potential targets is insufficient. Progress in characterizing lncRNAs and their targets from different species may hold the key to efficient identification of this class of ncRNAs from transcriptomic and potentially genomic resources. In wheat and barley, two of the most important crops, the knowledge about lncRNAs is very limited. However, recently published high-quality genomes of these crops are considered as promising resources for the identification of not only lncRNAs, but any class of molecules. Considering the increasing demand for food, these resources should be used efficiently to discover molecular mechanisms lying behind development and a/biotic stress responses. As our understanding of lncRNAs expands, interactions among ncRNA classes, as well as interactions with the coding sequences, will likely define novel functional networks that may be modulated for crop improvement.

Lysine, Lysine-Rich, Serine, and Serine-Rich Proteins: Link Between Metabolism, Development, and Abiotic Stress Tolerance and the Role of ncRNAs in Their Regulation

Lysine (Lys) is indispensable nutritionally, and its levels in plants are modulated by both transcriptional and post-transcriptional control during plant ontogeny. Animal glutamate receptor homologs have been detected in plants, which may participate in several plant processes through the Lys catabolic products. Interestingly, a connection between Lys and serotonin metabolism has been established recently in rice. 2-Aminoadipate, a catabolic product of Lys appears to play a critical role between serotonin accumulation and the color of rice endosperm/grain. It has also been shown that expression of some lysine-methylated proteins and genes encoding lysine-methyltransferases (KMTs) are regulated by cadmium even as it is known that Lys biosynthesis and its degradation are modulated by novel mechanisms. Three complex pathways co-exist in plants for serine (Ser) biosynthesis, and the relative preponderance of each pathway in relation to plant development or abiotic stress tolerance are being unfolded slowly. But the phosphorylated pathway of L-Ser biosynthesis (PPSB) appears to play critical roles and is essential in plant metabolism and development. Ser, which participates indirectly in purine and pyrimidine biosynthesis and plays a pivotal role in plant metabolism and signaling. Also, L-Ser has been implicated in plant responses to both biotic and abiotic stresses. A large body of information implicates Lys-rich and serine/arginine-rich (SR) proteins in a very wide array of abiotic stresses. Interestingly, a link exists between Lys-rich K-segment and stress tolerance levels. It is of interest to note that abiotic stresses largely influence the expression patterns of SR proteins and also the alternative splicing (AS) patterns. We have checked if any lncRNAs form a cohort of differentially expressed genes from the publicly available PPSB, sequence read archives of NCBI GenBank. Finally, we discuss the link between Lys and Ser synthesis, catabolism, Lys-proteins, and SR proteins during plant development and their myriad roles in response to abiotic stresses.

Transcriptome and Phytochemical Analyses Provide New Insights Into Long Non-Coding RNAs Modulating Characteristic Secondary Metabolites of Oolong Tea (Camellia sinensis) in Solar-Withering

Oolong tea is a popular and semi-fermented beverage. During the processing of tea leaves, withering is the first indispensable process for improving flavor. However, the roles of long non-coding RNAs (lncRNAs) and the characteristic secondary metabolites during the withering of oolong tea leaves remain unknown. In this study, phytochemical analyses indicated that total polyphenols, flavonoids, catechins, epigallocatechin (EGC), catechin gallate (CG), gallocatechin gallate (GCG), epicatechin gallate (ECG), and epigallocatechin gallate (EGCG) were all less abundant in the solar-withered leaves (SW) than in the fresh leaves (FL) and indoor-withered leaves (IW). In contrast, terpenoid, jasmonic acid (JA), and methyl jasmonate (MeJA) contents were higher in the SW than in the FL and IW. By analyzing the transcriptome data, we detected 32,036 lncRNAs. On the basis of the Kyoto Encyclopedia of Genes and Genomes analysis, the flavonoid metabolic pathway, the terpenoid metabolic pathway, and the JA/MeJA biosynthesis and signal transduction pathway were enriched pathways. Additionally, 63 differentially expressed lncRNAs (DE-lncRNAs) and 23 target genes were identified related to the three pathways. A comparison of the expression profiles of the DE-lncRNAs and their target genes between the SW and IW revealed four up-regulated genes (FLS, CCR, CAD, and HCT), seven up-regulated lncRNAs, four down-regulated genes (4CL, CHI, F3H, and F3'H), and three down-regulated lncRNAs related to flavonoid metabolism; nine up-regulated genes (DXS, CMK, HDS, HDR, AACT, MVK, PMK, GGPPS, and TPS), three up-regulated lncRNAs, and six down-regulated lncRNAs related to terpenoid metabolism; as well as six up-regulated genes (LOX, AOS, AOC, OPR, ACX, and MFP2), four up-regulated lncRNAs, and three down-regulated lncRNAs related to JA/MeJA biosynthesis and signal transduction. These results suggested that the expression of DE-lncRNAs and their targets involved in the three pathways may be related to the low abundance of the total polyphenols, flavonoids, and catechins (EGC, CG, GCG, ECG, and EGCG) and the high abundance of terpenoids in the SW. Moreover, solar irradiation, high JA and MeJA contents, and the endogenous target mimic (eTM)-related regulatory mechanism in the SW were also crucial for increasing the terpenoid levels. These findings provide new insights into the greater contribution of solar-withering to the high-quality flavor of oolong tea compared with the effects of indoor-withering.

Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales)

The limitations of RNA sequencing make it difficult to accurately predict alternative splicing (AS) and alternative polyadenylation (APA) events and long non-coding RNAs (lncRNAs), all of which reveal transcriptomic diversity and the complexity of gene regulation. Gnetum, a genus with ambiguous phylogenetic placement in seed plants, has a distinct stomatal structure and photosynthetic characteristics. In this study, a full-length transcriptome of Gnetum luofuense leaves at different developmental stages was sequenced with the latest PacBio Sequel platform. After correction by short reads generated by Illumina RNA-Seq, 80,496 full-length transcripts were obtained, of which 5269 reads were identified as isoforms of novel genes. Additionally, 1660 lncRNAs and 12,998 AS events were detected. In total, 5647 genes in the G. luofuense leaves had APA featured by at least one poly(A) site. Moreover, 67 and 30 genes from the bHLH gene family, which play an important role in stomatal development and photosynthesis, were identified from the G. luofuense genome and leaf transcripts, respectively. This leaf transcriptome supplements the reference genome of G. luofuense, and the AS events and lncRNAs detected provide valuable resources for future studies of investigating low photosynthetic capacity of Gnetum.

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.

Transcriptomic and evolutionary analysis of the mechanisms by which P. argentatum, a rubber producing perennial, responds to drought

BACKGROUND

Guayule (Parthenium argentatum Gray) is a drought tolerant, rubber producing perennial shrub native to northern Mexico and the US Southwest. Hevea brasiliensis, currently the world's only source of natural rubber, is grown as a monoculture, leaving it vulnerable to both biotic and abiotic stressors. Isolation of rubber from guayule occurs by mechanical harvesting of the entire plant. It has been reported that environmental conditions leading up to harvest have a profound impact on rubber yield. The link between rubber biosynthesis and drought, a common environmental condition in guayule's native habitat, is currently unclear.

RESULTS

We took a transcriptomic and comparative genomic approach to determine how drought impacts rubber biosynthesis in guayule. We compared transcriptional profiles of stem tissue, the location of guayule rubber biosynthesis, collected from field-grown plants subjected to water-deficit (drought) and well-watered (control) conditions. Plants subjected to the imposed drought conditions displayed an increase in production of transcripts associated with defense responses and water homeostasis, and a decrease in transcripts associated with rubber biosynthesis. An evolutionary and comparative analysis of stress-response transcripts suggests that more anciently duplicated transcripts shared among the Asteraceae, rather than recently derived duplicates, are contributing to the drought response observed in guayule. In addition, we identified several deeply conserved long non-coding RNAs (lncRNAs) containing microRNA binding motifs. One lncRNA in particular, with origins at the base of Asteraceae, may be regulating the vegetative to reproductive transition observed in water-stressed guayule by acting as a miRNA sponge for miR166.

CONCLUSIONS

These data represent the first genomic analyses of how guayule responds to drought like conditions in agricultural production settings. We identified an inverse relationship between stress-responsive transcripts and those associated with precursor pathways to rubber biosynthesis suggesting a physiological trade-off between maintaining homeostasis and plant productivity. We also identify a number of regulators of abiotic responses, including transcription factors and lncRNAs, that are strong candidates for future projects aimed at modulating rubber biosynthesis under water-limiting conditions common to guayules' native production environment.

Mechanisms and Functions of Long Non-Coding RNAs at Multiple Regulatory Levels

Long non-coding (lnc) RNAs are non-coding RNAs longer than 200 nt. lncRNAs primarily interact with mRNA, DNA, protein, and miRNA and consequently regulate gene expression at the epigenetic, transcriptional, post-transcriptional, translational, and post-translational levels in a variety of ways. They play important roles in biological processes such as chromatin remodeling, transcriptional activation, transcriptional interference, RNA processing, and mRNA translation. lncRNAs have important functions in plant growth and development; biotic and abiotic stress responses; and in regulation of cell differentiation, the cell cycle, and the occurrence of many diseases in humans and animals. In this review, we summarize the functions and mechanisms of lncRNAs in plants, humans, and animals at different regulatory levels.

MYC2 regulates ARR16, a component of cytokinin signaling pathways, in Arabidopsis seedling development

MYC2 is a basic helix-loop-helix transcription factor that acts as a repressor of blue light-mediated photomorphogenic growth; however, it promotes lateral root formation. MYC2 also regulates different phytohormone-signaling pathways in crucial manner. Arabidopsis response regulator 16 (ARR16) is a negative regulator of cytokinin signaling pathways. Here, we show that MYC2 directly binds to the E-box of ARR16 minimal promoter and negatively regulates its expression in a cytokinin-dependent manner. While ARR16 and MYC2 influence jasmonic acid and cytokinin signaling, the expression of ARR16 is regulated by cry1, GBF1, and HYH, the components of light signaling pathways. The transgenic studies show that the expression of ARR16 is regulated by MYC2 at various stages of development. The mutational studies reveal that ARR16 positively regulates the hypocotyl growth in blue light, and phenotypic analysis of atmyc2 arr16 double mutant further reveals that arr16 can suppress the short hypocotyl phenotype of atmyc2. Altogether, this work highlights MYC2-mediated transcriptional repression of ARR16 in Arabidopsis seedling development.

Systematic identification of long noncoding RNAs expressed during light-induced anthocyanin accumulation in apple fruit

Anthocyanin pigments contribute to the red color of apple (Malus × domestica) fruit and have a major influence on their ornamental, dietary and market value. In this study, we investigated the potential role of long noncoding RNAs (lncRNAs) in anthocyanin biosynthesis. RNA-seq analysis of apple peels from the 'Red Fuji' cultivar during light-induced rapid anthocyanin accumulation revealed 5297 putative lncRNAs. Differential expression analysis further showed that lncRNAs were induced during light treatment and were involved in photosynthesis. Using the miRNA-lncRNA-mRNA network and endogenous target mimic (eTM) analysis, we predicted that two differentially expressed lncRNAs, MLNC3.2 and MLNC4.6, were potential eTMs for miRNA156a and promoted the expression of the SPL2-like and SPL33 transcription factors. Transient expression in apple fruit and stable transformation of apple callus showed that overexpression of the eTMs and SPLs promoted anthocyanin accumulation, with the opposite results in eTM and SPL-silenced fruit. Silencing or overexpressing of miR156a also affected the expression of the identified eTMs and SPLs. These results indicated that MLNC3.2 and MLNC4.6 function as eTMs for miR156a and prevent cleavage of SPL2-like and SPL33 by miR156a during light-induced anthocyanin biosynthesis. Our study provides fundamental insights into lncRNA involvement in the anthocyanin biosynthetic pathway in apple fruit.

The long non-coding RNA lncRNA973 is involved in cotton response to salt stress

BACKGROUND

Long non-coding (lnc) RNAs are a class of functional RNA molecules greater than 200 nucleotides in length, and lncRNAs play important roles in various biological regulatory processes and response to the biotic and abiotic stresses. LncRNAs associated with salt stress in cotton have been identified through RNA sequencing, but the function of lncRNAs has not been reported. We previously identified salt stress-related lncRNAs in cotton (Gossypium spp.), and discovered the salt-related lncRNA-lncRNA973.

RESULTS

In this study, we identified the expression level, localization, function, and preliminary mechanism of action of lncRNA973. LncRNA973, which was localized in the nucleus, was expressed at a low level under nonstress conditions but can be significantly increased by salt treatments. Here lncRNA973 was transformed into Arabidopsis and overexpressed. Along with the increased expression compared with wild type under salt stress conditions in transgenic plants, the seed germination rate, fresh weights and root lengths of the transgenic plants increased. We also knocked down the expression of lncRNA973 using virus-induced gene silencing technology. The lncRNA973 knockdown plants wilted, and the leaves became yellowed and dropped under salt-stress conditions, indicating that the tolerance to salt stress had decreased compared with wild type. LncRNA973 may be involved in the regulation of reactive oxygen species-scavenging genes, transcription factors and genes involved in salt stress-related processes in response to cotton salt stress.

CONCLUSIONS

LncRNA973 was localized in the nucleus and its expression was increased by salt treatment. The lncRNA973-overexpression lines had increased salt tolerance compared with the wild type, while the lncRNA973 knockdown plants had reduced salt tolerance. LncRNA973 regulated cotton responses to salt stress by modulating the expression of a series of salt stress-related genes. The data provides a basis for further studies on the mechanisms of lncRNA973-associated responses to salt stress in cotton.

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.

Genome-Wide Identification of Long Noncoding RNAs and Their Responses to Salt Stress in Two Closely Related Poplars

Long noncoding RNAs (lncRNAs) are involved in various biological regulatory processes, but their roles in plants resistance to salt stress remain largely unknown. To systematically explore the characteristics of lncRNAs and their roles in plant salt responses, we conducted strand-specific RNA-sequencing of four tissue types with salt treatments in two closely related poplars (Populus euphratica and Populus alba var. pyramidalis), and a total of 10,646 and 10,531 lncRNAs were identified, respectively. These lncRNAs showed significantly lower values in terms of length, expression, and expression correction than with mRNA. We further found that about 40% and 60% of these identified lncRNAs responded to salt stress with tissue-specific expression patterns across the two poplars. Furthermore, lncRNAs showed weak evolutionary conservation in sequences and exhibited diverse regulatory styles; in particular, tissue- and species-specific responses to salt stress varied greatly in two poplars, for example, 322 lncRNAs were found highly expressed in P. euphratica but not in P. alba var. pyramidalis and 3,425 lncRNAs were identified to be species-specific in P. euphratica in response to salt stress. Moreover, tissue-specific expression of lncRNAs in two poplars were identified with predicted target genes included Aux/IAA, NAC, MYB, involved in regulating plant growth and the plant stress response. Taken together, the systematic analysis of lncRNAs between sister species enhances our understanding of the characteristics of lncRNAs and their roles in plant growth and salt response.

Long non-coding RNAs: Fine-tuning the developmental responses in plants

Plant developmental biology is associated with various gene regulatory pathways involved in different phases of their life cycle. In course of development, growth and differentiation of different organs in plants are regulated by specific sets of gene expression. With the advances in genomic and bioinformatic techniques, particularly high-throughput sequencing technology, many transcriptional units with no protein-coding potential have been discovered. Previously thought to be the dark matters of genome, long non-coding RNAs (lncRNAs) are gradually gaining importance as crucial players in gene regulation during different developmental phases. Some lncRNAs, showing complementarity to microRNAs (miRNAs), are used as endogenous target mimics of specific miRNA family. A number of lncRNAs can also act as natural antisense transcripts to attenuate the expression of coding genes. Although lncRNA-mediated regulations have extensively been studied in animals, plant lncRNA research is still in its initial phase. The present review highlights the regulatory mechanism and different physiological aspects of lncRNAs in plant development. In plants, lncRNAs are found to be associated with a number of plant developmental functions such as lateral root development, vernalization, photomorphogenesis, pollen development, fiber development and nodulation. Understanding these potent roles of lncRNAs in plant development can further provide novel tools for crop improvement programs in future.

Plant Noncoding RNAs: Hidden Players in Development and Stress Responses

A large and significant portion of eukaryotic transcriptomes consists of noncoding RNAs (ncRNAs) that have minimal or no protein-coding capacity but are functional. Diverse ncRNAs, including both small RNAs and long ncRNAs (lncRNAs), play essential regulatory roles in almost all biological processes by modulating gene expression at the transcriptional and posttranscriptional levels. In this review, we summarize the current knowledge of plant small RNAs and lncRNAs, with a focus on their biogenesis, modes of action, local and systemic movement, and functions at the nexus of plant development and environmental responses. The complex connections among small RNAs, lncRNAs, and small peptides in plants are also discussed, along with the challenges of identifying and investigating new classes of ncRNAs.

Subcellular Localization and Functions of Plant lncRNAs in Drought and Salt Stress Tolerance

Plant lncRNAs are expected to play important roles in various plant processes including response to abiotic stress, although the functions of most plant lncRNAs are still unknown. To investigate these potential functions, integrated approaches that employ genetic, transgenic, and molecular biology methods are required. Here, we describe general methods to study the function of lncRNAs in plant response to drought and salt stresses. First, the expression pattern and subcellular localization of lncRNA are analyzed by GUS staining and RNA fluorescence in situ hybridization (FISH). Then the responses of lncRNA mutant and overexpressing transgenic plants to drought and salt stress treatments are characterized, and their sensitivities to ABA are also assayed. To understand the molecular mechanism of lncRNAs' function in stress response, transcriptome sequencing (RNA-seq) and real-time quantitative PCR are performed to analyze altered expression of stress-related genes. Finally, proline content and ROS content are measured to reveal the accumulation of osmolytes and second messengers in these plants in response to drought and salt stress treatments.

Non-coding RNAs and transposable elements in plant genomes: emergence, regulatory mechanisms and roles in plant development and stress responses

This review will provide evidence for the indispensable function of these elements in regulating plant development and resistance to biotic and abiotic stresses, as well as their evolutionary role in facilitating plant adaptation. Over millions of years of evolution, plant genomes have acquired a complex constitution. Plant genomes consist not only of protein coding sequences, but also contain large proportions of non-coding sequences. These include introns of protein-coding genes, and intergenic sequences such as non-coding RNA, repeat sequences and transposable elements. These non-coding sequences help to regulate gene expression, and are increasingly being recognized as playing an important role in genome organization and function. In this review, we summarize the known molecular mechanisms by which gene expression is regulated by several species of non-coding RNAs (microRNAs, long non-coding RNAs, and circular RNAs) and by transposable elements. We further discuss how these non-coding RNAs and transposable elements evolve and emerge in the genome, and the potential influence and importance of these non-coding RNAs and transposable elements in plant development and in stress responses.

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.

Identification and functional prediction of cold-related long non-coding RNA (lncRNA) in grapevine

Plant long non-coding RNA (lncRNA) undergoes dynamic regulation and acts in developmental and stress regulation. In this study, we surveyed the expression dynamics of lncRNAs in grapevine (Vitis vinifera L.) under cold stress using high-throughput sequencing. Two-hundred and three known lncRNAs were significantly up-regulated and 144 known lncRNAs were significantly down-regulated in cold-treated grapevine. In addition, 2 088 novel lncRNA transcripts were identified in this study, with 284 novel lncRNAs significantly up-regulated and 182 novel lncRNAs significantly down-regulated in cold-treated grapevine. Two-hundred and forty-two differentially expressed grapevine lncRNAs were predicted to target 326 protein-coding genes in a cis-regulatory relationship. Many differentially expressed grapevine lncRNAs targeted stress response-related genes, such as CBF4 transcription factor genes, late embryogenesis abundant protein genes, peroxisome biogenesis protein genes, and WRKY transcription factor genes. Sixty-two differentially expressed grapevine lncRNAs were predicted to target 100 protein-coding genes in a trans-regulatory relationship. The expression of overall target genes in both cis and trans-regulatory relationships were positively related to the expression of lncRNAs in grapevines under cold stress. We identified 31 known lncRNAs as 34 grapevine micro RNA (miRNA) precursors and some miRNAs may be derived from multiple lncRNAs. We found 212 lncRNAs acting as targets of miRNAs in grapevines, involving 150 miRNAs; additionally, 120 grapevine genes were predicted as targets of grapevine miRNAs and lncRNAs. We found one gene cluster that was up-regulated and showed the same expression trend. In this cluster, many genes may be involved in abiotic stress response such as WRKY, Hsf, and NAC transcription factor genes.

Seasonal regulation of the lncRNA LDAIR modulates self-protective behaviours during the breeding season

To cope with seasonal environmental changes, animals adapt their physiology and behaviour in response to photoperiod. However, the molecular mechanisms underlying these adaptive changes are not completely understood. Here, using genome-wide expression analysis, we show that an uncharacterized long noncoding RNA (lncRNA), LDAIR, is strongly regulated by photoperiod in Japanese medaka fish (Oryzias latipes). Numerous transcripts and signalling pathways are activated during the transition from short- to long-day conditions; however, LDAIR is one of the first genes to be induced and its expression shows a robust daily rhythm under long-day conditions. Transcriptome analysis of LDAIR knockout fish reveals that the LDAIR locus regulates a gene neighbourhood, including corticotropin releasing hormone receptor 2, which is involved in the stress response. Behavioural analysis of LDAIR knockout fish demonstrates that LDAIR affects self-protective behaviours under long-day conditions. Therefore, we propose that photoperiodic regulation of corticotropin releasing hormone receptor 2 by LDAIR modulates adaptive behaviours to seasonal environmental changes.

Plant lncRNAs are enriched in and move systemically through the phloem in response to phosphate deficiency

In response to phosphate (Pi) deficiency, it has been shown that micro-RNAs (miRNAs) and mRNAs are transported through the phloem for delivery to sink tissues. Growing evidence also indicates that long non-coding RNAs (lncRNAs) are critical regulators of Pi homeostasis in plants. However, whether lncRNAs are present in and move through the phloem, in response to Pi deficiency, remains to be established. Here, using cucumber as a model plant, we show that lncRNAs are enriched in the phloem translocation stream and respond, systemically, to an imposed Pi-stress. A well-known lncRNA, IPS1, the target mimic (TM) of miRNA399, accumulates to a high level in the phloem, but is not responsive to early Pi deficiency. An additional 24 miRNA TMs were also detected in the phloem translocation stream; among them miRNA171 TMs and miR166 TMs were induced in response to an imposed Pi stress. Grafting studies identified 22 lncRNAs which move systemically into developing leaves and root tips. A CU-rich PTB motif was further identified in these mobile lncRNAs. Our findings revealed that lncRNAs respond to Pi deficiency, non-cell-autonomously, and may act as systemic signaling agents to coordinate early Pi deficiency signaling, at the whole-plant level.

Strand-specific RNA-seq based identification and functional prediction of drought-responsive lncRNAs in cassava

Background

Long noncoding RNAs (lncRNAs) have emerged as playing crucial roles in abiotic stress responsive regulation, however, the mechanism of lncRNAs underlying drought-tolerance remains largely unknown in cassava, an important tropical and sub-tropical root crop of remarkable drought tolerance.

Results

In this study, a total of 833 high-confidence lncRNAs, including 652 intergenic and 181 anti-sense lncRNAs, were identified in cassava leaves and root using strand-specific RNA-seq technology, of which 124 were drought-responsive. Trans-regulatory co-expression network revealed that lncRNAs exhibited tissue-specific expression patterns and they preferred to function differently in distinct tissues: e.g., cell-related metabolism, cell wall, and RNA regulation of transcription in folded leaf (FL); degradation of major carbohydrate (CHO) metabolism, calvin cycle and light reaction, light signaling, and tetrapyrrole synthesis in full expanded leaf (FEL); synthesis of major CHO metabolism, nitrogen-metabolism, photosynthesis, and redox in bottom leaf (BL); and hormone metabolism, secondary metabolism, calcium signaling, and abiotic stress in root (RT). In addition, 27 lncRNA-mRNA pairs referred to cis-acting regulation were identified, and these lncRNAs regulated the expression of their neighboring genes mainly through hormone metabolism, RNA regulation of transcription, and signaling of receptor kinase. Besides, 11 lncRNAs were identified acting as putative target mimics of known miRNAs in cassava. Finally, five drought-responsive lncRNAs and 13 co-expressed genes involved in trans-acting, cis-acting, or target mimic regulation were selected and confirmed by qRT-PCR.

Conclusions

These findings provide a comprehensive view of cassava lncRNAs in response to drought stress, which will enable in-depth functional analysis in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5585-5) contains supplementary material, which is available to authorized users.

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.

Genome-wide identification of long non-coding RNA targets of the tomato MADS box transcription factor RIN and function analysis

BACKGROUND AND AIMS

In recent years, increasing numbers of long non-coding RNAs (lncRNAs) have been identified in humans, animals and plants, and several of them have been shown to play important roles in diverse biological processes. However, little work has been performed on the regulation mechanism of lncRNA biogenesis and expression, especially in plants. Compared with studies of tomato MADS-box transcription factor RIPENING INHIBITOR (RIN) target coding genes, there are few reports on its relationship to non-coding RNAs. The aim of the present study was to identify and explore the specific role of RIN target lncRNAs in tomato fruit development and ripening.

METHODS

lncRNA targets of RIN were identified by chromatin immunoprecipitation sequencing (ChIP-seq) combined with RNA deep sequencing analysis. Six selected lncRNA targets were validated by quantitative real-time PCR, ChIP and electrophoretic mobility shift assays, and we further confirmed differential expression between wild-type and ripening-deficient mutant fruit, and RIN direct binding in the promoter regions. By means of virus-induced gene silencing (VIGS) assays and a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) genome editing strategy, the ripening-related function of a specific target lncRNA (lncRNA2155) was studied.

KEY RESULTS

We identified 187 lncRNAs as direct RIN targets, which exhibited RIN binding sites in their promoters and showed different expression between the wild-type and rin mutant. Six target lncRNAs were shown to bind with RIN directly in their promoters in vivo and in vitro. Moreover, using CRISPR/Cas9 technology to knock out the locus of the target lncRNA2155 indicated that it delayed fruit ripening in tomato.

CONCLUSIONS

Collectively, these findings provide new insight into RIN in the transcriptional regulation of lncRNAs and suggest that lncRNAs will contribute to a better understanding of the RIN regulatory network that controls fruit ripening.

Pol III-Dependent Cabbage BoNR8 Long ncRNA Affects Seed Germination and Growth in Arabidopsis

Long non-coding RNAs (lncRNAs) are non-protein-coding transcripts longer than 200 nt that are distributed widely in organisms and play many physiological roles. The BoNR8 lncRNA is a 272 nt long transcript yielded by RNA polymerase III in cabbage that was identified as the closest homolog of the AtR8 lncRNA in Arabidopsis. The BoNR8 lncRNA was expressed extensively in the epidermal tissue in the root elongation zone of germinated seeds, and its accumulation was induced by abiotic stresses, auxins and ABA. To investigate the correlation between the BoNR8 lncRNA and germination, BoNR8-overexpressing Arabidopsis plants (BoNR8-AtOX) were prepared. Three independent BoNR8-AtOX lines showed less primary root elongation, incomplete silique development and decreased germination rates. The germination efficiencies were affected strongly by ABA and slightly by salt stress, and ABA-related gene expression was changed in the BoNR8-AtOX lines.

Differential co-expression networks of long non-coding RNAs and mRNAs in Cleistogenes songorica under water stress and during recovery

BACKGROUND

Water stress seriously constrains plant growth and yield. Long non-coding RNAs (lncRNAs) serve as versatile regulators in various biological regulatory processes. To date, the systematic screening and potential functions of lncRNA have not yet been characterized in Cleistogenes songorica, especially under water stress conditions.

RESULTS

In this study, we obtained the root and shoot transcriptomes of young C. songorica plants subjected to different degrees of water stress and recovery treatments by Illumina-based RNA-seq. A total of 3397 lncRNAs were identified through bioinformatics analysis. LncRNA differential expression analysis indicated that the higher response of roots compared to shoots during water stress and recovery. We further identified the 1644 transcription factors, 189 of which were corresponded to 163 lncRNAs in C. songorica. Though comparative analyses with major Poaceae species based on blast, 81 water stress-related orthologues regulated to lncRNAs were identified as a core of evolutionary conserved genes important to regulate water stress responses in the family. Among these target genes, two genes were found to be involved in the abscisic acid (ABA) signalling pathway, and four genes were enriched for starch and sucrose metabolism. Additionally, the 52 lncRNAs were predicted as target mimics for microRNAs (miRNAs) in C. songorica. RT-qPCR results suggested that MSTRG.43964.1 and MSTRG.4400.2 may regulate the expression of miRNA397 and miRNA166, respectively, as target mimics under water stress and during recovery. Finally, a co-expression network was constructed based on the lncRNAs, miRNAs, protein-coding genes (PCgenes) and transcription factors under water stress and during recovery in C. songorica.

CONCLUSIONS

In C. songorica, lncRNAs, miRNAs, PCgenes and transcription factors constitute a complex transcriptional regulatory network which lncRNAs can regulate PCgenes and miRNAs under water stress and recovery. This study provides fundamental resources to deepen our knowledge on lncRNAs during ubiquitous water stress.

Target of Rapamycin (TOR) Regulates the Expression of lncRNAs in Response to Abiotic Stresses in Cotton

TOR (Target of Rapamycin) kinase is an evolutionarily conserved protein kinase, which integrates stress-related cues with growth and metabolic outputs. Long non-coding RNAs (lncRNAs) play a vital role in the regulation of eukaryotic genes. However, little is known about TOR's function in regulating the expression of lncRNAs in plants. In this study, four putative homologous genes encoding the TOR protein were identified by utilizing the recently completed cotton genome. Pharmacological experiments with TOR inhibitor AZD8055 and on silencing GhTOR genes resulted in obvious cotton growth retardation, indicating the conserved role of TOR in plant growth. The expression pattern analyses in different tissues reveal that TOR may play a role in root development, and the transcript levels of TOR genes were changed under different stress conditions. Importantly, we found TOR may be a key player in regulating the expression of long non-coding RNAs (lncRNAs). A total of 10,315 lncRNAs were discovered in cotton seedlings, 90.7% of which were long intergenic ncRNAs. Moreover, we identified the differentially expressed lncRNAs, of which 296 were significantly upregulated and 105 were downregulated in TOR inactivated plants. GO and KEGG analyses of differentially expressed lncRNA neighboring genes reveal that these differentially expressed lncRNA-targeted genes are involved in many life processes, including stress response, glutathione, and ribosomes in cotton. A series of differentially expressed lncRNAs potentially involved in plant stress response was identified under TOR inhibition. Collectively, these results suggest that cotton TOR proteins may directly modulate the expression of putative stress-related lncRNAs and eventually play a potential role in the cotton stress response.

Using Protein Interaction Profile Sequencing (PIP-seq) to Identify RNA Secondary Structure and RNA-Protein Interaction Sites of Long Noncoding RNAs in Plants

From the moment of transcription, RNA molecules are continuously bound by RNA-binding proteins (RBPs). While the majority of research has focused on how these RBPs regulate posttranscriptional gene regulation of messenger RNAs (mRNAs), the majority of cellular RNAs do not code for proteins, such as ribosomal RNAs, transfer RNAs, and microRNAs. Since these RNAs do not code for protein, their function is mainly determined by their interactions with RBPs as well as their intramolecular base pairing, or RNA secondary structure. One class of noncoding RNAs termed long noncoding RNAs (lncRNAs) have recently become the subject of intense research interest. To study the function of lncRNAs in eukaryotic cells, it is important to examine both their interactions with RBPs as well as their RNA secondary structure. Protein interaction profile sequencing (PIP-seq) is a genome-wide method that uses structure-specific ribonucleases (RNases) to identify regions of double-stranded and single-stranded RNA as well as regions that are protected from these RNases, which represent sites of RBP binding. This method is a very powerful way to examine RNA-protein interactions and RNA secondary structure of all lncRNAs expressed in cells and tissues and can reveal cell-type-specific or tissue-specific patterns of these RNA features. Here, we give a step-by-step account of performing this technique for comprehensively analyzing RNA-protein interactions and RNA secondary structure in plant transcriptomes.