Long noncoding RNA transcriptome of plants

Noncoding RNAs in regulation of plant secondary metabolism

Plant secondary metabolites (PSMs) are a large class of structurally diverse molecules, mainly consisting of terpenoids, phenolic compounds, and nitrogen-containing compounds, which play active roles in plant development and stress responses. The biosynthetic processes of PSMs are governed by a sophisticated regulatory network at multiple levels. Noncoding RNAs (ncRNAs) such as microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs) may serve as post-transcriptional regulators for plant secondary metabolism through acting on genes encoding either transcription factors or participating enzymes in relevant metabolic pathways. High-throughput sequencing technologies have facilitated the large-scale identifications of ncRNAs potentially involved in plant secondary metabolism in model plant species as well as certain species with enriched production of specific types of PSMs. Moreover, a series of miRNA-target modules have been functionally characterized to be responsible for regulating PSM biosynthesis and accumulation in plants under abiotic or biotic stresses. In this review, we will provide an overview of current findings on the ncRNA-mediated regulation of plant secondary metabolism with special attention to its participation in plant stress responses, and discuss possible issues to be addressed in future fundamental research and breeding practice.

The complex transcriptional regulation of heat stress response in maize

As one of the most important food and feed crops worldwide, maize suffers much more tremendous damages under heat stress compared to other plants, which seriously inhibits plant growth and reduces productivity. To mitigate the heat-induced damages and adapt to high temperature environment, plants have evolved a series of molecular mechanisms to sense, respond and adapt high temperatures and heat stress. In this review, we summarized recent advances in molecular regulations underlying high temperature sensing, heat stress response and memory in maize, especially focusing on several important pathways and signals in high temperature sensing, and the complex transcriptional regulation of ZmHSFs (Heat Shock Factors) in heat stress response. In addition, we highlighted interactions between ZmHSFs and several epigenetic regulation factors in coordinately regulating heat stress response and memory. Finally, we laid out strategies to systematically elucidate the regulatory network of maize heat stress response, and discussed approaches for breeding future heat-tolerance maize.

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Long noncoding RNA from Betula platyphylla, BplncSIR1, confers salt tolerance by regulating BpNAC2 to mediate reactive oxygen species scavenging and stomatal movement

Long noncoding RNAs (lncRNAs) play an important role in abiotic stress tolerance. However, their function in conferring abiotic stress tolerance is still unclear. Herein, we characterized the function of a salt-responsive nuclear lncRNA (BplncSIR1) from Betula platyphylla (birch). Birch plants overexpressing and knocking out for BplncSIR1 were generated. BplncSIR1 was found to improve salt tolerance by inducing antioxidant activity and stomatal closure, and also accelerate plant growth. Chromatin isolation by RNA purification (ChIRP) combined with RNA sequencing indicated that BplncSIR1 binds to the promoter of BpNAC2 (encoding NAC domain-containing protein 2) to activate its expression. Plants overexpressing and knocking out for BpNAC2 were generated. Consistent with that of BplncSIR1, overexpression of BpNAC2 also accelerated plant growth and conferred salt tolerance. In addition, BpNAC2 binds to different cis-acting elements, such as G-box and 'CCAAT' sequences, to regulate the genes involved in salt tolerance, resulting in reduced ROS accumulation and decreased water loss rate by stomatal closure. Taken together, BplncSIR1 serves as the regulator of BpNAC2 to induce its expression in response to salt stress, and activated BpNAC2 accelerates plant growth and improves salt tolerance. Therefore, BplncSIR1 might be a candidate gene for molecular breeding to cultivate plants with both a high growth rate and improved salt tolerance.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs Associated with Floral Scent Formation in Jasmine (Jasminum sambac)

Long non-coding RNAs (lncRNAs) have emerged as curial regulators of diverse biological processes in plants. Jasmine (Jasminum sambac) is a world-renowned ornamental plant for its attractive and exceptional flower fragrance. However, to date, no systematic screening of lncRNAs and their regulatory roles in the production of the floral fragrance of jasmine flowers has been reported. In this study, we identified a total of 31,079 novel lncRNAs based on an analysis of strand-specific RNA-Seq data from J. sambac flowers at different stages. The lncRNAs identified in jasmine flowers exhibited distinct characteristics compared with protein-coding genes (PCGs), including lower expression levels, shorter transcript lengths, and fewer exons. Certain jasmine lncRNAs possess detectable sequence conservation with other species. Expression analysis identified 2752 differentially expressed lncRNAs (DE_lncRNAs) and 8002 DE_PCGs in flowers at the full-blooming stage. DE_lncRNAs could potentially cis- and trans-regulate PCGs, among which DE_lincRNAs and their targets showed significant opposite expression patterns. The flowers at the full-blooming stage are specifically enriched with abundant phenylpropanoids and terpenoids potentially contributed by DE_lncRNA cis-regulated PCGs. Notably, we found that many cis-regulated DE_lncRNAs may be involved in terpenoid and phenylpropanoid/benzenoid biosynthesis pathways, which potentially contribute to the production of jasmine floral scents. Our study reports numerous jasmine lncRNAs and identifies floral-scent-biosynthesis-related lncRNAs, which highlights their potential functions in regulating the floral scent formation of jasmine and lays the foundations for future molecular breeding.

Insights into the roles of long noncoding RNAs in the communication between plants and the environment

In addition to coding proteins, RNA molecules, especially long noncoding RNAs (lncRNAs), have well-established functions in regulating gene expression. The number of studies focused on the roles played by different types of lncRNAs in a variety of plant biological processes has markedly increased. These lncRNA roles involve plant vegetative and reproductive growth and responses to biotic and abiotic stresses. In this review, we examine the classification, mechanisms, and functions of lncRNAs and then emphasize the roles played by these lncRNAs in the communication between plants and the environment mainly with respect to the following environmental factors: temperature, light, water, salt stress, and nutrient deficiencies. We also discuss the consensus among researchers and the remaining challenges and underscore the exciting ways lncRNAs may affect the biology of plants.

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.

Plant long non-coding RNAs: identification and analysis to unveil their physiological functions

Eukaryotic genomes encode thousands of RNA molecules; however, only a minimal fraction is translated into proteins. Among the non-coding elements, long non-coding RNAs (lncRNAs) play important roles in diverse biological processes. LncRNAs are associated mainly with the regulation of the expression of the genome; nonetheless, their study has just scratched the surface. This is somewhat due to the lack of widespread conservation at the sequence level, in addition to their relatively low and highly tissue-specific expression patterns, which makes their exploration challenging, especially in plant genomes where only a few of these molecules have been described completely. Recently published high-quality genomes of crop plants, along with new computational tools, are considered promising resources for studying these molecules in plants. This review briefly summarizes the characteristics of plant lncRNAs, their presence and conservation, the different protocols to find these elements, and the limitations of these protocols. Likewise, it describes their roles in different plant physiological phenomena. We believe that the study of lncRNAs can help to design strategies to reduce the negative effect of biotic and abiotic stresses on the yield of crop plants and, in the future, help create fruits and vegetables with improved nutritional content, higher amounts of compounds with positive effects on human health, better organoleptic characteristics, and fruits with a longer postharvest shelf life.

Systemic movement of long non-coding RNA ELENA1 attenuates leaf senescence under nitrogen deficiency

Nitrogen is an essential macronutrient that is absorbed by roots and stored in leaves, mainly as ribulose-1,5-bisphosphate carboxylase/oxygenase1,2. During nitrogen deficiency (-N), plants activate leaf senescence for source-to-sink nitrogen remobilization for adaptative growth3-6. However, how -N signals perceived by roots are propagated to shoots remains underexplored. We found that ELF18-INDUCED LONG NONCODING RNA 1 (ELENA1) is -N inducible and attenuates -N-induced leaf senescence in Arabidopsis. Analysis of plants expressing the ELENA1 promoter β-glucuronidase fusion gene showed that ELENA1 is transcribed specifically in roots under -N. Reciprocal grafting of the wild type and elena1 demonstrated that ELENA1 functions systemically. ELENA1 dissociates the MEDIATOR SUBUNIT 19a-ORESARA1 transcriptional complex, thereby calibrating senescence progression. Our observations establish the systemic regulation of leaf senescence by a root-derived long non-coding RNA under -N in Arabidopsis.

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.

