Functions of plants long non-coding RNAs

RNA-protein interactions reveals the pivotal role of lncRNA1840 in tomato fruit maturation

Long non-coding RNAs (lncRNAs) play crucial roles in various biological processes in plants. However, the functional mechanism of lncRNAs in fruit ripening, particularly the transition from unripe to ripe stages, remains elusive. One such lncRNA1840, reported by our group, was found to have important role in tomato fruit ripening. In the present study, we gain insight into its functional role in fruit ripening. CRISPR-Cas9 mediated lncRNA1840 mutants caused the delayed tomato fruit ripening. Notably, loss function of lncRNA1840 did not directly impact ethylene signaling but rather delay ethylene synthesis. Transcriptomic analysis revealed differences in the expression of ripening related genes in lncRNA1840 mutants, suggesting that it is involved in gene regulation of fruit ripening. We used Chromatin Isolation by RNA Purification (ChIRP)-Seq to identify lncRNA1840 binding sites on chromatin. ChIRP-seq suggested that lncRNA1840 had occupancy on 40 genes, but none of them is differentially expressed genes in transcriptomic analysis, which indicated lncRNA1840 might indirectly modulate the gene expression. ChIRP-mass spectrometry analysis identified potential protein interactors of lncRNA1840, Pre-mRNA processing splicing factor 8, highlighting its involvement in post-transcriptional regulatory pathways. In summary, lncRNA1840 is key player in tomato plant growth and fruit ripening, with multifaceted roles in gene expression and regulatory networks.

Exploration of GmDof11-lncRNA13082 Module Regulating Oil Synthesis in Plants

Soybean (Glycine max [Linn.] Merr.) is an important oilseed crop. Although transcription factors (TFs) can coordinate the expression of mRNA and lncRNA, their coordination in the soybean oil synthesis pathway remains unclear. This study examined the interaction between the TF GmDof11 and lncRNA13082 and found that overexpression of GmDof11 led to an increase in the number of Arabidopsis seeds, thousand seed weight, crude protein, hydrolysis amino acid, and soluble sugar. Additionally, it reduced the triglyceride and starch contents and affected the proportion of fatty acids, increasing the contents of palmitic acid, stearic acid, and linolenic acid. The yeast two-hybrid experiments revealed that GmDof11 interacts with GmBCCP1, GmLEC1b, and GmFAB2 proteins. In the RT-qPCR analysis of transgenic soybean roots, it was found that GmDof11 can activate the production of lncRNA13082 and work in conjunction with lncRNA13082 to oversee oil synthesis and nutrient storage. Our research provides robust theoretical evidence for a comprehensive resolution of TF-lncRNA regulation in the soybean oil synthesis network.

Transcription factors KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4 regulate long intergenic noncoding RNAs expressed in Arabidopsis roots

Thousands of long intergenic noncoding RNAs (lincRNAs) have been identified in plant genomes. While some lincRNAs have been characterized as important regulators in different biological processes, little is known about the transcriptional regulation for most plant lincRNAs. Through the integration of 8 annotation resources, we defined 6,599 high-confidence lincRNA loci in Arabidopsis (Arabidopsis thaliana). For lincRNAs belonging to different evolutionary age categories, we identified major differences in sequence and chromatin features, as well as in the level of conservation and purifying selection acting during evolution. Spatiotemporal gene expression profiles combined with transcription factor (TF) chromatin immunoprecipitation (ChIP) data were used to construct a TF-lincRNA regulatory network containing 2,659 lincRNAs and 15,686 interactions. We found that properties characterizing lincRNA expression, conservation, and regulation differ between plants and animals. Experimental validation confirmed the role of 3 TFs, KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4, as key regulators controlling root-specific lincRNA expression, demonstrating the predictive power of our network. Furthermore, we identified 58 lincRNAs, regulated by these TFs, showing strong root cell type-specific expression or chromatin accessibility, which are linked with genome-wide association studies genetic associations related to root system development and growth. The multilevel genome-wide characterization covering chromatin state information, promoter conservation, and chromatin immunoprecipitation-based TF binding, for all detectable lincRNAs across 769 expression samples, permits rapidly defining the biological context and relevance of Arabidopsis lincRNAs through regulatory networks.

Role of methylation in vernalization and photoperiod pathway: a potential flowering regulator?

Recognized as a pivotal developmental transition, flowering marks the continuation of a plant's life cycle. Vernalization and photoperiod are two major flowering pathways orchestrating numerous florigenic signals. Methylation, including histone, DNA and RNA methylation, is one of the recent foci in plant development. Considerable studies reveal that methylation seems to show an increasing potential regulatory role in plant flowering via altering relevant gene expression without altering the genetic basis. However, little has been reviewed about whether and how methylation acts on vernalization- and photoperiod-induced flowering before and after FLOWERING LOCUS C (FLC) reactivation, what role RNA methylation plays in vernalization- and photoperiod-induced flowering, how methylation participates simultaneously in both vernalization- and photoperiod-induced flowering, the heritability of methylation memory under the vernalization/photoperiod pathway, and whether and how methylation replaces vernalization/photoinduction to regulate flowering. Our review provides insight about the crosstalk among the genetic control of the flowering gene network, methylation (methyltransferases/demethylases) and external signals (cold, light, sRNA and phytohormones) in vernalization and photoperiod pathways. The existing evidence that RNA methylation may play a potential regulatory role in vernalization- and photoperiod-induced flowering has been gathered and represented for the first time. This review speculates about and discusses the possibility of substituting methylation for vernalization and photoinduction to promote flowering. Current evidence is utilized to discuss the possibility of future methylation reagents becoming flowering regulators at the molecular level.

Evaluation of efficacy of non-coding RNA in abiotic stress management of field crops: Current status and future prospective

