Survey of High Throughput RNA-Seq Data Reveals Potential Roles for lncRNAs during Development and Stress Response in Bread Wheat

Exploring the emerging role of long non-coding RNAs (lncRNAs) in plant biology: Functions, mechanisms of action, and future directions

Long non-coding RNAs (lncRNAs) are a class of RNA transcripts that surpass 200 nucleotides in length and lack discernible coding potential. LncRNAs that have been functionally characterized have pivotal functions in several plant processes, including the regulation of flowering, and development of lateral roots. It also plays a crucial role in the plant's response to abiotic stressors and exhibits vital activities in environmental adaptation. The progress in NGS (next-generation sequencing) and functional genomics technology has facilitated the discovery of lncRNA in plant species. This review is a brief explanation of lncRNA genomics, its molecular role, and the mechanism of action in plants. The review also addresses the challenges encountered in this field and highlights promising molecular and computational methodologies that can aid in the comparative and functional analysis of lncRNAs.

Exploring the Role of TaPLC1-2B in Heat Tolerance at Seedling and Adult Stages of Wheat through Transcriptome Analysis

Heat stress is a major abiotic stress that can cause serious losses of a crop. Our previous work identified a gene involved in heat stress tolerance in wheat, TaPLC1-2B. To further investigate its mechanisms, in the present study, TaPLC1-2B RNAi-silenced transgenic wheat and the wild type were comparatively analyzed at both the seedling and adult stages, with or without heat stress, using transcriptome sequencing. A total of 15,549 differentially expressed genes (DEGs) were identified at the adult stage and 20,535 DEGs were detected at the seedling stage. After heat stress, an enrichment of pathways such as phytohormones and mitogen-activated protein kinase signaling was mainly found in the seedling stage, and pathways related to metabolism, glycerophospholipid metabolism, circadian rhythms, and ABC transporter were enriched in the adult stage. Auxin and abscisic acid were downregulated in the seedling stage and vice versa in the adult stage; and the MYB, WRKY, and no apical meristem gene families were downregulated in the seedling stage in response to heat stress and upregulated in the adult stage in response to heat stress. This study deepens our understanding of the mechanisms of TaPLC1-2B in regard to heat stress in wheat at the seedling and adult stages.

Genome-wide identification of long non-coding RNAs and their potential functions in radish response to salt stress

Long non-coding RNAs (lncRNAs) are increasingly recognized as cis- and trans-acting regulators of protein-coding genes in plants, particularly in response to abiotic stressors. Among these stressors, high soil salinity poses a significant challenge to crop productivity. Radish (Raphanus sativus L.) is a prominent root vegetable crop that exhibits moderate susceptibility to salt stress, particularly during the seedling stage. Nevertheless, the precise regulatory mechanisms through which lncRNAs contribute to salt response in radish remain largely unexplored. In this study, we performed genome-wide identification of lncRNAs using strand-specific RNA sequencing on radish fleshy root samples subjected to varying time points of salinity treatment. A total of 7,709 novel lncRNAs were identified, with 363 of them displaying significant differential expression in response to salt application. Furthermore, through target gene prediction, 5,006 cis- and 5,983 trans-target genes were obtained for the differentially expressed lncRNAs. The predicted target genes of these salt-responsive lncRNAs exhibited strong associations with various plant defense mechanisms, including signal perception and transduction, transcription regulation, ion homeostasis, osmoregulation, reactive oxygen species scavenging, photosynthesis, phytohormone regulation, and kinase activity. Notably, this study represents the first comprehensive genome-wide analysis of salt-responsive lncRNAs in radish, to the best of our knowledge. These findings provide a basis for future functional analysis of lncRNAs implicated in the defense response of radish against high salinity, which will aid in further understanding the regulatory mechanisms underlying radish response to salt stress.

Genome-wide identification and characterization of Puccinia striiformis-responsive lncRNAs in Triticum aestivum

Wheat stripe rust (yellow rust) caused by Puccinia striiformis f. sp. tritici (Pst) is a serious biotic stress factor limiting wheat production worldwide. Emerging evidence demonstrates that long non-coding RNAs (lncRNAs) participate in various developmental processes in plants via post-transcription regulation. In this study, RNA sequencing (RNA-seq) was performed on a pair of near-isogenic lines-rust resistance line FLW29 and rust susceptible line PBW343-which differed only in the rust susceptibility trait. A total of 6,807 lncRNA transcripts were identified using bioinformatics analyses, among which 10 lncRNAs were found to be differentially expressed between resistance and susceptible lines. In order to find the target genes of the identified lncRNAs, their interactions with wheat microRNA (miRNAs) were predicted. A total of 199 lncRNAs showed interactions with 65 miRNAs, which further target 757 distinct mRNA transcripts. Moreover, detailed functional annotations of the target genes were used to identify the candidate genes, pathways, domains, families, and transcription factors that may be related to stripe rust resistance response in wheat plants. The NAC domain protein, disease resistance proteins RPP13 and RPM1, At1g58400, monodehydroascorbate reductase, NBS-LRR-like protein, rust resistance kinase Lr10-like, LRR receptor, serine/threonine-protein kinase, and cysteine proteinase are among the identified targets that are crucial for wheat stripe rust resistance. Semiquantitative PCR analysis of some of the differentially expressed lncRNAs revealed variations in expression profiles of two lncRNAs between the Pst-resistant and Pst-susceptible genotypes at least under one condition. Additionally, simple sequence repeats (SSRs) were also identified from wheat lncRNA sequences, which may be very useful for conducting targeted gene mapping studies of stripe rust resistance in wheat. These findings improved our understanding of the molecular mechanism responsible for the stripe rust disease that can be further utilized to develop wheat varieties with durable resistance to this disease.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Genome-wide identification and characterization of long non-coding RNA in barley roots in response to Piriformospora indica colonization

Currently, there is very limited information about long noncoding RNAs (lncRNAs) found in barley. It remains unclear whether barley lncRNAs are responsive to Piriformospora indica (P. indica) colonization.We found that barley roots exhibited fast development and that large roots branched after P. indica colonization. Genome-wide high-throughput RNA-seq and bioinformatic analysis showed that 4356 and 5154 differentially expressed LncRNAs (DELs) were found in response to P. indica at 3 and 7 days after colonization (dai), respectively, and 2456 DELs were found at 7 dai compared to 3 dai. Based on the coexpression correlation of lncRNAmRNA, we found that 98.6% of lncRNAs were positively correlated with 3430 mRNAs at 3 dai and 7 dai. Further GO analysis showed that 30 lncRNAs might be involved in the regulation of gene transcription; 23 lncRNAs might participate in cell cycle regulation. Moreover, the metabolite analysis indicated that chlorophyll a, sucrose, protein, gibberellin, and auxin were in accordance with the results of the transcriptome, and the respective lncRNAs were positively correlated with these target RNAs. Gene silencing suggested that lncRNA TCONS_00262342 is probably a key regulator of GA₃ synthesis pathway, which participates in P. indica and barley interactions. We concluded that acting as a molecular material basis and resource, lncRNAs respond to P. indica colonization by regulating metabolite content in barley and coordinate the complex regulatory process of higher life by constructing highly positive correlations with their target mRNAs.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Regulatory Networks of lncRNAs, miRNAs, and mRNAs in Response to Heat Stress in Wheat (Triticum Aestivum L.): An Integrated Analysis