Barley HOMOCYSTEINE METHYLTRANSFERASE 2 confers drought tolerance by improving polyamine metabolism

Drought stress poses a serious threat to crop production worldwide. Genes encoding homocysteine methyltransferase (HMT) have been identified in some plant species in response to abiotic stress, but its molecular mechanism in plant drought tolerance remains unclear. Here, transcriptional profiling, evolutionary bioinformatics, and population genetics were conducted to obtain insight into the involvement of HvHMT2 from Tibetan wild barley (Hordeum vulgare ssp. agriocrithon) in drought tolerance. We then performed genetic transformation coupled with physio-biochemical dissection and comparative multiomics approaches to determine the function of this protein and the underlying mechanism of HvHMT2-mediated drought tolerance. HvHMT2 expression was strongly induced by drought stress in tolerant genotypes in a natural Tibetan wild barley population and contributed to drought tolerance through S-adenosylmethionine (SAM) metabolism. Overexpression of HvHMT2 promoted HMT synthesis and efficiency of the SAM cycle, leading to enhanced drought tolerance in barley through increased endogenous spermine and less oxidative damage and growth inhibition, thus improving water status and final yield. Disruption of HvHMT2 expression led to hypersensitivity under drought treatment. Application of exogenous spermine reduced accumulation of reactive oxygen species (ROS), which was increased by exogenous mitoguazone (inhibitor of spermine biosynthesis), consistent with the association of HvHMT2-mediated spermine metabolism and ROS scavenging in drought adaptation. Our findings reveal the positive role and key molecular mechanism of HvHMT2 in drought tolerance in plants, providing a valuable gene not only for breeding drought-tolerant barley cultivars but also for facilitating breeding schemes in other crops in a changing global climate.

Comprehensive transcriptomic analysis of three varieties with different brown planthopper-resistance identifies leaf sheath lncRNAs in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) have been brought great attention for their crucial roles in diverse biological processes. However, systematic identification of lncRNAs associated with specialized rice pest, brown planthopper (BPH), defense in rice remains unexplored.

RESULTS

In this study, a genome-wide high throughput sequencing analysis was performed using leaf sheaths of susceptible rice Taichung Native 1 (TN1) and resistant rice IR36 and R476 with and without BPH feeding. A total of 2283 lncRNAs were identified, of which 649 lncRNAs were differentially expressed. During BPH infestation, 84 (120 in total), 52 (70 in total) and 63 (94 in total) of differentially expressed lncRNAs were found only in TN1, IR36 and R476, respectively. Through analyzing their cis-, trans-, and target mimic-activities, not only the lncRNAs targeting resistance genes (NBS-LRR and RLKs) and transcription factors, but also the lncRNAs acting as the targets of the well-studied stress-related miRNAs (miR2118, miR528, and miR1320) in each variety were identified. Before the BPH feeding, 238 and 312 lncRNAs were found to be differentially expressed in TN1 vs. IR36 and TN1 vs. R476, respectively. Among their putative targets, the plant-pathogen interaction pathway was significantly enriched. It is speculated that the resistant rice was in a priming state by the regulation of lncRNAs. Furthermore, the lncRNAs extensively involved in response to BPH feeding were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and the possible regulation networks of the key lncRNAs were constructed. These lncRNAs regulate different pathways that contribute to the basal defense and specific resistance of rice to the BPH.

CONCLUSION

In summary, we identified the specific lncRNAs targeting the well-studied stress-related miRNAs, resistance genes, and transcription factors in each variety during BPH infestation. Additionally, the possible regulating network of the lncRNAs extensively responding to BPH feeding revealed by WGCNA were constructed. These findings will provide further understanding of the regulatory roles of lncRNAs in BPH defense, and lay a foundation for functional research on the candidate lncRNAs.

Pokkali: A Naturally Evolved Salt-Tolerant Rice Shows a Distinguished Set of lncRNAs Possibly Contributing to the Tolerant Phenotype

Pokkali is a strong representation of how stress-tolerant genotypes have evolved due to natural selection pressure. Numerous omics-based investigations have indicated different categories of stress-related genes and proteins, possibly contributing to salinity tolerance in this wild rice. However, a comprehensive study towards understanding the role of long-noncoding RNAs (lncRNAs) in the salinity response of Pokkali has not been done to date. We have identified salt-responsive lncRNAs from contrasting rice genotypes IR64 and Pokkali. A total of 63 and 81 salinity-responsive lncRNAs were differentially expressed in IR64 and Pokkali, respectively. Molecular characterization of lncRNAs and lncRNA-miRNA-mRNA interaction networks helps to explore the role of lncRNAs in the stress response. Functional annotation revealed that identified lncRNAs modulate various cellular processes, including transcriptional regulation, ion homeostasis, and secondary metabolite production. Additionally, lncRNAs were predicted to bind stress-responsive transcription factors, namely ERF, DOF, and WRKY. In addition to salinity, expression profiling was also performed under other abiotic stresses and phytohormone treatments. A positive modulation in TCONS_00035411, TCONS_00059828, and TCONS_00096512 under both abiotic stress and phytohormone treatments could be considered as being of potential interest for the further functional characterization of IncRNA. Thus, extensive analysis of lncRNAs under various treatments helps to delineate stress tolerance mechanisms and possible cross-talk.

Identification of Long Intergenic Noncoding RNAs in Rhizoctonia cerealis following Inoculation of Wheat

Wheat sharp eyespot caused by Rhizoctonia cerealis is primarily a severe threat to worldwide wheat production. Currently, there are no resistant wheat cultivars, and the use of fungicides is the primary method for controlling this disease. Elucidating the mechanisms of R. cerealis pathogenicity can accelerate the pace of the control of this disease. Long intergenic noncoding RNAs (lincRNAs) that function in plant-pathogen interactions might provide a new perspective. We systematically analyzed lincRNAs and identified a total of 1,319 lincRNAs in R. cerealis. We found that lincRNAs are involved in various biological processes, as shown by differential expression analysis and weighted correlation network analysis (WGCNA). Next, one of nine hub lincRNAs in the blue module that was related to infection and growth processes, MSTRG.4380.1, was verified to reduce R. cerealis virulence on wheat by a host-induced gene silencing (HIGS) assay. Following that, RNA sequencing (RNA-Seq) analysis revealed that the significantly downregulated genes in the MSTRG.4380.1 knockdown lines were associated mainly with infection-related processes, including hydrolase, transmembrane transporter, and energy metabolism activities. Additionally, 23 novel microRNAs (miRNAs) were discovered during small RNA (sRNA) sequencing (sRNA-Seq) analysis of MSTRG.4380.1 knockdown, and target prediction of miRNAs suggested that MSTRG.4380.1 does not act as a competitive endogenous RNA (ceRNA). This study performed the first genome-wide identification of R. cerealis lincRNAs and miRNAs. It confirmed the involvement of a lincRNA in the infection process, providing new insights into the mechanism of R. cerealis infection and offering a new approach for protecting wheat from R. cerealis. IMPORTANCE Rhizoctonia cerealis, the primary causal agent of wheat sharp eyespot, has caused significant losses in worldwide wheat production. Since no resistant wheat cultivars exist, chemical control is the primary method. However, this approach is environmentally unfriendly and costly. RNA interference (RNAi)-mediated pathogenicity gene silencing has been proven to reduce the growth of Rhizoctonia and provides a new perspective for disease control. Recent studies have shown that lincRNAs are involved in various biological processes across species, such as biotic and abiotic stresses. Therefore, verifying the function of lincRNAs in R. cerealis is beneficial for understanding the infection mechanism. In this study, we reveal that lincRNAs could contribute to the virulence of R. cerealis, which provides new insights into controlling this pathogen.

Plant long noncoding RNAs: Recent progress in understanding their roles in growth, development, and stress responses

Long noncoding RNA (lncRNA) transcripts are longer than 200 nt and are not translated into proteins. LncRNAs function in a wide variety of processes in plants and animals, but, perhaps because of their lower expression and conservation levels, plant lncRNAs had attracted less attention than protein-coding mRNAs. Now, recent studies have made remarkable progress in identifying lncRNAs and understanding their functions. In this review, we discuss a number of lncRNAs that have important functions in growth, development, reproduction, responses to abiotic stresses, and regulation of disease and insect resistance in plants. Additionally, we describe the known mechanisms of action of plant lncRNAs according to their origins within the genome. This review thus provides a guide for identifying and functionally characterizing new lncRNAs in plants.

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.