Abiotic stresses are responsible for the major losses in crop yield all over the world. Stresses generate harmful ROS which can impair cellular processes in plants. Therefore, plants have evolved antioxidant systems in defence against the stress-induced damages. The frequency of occurrence of abiotic stressors has increased several-fold due to the climate change experienced in recent times and projected for the future. This had particularly aggravated the risk of yield losses and threatened global food security. Non-coding RNAs are the part of eukaryotic genome that does not code for any proteins. However, they have been recently found to have a crucial role in the responses of plants to both abiotic and biotic stresses. There are different types of ncRNAs, for example, miRNAs and lncRNAs, which have the potential to regulate the expression of stress-related genes at the levels of transcription, post-transcription, and translation of proteins. The lncRNAs are also able to impart their epigenetic effects on the target genes through the alteration of the status of histone modification and organization of the chromatins. The current review attempts to deliver a comprehensive account of the role of ncRNAs in the regulation of plants' abiotic stress responses through ROS homeostasis. The potential applications ncRNAs in amelioration of abiotic stresses in field crops also have been evaluated.

Regulatory mechanisms of fatty acids biosynthesis in Armeniaca sibirica seed kernel oil at different developmental stages

Background

Armeniaca sibirica seed kernel oil is rich in oleic acid and linoleic acid, thus holding potential value as a source of high-quality edible oils. However, some regulatory factors involved in fatty acids accumulation in A. sibirica seed kernels remain largely elusive. Thus, the aim of this study was to elucidate the regulatory mechanisms underlying fatty acids biosynthesis in A. sibirica developing seed kernels.

Methods

Seed kernels from six plants from a single A. sibirica clone were taken at five different developmental stages (days 30, 41, 52, 63, and 73 after anthesis). Fatty acid composition in seed kernel oil was determined by gas chromatography-mass spectrometry (GC-MS). In addition, transcriptome analysis was conducted using second-generation sequencing (SGS) and single-molecule real-time sequencing (SMRT).

Results

Rapid accumulation of fatty acids occurred throughout the different stages of seed kernels development, with oleic acid and linoleic acid as the main fatty acids. A total of 10,024, 9,803, 6,004, 6,719 and 9,688 unigenes were matched in the Nt, Nr, KOG, GO and KEGG databases, respectively. In the category lipid metabolism, 228 differentially expressed genes (DEGs) were annotated into 13 KEGG pathways. Specific unigenes encoding 12 key enzymes related to fatty acids biosynthesis were determined. Co-expression network analysis identified 11 transcription factors (TFs) and 13 long non-coding RNAs (lncRNAs) which putatively participate in the regulation of fatty acid biosynthesis. This study provides insights into the molecular regulatory mechanisms of fatty acids biosynthesis in A. sibirica developing seed kernels, and enabled the identification of novel candidate factors for future improvement of the production and quality of seed kernel oil by breeding.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

Unravelling the Role of Epigenetic Modifications in Development and Reproduction of Angiosperms: A Critical Appraisal

Epigenetics are the heritable changes in gene expression patterns which occur without altering DNA sequence. These changes are reversible and do not change the sequence of the DNA but can alter the way in which the DNA sequences are read. Epigenetic modifications are induced by DNA methylation, histone modification, and RNA-mediated mechanisms which alter the gene expression, primarily at the transcriptional level. Such alterations do control genome activity through transcriptional silencing of transposable elements thereby contributing toward genome stability. Plants being sessile in nature are highly susceptible to the extremes of changing environmental conditions. This increases the likelihood of epigenetic modifications within the composite network of genes that affect the developmental changes of a plant species. Genetic and epigenetic reprogramming enhances the growth and development, imparts phenotypic plasticity, and also ensures flowering under stress conditions without changing the genotype for several generations. Epigenetic modifications hold an immense significance during the development of male and female gametophytes, fertilization, embryogenesis, fruit formation, and seed germination. In this review, we focus on the mechanism of epigenetic modifications and their dynamic role in maintaining the genomic integrity during plant development and reproduction.

Long Non-Coding RNAs as Emerging Regulators of Pathogen Response in Plants

Long non-coding RNAs (lncRNAs) are transcripts without protein-coding potential that contain more than 200 nucleotides that play important roles in plant survival in response to different stresses. They interact with molecules such as DNA, RNA, and protein, and play roles in the regulation of chromatin remodeling, RNA metabolism, and protein modification activities. These lncRNAs regulate the expression of their downstream targets through epigenetic changes, at the level of transcription and post-transcription. Emerging information from computational biology and functional characterization of some of them has revealed their diverse mechanisms of action and possible roles in biological processes such as flowering time, reproductive organ development, as well as biotic and abiotic stress responses. In this review, we have mainly focused on the role of lncRNAs in biotic stress response due to the limited availability of knowledge in this domain. We have discussed the available molecular mechanisms of certain known lncRNAs against specific pathogens. Further, considering that fungal, viral, and bacterial diseases are major factors in the global food crisis, we have highlighted the importance of lncRNAs against pathogen responses and the progress in plant research to develop a better understanding of their functions and molecular mechanisms.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Genome-wide analysis uncovers tomato leaf lncRNAs transcriptionally active upon Pseudomonas syringae pv. tomato challenge