Climate change has become a major source of concern, particularly in agriculture, because it has a significant impact on the production of economically important crops such as wheat, rice, and maize. In the present study, an attempt has been made to identify differentially expressed heat stress-responsive long non-coding RNAs (lncRNAs) in the wheat genome using publicly available wheat transcriptome data (24 SRAs) representing two conditions, namely, control and heat-stressed. A total of 10,965 lncRNAs have been identified and, among them, 153, 143, and 211 differentially expressed transcripts have been found under 0 DAT, 1 DAT, and 4 DAT heat-stress conditions, respectively. Target prediction analysis revealed that 4098 lncRNAs were targeted by 119 different miRNA responses to a plethora of environmental stresses, including heat stress. A total of 171 hub genes had 204 SSRs (simple sequence repeats), and a set of target sequences had SNP potential as well. Furthermore, gene ontology analysis revealed that the majority of the discovered lncRNAs are engaged in a variety of cellular and biological processes related to heat stress responses. Furthermore, the modeled three-dimensional (3D) structures of hub genes encoding proteins, which had an appropriate range of similarity with solved structures, provided information on their structural roles. The current study reveals many elements of gene expression regulation in wheat under heat stress, paving the way for the development of improved climate-resilient wheat cultivars.

Identification of birch lncRNAs and mRNAs responding to salt stress and characterization of functions of lncRNA

Long noncoding RNAs (lncRNAs) are important in abiotic stress tolerance. Here, we identified salt-responsive genes and lncRNAs in the roots and leaves of Betula platyphylla Suk. (birch), and characterized their lncRNAs functions. In total, 2660 mRNAs and 539 lncRNAs responding to salt treatment were identified using RNA-seq. The salt-responsive genes were substantially enriched in 'cell wall biogenesis' and 'wood development' in the roots and were enriched in 'photosynthesis' and 'response to stimulus' in the leaves. Meanwhile, the potential target genes of the salt-responsive lncRNAs in roots and leaves were both enriched in 'nitrogen compound metabolic process' and 'response to stimulus'. We further built a method for quickly identifying abiotic stress tolerance of lncRNAs, which employed transient transformation for overexpression and knock-down of the lncRNA, enabling gain- and loss-of-function analysis. Using this method, 11 randomly selected salt-responsive lncRNAs were characterized. Among them, six lncRNAs confer salt tolerance, two lncRNAs confer salt sensitivity, and the other three lncRNAs are not involved in salt tolerance. In addition, a lncRNA, LncY1, was further characterized, which improves salt tolerance by regulating two transcription factors, BpMYB96 and BpCDF3. Taken together, our results suggested that lncRNAs play important roles in the salt response of birch plants.

Wheat Long Noncoding RNAs from Organelle and Nuclear Genomes Carry Conserved microRNA Precursors Which May Together Comprise Intricate Networks in Insect Responses

Long noncoding RNAs (lncRNAs) are a diverse class of noncoding RNAs that are typically longer than 200 nucleotides but lack coding potentials. Advances in deep sequencing technologies enabled a better exploration of this type of noncoding transcripts. The poor sequence conservation, however, complicates the identification and annotation of lncRNAs at a large scale. Wheat is among the leading food staples worldwide whose production is threatened by both biotic and abiotic stressors. Here, we identified putative lncRNAs from durum wheat varieties that differ in stem solidness, a major source of defense against wheat stem sawfly, a devastating insect pest. We also analyzed and annotated lncRNAs from two bread wheat varieties, resistant and susceptible to another destructive pest, orange wheat blossom midge, with and without infestation. Several putative lncRNAs contained potential precursor sequences and/or target regions for microRNAs, another type of regulatory noncoding RNAs, which may indicate functional networks. Interestingly, in contrast to lncRNAs themselves, microRNAs with potential precursors within the lncRNA sequences appeared to be highly conserved at the sequence and family levels. We also observed a few putative lncRNAs that have perfect to near-perfect matches to organellar genomes, supporting the recent observations that organellar genomes may contribute to the noncoding transcript pool of the cell.

Genome-wide identification and characterization of lncRNAs in sunflower endosperm

BACKGROUND

Long non-coding RNAs (lncRNAs), as important regulators, play important roles in plant growth and development. The expression and epigenetic regulation of lncRNAs remain uncharacterized generally in plant seeds, especially in the transient endosperm of the dicotyledons.

RESULTS

In this study, we identified 11,840 candidate lncRNAs in 12 day-after-pollination sunflower endosperm by analyzing RNA-seq data. These lncRNAs were evenly distributed in all chromosomes and had specific features that were distinct from mRNAs including tissue-specificity expression, shorter and fewer exons. By GO analysis of protein coding genes showing strong correlation with the lncRNAs, we revealed that these lncRNAs potential function in many biological processes of seed development. Additionally, genome-wide DNA methylation analyses revealed that the level of DNA methylation at the transcription start sites was negatively correlated with gene expression levels in lncRNAs. Finally, 36 imprinted lncRNAs were identified including 32 maternally expressed lncRNAs and four paternally expressed lncRNAs. In CG and CHG context, DNA methylation levels of imprinted lncRNAs in the upstream and gene body regions were slightly lower in the endosperm than that in embryo tissues, which indicated that the maternal demethylation potentially induce the paternally bias expression of imprinted lncRNAs in sunflower endosperm.

CONCLUSION

Our findings not only identified and characterized lncRNAs on a genome-wide scale in the development of sunflower endosperm, but also provide novel insights into the parental effects and epigenetic regulation of lncRNAs in dicotyledonous seeds.

Genome-wide identification of long intergenic non-coding RNAs for Ralstonia solanacearum resistance in tomato (Solanum lycopersicum)

There is growing evidences indicating that long intergenic ncRNAs (lincRNAs) play key roles in plant development and stress responses. To research tomato lincRNA functions during the interaction between tomato and Ralstonia solanacearum, RNA-seq data of tomato plants inoculated with R. solanacearum was analyzed. In this study, 315 possible lincRNAs were identified from RNA-seq data. Then 23 differentially expressed lincRNAs between tomato plants inoculated with R. solanacearum and control were identified and a total of 171 possible target genes for these differentially expressed lincRNAs were predicted. Through GO and KEGG analysis, we found that lincRNA might be involved in jasmonic acid and ethylene signaling pathways to respond to tomato bacterial wilt infection. Furthermore, lincRNA may also be involved in regulating the expression of AGO protein. Subsequently, analysis of expression patterns between differentially expressed lincRNAs and adjacent mRNAs by qRT-PCR revealed that part of lincRNAs and their possible target genes exhibited positive correlation. Taken together, these results suggest that lincRNAs play potential roles in tomato against R. solanacearum infection and will provide fundamental information about the lincRNA-based plant defense mechanisms.