Long Non-Coding RNA lncWOX11a Suppresses Adventitious Root Formation of Poplar by Regulating the Expression of PeWOX11a

Long non-coding RNAs (lncRNAs), a class of poorly conserved transcripts without protein-encoding ability, are widely involved in plant organogenesis and stress responses by mediating the transmission and expression of genetic information at the transcriptional, posttranscriptional, and epigenetic levels. Here, we cloned and characterized a novel lncRNA molecule through sequence alignment, Sanger sequencing, transient expression in protoplasts, and genetic transformation in poplar. lncWOX11a is a 215 bp transcript located on poplar chromosome 13, ~50 kbp upstream of PeWOX11a on the reverse strand, and the lncRNA may fold into a series of complex stem–loop structures. Despite the small open reading frame (sORF) of 51 bp within lncWOX11a, bioinformatics analysis and protoplast transfection revealed that lncWOX11a has no protein-coding ability. The overexpression of lncWOX11a led to a decrease in the quantity of adventitious roots on the cuttings of transgenic poplars. Further, cis-regulatory module prediction and CRISPR/Cas9 knockout experiments with poplar protoplasts demonstrated that lncWOX11a acts as a negative regulator of adventitious rooting by downregulating the WUSCHEL-related homeobox gene WOX11, which is supposed to activate adventitious root development in plants. Collectively, our findings imply that lncWOX11a is essential for modulating the formation and development of adventitious roots.

Comparative transcriptome-wide identification and differential expression of genes and lncRNAs in rice near-isogenic line (KW-Bph36-NIL) in response to BPH feeding

Brown planthopper (BPH) is the most devastating pest of rice in Asia, causing substantial yield losses and has become a challenging task to be controlled under field conditions. Although extensive measures have been taken over the past decades, which resulted in the evolution of new resistant BPH strains. Therefore, besides other possible approaches, equipping host plants with resistant genes is the most effective and environment-friendly technique for BPH control. Here, we systematically analyzed transcriptome changes in the susceptible rice variety Kangwenqingzhan (KW) and the resistant near-isogenic line (NIL) KW-Bph36-NIL, through RNA-seq, depicting the differential expression profiles of mRNAs and long non-coding RNAs (lncRNAs) in rice before and after BPH feeding. We observed a proportion of genes (1.48%) and (2.74%) were altered in KW and NIL, respectively, indicating different responses of rice strains against BPH feeding. Nevertheless, we characterized 384 differentially expressed long non-coding RNAs (DELs) that can be impacted by the two strains by alternatively changing the expression patterns of the respective coding genes, suggesting their certain involvement in response to BPH feeding. In BPH invasion, KW and NIL responded differently by modifying the synthesis, storage, and transformation of intracellular substances, adjusting the nutrient accumulation and utilization inside and outside the cells. In addition, NIL expressed stronger resistance by acutely up-regulating genes and other transcription factors related to stress resistance and plant immunity. Altogether, our study elaborates valuable insights into the genome-wide DEGs and DELs expression profiles of rice under BPH invasion by high throughput sequencing and further suggests that NILs can be utilized in BPH resistance breeding programs in developing high-resistance rice lines.

Construction of ceRNA and m6A-related lncRNA networks associated with anti-inflammation of AdipoAI

Background

Adiponectin (APN) is an endogenous adipokine secreted from adipocytes that exerts anti-inflammatory properties. AdipoAI is an orally active adiponectin receptor agonist identified by our group that can emulate APN's anti-inflammatory properties through mechanisms that are not fully understood. LncRNAs, a type of noncoding RNA more than 200 bp in length, have been demonstrated to have abundant biological functions, including in anti-inflammatory responses.

Materials and Result

In the current study, we performed a lncRNA microarray in LPS-induced Raw264.7 cells that were prestimulated with AdipoAI and screened 110 DElncRNAs and 190 DEmRNAs. Enrichment analyses were conducted on total mRNAs and DEmRNAs, including GSVA, ssGSEA, GO/KEGG, GSEA, and PPI analysis. Among all these processes, endocytosis was significantly enriched. A coexpression analysis was built based on DElncRNAs and DEmRNAs. Then, using TargetScan and miRwalk to predict related microRNAs of DElncRNAs and DEmRNAs, respectively, we established competing endogenous RNA (ceRNA) networks including 54 mRNAs from 8 GO items. Furthermore, 33 m6A methylation-related marker genes were obtained from a previous study and used for the construction of an m6A-related lncRNA network by coexpression analysis. We identified FTO as the hub gene of the network and 14 lncRNAs that interacted with it. The expression levels of 10 lncRNAs selected from ceRNA and FTO-related lncRNA networks were validated with qRT‒PCR. Finally, macrophage phenotype scores showed that AdipoAI could attenuate the M2b and M2c polarization of macrophages and correlate with the above lncRNAs.

Conclusion

Our work reveals that lncRNAs might be involved in the anti-inflammation process of AdipoAI in LPS-induced macrophages through the ceRNA network and the epigenetic regulation of m6A. Mechanistically, these lncRNAs associated with AdipoAI might be related to endocytosis and polarization in macrophages and provide new candidates for the anti-inflammatory application of APN and its receptor agonist.

Full-length transcriptome sequencing reveals the molecular mechanism of monoterpene and sesquiterpene biosynthesis in Cinnamomum burmannii

Essential oil of Cinnamomum burmannii is rich in monoterpenes and sesquiterpenes and is widely used in cosmetics and medicines. Knowledge about the enzymes that catalyze the formation of monoterpenes and sesquiterpenes in C. burmannii is insufficient. Therefore, anatomy observation of C. burmannii at the four developmental stages (7 days, CBS1; 14 days, CBS2; 21 days, CBS3, and 28 days, CBS4) were conducted to elucidate the origins of essential oil production. Twelve full-length transcriptomes of C. burmannii leaves at the four stages were generated using Oxford Nanopore Technologies. GC-MS analysis revealed 15 monoterpene and sesquiterpenes dramatically increased from CBS1 to CBS4. A weighted correlation network analysis (WGCNA) in association and differentially expressed genes across four developmental stages were performed. A total of 44 differentially expressed genes (DEGs) were involved in terpenoid syntheses during leaf development. Among them, the DEGs of the mevalonate acid (MVA) pathway were predominantly expressed at CBS1, while those of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway showed increased expression from CBS2 to CBS4. Besides, fourteen genes were associated with monoterpene synthesis and nine with sesquiterpene synthesis. Functions of these DEGs were further predicted with regard to gene expression profile and phylogenetic relationship with those characterized in previous studies. In addition, 922 long noncoding RNAs (lncRNAs) were detected, of which twelve were predicted to regulate monoterpene and sesquiterpene biosynthesis. The present study provided new insights the molecular mechanisms of monoterpenoid and sesquiterpenoid syntheses of C. burmannii.

A Medicago truncatula lncRNA MtCIR1 negatively regulates response to salt stress

A lncRNA MtCIR1 negatively regulates the response to salt stress in Medicago truncatula seed germination by modulating seedling growth and ABA metabolism and signaling by enhancing Na⁺ accumulation. Increasing evidence suggests that long non-coding RNAs (lncRNAs) are involved in the regulation of plant tolerance to varying abiotic stresses. A large number of lncRNAs that are responsive to abiotic stress have been identified in plants; however, the mechanisms underlying the regulation of plant responses to abiotic stress by lncRNAs are largely unclear. Here, we functionally characterized a salt stress-responsive lncRNA derived from the leguminous model plant M. truncatula, referred to as MtCIR1, by expressing MtCIR1 in Arabidopsis thaliana in which no such homologous sequence was observed. Expression of MtCIR1 rendered seed germination more sensitive to salt stress by enhanced accumulation of abscisic acid (ABA) due to suppressing the expression of the ABA catabolic enzyme CYP707A2. Expression of MtCIR1 also suppressed the expression of genes associated with ABA receptors and signaling. The ABA-responsive gene AtPGIP2 that was involved in degradation of cell wall during seed germination was up-regulated by expressing MtCIR1. On the other hand, expression of MtCIR1 in Arabidopsis thaliana enhanced foliar Na⁺ accumulation by down-regulating genes encoding Na⁺ transporters, thus rendering the transgenic plants more sensitive to salt stress. These results demonstrate that the M. truncatula lncRNA MtCIR1 negatively regulates salt stress response by targeting ABA metabolism and signaling during seed germination and foliar Na⁺ accumulation by affecting Na⁺ transport under salt stress during seedling growth. These novel findings would advance our knowledge on the regulatory roles of lncRNAs in response of plants to salt stress.