Plants rely on (in)direct detection of bacterial pathogens through plasma membrane-localized and intracellular receptor proteins. Surface pattern-recognition receptors (PRRs) participate in the detection of microbe-associated molecular patterns (MAMPs) and are required for the activation of pattern-triggered immunity (PTI). Pathogenic bacteria, such as Pseudomonas syringae pv. tomato (Pst) deploys ~ 30 effector proteins into the plant cell that contribute to pathogenicity. Resistant plants are capable of detecting the presence or activity of effectors and mount another response termed effector-triggered immunity (ETI). In order to investigate the involvement of tomato’s long non-coding RNAs (lncRNAs) in the immune response against Pst, we used RNA-seq data to predict and characterize those that are transcriptionally active in leaves challenged with a large set of treatments. Our prediction strategy was validated by sequence comparison with tomato lncRNAs described in previous works and by an alternative approach (RT-qPCR). Early PTI (30 min), late PTI (6 h) and ETI (6 h) differentially expressed (DE) lncRNAs were identified and used to perform a co-expression analysis including neighboring (± 100 kb) DE protein-coding genes. Some of the described networks could represent key regulatory mechanisms of photosynthesis, PRR abundance at the cell surface and mitigation of oxidative stress, associated to tomato-Pst pathosystem.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Systematic identification and characterization of circular RNAs involved in flag leaf senescence of rice

Circular RNAs (circRNAs) identification, expression profiles, and construction of circRNA-parental gene relationships and circRNA-miRNA-mRNA ceRNA networks indicate that circRNAs are involved in flag leaf senescence of rice. Circular RNAs (circRNAs) are a class of 3'-5' head-to-tail covalently closed non-coding RNAs which have been proved to play important roles in various biological processes. However, no systematic identification of circRNAs associated with leaf senescence in rice has been studied. In this study, a genome-wide high-throughput sequencing analysis was performed using rice flag leaves developing from normal to senescence. Here, a total of 6612 circRNAs were identified, among which, 113 circRNAs were differentially expressed (DE) during the leaf senescence process. Moreover, 4601 (69.59%) circRNAs were derived from the exons or introns of their parental genes, while 2110 (71%) of the parental genes produced only one circRNA. The sequence alignment analysis showed that hundreds of rice circRNAs were conserved among different plant species. Gene Ontology (GO) enrichment analysis revealed that parental genes of DE circRNAs were enriched in many biological processes closely related to leaf senescence. Through weighted gene co-expression network analysis (WGCNA), six continuously down-expressed circRNAs, 18 continuously up-expressed circRNAs and 15 turn-point high-expressed circRNAs were considered to be highly associated with leaf senescence. Additionally, a total of 17 senescence-associated circRNAs were predicted to have parental genes, in which, regulations of three circRNAs to their parental genes were validated by qRT-PCR. The competing endogenous RNA (ceRNA) networks were also constructed. And a total of 11 senescence-associated circRNAs were predicted to act as miRNA sponges to regulate mRNAs, in which, regulation of two circRNAs to eight mRNAs was validated by qRT-PCR. It is discussed that senescence-associated circRNAs were involved in flag leaf senescence probably through mediating their parental genes and ceRNA networks, to participate in several well-studied senescence-associated processes, mainly including the processes of transcription, translation, and posttranslational modification (especially protein glycosylation), oxidation-reduction process, involvement of senescence-associated genes, hormone signaling pathway, proteolysis, and DNA damage repair. This study not only showed the systematic identification of circRNAs involved in leaf senescence of rice, but also laid a foundation for functional research on candidate circRNAs.

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.

Single-Molecule Real-Time Transcript Sequencing of Turnips Unveiling the Complexity of the Turnip Transcriptome

To generate the full-length transcriptome of Xinjiang green and purple turnips, Brassica rapa var. Rapa, using single-molecule real-time (SMRT) sequencing. The samples of two varieties of Brassica rapa var. Rapa at five developmental stages were collected and combined to perform SMRT sequencing. Meanwhile, next generation sequencing was performed to correct SMRT sequencing data. A series of analyses were performed to investigate the transcript structure. Finally, the obtained transcripts were mapped to the genome of Brassica rapa ssp. pekinesis Chiifu to identify potential novel transcripts. For green turnip (F01), a total of 19.54 Gb clean data were obtained from 8 cells. The number of reads of insert (ROI) and full-length non-chimeric (FLNC) reads were 510,137 and 267,666. In addition, 82,640 consensus isoforms were obtained in the isoform sequences clustering, of which 69,480 were high-quality, and 13,160 low-quality sequences were corrected using Illumina RNA seq data. For purple turnip (F02), there were 20.41 Gb clean data, 552,829 ROIs, and 274,915 FLNC sequences. A total of 93,775 consensus isoforms were obtained, of which 78,798 were high-quality, and the 14,977 low-quality sequences were corrected. Following the removal of redundant sequences, there were 46,516 and 49,429 non-redundant transcripts for F01 and F02, respectively; 7,774 and 9,385 alternative splicing events were predicted for F01 and F02; 63,890 simple sequence repeats, 59,460 complete coding sequences, and 535 long-non coding RNAs were predicted. Moreover, 5,194 and 5,369 novel transcripts were identified by mapping to Brassica rapa ssp. pekinesis Chiifu. The obtained transcriptome data may improve turnip genome annotation and facilitate further study of the Brassica rapa var. Rapa genome and transcriptome.

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.

Intronic heterochromatin prevents cryptic transcription initiation in Arabidopsis

Intronic transposable elements (TEs) comprise a large proportion in eukaryotic genomes, but how they regulate the host genes remains to be explored. Our forward genetic screen disclosed the plant-specific RNA polymerases IV and V in suppressing intronic TE-mediated cryptic transcription initiation of a chimeric transcripts at FLC (FLC^TE ). Initiation of FLC^TE transcription is blocked by the locally formed intronic heterochromatin, which is directly associated with RNA Pol V to inhibit the entry of RNA Pol II and the occupancy of H3K4 methylation. Genome-wide Pol II Ser5p native elongation transcription sequencing revealed that a significant number of intronic heterochromatin-containing genes undergo this mechanism. This study sheds light on deeply understanding the function of intronic heterochromatin on host genes expression in eukaryotic genome.