Revealing the novel complexity of plant long non-coding RNA by strand-specific and whole transcriptome sequencing for evolutionarily representative plant species

Background

Previous studies on plant long noncoding RNAs (lncRNAs) lacked consistency and suffered from many factors like heterogeneous data sources and experimental protocols, different plant tissues, inconsistent bioinformatics pipelines, etc. For example, the sequencing of RNAs with poly(A) tails excluded a large portion of lncRNAs without poly(A), and use of regular RNA-sequencing technique did not distinguish transcripts’ direction for lncRNAs. The current study was designed to systematically discover and analyze lncRNAs across eight evolutionarily representative plant species, using strand-specific (directional) and whole transcriptome sequencing (RiboMinus) technique.

Results

A total of 39,945 lncRNAs (25,350 lincRNAs and 14,595 lncNATs) were identified, which showed molecular features of lncRNAs that are consistent across divergent plant species but different from those of mRNA. Further, transposable elements (TEs) were found to play key roles in the origination of lncRNA, as significantly large number of lncRNAs were found to contain TEs in gene body and promoter region, and transcription of many lncRNAs was driven by TE promoters. The lncRNA sequences were divergent even in closely related species, and most plant lncRNAs were genus/species-specific, amid rapid turnover in evolution. Evaluated with PhastCons scores, plant lncRNAs showed similar conservation level to that of intergenic sequences, suggesting that most lincRNAs were young and with short evolutionary age. INDUCED BY PHOSPHATE STARVATION (IPS) was found so far to be the only plant lncRNA group with conserved motifs, which may play important roles in the adaptation of terrestrial life during migration from aquatic to terrestrial. Most highly and specially expressed lncRNAs formed co-expression network with coding genes, and their functions were believed to be closely related to their co-expression genes.

Conclusion

The study revealed novel features and complexity of lncRNAs in plants through systematic analysis, providing important insights into the origination and evolution of plant lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08602-9.

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

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

Identification and characterization of heat-responsive lncRNAs in maize inbred line CM1

BACKGROUND

Frequent occurrence of extreme high temperature is a major threat to crop production. Increasing evidence demonstrates that long non-coding RNAs (lncRNAs) have important biological functions in the regulation of the response to heat stress. However, the regulatory mechanism of lncRNAs involved in heat response requires further exploration and the regulatory network remains poorly understood in maize.

RESULTS

In this research, high-throughput sequencing was adopted to systematically identify lncRNAs in maize inbred line CM1. In total, 53,249 lncRNAs (259 known lncRNAs and 52,990 novel lncRNAs) were detected, of which 993 lncRNAs showed significantly differential expression (DElncRNAs) under heat stress. By predicting the target genes, 953 common targets shared by cis- and trans-regulation of the DElncRNAs were identified, which exhibited differential expression between the control and the heat stress treatments. Functional annotation indicated that a number of important biological processes and pathways, including photosynthesis, metabolism, translation, stress response, hormone signal transduction, and spliceosome, were enriched for the common targets, suggesting that they play important roles in heat response. A lncRNA-mediated regulatory network was constructed to visualize the molecular response mechanism in response to heat stress, which represented the direct regulatory relationships of DElncRNAs, differentially expressed miRNAs, target genes, and functional annotations.

CONCLUSIONS

This study lays a foundation for further elucidation of the regulatory mechanism for the response to heat stress in the maize inbred line CM1. The findings provide important information for identification of heat-responsive genes, which will be beneficial for the molecular breeding in the cultivation of heat-tolerant maize germplasm.

Global identification and integrated analysis of heat-responsive long non-coding RNAs in contrasting rice cultivars

Identified 2743 rice lncRNAs LncRNAs in response to heat stress Function prediction of HRLs Network among HRLs, genes and miRNAs co-localization of HRLs with QTLs Significant motifs in HRL sequences Long non-coding RNAs (lncRNAs) play vital roles in plant responses to environmental challenges. A better understanding of the gene regulation mediated by lncRNAs and their systematic identification would provide great benefits for modern agriculture. In this study, we performed strand-specific RNA sequencing for two rice varieties, heat-tolerant ZS97B and heat-susceptible SYD2 under heat stress. In total, 2743 putative lncRNAs were identified, and their expression profiles in response to heat treatments were established. We identified 231 differentially expressed lncRNAs (DELs) under heat stress, including 31 DELs common to both varieties and 103 and 97 specific to ZS97B and SYD2, respectively, all defined as heat-responsive lncRNAs (HRLs). The target-coding genes of HRLs were predicted, and GO and KEGG annotations of HRL targets revealed functions in which HRLs might be involved. The interaction network between HRLs, target genes and relevant miRNAs was constructed. The HRLs and their targets were compared with publicly available QTLs for rice seedling growth under heat stimulus. Ten HRLs and twelve target genes were linked with five heat stress-relevant QTLs. Sequence analysis revealed several motifs significantly enriched within the 231 HRL sequences. Our findings provide a valuable resource for further characterization of lncRNAs in terms of heat response and plant heat tolerance improvement.

Identification and expression analysis of genome-wide long noncoding RNA responsive CO2 fluctuated environment in marine microalga Nannochloropsis oceanica

Long non-coding RNAs (lncRNAs) have been demonstrated to participate in plant growth and development as well as response to different biotic and abiotic stresses. However, the knowledge of lncRNA was limited in microalgae. In this study, by RNA deep sequencing, 134 lncRNAs were identified in marine Nannochloropsis oceanica in response to carbon dioxide fluctuation. Among them, there were 51 lncRNAs displayed differentially expressed between low and high CO₂ treatments, including 33 upregulation and 18 downregulation lncRNAs. Cellulose metabolic process, glucan metabolic process, polysaccharide metabolic process, and transmembrane transporter activity were functionally enriched. Multiple potential target genes of lncRNA and lncRNA-mRNA co-located gene network were analyzed. Subsequent analysis had demonstrated that lncRNAs would participate in many biological molecular processes, including gene expression, transcriptional regulation, protein expression and epigenetic regulation. In addition, alternative splicing events were firstly analyzed in response to CO₂ fluctuation. There were 2051 alternative splicing (AS events) identified, which might be associated with lncRNA. These observations will provide a novel insight into lncRNA function in Nannochloropsis and provide a series of targets for lncRNA-based gene editing in future.

Comprehensive Analysis of Long Non-coding RNA Modulates Axillary Bud Development in Tobacco (Nicotiana tabacum L.)