BPL3 binds the long non-coding RNA nalncFL7 to suppress FORKED-LIKE7 and modulate HAI1-mediated MPK3/6 dephosphorylation in plant immunity

RNA-binding proteins (RBPs) participate in a diverse set of biological processes in plants, but their functions and underlying mechanisms in plant-pathogen interactions are largely unknown. We previously showed that Arabidopsis thaliana BPA1-LIKE PROTEIN3 (BPL3) belongs to a conserved plant RBP family and negatively regulates reactive oxygen species (ROS) accumulation and cell death under biotic stress. In this study, we demonstrate that BPL3 suppresses FORKED-LIKE7 (FL7) transcript accumulation and raises levels of the cis-natural antisense long non-coding RNA (lncRNA) of FL7 (nalncFL7). FL7 positively regulated plant immunity to Phytophthora capsici while nalncFL7 negatively regulated resistance. We also showed that BPL3 directly binds to and stabilizes nalncFL7. Moreover, nalncFL7 suppressed accumulation of FL7 transcripts. Furthermore, FL7 interacted with HIGHLY ABA-INDUCED PP2C1 (HAI1), a type 2C protein phosphatase, and inhibited HAI1 phosphatase activity. By suppressing HAI1 activity, FL7 increased the phosphorylation levels of MITOGEN-ACTIVATED PROTEIN KINASE 3 (MPK3) and MPK6, thus enhancing immunity responses. BPL3 and FL7 are conserved in all plant species tested, but the BPL3-nalncFL7-FL7 cascade was specific to the Brassicaceae. Thus, we identified a conserved BPL3-nalncFL7-FL7 cascade that coordinates plant immunity.

Cucumber STACHYOSE SYNTHASE is regulated by its cis-antisense RNA asCsSTS to balance source-sink carbon partitioning

Photosynthate partitioning between source and sink is a key determinant of crop yield. In contrast to sucrose-transporting plants, cucumber (Cucumis sativus) plants mainly transport stachyose and stachyose synthase (CsSTS) synthesizes stachyose in the vasculature for loading. Therefore, CsSTS is considered a key regulator of carbon partitioning. We found that CsSTS expression and CsSTS enzyme activity were upregulated in the vasculature and downregulated in mesophyll tissues at fruiting. In situ hybridization and tissue enrichment experiments revealed that a cis-natural antisense noncoding transcript of CsSTS, named asCsSTS, is mainly expressed in mesophyll tissues. In vitro overexpression (OE), RNA interference (RNAi), and dual luciferase reporter experiments indicated that CsSTSs are negatively regulated by asCsSTS. Fluorescence in situ hybridization revealed that asCsSTS transcript localized in leaf cytoplasm, indicating that the regulation of CsSTS by asCsSTS is a posttranscriptional process. Further investigation revealed that this regulation occurred by reducing CsSTS transcript stability through a DICER-like protein-mediated pathway. Chemically induced OE and RNAi of asCsSTS led to promotion or inhibition, respectively, of assimilate export from leaves and altered fruit growth rates. Our results suggest that the regulation of CsSTSs between the mesophyll and vasculature reduces sugar storage in mesophyll tissue and promotes assimilate export from the leaf when the plant carries fruit.

A long non-coding RNA PelncRNA1 is involved in Phyllostachys edulis response to UV-B stress

Phyllostachys edulis (moso bamboo) is China's most widespread bamboo species, with significant economic and ecological values. Long non-coding RNA (lncRNA) is a type of regulatory RNA that is longer than 200 nucleotides and incapable of encoding proteins, and is frequently involved in regulating biotic and abiotic stress and plant development. However, the biological functions of lncRNA in moso bamboo are unknown. In this study, a lncRNA (named PelncRNA1) differentially expressed following UV-B treatment was discovered in the whole transcriptome sequencing database of moso bamboo. The target genes were filtered and defined by correlation analysis of PelncRNA1 and gene expression pattern. The expression levels of PelncRNA1 and its target genes were verified using qRT-PCR. The results demonstrated that the expression levels of PelncRNA1 and its target genes increased during UV-B treatment. In Arabidopsis transgenic seedlings and moso bamboo protoplasts, PelncRNA1 was discovered to influence the expression of its target genes when overexpressed. In addition, transgenic Arabidopsis showed higher tolerance to UV-B stress. These results suggest that PelncRNA1 and its target genes are involved in the response of moso bamboo to UV-B stress. The novel findings would contribute to our understanding of how lncRNAs regulate the response to abiotic stresses in moso bamboo.

Epitranscriptome changes triggered by ammonium nutrition regulate the proteome response of maritime pine roots

Epitranscriptome constitutes a gene expression checkpoint in all living organisms. Nitrogen is an essential element for plant growth and development that influences gene expression at different levels such as epigenome, transcriptome, proteome, and metabolome. Therefore, our hypothesis is that changes in the epitranscriptome may regulate nitrogen metabolism. In this study, epitranscriptomic modifications caused by ammonium nutrition were monitored in maritime pine roots using Oxford Nanopore Technology. Transcriptomic responses mainly affected transcripts involved in nitrogen and carbon metabolism, defense, hormone synthesis/signaling, and translation. Global detection of epitranscriptomic marks was performed to evaluate this posttranscriptional mechanism in un/treated seedlings. Increased N⁶-methyladenosine (m⁶A) deposition in the 3'-UTR was observed in response to ammonium, which seems to be correlated with poly(A) lengths and changes in the relative abundance of the corresponding proteins. The results showed that m⁶A deposition and its dynamics seem to be important regulators of translation under ammonium nutrition. These findings suggest that protein translation is finely regulated through epitranscriptomic marks likely by changes in mRNA poly(A) length, transcript abundance and ribosome protein composition. An integration of multiomics data suggests that the epitranscriptome modulates responses to nutritional, developmental and environmental changes through buffering, filtering, and focusing the final products of gene expression.

A survey of the full-length transcriptome of Gracilariopsis lemaneiformis using single-molecule long-read sequencing

BACKGROUND

Posttranscriptional processing of precursor mRNAs contributes to transcriptome and protein diversity and gene regulatory mechanisms in eukaryotes. However, this posttranscriptional mechanism has not been studied in the marine macroalgae Gracilariopsis lemaneiformis, which is the most cultivated red seaweed species in China.

RESULTS

In the present study, third-generation sequencing (Pacific Biosciences single-molecule real-time long-read sequencing, SMRT-Seq) was used to sequence the full-length transcriptome of G. lemaneiformis to identify alternatively spliced transcripts and alternative polyadenylation (APA) sites in this species. RNAs were isolated from G. lemaneiformis under various treatments including abiotic stresses and exogenous phytohormones, and then equally pooled for SMRT-Seq. In summary, 346,544 full-length nonchimeric reads were generated, from which 13,630 unique full-length transcripts were obtained in G. lemaneiformis. Compared with the known splicing events in the gene models, more than 3000 new alternative splicing (AS) events were identified in the SMRT-Seq reads. Additionally, 810 genes were found to have poly (A) sites and 91 microRNAs (miRNAs), 961 long noncoding RNAs and 1721 novel genes were identified in G. lemaneiformis. Moreover, validation experiments showed that abiotic stresses and phytohormones could induce some specific AS events, especially intron retain isoforms, cause some alterations to the relative ratios of transcripts annotated to the same gene, and generate novel 3' ends because of differential APA. The growth of G. lemaneiformis was inhibited by Cu stress, while this inhibition was alleviated by ACC treatment. RNA-Seq analysis further revealed that 211 differential alternative splicing (DAS) events and 142 DAS events was obtained in CK vs Cu and Cu vs Cu + ACC, respectively, suggesting that AS of functional genes could be regulated by Cu stress and ACC. Compared with Cu stress, the expression of transcripts with DAS events mainly involved in the carbon fixation in photosynthetic organisms and oxidative phosphorylation pathway was upregulated in Cu + ACC treatment, revealing that ACC alleviated the growth inhibition by Cu stress by increasing carbon fixation and oxidative phosphorylation.