Genome-Wide Identification of lncRNAs During Rice Seed Development

Rice seed is a pivotal reproductive organ that directly determines yield and quality. Long non-coding RNAs (lncRNAs) have been recognized as key regulators in plant development, but the roles of lncRNAs in rice seed development remain unclear. In this study, we performed a paired-end RNA sequencing in samples of rice pistils and seeds at three and seven days after pollination (DAP) respectively. A total of 540 lncRNAs were obtained, among which 482 lncRNAs had significantly different expression patterns during seed development. Results from semi-qPCR conducted on 15 randomly selected differentially expressed lncRNAs suggested high reliability of the transcriptomic data. RNA interference of TCONS_00023703, which is predominantly transcribed in developing seeds, significantly reduced grain length and thousand-grain weight. These results expanded the dataset of lncRNA in rice and enhanced our understanding of the biological functions of lncRNAs in rice seed development.

The Arabidopsis Hypoxia Inducible AtR8 Long Non-Coding RNA also Contributes to Plant Defense and Root Elongation Coordinating with WRKY Genes under Low Levels of Salicylic Acid

AtR8 lncRNA was previously identified in the flowering plant Arabidopsis thaliana as an abundant Pol III-transcribed long non-coding RNA (lncRNA) of approximately 260 nt. AtR8 lncRNA accumulation is responsive to hypoxic stress and salicylic acid (SA) treatment in roots, but its function has not yet been identified. In this study, microarray analysis of an atr8 mutant and wild-type Arabidopsis indicated a strong association of AtR8 lncRNA with the defense response. AtR8 accumulation exhibited an inverse correlation with an accumulation of two WRKY genes (WRKY53/WRKY70) when plants were exposed to exogenous low SA concentrations (20 µM), infected with Pseudomonas syringae, or in the early stage of development. The highest AtR8 accumulation was observed 5 days after germination, at which time no WRKY53 or WRKY70 mRNA was detectable. The presence of low levels of SA resulted in a significant reduction of root length in atr8 seedlings, whereas wrky53 and wrky70 mutants exhibited the opposite phenotype. Taken together, AtR8 lncRNA participates in Pathogenesis-Related Proteins 1 (PR-1)-independent defense and root elongation, which are related to the SA response. The mutual regulation of AtR8 lncRNA and WRKY53/WRKY70 is mediated by Nonexpressor of Pathogenesis-Related Gene 1 (NPR1).

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.

Screening and characterization of long noncoding RNAs involved in the albinism of Ananas comosus var. bracteatus leaves

Long noncoding RNAs (lncRNAs) have been reported to play key regulatory roles in plant growth, development, and biotic and abiotic stress physiology. Revealing the mechanism of lncRNA regulation in the albino portions of leaves is important for understanding the development of chimeric leaves in Ananas comosus var. bracteatus. In this study, a total of 3,543 candidate lncRNAs were identified, among which 1,451 were differentially expressed between completely green (CGr) and completely white (CWh) leaves. LncRNAs tend to have shorter transcripts, lower expression levels, and greater expression specificity than protein-coding genes. Predicted lncRNA targets were functionally annotated by the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A lncRNA-mRNA interaction network was constructed, and 36 target mRNAs related to chlorophyll metabolism were predicted to interact with 86 lncRNAs. Among these, 25 significantly differentially expressed lncRNAs putatively interacted with 16 target mRNAs. Based on an expression pattern analysis of the lncRNAs and their target mRNAs, the lncRNAs targeting magnesium chelatase subunit H (ChlH), protochlorophyllide oxidoreductase (POR), and heme o synthase (COX10) were suggested as key regulators of chlorophyll metabolism. This study provides the first lncRNA database for A. comosus var. bracteatus and contributes greatly to understanding the mechanism of epigenetic regulation of leaf albinism.

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.

Barley long non-coding RNAs (lncRNA) responsive to excess boron

Long non-coding RNA (lncRNA) has a misleading name, since although they do not encode proteins, they may encode small peptides. Such transcripts are emerging as regulatory molecules. With the advent of next-generation sequencing technologies and novel bioinformatics tools, a tremendous amount of lncRNAs have been identified in several plant species. Recent reports demonstrated roles of plant lncRNAs such as development and environmental response. Here, we reported a genome-wide discovery of ~8000 barley lncRNAs and measured their expression pattern upon excessive boron (B) treatment. According to the tissue-based comparison, leaves have a greater number of B-responsive differentially expressed lncRNAs than the root. Functional annotation of the coding transcripts, which were co-expressed with lncRNAs, revealed that molecular function of the ion transport, establishment of localization, and response to stimulus significantly enriched only in the leaf. On the other hand, 32 barley endogenous target mimics (eTM) as lncRNAs, which potentially decoy the transcriptional suppression activity of 18 miRNAs, were obtained. Also, six lncRNAs, differentially expressed upon B-treatment, were selected and quantitatively analyzed in both B-sensitive and B-tolerant cultivars treated by excess B-level. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis confirmed the B-responsive expressional changes obtained by RNA sequencing. Notably, some lncRNAs (i.e., TCONS_00045190 and TCONS_00056415) over-expressed only in B-tolerant cultivar upon excess B treatment. Presented data including identification, expression measurement, and functional characterization of barley lncRNAs suggest that B-stress response might also be regulated by lncRNA expression, via cooperative interaction of miRNA-eTM-coding target transcript modules.

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.

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.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

BACKGROUND

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

RESULTS

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

CONCLUSION

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Distribution, Characteristics, and Regulatory Potential of Long Noncoding RNAs in Brown-Rot Fungi