Long non-coding RNAs (lncRNAs) regulate gene expression and are crucial for plant growth and development. However, the mechanisms underlying the effects of activated lncRNAs on axillary bud development remain largely unknown. By lncRNA transcriptomes of axillary buds in topped and untopped tobacco plants, we identified a total of 13,694 lncRNAs. LncRNA analysis indicated that the promoted growth of axillary bud by topping might be partially ascribed to the genes related to hormone signal transduction and glycometabolism, trans-regulated by differentially expressed lncRNAs, such as MSTRG.52498.1, MSTRG.60026.1, MSTRG.17770.1, and MSTRG.32431.1. Metabolite profiling indicated that auxin, abscisic acid and gibberellin were decreased in axillary buds of topped tobacco lines, while cytokinin was increased, consistent with the expression levels of related lncRNAs. MSTRG.52498.1, MSTRG.60026.1, MSTRG.17770.1, and MSTRG.32431.1 were shown to be influenced by hormones and sucrose treatments, and were associated with changes of axillary bud growth in the overexpression of NtCCD8 plants (with reduced axillary buds) and RNA interference of NtTB1 plants (with increased axillary buds). Moreover, MSTRG.28151.1 was identified as the antisense lncRNA of NtTB1. Silencing of MSTRG.28151.1 in tobacco significantly attenuated the expression of NtTB1 and resulted in larger axillary buds, suggesting the vital function of MSTRG.28151.1 axillary bud developmen by NtTB1. Our findings shed light on lncRNA-mRNA interactions and their functional roles in axillary bud growth, which would improve our understanding of lncRNAs as important regulators of axillary bud development and plant architecture.

Conservation and Divergence of SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family between Wheat and Rice

The SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family affects plant architecture, panicle structure, and grain development, representing key genes for crop improvements. The objective of the present study is to utilize the well characterized SPLs' functions in rice to facilitate the functional genomics of TaSPL genes. To achieve these goals, we combined several approaches, including genome-wide analysis of TaSPLs, comparative genomic analysis, expression profiling, and functional study of TaSPL3 in rice. We established the orthologous relationships of 56 TaSPL genes with the corresponding OsSPLs, laying a foundation for the comparison of known SPL functions between wheat and rice. Some TaSPLs exhibited different spatial-temporal expression patterns when compared to their rice orthologs, thus implicating functional divergence. TaSPL2/6/8/10 were identified to respond to different abiotic stresses through the combination of RNA-seq and qPCR expression analysis. Additionally, ectopic expression of TaSPL3 in rice promotes heading dates, affects leaf and stem development, and leads to smaller panicles and decreased yields per panicle. In conclusion, our work provides useful information toward cataloging of the functions of TaSPLs, emphasized the conservation and divergence between TaSPLs and OsSPLs, and identified the important SPL genes for wheat improvement.

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.

Expression of lncRNAs in response to bacterial infections of goat mammary epithelial cells reveals insights into mammary gland diseases

Mastitis or inflammation of the mammary gland is a highly economic and deadly alarming disease for the dairy sector as well as policymakers caused by microbial infection. Transcriptomic and proteomic approaches have been widely employed to identify the underlying molecular mechanisms of bacterial infections in the mammary gland. Numerous differentially expressed mRNAs, miRNAs, and proteins together with their associated signaling pathways have been identified during bacterial infection, paving the way for analysis of their biological functions. Long noncoding RNAs (lncRNAs) are important regulators of multiple biological processes. However, little is known regarding their role in bacterial infection in mammary epithelial cells. Hence, RNA-sequencing was performed by infecting primary mammary epithelial cells (pMECs) with both gram-negative (E. coli) and gram-positive bacteria (S. aureus). Using stringent pipeline, a set of 1957 known and 1175 novel lncRNAs were identified, among which, 112 lncRNAs were found differentially expressed in bacteria challenged PMECs compared with the control. Additionally, potential targets of the lncRNAs were predicted in cis- and trans-configuration. KEGG analysis revealed that DE lncRNAs were associated with at least 15 immune-related pathways. Therefore, our study revealed that bacterial challenge triggers the expression of lncRNAs associated with immune response and defense mechanisms in goat mammary epithelial cells.

Identification and functional deciphering suggested the regulatory roles of long intergenic ncRNAs (lincRNAs) in increasing grafting pepper resistance to Phytophthora capsici

BACKGROUND

As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown.

RESULTS

In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks.

CONCLUSIONS

This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Regulatory short RNAs: A decade's tale for manipulating salt tolerance in plants

Salt stress is a globally increasing environmental detriment to crop growth and productivity. Exposure to salt stress evokes a complex medley of cellular signals, which rapidly reprogram transcriptional and metabolic networks to shape plant phenotype. To date, genetic engineering approaches were used with success to enhance salt tolerance; however, their performance is yet to be evaluated under realistic field conditions. Regulatory short non-coding RNAs (rsRNAs) are emerging as next-generation candidates for engineering salt tolerance in crops. In view of this, the present review provides a comprehensive analysis of a decade's worth of functional studies on non-coding RNAs involved in salt tolerance. Further, we have integrated this knowledge of rsRNA-mediated regulation with the current paradigm of salt tolerance to highlight two regulatory complexes (RCs) for regulating salt tolerance in plants. Finally, a knowledge-driven roadmap is proposed to judiciously utilize RC component(s) for enhancing salt tolerance in crops.

Identification of Long Non-Coding RNAs Associated with Tomato Fruit Expansion and Ripening by Strand-Specific Paired-End RNA Sequencing

As emerging essential regulators in plant development, long non-coding RNAs (lncRNAs) have been extensively investigated in multiple horticultural crops, as well as in different tissues of plants. Tomato fruits are an indispensable part of people’s diet and are consumed as fruits and vegetables. Meanwhile, tomato is widely used as a model to study the ripening mechanism in fleshy fruit. Although increasing evidence shows that lncRNAs are involved in lots of biological processes in tomato plants, the comprehensive identification of lncRNAs in tomato fruit during its expansion and ripening and their functions are partially known. Here, we performed strand-specific paired-end RNA sequencing (ssRNA-seq) of tomato Heinz1706 fruits at five different developmental stages, as well as flowers and leaves. We identified 17,674 putative lncRNAs by referencing the recently released SL4.0 and annotation ITAG4.0 in tomato plants. Many lncRNAs show different expression patterns in fleshy fruit at different developmental stages compared with leaves or flowers. Our results indicate that lncRNAs play an important role in the regulation of tomato fruit expansion and ripening, providing informative lncRNA candidates for further studies in tomato fruits. In addition, we also summarize the recent advanced progress in lncRNAs mediated regulation on horticultural fruits. Hence, our study updates the understanding of lncRNAs in horticultural plants and provides resources for future studies relating to the expansion and ripening of tomato fruits.

Epigenetic regulation of salinity stress responses in cereals