CONCLUSIONS

Our results provide the first comprehensive picture of the full-length transcriptome and posttranscriptional mechanism in red macroalgae, including transcripts that appeared in the presence of common abiotic stresses and phytohormones, which will improve the gene annotations of Gracilariopsis and contribute to the study of gene regulation in this important cultivated seaweed.

Long non-coding RNAs as the regulatory hubs in rice response to salt stress

Salinity seriously constrains growth and fertility of rice worldwide. Long non-coding RNAs (lncRNAs) play crucial roles in plant abiotic stress response. However, salt responsive lncRNAs are poorly understood in rice. Herein, salt responsive lncRNAs (DE-lncRNAs) were identified in FL478 (salt tolerant) compared to its susceptible parent (IR29) using RNA-seq in root tissues at seedling stage. In FL478 and IR29, 8724 and 9235 transcripts with length of > 200 bp were nominated as potential lncRNAs, respectively. Rigorous filtering left four (in FL478) and nine (in IR29) DE-lncRNAs with only 2 DE-lncRNAs in common. ATAC-seq data showed that the genomic regions of all four lncRNAs in FL478 and 6/9 in IR29 are significantly accessible for transcription. Weighted correlation network analysis (WGCNA) revealed that lncRNA.2-FL was highly correlated with 173 mRNAs as trans-targets and a gene encoding pentatricopeptide repeat (PPR) protein was predicted as cis-target of lncRNA.2-FL. In silico mutagenesis analysis proposed the same transcription factor binding sites (TFBSs) in vicinity of the trans- and cis-regulatory target genes of lncRNA.2-FL, which significantly affect their transcription start site (TSS). This study provides new insights into involvement of the DE-lncRNAs in rice response to salt stress. Among them, lncRNA.2-FL may play a significant regulatory role in the salt stress tolerance of FL478.

Whole transcriptome analysis of HCT-8 cells infected by Cryptosporidium parvum

BACKGROUND

Cryptosporidium species are zoonotic protozoans that are important causes of diarrhoeal disease in both humans and animals. Non-coding RNAs (ncRNAs) play an important role in the innate immune defense against Cryptosporidium infection, but the underlying molecular mechanisms in the interaction between human ileocecal adenocarcinoma (HCT-8) cells and Cryptosporidium species have not been entirely revealed.

METHODS

The expression profiles of messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs) in the early phase of infection of HCT-8 cells with Cryptosporidium parvum and at 3 and 12 h post infection were analyzed using the RNA-sequencing technique. The biological functions of differentially expressed RNAs (dif-RNAs) were discovered through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The targeting relationships between three ncRNAs and mRNAs were analyzed using bioinformatics methods, followed by building a competing endogenous RNA (ceRNA) regulatory network centered on miRNAs.

RESULTS

After strictly filtering the raw data, our analysis revealed 393 dif-lncRNAs, 69 dif-miRNAs and 115 dif-mRNAs at 3 hpi, and 450 dif-lncRNAs, 129 dif-miRNAs, 117 dif-mRNAs and one dif-circRNA at 12 hpi. Of these, 94 dif-lncRNAs, 24 dif-miRNAs and 22 dif-mRNAs were detected at both post-infection time points. Eleven dif-lncRNAs, seven dif-miRNAs, eight dif-mRNAs and one circRNA were randomly selected and confirmed using the quantitative real-time PCR. Bioinformatics analyses showed that the dif-mRNAs were significantly enriched in nutritional absorption, metabolic processes and metabolism-related pathways, while the dif-lncRNAs were mainly involved in the pathways related to the infection and pathogenicity of C. parvum (e.g. tight junction protein) and immune-related pathways (e.g. cell adhesion molecules). In contrast, dif-miRNAs and dif-circRNA were significantly enriched in apoptosis and apoptosis-related pathways. Among the downregulated RNAs, the miRNAs has-miR-324-3p and hsa-miR-3127-5p appear to be crucial miRNAs which could negatively regulate circRNA, lncRNA and mRNA.

CONCLUSIONS

The whole transcriptome profiles of HCT-8 cells infected with C. parvum were obtained in this study. The results of the GO and KEGG pathway analyses suggest significant roles for these dif-RNAs during the course of C. parvum infection. A ceRNA regulation network containing miRNA at its center was constructed for the first time, with hsa-miR-324-3p and hsa-miR-3127-5p being the crucial miRNAs. These findings provide novel insights into the responses of human intestinal epithelial cells to C. parvum infection.

Analysis of Long Non-Coding RNA-Mediated Regulatory Networks of Plutella xylostella in Response to Metarhizium anisopliae Infection

Long non-coding RNAs (lncRNAs) represent a diverse class of RNAs that are structurally similar to messenger RNAs (mRNAs) but do not encode proteins. Growing evidence suggests that in response to biotic and abiotic stresses, the lncRNAs play crucial regulatory roles in plants and animals. However, the potential role of lncRNAs during fungal infection has yet to be characterized in Plutella xylostella, a devastating pest of cruciferous crops. In the current study, we performed a strand-specific RNA sequencing of Metarhizium anisopliae-infected (Px36hT, Px72hT) and uninfected (Px36hCK, Px72hCK) P. xylostella fat body tissues. Comprehensive bioinformatic analysis revealed a total of 5665 and 4941 lncRNAs at 36 and 72-h post-infection (hpi), including 563 (Px36hT), 532 (Px72hT) known and 5102 (Px36hT), 4409 (Px72hT) novel lncRNA transcripts. These lncRNAs shared structural similarities with their counterparts in other species, including shorter exon and intron length, fewer exon numbers, and a lower expression profile than mRNAs. LncRNAs regulate the expression of neighboring protein-coding genes by acting in a cis and trans manner. Functional annotation and pathway analysis of cis-acting lncRNAs revealed their role in several immune-related genes, including Toll, serpin, transferrin, βGRP etc. Furthermore, we identified multiple lncRNAs acting as microRNA (miRNA) precursors. These miRNAs can potentially regulate the expression of mRNAs involved in immunity and development, suggesting a crucial lncRNA-miRNA-mRNA complex. Our findings will provide a genetic resource for future functional studies of lncRNAs involved in P. xylostella immune responses to M. anisopliae infection and shed light on understanding insect host-pathogen interactions.

Noncoding RNAs responsive to nitric oxide and their protein-coding gene targets shed light on root hair formation in Arabidopsis thaliana

An overview of the total Arabidopsis thaliana transcriptome, described previously by our research group, pointed some noncoding RNA (ncRNA) as participants in the restoration of hair-root phenotype in A. thaliana rhd6 mutants, leading us to a deeper investigation. A transcriptional gene expression profiling of seedling roots was performed aiming to identify ncRNA responsive to nitric oxide (GSNO) and auxin (IAA), and their involvement in root hair formation in the rhd6 null mutant. We identified 3,631 ncRNAs, including new ones, in A. thaliana and differential expression (DE) analysis between the following: 1) GSNO-treated rhd6 vs. untreated rhd6, 2) IAA-treated rhd6 vs. untreated rhd6, 3) GSNO-treated rhd6 vs. IAA-treated rhd6, and 4) WS-2 vs. untreated rhd6 detected the greatest number of DE genes in GSNO-treated rhd6. We detected hundreds of in silico interactions among ncRNA and protein-coding genes (PCGs), highlighting MIR5658 and MIR171 precursors highly upregulated in GSNO-treated rhd6 and wild type, respectively. Those ncRNA interact with many DE PCGs involved in hormone signaling, cell wall development, transcription factors, and root hair formation, becoming candidate genes in cell wall modulation and restoration of root hair phenotype by GSNO treatment. Our data shed light on how GSNO modulates ncRNA and their PCG targets in A. thaliana root hair formation.

Beat the heat: plant- and microbe-mediated strategies for crop thermotolerance

Heat stress (HS) affects plant growth and development, and reduces crop yield. To combat HS, plants have evolved several sophisticated strategies. The primary HS response in plants involves the activation of heat-shock transcription factors and heat-shock proteins (HSPs). Plants also deploy more advanced epigenetic mechanisms in response to recurring HS conditions. In addition, beneficial microbes can reprogram the plant epitranscriptome to induce thermotolerance, and have the potential to improve crop yield productivity by mitigating HS-induced inhibition of growth and development. We summarize the latest advances in plant epigenetic regulation and highlight microbe-mediated thermotolerance in plants.