Long noncoding RNAs have been thoroughly studied in plants, animals, and yeasts, where they play important roles as regulators of transcription. Nevertheless, almost nothing is known about their presence and characteristics in filamentous fungi, especially in basidiomycetes. In the present study, we have carried out an exhaustive annotation and characterization of lncRNAs in two lignin degrader basidiomycetes, Coniophora puteana and Serpula lacrymans. We identified 2,712 putative lncRNAs in the former and 2,242 in the latter, mainly originating from intergenic locations of transposon-sparse genomic regions. The lncRNA length, GC content, expression levels, and stability of the secondary structure differ from coding transcripts but are similar in these two species and resemble that of other eukaryotes. Nevertheless, they lack sequence conservation. Also, we found that lncRNAs are transcriptionally regulated in the same proportion as genes when the fungus actively decomposes soil organic matter. Finally, up to 7% of the upstream gene regions of Coniophora puteana and Serpula lacrymans are transcribed and produce lncRNAs. The study of expression trends in these gene-lncRNA pairs uncovered groups with similar and opposite transcriptional profiles which may be the result of cis-transcriptional regulation.

Genome-Wide Identification of Long Non-coding RNAs Responsive to Lasiodiplodia theobromae Infection in Grapevine

Long non-coding RNAs (lncRNAs) refer to a class of RNA molecules that are longer than 200 nucleotides and do not encode proteins. Numerous lncRNAs have recently emerged as important regulators of many biological processes in animals and plants, including responses to environmental stress and pathogens. Botryosphaeria dieback is one of the more severe grapevine trunk diseases worldwide. However, how lncRNAs function during Botryosphaeriaceae infection is largely unknown. We performed high-throughput RNA-sequencing (RNA-seq) of susceptible and more tolerant grapevine cultivars infected with Lasiodiplodia theobromae. Overall, we predicted 1826 novel candidate lncRNAs, including long intergenic non-coding RNAs (lincRNAs) and natural antisense transcripts (lncNATs). The data reveal the functions of a set of lncRNAs that were differentially expressed between the resistant cultivar Merlot and the susceptible cultivar Cabernet Franc. Several lncRNAs were predicted to be precursors for grape microRNAs involved in the L theobromae infection. These results provide new insight into the lncRNAs of grapevine that are involved in the response to L theobromae infection.

Comprehensive identification of the full-length transcripts and alternative splicing related to the secondary metabolism pathways in the tea plant (Camellia sinensis)

Flavonoids, theanine and caffeine are the main secondary metabolites of the tea plant (Camellia sinensis), which account for the tea's unique flavor quality and health benefits. The biosynthesis pathways of these metabolites have been extensively studied at the transcriptional level, but the regulatory mechanisms are still unclear. In this study, to explore the transcriptome diversity and complexity of tea plant, PacBio Iso-Seq and RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression during the leaf development. A total of 1,388,066 reads of insert (ROI) were generated with an average length of 1,762 bp, and more than 54% (755,716) of the ROIs were full-length non-chimeric (FLNC) reads. The Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness was 92.7%. A total of 93,883 non-redundant transcripts were obtained, and 87,395 (93.1%) were new alternatively spliced isoforms. Meanwhile, 7,650 differential expression transcripts (DETs) were identified. A total of 28,980 alternative splicing (AS) events were predicted, including 1,297 differential AS (DAS) events. The transcript isoforms of the key genes involved in the flavonoid, theanine and caffeine biosynthesis pathways were characterized. Additionally, 5,777 fusion transcripts and 9,052 long non-coding RNAs (lncRNAs) were also predicted. Our results revealed that AS potentially plays a crucial role in the regulation of the secondary metabolism of the tea plant. These findings enhanced our understanding of the complexity of the secondary metabolic regulation of tea plants and provided a basis for the subsequent exploration of the regulatory mechanisms of flavonoid, theanine and caffeine biosynthesis in tea plants.

Spatio-Temporal Transcriptional Dynamics of Maize Long Non-Coding RNAs Responsive to Drought Stress

Long non-coding RNAs (lncRNAs) have emerged as important regulators in plant stress response. Here, we report a genome-wide lncRNA transcriptional analysis in response to drought stress using an expanded series of maize samples collected from three distinct tissues spanning four developmental stages. In total, 3488 high-confidence lncRNAs were identified, among which 1535 were characterized as drought responsive. By characterizing the genomic structure and expression pattern, we found that lncRNA structures were less complex than protein-coding genes, showing shorter transcripts and fewer exons. Moreover, drought-responsive lncRNAs exhibited higher tissue- and development-specificity than protein-coding genes. By exploring the temporal expression patterns of drought-responsive lncRNAs at different developmental stages, we discovered that the reproductive stage R1 was the most sensitive growth stage with more lncRNAs showing altered expression upon drought stress. Furthermore, lncRNA target prediction revealed 653 potential lncRNA-messenger RNA (mRNA) pairs, among which 124 pairs function in cis-acting mode and 529 in trans. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to oxidoreductase activity, water binding, and electron carrier activity. Multiple promising targets of drought-responsive lncRNAs were discovered, including the V-ATPase encoding gene, vpp4. These findings extend our knowledge of lncRNAs as important regulators in maize drought response.

Identification and characterization of circular RNAs during the sea buckthorn fruit development