Cereals are important crops and are exposed to various types of environmental stresses that affect the overall growth and yield. Among the various abiotic stresses, salt stress is a major environmental factor that influences the genetic, physiological, and biochemical responses of cereal crops. Epigenetic regulation which includes DNA methylation, histone modification, and chromatin remodelling plays an important role in salt stress tolerance. Recent studies in rice genomics have highlighted that the epigenetic changes are heritable and therefore can be considered as molecular signatures. An epigenetic mechanism under salinity induces phenotypic responses involving modulations in gene expression. Association between histone modification and altered DNA methylation patterns and differential gene expression has been evidenced for salt sensitivity in rice and other cereal crops. In addition, epigenetics also creates stress memory that helps the plant to better combat future stress exposure. In the present review, we have discussed epigenetic influences in stress tolerance, adaptation, and evolution processes. Understanding the epigenetic regulation of salinity could help for designing salt-tolerant varieties leading to improved crop productivity.

Genome-wide identification and functional prediction of salt- stress related long non-coding RNAs (lncRNAs) in chickpea (Cicer arietinum L.)

LncRNAs (long noncoding RNAs) are 200 bp length crucial RNA molecules, lacking coding potential and having important roles in regulating gene expression, particularly in response to abiotic stresses. In this study, we identified salt stress-induced lncRNAs in chickpea roots and predicted their intricate regulatory roles. A total of 3452 novel lncRNAs were identified to be distributed across all 08 chickpea chromosomes. On comparing salt-tolerant (ICCV 10, JG 11) and salt-sensitive cultivars (DCP 92-3, Pusa 256), 4446 differentially expressed lncRNAs were detected under various salt  treatments. We predicted 3373 lncRNAs to be regulating their target genes in cis regulating manner and 80 unique lncRNAs were observed as interacting with 136 different miRNAs, as eTMs (endogenous target mimic) targets of miRNAs and implicated them in the regulatory network of salt stress response. Functional analysis of these lncRNA revealed their association in targeting salt stress response-related genes like potassium transporter, transporter family genes, serine/threonine-protein kinase, aquaporins like TIP1-2, PIP2-5 and transcription factors like, AP2, NAC, bZIP, ERF, MYB and WRKY. Furthermore, about 614 lncRNA-SSRs (simple sequence repeats) were identified as a new generation of molecular markers with higher efficiency and specificity in chickpea. Overall, these findings will pave the understanding of comprehensive functional role of potential lncRNAs, which can help in providing insight into the molecular mechanism of salt tolerance in chickpea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01093-0.

The Roles of Temperature-Related Post-Transcriptional Regulation in Cereal Floral Development

Temperature is a critical environmental signal in the regulation of plant growth and development. The temperature signal varies across a daily 24 h period, between seasons and stochastically depending on local environmental events. Extracting important information from these complex signals has led plants to evolve multiple temperature responsive regulatory mechanisms at the molecular level. In temperate cereals, we are starting to identify and understand these molecular mechanisms. In addition, we are developing an understanding of how this knowledge can be used to increase the robustness of crop yield in response to significant changes in local and global temperature patterns. To enable this, it is becoming apparent that gene regulation, regarding expression and post-transcriptional regulation, is crucial. Large transcriptomic studies are identifying global changes in spliced transcript variants and regulatory non-coding RNAs in response to seasonal and stress temperature signals in many of the cereal crops. Understanding the functions of these variants and targets of the non-coding RNAs will greatly increase how we enable the adaptation of crops. This review considers our current understanding and areas for future development.

Identification and characterization of long non-coding RNA and their response against citrus bark cracking viroid infection in Humulus lupulus

Long non-coding RNAs (lncRNAs) are a highly heterogeneous class of non-protein-encoding transcripts that play an essential regulatory role in diverse biological processes, including stress responses. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) poses a major threat to the production of Humulus lupulus (hop) plants. In this study, we systematically investigate the characteristics of the lncRNAs in hop and their role in CBCVd-infection using RNA-sequencing data. Following a stringent filtration criterion, a total of 3598 putative lncRNAs were identified with a high degree of certainty, of which 19% (684) of the lncRNAs were significantly differentially expressed (DE) in CBCVd-infected hop, which were predicted to be mainly involved in plant-pathogen interactions, kinase cascades, secondary metabolism and phytohormone signal transduction. Besides, several lncRNAs and CBCVd-responsive lncRNAs were identified as the precursor of microRNAs and predicted as endogenous target mimics (eTMs) for hop microRNAs involved in CBCVd-infection.

Genome wide identification and characterization of abiotic stress responsive lncRNAs in Capsicum annuum

Long non-coding RNAs (lncRNAs) are a type of non-coding transcripts having length of more than 200 nucleotides lacking protein-coding ability. In the present study, 12807 lncRNAs were identified in Capsicum annuum tissues exposed to abiotic stress conditions viz. heat, cold, osmotic and salinity stress. Expression analysis of lncRNAs in different treatment conditions demonstrates their stress-specific expression. Thirty lncRNAs were found to act as precursors for 10 microRNAs (miRNAs) of C. annuum. Additionally, a total of 1807 lncRNAs were found to interact with 194 miRNAs which targeted 621 mRNAs of C. annuum. Among these, 344 lncRNAs were found to act as target mimics for 621 genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that out of those 621 gene sequences, 546 were tagged with GO terms, 105 Enzyme Code (EC) numbers were assigned to 246 genes and 223 genes are found to be involved in 63 biological pathways. In this report, we have highlighted the prospective role of lncRNAs in different abiotic stress conditions by interacting with miRNAs and regulating stress responsive transcription factors (TFs) such as DoF, WRKY, MYB, bZIP and ERF in C. annuum.

Comparative analysis of long noncoding RNAs in angiosperms and characterization of long noncoding RNAs in response to heat stress in Chinese cabbage

Long noncoding RNAs (lncRNAs) are widely present in different species and play critical roles in response to abiotic stresses. However, the functions of lncRNAs in Chinese cabbage under heat stress remain unknown. Here, we first conducted a global comparative analysis of 247,242 lncRNAs among 37 species. The results indicated that lncRNAs were poorly conserved among different species, and only 960 lncRNAs were homologous to 524 miRNA precursors. We then carried out lncRNA sequencing for a genome-wide analysis of lncRNAs and their target genes in Chinese cabbage at different stages of heat treatment. In total, 18,253 lncRNAs were identified, of which 1229 differentially expressed (DE) lncRNAs were characterized as being heat-responsive. The ceRNA network revealed that 38 lncRNAs, 16 miRNAs, and 167 mRNAs were involved in the heat response in Chinese cabbage. Combined analysis of the cis- and trans-regulated genes indicated that the targets of DE lncRNAs were significantly enriched in the "protein processing in endoplasmic reticulum" and "plant hormone signal transduction" pathways. Furthermore, the majority of HSP and PYL genes involved in these two pathways exhibited similar expression patterns and responded to heat stress rapidly. Based on the networks of DE lncRNA-mRNAs, 29 and 22 lncRNAs were found to interact with HSP and PYL genes, respectively. Finally, the expression of several critical lncRNAs and their targets was verified by qRT-PCR. Overall, we conducted a comparative analysis of lncRNAs among 37 species and performed a comprehensive analysis of lncRNAs in Chinese cabbage. Our findings expand the knowledge of lncRNAs involved in the heat stress response in Chinese cabbage, and the identified lncRNAs provide an abundance of resources for future comparative and functional studies.