Identification and Characterization of Long Non-coding RNA in Tomato Roots Under Salt Stress

As one of the most important vegetable crops in the world, the production of tomatoes was restricted by salt stress. Therefore, it is of great interest to analyze the salt stress tolerance genes. As the non-coding RNAs (ncRNAs) with a length of more than 200 nucleotides, long non-coding RNAs (lncRNAs) lack the ability of protein-coding, but they can play crucial roles in plant development and response to abiotic stresses by regulating gene expression. Nevertheless, there are few studies on the roles of salt-induced lncRNAs in tomatoes. Therefore, we selected wild tomato Solanum pennellii (S. pennellii) and cultivated tomato M82 to be materials. By high-throughput sequencing, 1,044 putative lncRNAs were identified here. Among them, 154 and 137 lncRNAs were differentially expressed in M82 and S. pennellii, respectively. Through functional analysis of target genes of differentially expressed lncRNAs (DE-lncRNAs), some genes were found to respond positively to salt stress by participating in abscisic acid (ABA) signaling pathway, brassinosteroid (BR) signaling pathway, ethylene (ETH) signaling pathway, and anti-oxidation process. We also construct a salt-induced lncRNA-mRNA co-expression network to dissect the putative mechanisms of high salt tolerance in S. pennellii. We analyze the function of salt-induced lncRNAs in tomato roots at the genome-wide levels for the first time. These results will contribute to understanding the molecular mechanisms of salt tolerance in tomatoes from the perspective of lncRNAs.

Genome-Wide Identification of Powdery Mildew Responsive Long Non-Coding RNAs in Cucurbita pepo

Cucurbita pepo L. is an essential economic vegetable crop worldwide, and its production is severely affected by powdery mildew (PM). However, our understanding of the molecular mechanism of PM resistance in C. pepo is very limited. Long non-coding RNAs (lncRNAs) play an important role in regulating plant responses to biotic stress. Here, we systematically identified 2,363 reliably expressed lncRNAs from the leaves of PM-susceptible (PS) and PM-resistant (PR) C. pepo. The C. pepo lncRNAs are shorter in length and expressed at a lower level than the protein-coding transcripts. Among the 2,363 lncRNAs, a total of 113 and 146 PM-responsive lncRNAs were identified in PS and PR, respectively. Six PM-responsive lncRNAs were predicted as potential precursors of microRNAs (miRNAs). In addition, 58 PM-responsive lncRNAs were predicted as targets of miRNAs and one PM-responsive lncRNA was predicted as an endogenous target mimic (eTM). Furthermore, a total of 5,200 potential cis target genes and 5,625 potential trans target genes were predicted for PM-responsive lncRNAs. Functional enrichment analysis showed that these potential target genes are involved in different biological processes, such as the plant-pathogen interaction pathway, MAPK signaling pathway, and plant hormone signal transduction pathway. Taken together, this study provides a comprehensive view of C. pepo lncRNAs and explores the putative functions of PM-responsive lncRNAs, thus laying the foundation for further study of the regulatory mechanisms of lncRNAs responding to PM.

Genome-wide identification and characterization of long noncoding RNAs during peach (Prunus persica) fruit development and ripening

LncRNAs represent a class of RNA transcripts of more than 200 nucleotides (nt) in length without discernible protein-coding potential. The expression levels of lncRNAs are significantly affected by stress or developmental cues. Recent studies have shown that lncRNAs participate in fruit development and ripening processes in tomato and strawberry; however, in other fleshy fruits, the association between lncRNAs and fruit ripening remains largely elusive. Here, we constructed 9 ssRNA-Seq libraries from three different peach (Prunus persica) fruit developmental stages comprising the first and second exponential stages and the fruit-ripening stage. In total, 1500 confident lncRNAs from 887 loci were obtained according to the bioinformatics analysis. The lncRNAs identified in peach fruits showed distinct characteristics compared with protein-coding mRNAs, including lower expression levels, lower complexity of alternative splicing, shorter isoforms and smaller numbers of exons. Expression analysis identified 575 differentially expressed lncRNAs (DELs) classified into 6 clusters, among which members of Clusters 1, 2, 4 and 5 were putatively associated with fruit development and ripening processes. Quantitative real-time PCR revealed that the DELs indeed had stage-specific expression patterns in peach fruits. GO and KEGG enrichment analysis revealed that DELs might be associated with fruit-ripening-related physiological and metabolic changes, such as flavonoid biosynthesis, fruit texture softening, chlorophyll breakdown and aroma compound accumulation. Finally, the similarity analysis of lncRNAs within different plant species indicated the low sequence conservation of lncRNAs. Our study reports a large number of fruit-expressed lncRNAs and identifies fruit development phase-specific expressed lncRNA members, which highlights their potential functions in fruit development and ripening processes and lays the foundations for future functional research.

Whole-Transcriptome Analysis Reveals Autophagy Is Involved in Early Senescence of zj-es Mutant Rice

Senescence is a necessary stage of plant growth and development, and the early senescence of rice will lead to yield reduction and quality decline. However, the mechanisms of rice senescence remain obscure. In this study, we characterized an early-senescence rice mutant, designated zj-es (ZheJing-early senescence), which was derived from the japonica rice cultivar Zhejing22. The mutant zj-es exhibited obvious early-senescence phenotype, such as collapsed chloroplast, lesions in leaves, declined fertility, plant dwarf, and decreased agronomic traits. The ZJ-ES gene was mapped in a 458 kb-interval between the molecular markers RM5992 and RM5813 on Chromosome 3, and analysis suggested that ZJ-ES is a novel gene controlling rice early senescence. Subsequently, whole-transcriptome RNA sequencing was performed on zj-es and its wild-type rice to dissect the underlying molecular mechanism for early senescence. Totally, 10,085 differentially expressed mRNAs (DEmRNAs), 1,253 differentially expressed lncRNAs (DElncRNAs), and 614 differentially expressed miRNAs (DEmiRNAs) were identified, respectively, in different comparison groups. Based on the weighted gene co-expression network analysis (WGCNA), the co-expression turquoise module was found to be the key for the occurrence of rice early senescence. Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 14 lncRNAs possibly regulated 16 co-expressed mRNAs through 8 miRNAs, and enrichment analysis showed that most of the DEmRNAs and the targets of DElncRNAs and DEmiRNAs were involved in reactive oxygen species (ROS)-triggered autophagy-related pathways. Further analysis showed that, in zj-es, ROS-related enzyme activities were markedly changed, ROS were largely accumulated, autophagosomes were obviously observed, cell death was significantly detected, and lesions were notably appeared in leaves. Totally, combining our results here and the remaining research, we infer that ROS-triggered autophagy induces the programmed cell death (PCD) and its coupled early senescence in zj-es mutant rice.

Plasmodium manipulates the expression of host long non-coding RNA during red blood cell intracellular infection

BACKGROUND

Parasites interact with their host through "direct" and/or "indirect" mechanisms. Plasmodium, for example, either mediates direct physical interactions with host factors or triggers the immune system of the host indirectly, leading to changes in infectious outcomes. Long non-coding RNAs (lncRNAs) participate in regulating biological processes, especially host-pathogen interactions. However, research on the role of host lncRNAs during Plasmodium infection is limited.

METHODS

A RNA sequencing method (RNA-seq) was used to confirm the differential expression profiles of lncRNAs in Plasmodium yeolii 17XL (P.y17XL)-infected BALB/c mice. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to elucidate the potential functions of Plasmodium-induced genes. Subsequently, the effect of specific lncRNAs on the modulation of immune-related signaling pathways in malaria was determined by fluorescence-activated cell sorting, western blot and enzyme-linked immunosorbent assay.

RESULTS

The data showed that in P.y17XL-infected BALB/c mice, Plasmodium upregulated the expression of 132 lncRNAs and downregulated the expression of 159 lncRNAs. Differentially expressed lncRNAs clearly associated with malaria infection were annotated, including four novel dominant lncRNAs: ENMSUSG00000111521.1, XLOC_038009, XLOC_058629 and XLOC_065676. GO and KEGG pathway analyses demonstrated that these four differentially expressed lncRNAs were associated with co-localized/co-expressed protein-coding genes that were totally enriched in malaria and with the transforming growth factor beta (TGF-β) signaling pathway. Using the models of P.y17XL-infected BALB/c mice, data certified that the level of TGF-β production and activation of TGF-β/Smad_2/3 signaling pathway were obviously changed in malaria infection.