As a rising star of noncoding RNA, circular RNAs (circRNAs) have a covalently closed loop structure, which formed by 3'-5' ligation during splicing. A few circRNAs were identified and thought to be transcriptional noise due to cognitive defect over the past 40 years. Recently, with the development of high-throughput RNA sequencing techniques and specific algorithms for circRNA detection and quantification, plenty of potential circRNAs were identified in many species which play important roles in various biological processes. However, researches on circRNAs in fruit ripening process were lacking. Here, we totally identified 2616 circRNAs in sea buckthorn fruit development process, which uniformly distributed in sea buckthorn chromosome. Among them, 1721 (65.8%) circRNAs were arising from the exons of their host genes, 252 circRNAs were identified as the differentially expressed circRNAs (DEcircRNAs) between three different development stages, and 181 (71.8%) DEcircRNAs had sequence similarity with 235 identified circRNAs from five know plant species. Functional annotation revealed that host genes of DEcircRNAs were predicted to be involved in carotenoid biosynthesis, lipid synthesis and plant hormone signal transduction. Additionally, 53 DEcircRNAs were predicted as the corresponding nine miRNAs sponges in sea buckthorn. Divergent reverse-transcription PCR and RT-qPCR were used for validate the differential expression and back-splicing sites of six DEcircRNAs. These results revealed the role of circRNAs in sea buckthorn fruit ripening process and promoted the noncoding RNA researches in plants.

Long Non-Coding RNAs: Rising Regulators of Plant Reproductive Development

Long non-coding RNAs (lncRNAs) exert a multitude of functions in regulating numerous biological processes. Recent studies have uncovered a growing number of lncRNAs within the plant genome. These molecules show striking tissue-specific expression patterns, suggesting that they exert regulatory functions in the growth and development processes of plants. Plant reproductive development is tightly regulated by both environmental and endogenous factors. As plant reproductive development is a crucial aspect of crop breeding, lncRNAs that modulate reproductive development are now particularly worth regarding. Here, we summarize findings that implicate lncRNAs in the control of plant reproductive development, especially in flowering control. Additionally, we elaborate on the regulation mechanisms of lncRNAs, tools for research on their function and mechanism, and potential directions of future research.

Identification of genes associated with ricinoleic acid accumulation in Hiptage benghalensis via transcriptome analysis

BACKGROUND

Ricinoleic acid is a high-value hydroxy fatty acid with broad industrial applications. Hiptage benghalensis seed oil contains a high amount of ricinoleic acid (~ 80%) and represents an emerging source of this unusual fatty acid. However, the mechanism of ricinoleic acid accumulation in H. benghalensis is yet to be explored at the molecular level, which hampers the exploration of its potential in ricinoleic acid production.

RESULTS

To explore the molecular mechanism of ricinoleic acid biosynthesis and regulation, H. benghalensis seeds were harvested at five developing stages (13, 16, 19, 22, and 25 days after pollination) for lipid analysis. The results revealed that the rapid accumulation of ricinoleic acid occurred at the early-mid-seed development stages (16-22 days after pollination). Subsequently, the gene transcription profiles of the developing seeds were characterized via a comprehensive transcriptome analysis with second-generation sequencing and single-molecule real-time sequencing. Differential expression patterns were identified in 12,555 transcripts, including 71 enzymes in lipid metabolic pathways, 246 putative transcription factors (TFs) and 124 long noncoding RNAs (lncRNAs). Twelve genes involved in diverse lipid metabolism pathways, including fatty acid biosynthesis and modification (hydroxylation), lipid traffic, triacylglycerol assembly, acyl editing and oil-body formation, displayed high expression levels and consistent expression patterns with ricinoleic acid accumulation in the developing seeds, suggesting their primary roles in ricinoleic acid production. Subsequent co-expression network analysis identified 57 TFs and 35 lncRNAs, which are putatively involved in the regulation of ricinoleic acid biosynthesis. The transcriptome data were further validated by analyzing the expression profiles of key enzyme-encoding genes, TFs and lncRNAs with quantitative real-time PCR. Finally, a network of genes associated with ricinoleic acid accumulation in H. benghalensis was established.

CONCLUSIONS

This study was the first step toward the understating of the molecular mechanisms of ricinoleic acid biosynthesis and oil accumulation in H. benghalensis seeds and identified a pool of novel genes regulating ricinoleic acid accumulation. The results set a foundation for developing H. benghalensis into a novel ricinoleic acid feedstock at the transcriptomic level and provided valuable candidate genes for improving ricinoleic acid production in other plants.

Identification of genes associated with ricinoleic acid accumulation in Hiptage benghalensis via transcriptome analysis

BACKGROUND

Ricinoleic acid is a high-value hydroxy fatty acid with broad industrial applications. Hiptage benghalensis seed oil contains a high amount of ricinoleic acid (~ 80%) and represents an emerging source of this unusual fatty acid. However, the mechanism of ricinoleic acid accumulation in H. benghalensis is yet to be explored at the molecular level, which hampers the exploration of its potential in ricinoleic acid production.

RESULTS

To explore the molecular mechanism of ricinoleic acid biosynthesis and regulation, H. benghalensis seeds were harvested at five developing stages (13, 16, 19, 22, and 25 days after pollination) for lipid analysis. The results revealed that the rapid accumulation of ricinoleic acid occurred at the early-mid-seed development stages (16-22 days after pollination). Subsequently, the gene transcription profiles of the developing seeds were characterized via a comprehensive transcriptome analysis with second-generation sequencing and single-molecule real-time sequencing. Differential expression patterns were identified in 12,555 transcripts, including 71 enzymes in lipid metabolic pathways, 246 putative transcription factors (TFs) and 124 long noncoding RNAs (lncRNAs). Twelve genes involved in diverse lipid metabolism pathways, including fatty acid biosynthesis and modification (hydroxylation), lipid traffic, triacylglycerol assembly, acyl editing and oil-body formation, displayed high expression levels and consistent expression patterns with ricinoleic acid accumulation in the developing seeds, suggesting their primary roles in ricinoleic acid production. Subsequent co-expression network analysis identified 57 TFs and 35 lncRNAs, which are putatively involved in the regulation of ricinoleic acid biosynthesis. The transcriptome data were further validated by analyzing the expression profiles of key enzyme-encoding genes, TFs and lncRNAs with quantitative real-time PCR. Finally, a network of genes associated with ricinoleic acid accumulation in H. benghalensis was established.