High-Resolution Transcriptome Atlas and Improved Genome Assembly of Common Buckwheat, Fagopyrum esculentum

Common buckwheat (Fagopyrum esculentum) is an important non-cereal grain crop and a prospective component of functional food. Despite this, the genomic resources for this species and for the whole family Polygonaceae, to which it belongs, are scarce. Here, we report the assembly of the buckwheat genome using long-read technology and a high-resolution expression atlas including 46 organs and developmental stages. We found that the buckwheat genome has an extremely high content of transposable elements, including several classes of recently (0.5-1 Mya) multiplied TEs ("transposon burst") and gradually accumulated TEs. The difference in TE content is a major factor contributing to the three-fold increase in the genome size of F. esculentum compared with its sister species F. tataricum. Moreover, we detected the differences in TE content between the wild ancestral subspecies F. esculentum ssp. ancestrale and buckwheat cultivars, suggesting that TE activity accompanied buckwheat domestication. Expression profiling allowed us to test a hypothesis about the genetic control of petaloidy of tepals in buckwheat. We showed that it is not mediated by B-class gene activity, in contrast to the prediction from the ABC model. Based on a survey of expression profiles and phylogenetic analysis, we identified the MYB family transcription factor gene tr_18111 as a potential candidate for the determination of conical cells in buckwheat petaloid tepals. The information on expression patterns has been integrated into the publicly available database TraVA: http://travadb.org/browse/Species=Fesc/. The improved genome assembly and transcriptomic resources will enable research on buckwheat, including practical applications.

Global Profiling of lncRNAs Expression Responsive to Allopolyploidization in Cucumis

Long non-coding RNAs (lncRNAs) play critical regulatory roles in various biological processes. However, the presence of lncRNAs and how they function in plant polyploidy are still largely unknown. Hence, we examined the profile of lncRNAs in a nascent allotetraploid Cucumis hytivus (S₁₄), its diploid parents, and the F₁ hybrid, to reveal the function of lncRNAs in plant-interspecific hybridization and whole genome duplication. Results showed that 2206 lncRNAs evenly transcribed from all 19 chromosomes were identified in C. hytivus, 44.6% of which were from intergenic regions. Based on the expression trend in allopolyploidization, we found that a high proportion of lncRNAs (94.6%) showed up-regulated expression to varying degrees following hybridization. However, few lncRNAs (33, 2.1%) were non-additively expressed after genome duplication, suggesting the significant effect of hybridization on lncRNAs, rather than genome duplication. Furthermore, 253 cis-regulated target genes were predicted for these differentially expressed lncRNAs in S₁₄, which mainly participated in chloroplast biological regulation (e.g., chlorophyll synthesis and light harvesting system). Overall, this study provides new insight into the function of lncRNAs during the processes of hybridization and polyploidization in plant evolution.

Genome-wide identification and characterization of long non-coding RNAs related to grain yield in foxtail millet [Setaria italica (L.) P. Beauv.]

Background

Long noncoding RNAs (lncRNAs) have been reported to play critical roles in diverse growth and development processes in plants. However, the systematic identification and characterization of lncRNAs in foxtail millet is nearly blank.

Results

In this study, we performed high-throughput sequencing of young spikelets from four foxtail millet varieties in different yield levels at booting stage. As a result, a total of 12,378 novel lncRNAs were identified, and 70 were commonly significantly differentially expressed in comparisons between high-yield varieties and conventional varieties, suggesting that they involved in yield formation and regulation in foxtail millet. Functional analysis revealed that among the 70 significantly differentially expressed lncRNAs, 67 could transcriptionally modulate target genes in cis and in trans. Moreover, 18 lncRNAs related to grain yield in foxtail millet were predicted to function as miRNA target mimics and regulate gene expression by competing for the interaction between miRNAs and their target mRNAs.

Conclusion

Our results will provide materials for elucidation of the molecular mechanisms of lncRNAs participate in yield regulation, and will contribute to high yield foxtail millet breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07272-9.

Identification of conserved genes linked to responses to abiotic stresses in leaves among different plant species

As a consequence of global climate change, certain stress factors that have a negative impact on crop productivity such as heat, cold, drought and salinity are becoming increasingly prevalent. We conducted a meta-analysis to identify genes conserved across plant species involved in (1) general abiotic stress conditions, and (2) specific and unique abiotic stress factors (drought, salinity, extreme temperature) in leaf tissues. We collected raw data and re-analysed eight RNA-Seq studies using our previously published bioinformatic pipeline. A total of 68 samples were analysed. Gene set enrichment analysis was performed using MapMan and PageMan whereas DAVID (Database for Annotation, Visualisation and Integrated Discovery) was used for metabolic process enrichment analysis. We identified of a total of 5122 differentially expressed genes when considering all abiotic stresses (3895 were upregulated and 1227 were downregulated). Jasmonate-related genes were more commonly upregulated by drought, whereas gibberellin downregulation was a key signal for drought and heat. In contrast, cold stress clearly upregulated genes involved in ABA (abscisic acid), cytokinin and gibberellins. A gene (non-phototrophic hypocotyl) involved in IAA (indoleacetic acid) response was induced by heat. Regarding secondary metabolism, as expected, MVA pathway (mevalonate pathway), terpenoids and alkaloids were generally upregulated by all different stresses. However, flavonoids, lignin and lignans were more repressed by heat (cinnamoyl coA reductase 1 and isopentenyl pyrophosphatase). Cold stress drastically modulated genes involved in terpenoid and alkaloids. Relating to transcription factors, AP2-EREBP, MADS-box, WRKY22, MYB, homoebox genes members were significantly modulated by drought stress whereas cold stress enhanced AP2-EREBPs, bZIP members, MYB7, BELL 1 and one bHLH member. C2C2-CO-LIKE, MADS-box and a homeobox (HOMEOBOX3) were mostly repressed in response to heat. Gene set enrichment analysis showed that ubiquitin-mediated protein degradation was enhanced by heat, which unexpectedly repressed glutaredoxin genes. Cold stress mostly upregulated MAP kinases (mitogen-activated protein kinase). Findings of this work will allow the identification of new molecular markers conserved across crops linked to major genes involved in quantitative agronomic traits affected by different abiotic stress.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-Wide Identification and Characterization of Fusarium graminearum-Responsive lncRNAs in Triticum aestivum