CONCLUSIONS

These differentially expressed immune-related genes were deemed to have a role in the process of Plasmodium infection in the host via dendritic/T regulatory cells and the TGF-β/Smad_2/3 signaling pathway. The results of the present study confirmed that Plasmodium infection-induced lncRNA expression is a novel mechanism used by Plasmodium parasites to modify host immune signaling. These results further enhance current understanding of the interaction between Plasmodium and host cells.

Expression Profiles and Characteristics of Apple lncRNAs in Roots, Phloem, Leaves, Flowers, and Fruit

LncRNAs impart crucial effects on various biological processes, including biotic stress responses, abiotic stress responses, fertility and development. The apple tree is one of the four major fruit trees in the world. However, lncRNAs's roles in different tissues of apple are unknown. We identified the lncRNAs in five tissues of apples including the roots, phloem, leaves, flowers, and fruit, and predicted the intricate regulatory networks. A total of 9440 lncRNAs were obtained. LncRNA target prediction revealed 10,628 potential lncRNA-messenger RNA (mRNA) pairs, 9410 pairs functioning in a cis-acting fashion, and 1218 acting in a trans-acting fashion. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to photosynthesis-antenna proteins, single-organism metabolic process and glutathione metabolism. Additionally, a total of 88 lncRNAs have various functions related to microRNAs (miRNAs) as miRNA precursors. Interactions between lncRNAs and miRNAs were predicted, 1341 possible interrelations between 187 mdm-miRNAs and 174 lncRNAs (1.84%) were identified. MSTRG.121644.5, MSTRG.121644.8, MSTRG.2929.2, MSTRG.3953.2, MSTRG.63448.2, MSTRG.9870.2, and MSTRG.9870.3 could participate in the functions in roots as competing endogenous RNAs (ceRNAs). MSTRG.11457.2, MSTRG.138614.2, and MSTRG.60895.2 could adopt special functions in the fruit by working with miRNAs. A further analysis showed that different tissues formed special lncRNA-miRNA-mRNA networks. MSTRG.60895.2-mdm-miR393-MD17G1009000 may participate in the anthocyanin metabolism in the fruit. These findings provide a comprehensive view of potential functions for lncRNAs, corresponding target genes, and related lncRNA-miRNA-mRNA networks, which will increase our knowledge of the underlying development mechanism in apple.

Comprehensive Transcriptome Analysis of Stem-Differentiating Xylem Upon Compression Stress in Cunninghamia Lanceolata

Compression wood (CW) in gymnosperm brings great difficulties to wood industry using wood as raw materials since CW presents special wood structure and have different physical and chemical properties from those of normal wood (NW). Chinese fir (Cunninghamia lanceolata) is widely distributed in China. However, global transcriptome profiling of coding and long non-coding RNA in response to compression stress has not been reported in the gymnosperm species. In this study, we revealed that CW in Chinese fir exhibited distinct morphology and cytology properties compared with those of NW, including high lignin content, thick and round tracheid cells. Furthermore, we combined both PacBio long-read SMRT sequencing (Iso-Seq) and Illumina short-read RNA-Seq to reveal the transcriptome in stem-differentiating xylem (SDX) under different time points (2, 26, and 74 h) upon compression stress in NW, CW, and OW (opposite wood), respectively. Iso-Seq was successfully assembled into 41,253 de-novo full-length transcriptome reference (average length 2,245 bp). Moreover, there were striking differences in expression upon compression stress, which were involved 13 and 7 key enzyme genes in the lignin and cellulose synthesis, respectively. Especially, we revealed 11 secondary growth-related transcription factors show differential expression under compression stress, which was further validated by qRT-PCR. Finally, the correlation between 6,533 differentially expressed coding genes and 372 differentially expressed long non-coding RNAs (lncRNAs) indicates that these lncRNAs may affect cell wall biogenesis and xyloglucan metabolism. In conclusion, our results provided comprehensive cytology properties and full-length transcriptome profiling of wood species upon compression stress. Especially we explored candidate genes, including both coding and long non-coding genes, and provided a theoretical basis for further research on the formation mechanism of CW in gymnosperm Chinese fir.

Non-Coding RNA Analyses of Seasonal Cambium Activity in Populus tomentosa

Non-coding RNA, known as long non-coding RNA (lncRNA), circular RNA (circRNA) and microRNA (miRNA), are taking part in the multiple developmental processes in plants. However, the roles of which played during the cambium activity periodicity of woody plants remain poorly understood. Here, lncRNA/circRNA-miRNA-mRNA regulatory networks of the cambium activity periodicity in Populus tomentosa was constructed, combined with morphologic observation and transcriptome profiling. Light microscopy and Periodic Acid Schiff (PAS) staining revealed that cell walls were much thicker and number of cell layers was increased during the active-dormant stage, accompanied by abundant change of polysaccharides. The novel lncRNAs and circRNAs were investigated, and we found that 2037 lncRNAs and 299 circRNAs were differentially expression during the vascular cambium period, respectively. Moreover, 1046 genes were identified as a target gene of 2037 novel lncRNAs, and 89 of which were the miRNA precursors or targets. By aligning miRNA precursors to the 7655 lncRNAs, 21 lncRNAs were identified as precursors tof 19 known miRNAs. Furthermore, the target mRNA of lncRNA/circRNA-miRNA network mainly participated in phytohormone, cell wall alteration and chlorophyll metabolism were analyzed by GO enrichment and KEGG pathway. Especially, circRNA33 and circRNA190 taking part in the phytohormone signal pathway were down-regulated during the active-dormant transition. Xyloglucan endotransglucosylase/hydrolase protein 24-like and UDP-glycosyltransferase 85A1 involved in the cell wall modification were the targets of lncRNA MSTRG.11198.1 and MSTRG.1050.1. Notably, circRNA103 and MSTRG.10851.1 regulate the cambium periodicity may interact with the miR482. These results give a new light into activity–dormancy regulation, associated with transcriptional dynamics and non-coding RNA networks of potential targets identification.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Biologia Futura: progress and future perspectives of long non-coding RNAs in forest trees

Forest trees are affected by climate change, anthropogenic pressure, as well as abiotic and biotic stresses. Conventional tree breeding has so far been limited to enhance overall productivity, and our understanding of the genetic basis of quantitative traits is still inadequate. Quantum leaps in next-generation sequencing technologies and bioinformatics have permitted the exploration and identification of various non-coding regions of the genome other than protein coding genes. These genomic regions produce various types of non-coding RNAs and regulate myriads of biological functions at epigenetic, transcriptional and translational levels. Recently, long non-coding RNAs (lncRNAs) which act as molecular switch have been identified to be pivotal molecules in forest trees. This review focuses on progress made in regulatory mechanisms in various developmental phases like wood formation, adventitious rooting and flowering and stress responses. It was predicted that complex regulatory interactions among lncRNA, miRNA and gene exist. LncRNAs can function as a sponge for miRNAs, reducing the suppressive effect of miRNAs on target mRNAs and perhaps adding a new layer of regulatory interactions among non-coding RNA classes in trees. Furthermore, network analysis revealed the interactions of lncRNA and genes during the expression of several important genes. The insights generated about lncRNAs in forest trees would enable improvement of economically important traits including the devastating abiotic and biotic stresses. In addition, solid understanding on the wide range of regulatory functions of lncRNAs on traits influencing biomass productivity and adaptation would aid the applications of biotechnology in genetic improvement of forest trees.

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.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Characterization of lncRNAs and mRNAs Involved in Powdery Mildew Resistance in Cucumber