CONCLUSIONS

This study was the first step toward the understating of the molecular mechanisms of ricinoleic acid biosynthesis and oil accumulation in H. benghalensis seeds and identified a pool of novel genes regulating ricinoleic acid accumulation. The results set a foundation for developing H. benghalensis into a novel ricinoleic acid feedstock at the transcriptomic level and provided valuable candidate genes for improving ricinoleic acid production in other plants.

LncRNA expression profile and ceRNA analysis in tomato during flowering

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that play essential regulatory roles in various developmental processes and stress responses. However, the functions of lncRNAs during the flowering period of tomato are largely unknown. To explore the lncRNA profiles and functions during flowering in tomato, we performed strand-specific paired-end RNA sequencing of tomato leaves, flowers and roots, with three biological replicates. We identified 10919 lncRNAs including 248 novel lncRNAs, of which 65 novel lncRNAs were significantly differentially expressed (DE) in the flowers, leaves, and roots. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out to identify the cis target gene of DE lncRNAs. The results showed that the lncRNAs might play an important role in the growth, development, and apoptosis of flowering tomato plant by regulating the formation of intima in flower tissues, binding to various molecules, influencing metabolic pathways, and inducing apoptosis. Moreover, we identified the interaction between 32, 78, and 397 kinds of miRNAs, lncRNAs, and mRNAs. The results suggest that the lncRNAs can regulate the expression of mRNA during flowering period in tomato by forming competitive endogenous RNA, and further regulate various biological metabolism pathways 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.

Genome-wide identification of oil biosynthesis-related long non-coding RNAs in allopolyploid Brassica napus

Background

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 bp that do not encode proteins but nonetheless have been shown to play important roles in various biological processes in plants. Brassica napus is an important seed oil crop worldwide and the target of many genetic improvement activities. To understand better the function of lncRNAs in regulating plant metabolic activities, we carried out a genome-wide lncRNA identification of lncRNAs in Brassica napus with a focus on lncRNAs involved in lipid metabolism. Twenty ribosomal RNA depleted strand specific RNA-seq (ssRNA-seq) datasets were generatred using RNAs isolated from B. napus seeds at four developmental stages. For comparison we also included 30 publically available RNA-seq datasets generated from poly(A) enriched mRNAs isolated from from various Brassica napus tissues in our analysis.

Results

A total of 8905 lncRNA loci were identified, including 7100 long intergenic noncoding RNA (lincRNA) loci and 1805 loci generating long noncoding natural antisense transcript (lncNAT). Many lncRNAs were identified only in the ssRNA-seq and poly(A) RNA-seq dataset, suggesting that B. napus has a large lncRNA repertoire and it is necessary to use libraries prepared from different tissues and developmental stages as well as different library preparation approaches to capture the whole spectrum of lncRNAs. Analysis of coexpression networks revealed that among the regulatory modules are networks containing lncRNAs and protein-coding genes related to oil biosynthesis indicating a possible role of lncRNAs in the control of lipid metabolism. One such example is that several lncRNAs are potential regulators of BnaC08g11970D that encodes oleosin1, a protein found in oil bodies and involved in seed lipid accumulation. We also observed that the expression levels of B. napus lncRNAs is positively correlated with their conservation levels.

Conclusions

We demonstrated that the B. napus genome has a large number of lncRNA and that these lncRNAs are expressed broadly across many developmental times and in different tissue types. We also provide evidence indicating that specific lncRNAs appear to be important regulators of lipid biosynthesis forming regulatory networks with transcripts involved in lipid biosynthesis. We also provide evidence that these lncRNAs are conserved in other species of the Brassicaceae family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5117-8) contains supplementary material, which is available to authorized users.

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.

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.

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.

Arabidopsis HEAT SHOCK TRANSCRIPTION FACTORA1b regulates multiple developmental genes under benign and stress conditions

In Arabidopsis thaliana, HEAT SHOCK TRANSCRIPTION FACTORA1b (HSFA1b) controls resistance to environmental stress and is a determinant of reproductive fitness by influencing seed yield. To understand how HSFA1b achieves this, we surveyed its genome-wide targets (ChIP-seq) and its impact on the transcriptome (RNA-seq) under non-stress (NS), heat stress (HS) in the wild type, and in HSFA1b-overexpressing plants under NS. A total of 952 differentially expressed HSFA1b-targeted genes were identified, of which at least 85 are development associated and were bound predominantly under NS. A further 1780 genes were differentially expressed but not bound by HSFA1b, of which 281 were classified as having development-associated functions. These genes are indirectly regulated through a hierarchical network of 27 transcription factors (TFs). Furthermore, we identified 480 natural antisense non-coding RNA (cisNAT) genes bound by HSFA1b, defining a further mode of indirect regulation. Finally, HSFA1b-targeted genomic features not only harboured heat shock elements, but also MADS box, LEAFY, and G-Box promoter motifs. This revealed that HSFA1b is one of eight TFs that target a common group of stress defence and developmental genes. We propose that HSFA1b transduces environmental cues to many stress tolerance and developmental genes to allow plants to adjust their growth and development continually in a varying environment.

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.

Photoperiodism dynamics during the domestication and improvement of soybean