Although the war between wheat and Fusarium has been widely investigated for years, long noncoding RNAs (lncRNAs), which have been proven to regulate important processes in the development and stress responses of plants, are still poorly known in wheat against Fusarium. Herein, we systematically reveal the roles of wheat lncRNAs in the process of Fusarium graminearum infection by high-throughput RNA sequencing. Well over 4130 of the total 4276 differentially expressed lncRNAs were already specifically expressed at 12 h postinoculation (hpi), but only 89 of these were specifically expressed at 24 hpi, indicating that the initial stage was the crucial stage for lncRNA-mediated gene regulation of wheat defense against F. graminearum. Target analysis showed the lncRNAs participated in various biological stress processes and had exclusive regulation models at different infection stages. Further H₂O₂ accumulation and protein ubiquitination assays supported this idea. Moreover, two lncRNAs (XLOC_302848 and XLOC_321638) were identified as Fusarium seedling blight resistance candidates by lncRNA-target expression pattern validation, and two lncRNAs (XLOC_113815, XLOC_123624) were Fusarium head blight resistance potential regulators by cross-validating the RNAseq data with the refined meta-QTL of wheat FHB resistance. These findings extend our knowledge on wheat lncRNAs response to F. graminearum attack and provide new insights for the functional and molecular research of future interactions between wheat and Fusarium.

Long non-coding RNA HULC exerts oncogenic activity on papillary thyroid cancer in vitro and in vivo

Thyroid cancer is a frequently happened malignancy in human endocrine system. Papillary thyroid cancer (PTC) presents 70-80% of all thyroid cancer cases. Herein, we probed the possible oncogenic function of long non-coding RNA (lncRNA) highly up-regulated in liver cancer (HULC) in PTC. First, the HULC and microRNA-106a (miR-106a) expressions in PTC tissues and cells were tested. Plasmids or miRNAs transfections were done for altering HULC and miR-106a expressions. Then, cells viability and apoptosis, along with cell proliferative, migratory and invasive abilities, were tested, respectively. The PI3K/AKT and Wnt/β-catenin pathways activities were measured. Finally, the animal model of PTC was constructed and the tumour volumes and weights were gauged. We discovered that HULC and miR-106a had relative high expression levels in PTC tissues and cells. HULC overexpression enhanced TPC-1 cells viability and cell proliferative, migratory and invasive abilities. Silencing HULC induced TPC-1 cell apoptosis. miR-106a engaged in the oncogenic impacts of HULC. Moreover, HULC overexpression boosted PI3K/AKT and Wnt/β-catenin pathways activities via raising miR-106a expression. Besides, HULC overexpression enhanced the volumes and weights of PTC tumours. To sum up, HULC exhibited oncogenic function on PTC in vitro and in vivo.

Identification and characterization of long noncoding RNAs provide insight into the regulation of gene expression in response to heat stress in rainbow trout (Oncorhynchus mykiss)

Rainbow trout are typical cold-water fish species. However, with the intensification of global warming, high temperatures have severely restricted the development of aquaculture during the summer. Understanding the molecular regulatory mechanisms of rainbow trout responses to heat stress will be beneficial for alleviating heat stress-related damage. In this study, we performed RNA-seq of liver tissues from rainbow trout under heat stress (24 °C) and control conditions (18 °C) to identify lncRNAs and target genes by strand-specific library. Changes in nonspecific immune parameters revealed that a strong stress response occurred in rainbow trout at 24 °C. More than 658 million filtered reads and 5916 lncRNAs were identified from six libraries. A total of 927 novel lncRNAs were identified, and 428 differentially expressed lncRNAs were screened with stringent thresholds. The RNA-seq results were verified by RT-qPCR. In addition, a regulatory network of lncRNA-mRNA functional interactions was constructed, and the potential antisense, cis and trans targets of lncRNAs were predicted. GO and KEGG enrichment analyses showed that many target genes involved in maintenance of homeostasis or adaptation to stress and stimuli were highly induced under heat stress. Several regulatory pathways were also found to be involved in heat stress, including the thyroid hormone signaling pathway, the PI3K-Akt signaling pathway, and the estrogen signaling pathway, among others. These results broaden our understanding of lncRNAs associated with heat stress and provide new insights into the lncRNA mediated regulation of the rainbow trout heat stress response.

Identification and characterization of long non-coding RNAs as competing endogenous RNAs in the cold stress response of Triticum aestivum

Long non-coding RNAs (lncRNAs) play important roles in plant development and stress responses. MicroRNAs (miRNAs) are involved in transcriptional and post-transcriptional gene regulation. It is not clear how lncRNA-mediated plant responses to cold stress and how lncRNAs, miRNAs and target mRNAs cooperate subject to the competing endogenous RNA (ceRNA). We interpreted the function of lncRNAs in the winter wheat cultivar Dongnongdongmai 1 (Dn1). A total of 9970 putative lncRNAs were initially identified from three Dn1 lncRNA libraries (5 °C, -10 °C and -25 °C) using high-throughput sequencing. Among the 14,626 genes detected via weighted gene co-expression network analysis, 7435 lncRNAs were co-expressed with 7191 mRNAs. We found six modules related to cold resistance in the lncRNA-mRNA weighted co-expression network, and the functions of mRNAs were similar in each module. Antioxidant systems and hormones played important roles in low-temperature responses. RNA sequencing analysis revealed that interactions between the 384 lncRNAs and 70 miRNAs were required for ceRNA activity. According to ceRNA activity, 225 lncRNAs, 60 miRNAs and 621 target mRNAs were involved in the regulatory networks of the cold stress response. Notably, a conserved region was found in the complementary regions of lncRNAs and miR164/408 but had reverse expression trends in the ceRNA network. Our results reveal possible roles of lncRNAs-mRNAs in the regulatory networks associated with tolerance to low temperature and provide useful information for more strategic use of genomic resources in wheat breeding.

Genome-wide identification and integrated analysis of lncRNAs in rice backcross introgression lines (BC2F12)

BACKGROUND

Distant hybridization is an important way to create interspecific genetic variation and breed new varieties in rice. A lot of backcross introgression lines (BILs) had been constructed for the scientific issues in rice. However, studies on the critical regulatory factor lncRNA in cultivated rice, wild rice and their BIL progenies were poorly reported.

RESULTS

Here, high-throughput RNA sequencing technology was used to explore the functional characteristics and differences of lncRNAs in O. sativa, O. longistaminata and their three BC₂F₁₂ progenies. A total of 1254 lncRNAs were screened out, and the number of differentially expressed lncRNAs between progenies and O. sativa were significantly less than that between progenies and O. longistaminata. Some lncRNAs regulated more than one mRNA, and 89.5% of lncRNAs regulated the expression of target genes through cis-acting. A total of 78 lncRNAs and 271 mRNAs were targeted by 280 miRNAs, and 22 lncRNAs were predicted to be the precursor of 20 microRNAs. Some miRNAs were found to target their own potential precursor lncRNAs. Over 50% of lncRNAs showed parental expression level dominance (ELD) in all three progenies, and most lncRNAs showed ELD-O. sativa rather than ELD-O. longistaminata. Further analysis showed that lncRNAs might regulate the expression of plant hormone-related genes and the adaptability of O. sativa, O. longistaminata and their progenies.