Powdery mildew (PM) is a severe fungal disease of cucumber worldwide. Identification of genetic factors resistant to PM is of great importance for marker-assisted breeding to ensure cucumber production. Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to play important roles in plant development and immunity; however, whether they have a role in PM response in cucurbit crops remains unknown. We performed strand-specific RNA sequencing and miRNA sequencing using RNA from cucumber leaves of two near-isogenic lines (NILs), S1003 and NIL (Pm5.1) infected with PM, and systematically characterized the profiles of cucumber lncRNAs and messenger RNA (mRNAs) responsive to PM. In total, we identified 12,903 lncRNAs and 25,598 mRNAs responsive to PM. Differential expression (DE) analysis showed that 119 lncRNAs and 136 mRNAs correlated with PM resistance. Functional analysis of these DE lncRNAs and DE mRNAs revealed that they are significantly associated with phenylpropanoid biosynthesis, phenylalanine metabolism, ubiquinone and other terpenoid-quinone biosynthesis, and endocytosis. Particularly, two lncRNAs, LNC_006805 and LNC_012667, might play important roles in PM resistance. In addition, we also predicted mature miRNAs and competing endogenous RNA (ceRNA) networks of lncRNA-miRNA-mRNA involved in PM resistance. A total of 49 DE lncRNAs could potentially act as target mimics for 106 miRNAs. Taken together, our results provide an abundant resource for further exploration of cucumber lncRNAs, mRNAs, miRNAs, and ceRNAs in PM resistance, and will facilitate the molecular breeding for PM-resistant varieties to control this severe disease in cucumber.

RNA Sequencing Reveals Dynamic Carbohydrate Metabolism and Phytohormone Signaling Accompanying Post-mowing Regeneration of Forage Winter Wheat (Triticum aestivum L.)

Winter wheat (Triticum aestivum L.) is used as fresh green winter forage worldwide, and its ability to regenerate after mowing determines whether it can be used for forage production; however, the molecular mechanism of regeneration is poorly understood. This study identified long-chain coding and non-coding RNAs in the wheat cultivar "XN9106," which is cultivated for forage and grain production separately in winter and summer, and analyzed their function during post-mowing regeneration. The results showed that the degradation of carbohydrate plays an important role in regeneration, as demonstrated by decreased carbohydrate content. The increased gene expression of enzymes including β-amylase, β-fructofuranosidase, sucrose synthase, sucrose-6-phosphate synthase, trehalose-6-phosphate synthase, and trehalose-6-phosphate phosphatase in mowed seedlings suggests regeneration is fueled by degraded carbohydrates that provide energy and carbon skeletons for the Krebs cycle and amino acid synthesis. The decreased auxin content relieved the inhibition of cytokinin synthesis, that controls the transition from cell division to cell expansion and stimulates cell expansion and differentiation during the cell expansion phase, and eventually accelerate post-mowing regeneration of seedlings. Additionally, differentially expressed long-chain non-coding RNAs (lncRNAs) might participate in the regulation of gene expression related to carbohydrate metabolism and hormone signal transduction. This study demonstrated the responses of key mRNAs and lncRNAs during post-mowing regeneration of winter wheat and revealed the importance of carbohydrate and hormone during regeneration, providing valuable information for genetic improvement of forage wheat.

Investigation of long non-coding RNAs as regulatory players of grapevine response to powdery and downy mildew infection

BACKGROUND

Long non-coding RNAs (lncRNAs) are regulatory transcripts of length > 200 nt. Owing to the rapidly progressing RNA-sequencing technologies, lncRNAs are emerging as considerable nodes in the plant antifungal defense networks. Therefore, we investigated their role in Vitis vinifera (grapevine) in response to obligate biotrophic fungal phytopathogens, Erysiphe necator (powdery mildew, PM) and Plasmopara viticola (downy mildew, DM), which impose huge agro-economic burden on grape-growers worldwide.

RESULTS

Using computational approach based on RNA-seq data, 71 PM- and 83 DM-responsive V. vinifera lncRNAs were identified and comprehensively examined for their putative functional roles in plant defense response. V. vinifera protein coding sequences (CDS) were also profiled based on expression levels, and 1037 PM-responsive and 670 DM-responsive CDS were identified. Next, co-expression analysis-based functional annotation revealed their association with gene ontology (GO) terms for 'response to stress', 'response to biotic stimulus', 'immune system process', etc. Further investigation based on analysis of domains, enzyme classification, pathways enrichment, transcription factors (TFs), interactions with microRNAs (miRNAs), and real-time quantitative PCR of lncRNAs and co-expressing CDS pairs suggested their involvement in modulation of basal and specific defense responses such as: Ca²⁺-dependent signaling, cell wall reinforcement, reactive oxygen species metabolism, pathogenesis related proteins accumulation, phytohormonal signal transduction, and secondary metabolism.

CONCLUSIONS

Overall, the identified lncRNAs provide insights into the underlying intricacy of grapevine transcriptional reprogramming/post-transcriptional regulation to delay or seize the living cell-dependent pathogen growth. Therefore, in addition to defense-responsive genes such as TFs, the identified lncRNAs can be further examined and leveraged to candidates for biotechnological improvement/breeding to enhance fungal stress resistance in this susceptible fruit crop of economic and nutritional importance.

From "Dark Matter" to "Star": Insight Into the Regulation Mechanisms of Plant Functional Long Non-Coding RNAs

Long non-coding RNAs (lncRNAs) play a vital role in a variety of biological functions in plant growth and development. In this study, we provided an overview of the molecular mechanisms of lncRNAs in interacting with other biomolecules with an emphasis on those lncRNAs validated only by low-throughput experiments. LncRNAs function through playing multiple roles, including sponger for sequestering RNA or DNA, guider or decoy for recruiting or hijacking transcription factors or peptides, and scaffold for binding with chromatin modification complexes, as well as precursor of microRNAs or small interfering RNAs. These regulatory roles have been validated in several plant species with a comprehensive list of 73 lncRNA-molecule interaction pairs in 16 plant species found so far, suggesting their commonality in the plant kingdom. Such initial findings of a small number of functional plant lncRNAs represent the beginning of what is to come as lncRNAs with unknown functions were found in orders of magnitude more than proteins.

The transcriptomic response to heat stress of a jujube (Ziziphus jujuba Mill.) cultivar is featured with changed expression of long noncoding RNAs

Long non-coding RNA (lncRNA) of plant species undergoes dynamic regulation and acts in developmental and stress regulation. Presently, there is little information regarding the identification of lncRNAs in jujube (Ziziphus jujuba Mill.), and it is uncertain whether the lncRNAs could respond to heat stress (HS) or not. In our previous study, a cultivar (Hqing1-HR) of Z. jujuba were treated by HS (45°C) for 0, 1, 3, 5 and 7 days, and it was found that HS globally changed the gene expression by RNA sequencing (RNA-seq) experiments and informatics analyses. In the current study, 8260 lncRNAs were identified successfully from the previous RNA-seq data, and it indicated that lncRNAs expression was also altered globally, suggesting that the lncRNAs might play vital roles in response to HS. Furthermore, bioinformatics analyses of potential target mRNAs of lncRNAs with cis-acting mechanism were performed, and it showed that multiple differentially expressed (DE) mRNAs co-located with DElncRNAs were highly enriched in pathways associated with response to stress and regulation of metabolic process. Taken together, these findings not only provide a comprehensive identification of lncRNAs but also useful clues for molecular mechanism response to HS in jujube.

Revealing the developmental dynamics in male strobilus transcriptome of Gnetum luofuense using nanopore sequencing technology

Gnetum is a pantropical distributed gymnosperm genus. As being dioecious, Gnetum species apply female and male strobili to attract and provide nutrition to insect pollinators. Due to its unique gross morphology, a Gnetum male strobilus receives much attention in previous taxonomic and evolutionary studies. However, underlying molecular mechanisms that control male strobilus development and pollination adaptation have not been well studied. In the present study, nine full-length transcriptomes were sequenced from three developmental stages of the G. luofuense male strobili using Oxford Nanopore Technologies. In addition, weighted gene co-expression network analysis (WGCNA), and RT-qPCR analysis were performed. Our results show that a total of 3138 transcription factors and 466 long non-coding RNAs (lncRNAs) were identified, and differentially expressed lncRNAs and TFs reveal a dynamic pattern during the male strobilus development. Our results show that MADS-box and Aux/IAA TFs were differentially expressed at the three developmental stages, suggesting their important roles in the regulation of male strobilus development of G. luofuense. Results of WGCNA analysis and annotation of differentially expressed transcripts corroborate that the male strobilus development of G. luofuense is closely linked to plant hormone changes, photosynthesis, pollination drop secretion and reproductive organ defense. Our results provide a valuable resource for understanding the molecular mechanisms that drive organ evolution and pollination biology in Gnetum.