Soybean (Glycine max) is a facultative short-day plant with a sensitive photoperiod perception and reaction system, which allows it to adjust its physiological state and gene regulatory networks to seasonal and diurnal changes in environmental conditions. In the past few decades, soybean cultivation has spread from East Asia to areas throughout the world. Biologists and breeders must now confront the challenge of understanding the molecular mechanism of soybean photoperiodism and improving agronomic traits to enable this important crop to adapt to geographical and environmental changes. In this review, we summarize the genetic regulatory network underlying photoperiodic responses in soybean. Genomic and genetic studies have revealed that the circadian clock, in conjunction with the light perception pathways, regulates photoperiodic flowering. Here, we provide an annotated list of 844 candidate flowering genes in soybean, with their putative biological functions. Many photoperiod-related genes have been intensively selected during domestication and crop improvement. Finally, we describe recent progress in engineering photoperiod-responsive genes for improving agronomic traits to enhance geographic adaptation in soybean, as well as future prospects for research on soybean photoperiodic responses.

The Role of Non-Coding RNAs in Cytoplasmic Male Sterility in Flowering Plants

The interactions between mitochondria and nucleus substantially influence plant development, stress response and morphological features. The prominent example of a mitochondrial-nuclear interaction is cytoplasmic male sterility (CMS), when plants produce aborted anthers or inviable pollen. The genes responsible for CMS are located in mitochondrial genome, but their expression is controlled by nuclear genes, called fertility restorers. Recent explosion of high-throughput sequencing methods enabled to study transcriptomic alterations in the level of non-coding RNAs under CMS biogenesis. We summarize current knowledge of the role of nucleus encoded regulatory non-coding RNAs (long non-coding RNA, microRNA as well as small interfering RNA) in CMS. We also focus on the emerging data of non-coding RNAs encoded by mitochondrial genome and their possible involvement in mitochondrial-nuclear interactions and CMS development.

The antiphasic regulatory module comprising CDF5 and its antisense RNA FLORE links the circadian clock to photoperiodic flowering

Circadian rhythms of gene expression are generated by the combinatorial action of transcriptional and translational feedback loops as well as chromatin remodelling events. Recently, long noncoding RNAs (lncRNAs) that are natural antisense transcripts (NATs) to transcripts encoding central oscillator components were proposed as modulators of core clock function in mammals (Per) and fungi (frq/qrf). Although oscillating lncRNAs exist in plants, their functional characterization is at an initial stage. By screening an Arabidopsis thaliana lncRNA custom-made array we identified CDF5 LONG NONCODING RNA (FLORE), a circadian-regulated lncRNA that is a NAT of CDF5. Quantitative real-time RT-PCR confirmed the circadian regulation of FLORE, whereas GUS-staining and flowering time evaluation were used to determine its biological function. FLORE and CDF5 antiphasic expression reflects mutual inhibition in a similar way to frq/qrf. Moreover, whereas the CDF5 protein delays flowering by directly repressing FT transcription, FLORE promotes it by repressing several CDFs (CDF1, CDF3, CDF5) and increasing FT transcript levels, indicating both cis and trans function. We propose that the CDF5/FLORE NAT pair constitutes an additional circadian regulatory module with conserved (mutual inhibition) and unique (function in trans) features, able to fine-tune its own circadian oscillation, and consequently, adjust the onset of flowering to favourable environmental conditions.

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.

Whole transcriptome sequencing of Pseudomonas syringae pv. actinidiae-infected kiwifruit plants reveals species-specific interaction between long non-coding RNA and coding genes

An outbreak of kiwifruit bacterial canker disease caused by Pseudomonas syringae pv. actinidiae (Psa) beginning in 2008 caused disaster to the kiwifruit industry. However the mechanisms of interaction between kiwifruit and Psa are unknown. Long noncoding RNAs (lncRNAs) are known to regulate many biological processes, but comprehensive repertoires of kiwifruit lncRNAs and their effects on the interaction between kiwifruit and Psa are unknown. Here, based on in-depth transcriptomic analysis of four kiwifruit materials at three stages of infection with Psa, we identified 14,845 transcripts from 12,280 loci as putative lncRNAs. Hierarchical clustering analysis of differentially-expressed transcripts reveals that both protein-coding and lncRNA transcripts are expressed species-specifically. Comparing differentially-expressed transcripts from different species, variations in pattern-triggered immunity (PTI) were the main causes of species-specific responses to infection by Psa. Using weighted gene co-expression network analysis, we identified species-specific expressed key lncRNAs which were closely related to plant immune response and signal transduction. Our results illustrate that different kiwifruit species employ multiple different plant immunity layers to fight against Psa infection, which causes distinct responses. We also discovered that lncRNAs might affect kiwifruit responses to Psa infection, indicating that both protein-coding regions and noncoding regions can affect kiwifruit response to Psa infection.

Plant hormone signaling in flowering: An epigenetic point of view

Reproduction is one of the most important phases in an organism's lifecycle. In the case of angiosperm plants, flowering provides the major developmental transition from the vegetative to the reproductive stage, and requires genetic and epigenetic reprogramming to ensure the success of seed production. Flowering is regulated by a complex network of genes that integrate multiple environmental cues and endogenous signals so that flowering occurs at the right time; hormone regulation, signaling and homeostasis are very important in this process. Working alone or in combination, hormones are able to promote flowering by epigenetic regulation. Some plant hormones, such as gibberellins, jasmonic acid, abscisic acid and auxins, have important effects on chromatin compaction mediated by DNA methylation and histone posttranslational modifications, which hints at the role that epigenetic regulation may play in flowering through hormone action. miRNAs have been viewed as acting independently from DNA methylation and histone modification, ignoring their potential to interact with hormone signaling - including the signaling of auxins, gibberellins, ethylene, jasmonic acid, salicylic acid and others - to regulate flowering. Therefore, in this review we examine new findings about interactions between epigenetic mechanisms and key players in hormone signaling to coordinate flowering.