CONCLUSIONS

Taken together, the above results provided valuable clues for elucidating the functional features and expression differences of lncRNAs between O. sativa, O. longistaminata and their BIL progenies, and expanded our understanding about the biological functions of lncRNAs in rice.

Identification of lncRNAs involved in wheat tillering development in two pairs of near-isogenic lines

Emerging evidence demonstrates that lncRNAs participate in various developmental processes in plants via post-transcription regulation. However, few lncRNAs have been identified as regulators of tiller development in wheat (Triticum aestivum L.). In this study, high-throughput ribosomal depleted RNA sequencing was performed on the tillering nodes of two pairs of near-isogenic lines that differed only in the tillering trait. We identified 5399 lncRNA transcripts using bioinformational analyses. KEGG pathway analysis revealed 74 common differentially expressed lncRNAs substantially enriched in photosynthesis-related, phenylpropanoid biosynthesis, phosphatidylinositol signaling, brassinosteroid biosynthesis, zeatin biosynthesis, and carotenoid biosynthesis pathways. Detailed functional annotations of target genes were used to identify 27 tillering-associated lncRNAs. Among these, 10 were in photosynthesis-related pathways; 15 were in secondary metabolite pathways; and 8 were in plant hormone pathways, with 6 enriched in two kinds of pathways. These findings contribute to identifying tillering-associated lncRNAs in wheat and enable further investigation into the functions and roles of key candidate lncRNAs, and more experimental evidence was also needed if breeders wanted to utilize these candidate lncRNAs in wheat crop yield improvement in the future.

Bacillus amyloliquefaciens LZ04 improves the resistance of Arabidopsis thaliana to high calcium stress and the potential role of lncRNA-miRNA-mRNA regulatory network in the resistance

Bacillus amyloliquefaciens is a non-pathogenic and plant growth-promoting rhizobacterium that enhances plant resistance to drought and diseases. Arabidopsis thaliana is a multipurpose model plant for exploring microorganism-plant interactions and a crucial vegetal tool for molecular research. Non-coding RNAs are RNA molecules involved in the regulation of various biological functions and constitute a research hotspot in the field of plant biology. In this study, the effect of B. amyloliquefaciens treatment on the resistance of A. thaliana to high calcium stress was analyzed. The transcriptome sequencing of A. thaliana roots under four treatment conditions was performed to screen differentially expressed lncRNAs, mRNAs and miRNAs. Functional analysis was also performed to understand the potential mechanism by which B. amyloliquefaciens-regulated lncRNAs, miRNAs and mRNAs affect the resistance of A. thaliana to high calcium stress. The results indicated that B. amyloliquefaciens treatment increased the resistance of A. thaliana to high calcium stress. A set of differentially expressed lncRNAs, mRNAs and miRNAs were screened between the high calcium and control group on one hand, and high calcium and high calcium + B. amyloliquefaciens groups on the other hand. Functional analysis indicated that the differentially expressed mRNAs and miRNA were involved in various biological functions and that transcriptional dysregulation caused by high calcium stress involves metabolic processes rather than defense responses. Conclusively, B. amyloliquefaciens may improve the resistance of A. thaliana to high calcium stress via a lncRNA-miRNA-mRNA regulatory network. These findings will contribute to the development of agriculture in karst regions with high calcium content.

Identification and characterization of long non-coding RNAs regulating resistant starch biosynthesis in bread wheat (Triticum aestivum L.)

Resistant starch (RS) also known as healthy starch has shown several health benefits. Enhancing the RS through genetic modification approaches has huge commercial importance. Regulatory RNA like long non-coding RNA (lncRNA) plays an important role in gene regulation. In this study, we mined 63 transcriptome datasets of wheat belonging to 35 genotypes representing two seed developmental stages. Contrasting expression of a subset of lncRNAs in RS mutant lines compared to parent wheat variety 'C 306' signifies their probable role in RS biosynthesis. Further, lncRNA- TCONS_00130663 showed strong positive correlation (r² = 1) with LYPL gene and strong negative correlation with SBEIIb (r² = -0.94). We found TCONS_00130663 as positive regulator of LYPL gene through interaction with miR1128. Based on relative expression, in silico interaction and DSC analysis we hypothesize the dual role of TCONS_00130663 in RS type 2 and type 5. The study provides a useful resource for functional mechanism of lncRNAs.

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.

A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

Genome-wide discovery and functional prediction of salt-responsive lncRNAs in duckweed

Background

Salt significantly depresses the growth and development of the greater duckweed, Spirodela polyrhiza, a model species of floating aquatic plants. Physiological responses of this plant to salt stress have been characterized, however, the roles of long noncoding RNAs (lncRNAs) remain unknown.

Results

In this work, totally 2815 novel lncRNAs were discovered in S. polyrhiza by strand-specific RNA sequencing, of which 185 (6.6%) were expressed differentially under salinity condition. Co-expression analysis indicated that the trans-acting lncRNAs regulated their co-expressed genes functioning in amino acid metabolism, cell- and cell wall-related metabolism, hormone metabolism, photosynthesis, RNA transcription, secondary metabolism, and transport. In total, 42 lncRNA-mRNA pairs that might participate in cis-acting regulation were found, and these adjacent genes were involved in cell wall, cell cycle, carbon metabolism, ROS regulation, hormone metabolism, and transcription factor. In addition, the lncRNAs probably functioning as miRNA targets were also investigated. Specifically, TCONS_00033722, TCONS_00044328, and TCONS_00059333 were targeted by a few well-studied salt-responsive miRNAs, supporting the involvement of miRNA and lncRNA interactions in the regulation of salt stress responses. Finally, a representative network of lncRNA-miRNA-mRNA was proposed and discussed to participate in duckweed salt stress via auxin signaling.

Conclusions

This study is the first report on salt-responsive lncRNAs in duckweed, and the findings will provide a solid foundation for in-depth functional characterization of duckweed lncRNAs in response to salt stress.

Insight into the Role of Epigenetic Processes in Abiotic and Biotic Stress Response in Wheat and Barley

Environmental stresses such as salinity, drought, heat, freezing, heavy metal and even pathogen infections seriously threaten the growth and yield of important cereal crops including wheat and barley. There is growing evidence indicating that plants employ sophisticated epigenetic mechanisms to fine-tune their responses to environmental stresses. Here, we provide an overview of recent developments in understanding the epigenetic processes and elements—such as DNA methylation, histone modification, chromatin remodeling, and non-coding RNAs—involved in plant responses to abiotic and biotic stresses in wheat and barley. Potentials of exploiting epigenetic variation for the improvement of wheat and barley are discussed.