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

Transcriptome analysis suggested that lncRNAs regulate rapeseed seedlings in responding to drought stress by coordinating the phytohormone signal transduction pathways

The growth, yield, and seed quality of rapeseed are negatively affected by drought stress. Therefore, it is of great value to understand the molecular mechanism behind this phenomenon. In a previous study, long non-coding RNAs (lncRNAs) were found to play a key role in the response of rapeseed seedlings to drought stress. However, many questions remained unanswered. This study was the first to investigate the expression profile of lncRNAs not only under control and drought treatment, but also under the rehydration treatment. A total of 381 differentially expressed lncRNA and 10,253 differentially expressed mRNAs were identified in the comparison between drought stress and control condition. In the transition from drought stress to rehydration, 477 differentially expressed lncRNAs and 12,543 differentially expressed mRNAs were detected. After identifying the differentially expressed (DE) lncRNAs, the comprehensive lncRNAs-engaged network with the co-expressed mRNAs in leaves under control, drought and rehydration was investigated. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of co-expressed mRNAs identified the most significant pathways related with plant hormones (expecially abscisic acid, auxin, cytokinins, and gibberellins) in the signal transduction. The genes, co-expressed with the most-enriched DE-lncRNAs, were considered as the most effective candidates in the water-loss and water-recovery processes, including protein phosphatase 2 C (PP2C), ABRE-binding factors (ABFs), and SMALL AUXIN UP-REGULATED RNAs (SAURs). In summary, these analyses clearly demonstrated that DE-lncRNAs can act as a regulatory hub in plant-water interaction by controlling phytohormone signaling pathways and provided an alternative way to explore the complex mechanisms of drought tolerance in rapeseed.

Plant non-coding RNAs: The new frontier for the regulation of plant development and adaptation to stress

Most plant transcriptomes constitute functional non-coding RNAs (ncRNAs) that lack the ability to encode proteins. In recent years, more research has demonstrated that ncRNAs play important regulatory roles in almost all plant biological processes by modulating gene expression. Thus, it is important to study the biogenesis and function of ncRNAs, particularly in plant growth and development and stress tolerance. In this review, we systematically explore the process of formation and regulatory mechanisms of ncRNAs, particularly those of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, we provide a comprehensive overview of the recent advancements in ncRNAs research, including their regulation of plant growth and development (seed germination, root growth, leaf morphogenesis, floral development, and fruit and seed development) and responses to abiotic and biotic stress (drought, heat, cold, salinity, pathogens and insects). We also discuss research challenges and provide recommendations to advance the understanding of the roles of ncRNAs in agronomic applications.

Exploiting DNA methylation in cassava under water deficit for crop improvement

DNA methylation plays a key role in the development and plant responses to biotic and abiotic stresses. This work aimed to evaluate the DNA methylation in contrasting cassava genotypes for water deficit tolerance. The varieties BRS Formosa (bitter) and BRS Dourada (sweet) were grown under greenhouse conditions for 50 days, and afterwards, irrigation was suspended. The stressed (water deficit) and non-stressed plants (negative control) consisted the treatments with five plants per variety. The DNA samples of each variety and treatment provided 12 MethylRAD-Seq libraries (two cassava varieties, two treatments, and three replicates). The sequenced data revealed methylated sites covering 18 to 21% of the Manihot esculenta Crantz genome, depending on the variety and the treatment. The CCGG methylated sites mapped mostly in intergenic regions, exons, and introns, while the CCNGG sites mapped mostly intergenic, upstream, introns, and exons regions. In both cases, methylated sites in UTRs were less detected. The differentially methylated sites analysis indicated distinct methylation profiles since only 12% of the sites (CCGG and CCNGG) were methylated in both varieties. Enriched gene ontology terms highlighted the immediate response of the bitter variety to stress, while the sweet variety appears to suffer more potential stress-damages. The predicted protein-protein interaction networks reinforced such profiles. Additionally, the genomes of the BRS varieties uncovered SNPs/INDELs events covering genes stood out by the interactomes. Our data can be useful in deciphering the roles of DNA methylation in cassava drought-tolerance responses and adaptation to abiotic stresses.

Exploring the Regulatory Dynamics of BrFLC-Associated lncRNA in Modulating the Flowering Response of Chinese Cabbage

Vernalization plays a crucial role in the flowering and yield of Chinese cabbage, a process intricately influenced by long non-coding RNAs (lncRNAs). Our research focused on lncFLC1, lncFLC2a, and lncFLC2b, which emerged as key players in this process. These lncRNAs exhibited an inverse expression pattern to the flowering repressor genes FLOWERING LOCUS C 1 (BrFLC1) and FLOWERING LOCUS C 2 (BrFLC2) during vernalization, suggesting a complex regulatory mechanism. Notably, their expression in the shoot apex and leaves was confirmed through in fluorescent in situ hybridization (FISH). Furthermore, when these lncRNAs were overexpressed in Arabidopsis, a noticeable acceleration in flowering was observed, unveiling functional similarities to Arabidopsis's COLD ASSISTED INTRONIC NONCODING RNA (COOLAIR). This resemblance suggests a potentially conserved regulatory mechanism across species. This study not only enhances our understanding of lncRNAs in flowering regulation, but also opens up new possibilities for their application in agricultural practices.

Different transcription of novel, functional long non-coding RNA genes by UV-B in green algae, Volvox carteri

Long non-coding RNAs (lncRNAs) are identified as important regulatory molecules related to diverse biological processes. In recent years, benefiting from the rapid development of high-throughput sequencing technology, RNA-seq, and analysis methods, more lncRNAs have been identified and discovered in various plant and algal species. However, so far, only limited studies related to algal lncRNAs are available. Volvox carteri f. nagariensis is the best multicellular model organism to study in developmental and evolutionary biology; therefore, studying and increasing information about this species is important. This study identified lncRNAs in the multicellular green algae Volvox carteri and 1457 lncRNAs were reported, using RNA-seq data and with the help of bioinformatics tools and software. This study investigated the effect of low-dose UV-B radiation on changes in the expression profile of lncRNAs in gonidial and somatic cells. The differential expression of lncRNAs was analyzed between the treatment (UV-B) and the control (WL) groups in gonidial and somatic cells. A total of 37 and 26 lncRNAs with significant differential expression in gonidial and somatic cells, respectively, were reported. Co-expression analysis between the lncRNAs and their neighbor protein-coding genes (in the interval of ± 10 Kb) was accomplished. In gonidial cells, 184 genes with a positive correlation and 13 genes with a negative correlation (greater than 0.95), and in somatic cells, 174 genes with a positive correlation, and 18 genes with a negative correlation were detected. Functional analysis of neighboring coding genes was also performed based on gene ontology. The results of the current work may help gain deeper insight into the regulation of gene expression in the studied model organism, Volvox carteri.

Long non-coding RNA-mediated epigenetic response for abiotic stress tolerance in plants

Plants perceive environmental fluctuations as stress and confront several stresses throughout their life cycle individually or in combination. Plants have evolved their sensing and signaling mechanisms to perceive and respond to a variety of stresses. Epigenetic regulation plays a critical role in the regulation of genes, spatiotemporal expression of genes under stress conditions and imparts a stress memory to encounter future stress responses. It is quintessential to integrate our understanding of genetics and epigenetics to maintain plant fitness, achieve desired genetic gains with no trade-offs, and durable long-term stress tolerance. The long non-coding RNA >200 nts having no coding potential (or very low) play several roles in epigenetic memory, contributing to the regulation of gene expression and the maintenance of cellular identity which include chromatin remodeling, imprinting (dosage compensation), stable silencing, facilitating nuclear organization, regulation of enhancer-promoter interactions, response to environmental signals and epigenetic switching. The lncRNAs are involved in a myriad of stress responses by activation or repression of target genes and hence are potential candidates for deploying in climate-resilient breeding programs. This review puts forward the significant roles of long non-coding RNA as an epigenetic response during abiotic stresses in plants and the prospects of deploying lncRNAs for designing climate-resilient plants.

Metabolic GWAS-based dissection of genetic basis underlying nutrient quality variation and domestication of cassava storage root

BACKGROUND

Metabolites play critical roles in regulating nutritional qualities of plants, thereby influencing their consumption and human health. However, the genetic basis underlying the metabolite-based nutrient quality and domestication of root and tuber crops remain largely unknown.

RESULTS

We report a comprehensive study combining metabolic and phenotypic genome-wide association studies to dissect the genetic basis of metabolites in the storage root (SR) of cassava. We quantify 2,980 metabolic features in 299 cultivated cassava accessions. We detect 18,218 significant marker-metabolite associations via metabolic genome-wide association mapping and identify 12 candidate genes responsible for the levels of metabolites that are of potential nutritional importance. Me3GT, MeMYB4, and UGT85K4/UGT85K5, which are involved in flavone, anthocyanin, and cyanogenic glucoside metabolism, respectively, are functionally validated through in vitro enzyme assays and in vivo gene silencing analyses. We identify a cluster of cyanogenic glucoside biosynthesis genes, among which CYP79D1, CYP71E7b, and UGT85K5 are highly co-expressed and their allelic combination contributes to low linamarin content. We find MeMYB4 is responsible for variations in cyanidin 3-O-glucoside and delphinidin 3-O-rutinoside contents, thus controlling SR endothelium color. We find human selection affects quercetin 3-O-glucoside content and SR weight per plant. The candidate gene MeFLS1 is subject to selection during cassava domestication, leading to decreased quercetin 3-O-glucoside content and thus increased SR weight per plant.

CONCLUSIONS

These findings reveal the genetic basis of cassava SR metabolome variation, establish a linkage between metabolites and agronomic traits, and offer useful resources for genetically improving the nutrition of cassava and other root crops.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

Long non-coding RNA-mediated competing endogenous RNA regulatory network during flower development and color formation in Melastoma candidum

M. candidum, an evergreen shrubby flower known for its superior adaptation ability in South China, has gained increased attention in garden applications. However, scant attention has been paid to its flower development and color formation process at the non-coding RNA level. To fill this gap, we conducted a comprehensive analysis based on long non-coding RNA sequencing (lncRNA-seq), RNA-seq, small RNA sequencing (sRNA-seq), and widely targeted metabolome detection of three different flower developmental stages of M. candidum. After differentially expressed lncRNAs (DElncRNAs), differentially expressed mRNAs (DEmRNAs), differentially expressed microRNAs (DEmiRNAs), and differentially synthesized metabolites (DSmets) analyses between the different flower developmental stages, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to identify some key genes and metabolites in flavonoid, flavone, anthocyanin, carotenoid, and alkaloid-related GO terms and biosynthetic pathways. Three direct-acting models, including antisense-acting, cis-acting, and trans-acting between lncRNAs and mRNAs, were detected to illustrate the direct function of lncRNAs on target genes during flower development and color formation. Based on the competitive endogenous RNA (ceRNA) regulatory theory, we constructed a lncRNA-mediated regulatory network composed of DElncRNAs, DEmiRNAs, DEmRNAs, and DSmets to elucidate the indirect role of lncRNAs in the flower development and color formation of M. candidum. By utilizing correlation analyses between DERNAs and DSmets within the ceRNA regulatory network, alongside verification trials of the ceRNA regulatory mechanism, the study successfully illustrated the significance of lncRNAs in flower development and color formation process. This research provides a foundation for improving and regulating flower color at the lncRNA level in M. candidum, and sheds light on the potential applications of non-coding RNA in studies of flower development.

Analysis of Alternative Splicing and Long Noncoding RNAs After the Edwardsiella anguillarum Infected the Immunized European Eels (Anguilla anguilla) Revealed the Role of Outer Membrane Protein A in OmpA Subunit Vaccine

Edwardsiella anguillarum is a bacterium that commonly infects cultivated eels. Outer membrane protein A (OmpA) emulsified with Freund's adjuvant has been shown to be an effective fishery vaccine against this pathogen. However, the specific roles of OmpA in the vaccine have not been fully explored. In this study, we performed RNA-seq in the liver of a European eel (Anguilla anguilla) after challenge with E. anguillarum in eels previously immunized with an OmpA subunit vaccine. Our aim was to elucidate the differentially alternative splicing (DAS) and differentially expressed long noncoding RNAs (DE-lncRNAs) using a genome-wide transcriptome. The results showed after that at 28 days post-immunization, eels challenged with E. anguillarum (Con_inf) exhibited severe pathological changes in the liver. In contrast, the OmpA infused eels (OmpA_inf group) showed infiltrated lymphocytes, while Freund's adjuvant-inoculated eels (FCIA_inf group) showed edema of hepatocytes and blood coagulation. The relative percent survival (RPS) was 77.7% and 44.4% for OmpA_inf and FCIA_inf compared to the Con_inf group. We identified 37 DE-lncRNAs and 293 DAS genes between OmpA_inf and FCIA_inf. Interactions between DAS gene-expressed proteins indicated that 66 expressed proteins formed 20 networks. Additionally, 33 DE-lncRNAs interacted with 194 target genes formed 246 and 41 networks in co-expression and co-location. Taken together, our findings demonstrate that the OmpA subunit vaccine elicits a higher RPS and provides novel insights into the role of OmpA through DAS genes and DE-lncRNAs perspective. These results are significant for the development of fishery subunit vaccines.

ASLncR: a novel computational tool for prediction of abiotic stress-responsive long non-coding RNAs in plants

Abiotic stresses are detrimental to plant growth and development and have a major negative impact on crop yields. A growing body of evidence indicates that a large number of long non-coding RNAs (lncRNAs) are key to many abiotic stress responses. Thus, identifying abiotic stress-responsive lncRNAs is essential in crop breeding programs in order to develop crop cultivars resistant to abiotic stresses. In this study, we have developed the first machine learning-based computational model for predicting abiotic stress-responsive lncRNAs. The lncRNA sequences which were responsive and non-responsive to abiotic stresses served as the two classes of the dataset for binary classification using the machine learning algorithms. The training dataset was created using 263 stress-responsive and 263 non-stress-responsive sequences, whereas the independent test set consists of 101 sequences from both classes. As the machine learning model can adopt only the numeric data, the Kmer features ranging from sizes 1 to 6 were utilized to represent lncRNAs in numeric form. To select important features, four different feature selection strategies were utilized. Among the seven learning algorithms, the support vector machine (SVM) achieved the highest cross-validation accuracy with the selected feature sets. The observed 5-fold cross-validation accuracy, AU-ROC, and AU-PRC were found to be 68.84, 72.78, and 75.86%, respectively. Furthermore, the robustness of the developed model (SVM with the selected feature) was evaluated using an independent test dataset, where the overall accuracy, AU-ROC, and AU-PRC were found to be 76.23, 87.71, and 88.49%, respectively. The developed computational approach was also implemented in an online prediction tool ASLncR accessible at https://iasri-sg.icar.gov.in/aslncr/ . The proposed computational model and the developed prediction tool are believed to supplement the existing effort for the identification of abiotic stress-responsive lncRNAs in plants.

Long Non-Coding RNAs of Plants in Response to Abiotic Stresses and Their Regulating Roles in Promoting Environmental Adaption

Abiotic stresses triggered by climate change and human activity cause substantial agricultural and environmental problems which hamper plant growth. Plants have evolved sophisticated mechanisms in response to abiotic stresses, such as stress perception, epigenetic modification, and regulation of transcription and translation. Over the past decade, a large body of literature has revealed the various regulatory roles of long non-coding RNAs (lncRNAs) in the plant response to abiotic stresses and their irreplaceable functions in environmental adaptation. LncRNAs are recognized as a class of ncRNAs that are longer than 200 nucleotides, influencing a variety of biological processes. In this review, we mainly focused on the recent progress of plant lncRNAs, outlining their features, evolution, and functions of plant lncRNAs in response to drought, low or high temperature, salt, and heavy metal stress. The approaches to characterize the function of lncRNAs and the mechanisms of how they regulate plant responses to abiotic stresses were further reviewed. Moreover, we discuss the accumulating discoveries regarding the biological functions of lncRNAs on plant stress memory as well. The present review provides updated information and directions for us to characterize the potential functions of lncRNAs in abiotic stresses in the future.

Using Exogenous Melatonin, Glutathione, Proline, and Glycine Betaine Treatments to Combat Abiotic Stresses in Crops

Abiotic stresses, such as drought, salinity, heat, cold, and heavy metals, are associated with global climate change and hamper plant growth and development, affecting crop yields and quality. However, the negative effects of abiotic stresses can be mitigated through exogenous treatments using small biomolecules. For example, the foliar application of melatonin provides the following: it protects the photosynthetic apparatus; it increases the antioxidant defenses, osmoprotectant, and soluble sugar levels; it prevents tissue damage and reduces electrolyte leakage; it improves reactive oxygen species (ROS) scavenging; and it increases biomass, maintains the redox and ion homeostasis, and improves gaseous exchange. Glutathione spray upregulates the glyoxalase system, reduces methylglyoxal (MG) toxicity and oxidative stress, decreases hydrogen peroxide and malondialdehyde accumulation, improves the defense mechanisms, tissue repairs, and nitrogen fixation, and upregulates the phytochelatins. The exogenous application of proline enhances growth and other physiological characteristics, upregulates osmoprotection, protects the integrity of the plasma lemma, reduces lipid peroxidation, increases photosynthetic pigments, phenolic acids, flavonoids, and amino acids, and enhances stress tolerance, carbon fixation, and leaf nitrogen content. The foliar application of glycine betaine improves growth, upregulates osmoprotection and osmoregulation, increases relative water content, net photosynthetic rate, and catalase activity, decreases photorespiration, ion leakage, and lipid peroxidation, protects the oxygen-evolving complex, and prevents chlorosis. Chemical priming has various important advantages over transgenic technology as it is typically more affordable for farmers and safe for plants, people, and animals, while being considered environmentally acceptable. Chemical priming helps to improve the quality and quantity of the yield. This review summarizes and discusses how exogenous melatonin, glutathione, proline, and glycine betaine can help crops combat abiotic stresses.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Physiological and Proteomic Responses of Cassava to Short-Term Extreme Cool and Hot Temperature

Temperature is one of the most critical factors affecting cassava metabolism and growth. This research was conducted to investigate the effects of short-term exposure to extreme cool (15 °C) and hot (45 °C) temperature on photosynthesis, biochemical and proteomics changes in potted plants of two cassava cultivars, namely Rayong 9 and Kasetsart 50. One-month-old plants were exposed to 15, 30, and 45 °C for 60 min in a temperature chamber under light intensity of 700 μmol m⁻² s⁻¹. Compared to the optimum temperature (30 °C), exposure to 15 °C resulted in 28% reduction in stomatal conductance (gs) and 62% reduction in net photosynthesis rate (P_n). In contrast, gs under 45 °C increased 2.61 folds, while P_n was reduced by 50%. The lower P_n but higher electron transport rate (ETR) of the cold-stressed plants indicated that a greater proportion of electrons was transported via alternative pathways to protect chloroplast from being damaged by reactive oxygen species (ROS). Moreover, malondialdehyde (MDA) contents, a marker related to the amount of ROS, were significantly higher at low temperature. Proteomics analysis revealed some interesting differentially expressed proteins (DEPs) including annexin, a multi-functional protein functioning in early events of heat stress signaling. In response to low-temperature stress, AP2/ERF domain-containing protein (a cold-related transcription factor) and glutaredoxin domain-containing protein (a component of redox signaling network under cold stress) were detected. Taken together, both cultivars were more sensitive to low than high temperature. Moreover, Rayong 9 displayed higher P_n under both temperature stresses, and was more efficient in controlling ROS under cold stress than Kasetsart 50.

A Long Noncoding RNA Derived from lncRNA-mRNA Networks Modulates Seed Vigor

The discovery of long noncoding RNAs (lncRNAs) has filled a great gap in our understanding of posttranscriptional gene regulation in a variety of biological processes related to plant stress responses. However, systematic analyses of the lncRNAs expressed in rice seeds that germinate under cold stress have been elusive. In this study, we performed strand-specific whole transcriptome sequencing in germinated rice seeds under cold stress and normal temperature. A total of 6258 putative lncRNAs were identified and expressed in a stage-specific manner compared to mRNA. By investigating the targets of differentially expressed (DE) lncRNAs of LT-I (phase I of low temperature)/NT-I (phase I of normal temperature), it was shown that the auxin-activated signaling pathway was significantly enriched, and twenty-three protein-coding genes with most of the members of the SAUR family located in chromosome 9 were identified as the candidate target genes that may interact with five lncRNAs. A seed vigor-related lncRNA, SVR, which interplays with the members of the SAUR gene family in cis was eventually identified. The CRISPR/Cas 9 engineered mutations in SVR cause delay of germination. The findings provided new insights into the connection between lncRNAs and the auxin-activated signaling pathway in the regulation of rice seed vigor.

Epigenetic marks for mitigating abiotic stresses in plants

Abiotic stressors are one of the major factors affecting agricultural output. Plants have evolved adaptive systems to respond appropriately to various environmental cues. These responses can be accomplished by modulating or fine-tuning genetic and epigenetic regulatory mechanisms. Understanding the response of plants' molecular features to abiotic stress is a priority in the current period of continued environmental changes. Epigenetic modifications are necessary that control gene expression by changing chromatin status and recruiting various transcription regulators. The present study summarized the current knowledge on epigenetic modifications concerning plant responses to various environmental stressors. The functional relevance of epigenetic marks in regulating stress tolerance has been revealed, and epigenetic changes impact the effector genes. This study looks at the epigenetic mechanisms that govern plant abiotic stress responses, especially DNA methylation, histone methylation/acetylation, chromatin remodeling, and various metabolites. Plant breeders will benefit from a thorough understanding of these processes to create alternative crop improvement approaches. Genome editing with clustered regularly interspaced short palindromic repeat/CRISPR-associated proteins (CRISPR/Cas) provides genetic tools to make agricultural genetic engineering more sustainable and publicly acceptable.

Transcriptome and DNA methylome provide insights into the molecular regulation of drought stress in sea buckthorn

Sea buckthorn is a typical drought-resistant tree species. However, there is a general lack of understanding of the pattern of DNA methylation linked with sea buckthorn responses to drought, and its relationship with drought tolerance mechanisms. In this study, we performed whole-transcriptome RNA sequencing and methylome sequencing in response to drought stress to explore differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs in sea buckthorn leaves. Based on predicted DE pairs, we constructed a competitive endogenous RNA network, which revealed potential transcriptional regulatory roles in response to drought stress. The results of methylome sequencing revealed that the DNA methylation level was increased in sea buckthorn leaves under drought stress. We identified 13,405 differentially methylated regions between CK and TR. We found one DMR-associated DEG (Vacuolar-sorting receptor 6) involved in the ABA accumulation pathway. In addition, two DNA methyltransferases (HrMET1 and HrDRM1) were closely associated with drought-induced hypermethylation in sea buckthorn. Together, we firstly conducted a comprehensive transcriptomic and epigenetic analysis of sea buckthorn under drought stress, providing a resource for further study of the potential functions of genes, miRNAs, lncRNAs, circRNAs and DNA methyltransferases.

Global identification of full-length cassava lncRNAs unveils the role of cold-responsive intergenic lncRNA 1 in cold stress response

Long noncoding RNAs (lncRNAs) have been considered to be important regulators of gene expression in a range of biological processes in plants. A large number of lncRNAs have been identified in plants. However, most of their biological functions still remain to be determined. Here, we identified a total of 3004 lncRNAs in cassava under normal or cold-treated conditions from Iso-seq data. We further characterized a cold-responsive intergenic lncRNA 1 (CRIR1) as a novel positive regulator of the plant response to cold stress. CRIR1 can be significantly induced by cold treatment. Ectopic expression of CRIR1 in cassava enhanced the cold tolerance of transgenic plants. Transcriptome analysis demonstrated that CRIR1 regulated a range of cold stress-related genes in a CBF-independent pathway. We further found that CRIR1 RNA can interact with cassava cold shock protein 5 (MeCSP5), which acts as an RNA chaperone, indicating that CRIR1 may recruit MeCSP5 to improve the translation efficiency of messenger RNA. In summary, our study extends the repertoire of lncRNAs in plants as well as their role in cold stress responses. Moreover, it reveals a mechanism by which CRIR1 affected cold stress response by modulating the expression of stress-responsive genes and increasing their translational yield.

Differential Expression of lncRNAs and miRNAs Between Self-Rooting Juvenile and Donor Clones Unveils Novel Insight Into the Molecular Regulation of Rubber Biosynthesis in Hevea brasiliensis

The rubber tree (Hevea brasiliensis Muell. Arg.) is a tropical tree species that produce natural rubber. Self-rooted juvenile clones (SRJCs) are novel rubber tree planting materials developed through primary somatic embryogenesis. SRJCs have a higher rubber yield compared with donor clones (DCs). The molecular basis underlying increased rubber yield in SRJCs remains largely unknown. Here, the latex from SRJCs and DCs were collected for strand-specific and small RNA-seq methods. A total of 196 differentially expressed long noncoding RNAs (DELs), and 11 differentially expressed microRNAs were identified in latex between SRJCs and DCs. Targeted genes of DELs were markedly enriched for various biological pathways related to plant hormone signal transduction, photosynthesis, glutathione metabolism, and amino acids biosynthesis. DELs probably acted as cis-acting regulation was calculated, and these DELs relevant to potentially regulate rubber biosynthesis, reactive oxygen species metabolism, and epigenetic modification. Furthermore, the DELs acting as microRNA targets were studied. The interaction of microRNA and DELs might involve in the regulation of natural rubber biosynthesis.

Bidirectional lncRNA Transfer between Cuscuta Parasites and Their Host Plant

Dodder species (Cuscuta spp.) are holoparasites that have extensive material exchange with their host plants through vascular connections. Recent studies on cross-species transfer have provided breakthrough insights, but little is known about the interaction mechanisms of the inter-plant mobile substances in parasitic systems. We sequenced the transcriptomes of dodder growing on soybean hosts to characterize the long non-coding RNA (lncRNA) transfer between the two species, and found that lncRNAs can move in high numbers (365 dodder lncRNAs and 14 soybean lncRNAs) in a bidirectional manner. Reverse transcription-polymerase chain reaction further confirmed that individual lncRNAs were trafficked in the dodder-soybean parasitic system. To reveal the potential functions of mobile transcripts, the Gene Ontology terms of mobile lncRNA target genes were predicted, and mobile dodder target genes were found to be mainly enriched in "metabolic process", "catalytic activity", "signaling", and "response to stimulus" categories, whereas mobile soybean target genes were enriched in organelle-related categories, indicating that specific mobile lncRNAs may be important in regulating dodder parasitism. Our findings reveal that lncRNAs are transferred between dodder and its host soybean plants, which may act as critical regulators to coordinate the host-dodder interaction at the whole parasitic level.

Overexpression of cassava RSZ21b enhances drought tolerance in Arabidopsis

Drought is one of the major environmental constraints affecting crop productivity. Plants have to adjust their developmental and physiological processes to cope with drought. We previously identified 18 cassava serine/arginine-rich (SR) proteins that had a pivotal role in alternative splicing in response to environmental stress. However, functional characterization of SR proteins is rarely explored. Here, we characterized the RSZ subfamily gene MeRSZ21b in cassava. The RSZ21b belongs to the RSZ subfamily, which was widely distributed in major crops and was highly conserved. Quantitative RT-PCR assay showed that the expression of MeRSZ21b was significantly induced by drought. Moreover, overexpression of MeRSZ21b in Arabidopsis was hypersensitive to abscisic acid (ABA) in the phases of seed germination and post-germination seedling growth. Meantime, MeRSZ21b overexpression lines were resistant to sorbitol treatment, and quickly closed the stomata when compared with Col-0 under drought condition. Importantly, overexpression of MeRSZ21b resulted in improved drought tolerance through modulating ABA-dependent signaling. Therefore, our findings refine our knowledge of the SR protein-coding genes and provide novel insights for enhancing plant resistance to environmental stress.

Expression Patterns and Regulation of Non-Coding RNAs during Synthesis of Cellulose in Eucalyptus grandis Hill

Cellulose, an essential structural component in the plant cell wall and a renewable biomass resource, plays a significant role in nature. Eucalyptus’s excellent timber tree species (including Eucalyptus grandis Hill) provide many raw materials for the paper and wood industries. The synthesis of cellulose is a very complex process involving multiple genes and regulated by various biological networks. However, research on regulating associated genes and non-coding RNAs during cellulose synthesis in E. grandis remains lacking. In this study, the wood anatomical characteristics and chemical indexes of E. grandis were analyzed by taking three different parts (diameter at breast height (DBH), middle and upper part of the trunk) from the main stem of E. grandis as raw materials. The role of non-coding RNAs (Long non-coding RNA, lncRNA; Micro RNA, miRNA; Circle RNA, circRNA) on regulating candidate genes was presented, and the network map of ceRNA (Competing endogenous RNA) regulation during wood cellulose biosynthesis of E. grandis was constructed. The transcriptome sequencing of nine samples obtained from the trunk of the immature xylem in E. grandis at DBH, middle and upper parts had a 95.81 G clean reading, 57,480 transcripts, 7365 lncRNAs, and 5180 circRNAs. Each sample had 172–306 known miRNAs and 1644–3508 new miRNAs. A total of 190 DE-lncRNAs (Differentially expressed long non-coding RNAs), 174 DE-miRNAs (Differentially expressed micro RNAs), and 270 DE-circRNAs (Differentially expressed circle RNAs) were obtained by comparing transcript expression levels. Four lncRNAs and nine miRNAs were screened out, and the ceRNA regulatory network was constructed. LncRNA1 and lncRNA4 regulated the genes responsible for cellulose synthesis in E. grandis, which were overexpressed in 84K (Populus Alba × Populus glandulosa) poplar. The cellulose and lignin content in lncRNA4-oe were significantly higher than wild type 84K poplar and lncRNA1-oe. The average plant height, middle and basal part of the stem diameter in lncRNA4-oe were significantly higher than the wild type. However, there was no significant difference between the growth of lncRNA1-oe and the wild type. Further studies are warranted to explore the molecular regulatory mechanism of cellulose biosynthesis in Eucalyptus species.

Dysregulated lncRNAs are Involved in the Progress of Sepsis by Constructing Regulatory Networks in Whole Blood Cells

Sepsis is a highly heterogeneous syndrome that is caused by an unbalanced host response to an infection. Long noncoding RNAs (lncRNAs) have been reported to exert regulatory roles in a variety of biological processes, and became potential biomarkers and therapeutic targets for diverse diseases. However, current understanding on the roles of lncRNAs in sepsis is extremely limited. Herein, to decipher the underlying functions of lncRNAs, we reexplored the 83 transcriptome datasets from specimens with sepsis, no_sepsis by final diagnosis, and control. The results of differentially expressed genes (DEGs), differentially expressed lncRNA (DElncRNA) analysis, and co-expression analysis of lncRNA–mRNA pairs were obtained. We found that the expression pattern of lncRNAs was significantly activated in sepsis specimens, which was clearly distinguished in sepsis from no_sepsis and control specimens. By performing co-expression analysis, we found DElncRNAs were closely related to T-cell activation and immune response–related terms in sepsis by regulating mRNA expression in the trans manner. The lncRNA–mRNA network and the qRT-PCR test revealed that lncRNAs LINC00861, RP11-284N8.3, and CTB-61M7.2 were significantly correlated with the pathogenesis of sepsis. In addition, weighted gene co-expression analysis (WGCNA) and cis-regulation analysis also revealed sepsis-specific lncRNAs were highly associated with important biological processes correlated with sepsis. In summary, the systematic dysregulation of lncRNAs is tightly involved in the remodeling of gene expression regulatory network in sepsis, and the lncRNA–mRNA expression network may be used to refine biomarker predictions for developing novel therapeutic approaches in sepsis.

Integrative lncRNA landscape reveals lncRNA-coding gene networks in the secondary cell wall biosynthesis pathway of moso bamboo (Phyllostachys edulis)

Background

LncRNAs are extensively involved in plant biological processes. However, the lack of a comprehensive lncRNA landscape in moso bamboo has hindered the molecular study of lncRNAs. Moreover, the role of lncRNAs in secondary cell wall (SCW) biosynthesis of moso bamboo is elusive.

Results

For comprehensively identifying lncRNA throughout moso bamboo genome, we collected 231 RNA-Seq datasets, 1 Iso-Seq dataset, and 1 full-length cDNA dataset. We used a machine learning approach to improve the pipeline of lncRNA identification and functional annotation based on previous studies and identified 37,009 lncRNAs in moso bamboo. Then, we established a network of potential lncRNA-coding gene for SCW biosynthesis and identified SCW-related lncRNAs. We also proposed that a mechanism exists in bamboo to direct phenylpropanoid intermediates to lignin or flavonoids biosynthesis through the PAL/4CL/C4H genes. In addition, we identified 4 flavonoids and 1 lignin-preferred genes in the PAL/4CL/C4H gene families, which gained implications in molecular breeding.

Conclusions

We provided a comprehensive landscape of lncRNAs in moso bamboo. Through analyses, we identified SCW-related lncRNAs and improved our understanding of lignin and flavonoids biosynthesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07953-z.

Transcriptome Profile Analysis of Arabidopsis Reveals the Drought Stress-Induced Long Non-coding RNAs Associated With Photosynthesis, Chlorophyll Synthesis, Fatty Acid Synthesis and Degradation

Long non-coding RNAs (lncRNAs) play an important role in the response of plants to drought stress. The previous studies have reported that overexpression of LEA3 and VOC could enhance drought tolerance and improve the oil content in Brassica napus and Arabidopsis thaliana, and most of the efforts have been invested in the gene function analysis, there is little understanding of how genes that involved in these important pathways are regulated. In the present study, the transcriptomic results of LEA3 and VOC over-expressed (OE) lines were compared with the RNAi lines, mutant lines and control lines under long-term and short-term drought treatment, a series of differentially expressed lncRNAs were identified, and their regulation patterns in mRNA were also investigated in above mentioned materials. The regulation of the target genes of differentially expressed lncRNAs on plant biological functions was studied. It was revealed that the mutant lines had less drought-response related lncRNAs than that of the OE lines. Functional analysis demonstrated that multiple genes were involved in the carbon-fixing and chlorophyll metabolism, such as CDR1, CHLM, and CH1, were regulated by the upregulated lncRNA in OE lines. In LEA-OE, AT4G13180 that promotes the fatty acid synthesis was regulated by five lncRNAs that were upregulated under both long-term and short-term drought treatments. The key genes, including of SHM1, GOX2, and GS2, in the methylglyoxal synthesis pathway were all regulated by a number of down-regulated lncRNAs in OE lines, thereby reducing the content of such harmful compounds produced under stress in plants. This study identified a series of lncRNAs related to the pathways that affect photosynthesis, chlorophyll synthesis, fatty acid synthesis, degradation, and other important effects on drought resistance and oil content. The present study provided a series of lncRNAs for further improvement of crop varieties, especially drought resistant and oil content traits.

Regulatory non-coding RNAs: a new frontier in regulation of plant biology

Beyond the most crucial roles of RNA molecules as a messenger, ribosomal, and transfer RNAs, the regulatory role of many non-coding RNAs (ncRNAs) in plant biology has been recognized. ncRNAs act as riboregulators by recognizing specific nucleic acid targets through homologous sequence interactions to regulate plant growth, development, and stress responses. Regulatory ncRNAs, ranging from small to long ncRNAs (lncRNAs), exert their control over a vast array of biological processes. Based on the mode of biogenesis and their function, ncRNAs evolved into different forms that include microRNAs (miRNAs), small interfering RNAs (siRNAs), miRNA variants (isomiRs), lncRNAs, circular RNAs (circRNAs), and derived ncRNAs. This article explains the different classes of ncRNAs and their role in plant development and stress responses. Furthermore, the applications of regulatory ncRNAs in crop improvement, targeting agriculturally important traits, have been discussed.

Identification of the specific long-noncoding RNAs involved in night-break mediated flowering retardation in Chenopodium quinoa

Background

Night-break (NB) has been proven to repress flowering of short-day plants (SDPs). Long-noncoding RNAs (lncRNAs) play key roles in plant flowering. However, investigation of the relationship between lncRNAs and NB responses is still limited, especially in Chenopodium quinoa, an important short-day coarse cereal.

Results

In this study, we performed strand-specific RNA-seq of leaf samples collected from quinoa seedlings treated by SD and NB. A total of 4914 high-confidence lncRNAs were identified, out of which 91 lncRNAs showed specific responses to SD and NB. Based on the expression profiles, we identified 17 positive- and 7 negative-flowering lncRNAs. Co-expression network analysis indicated that 1653 mRNAs were the common targets of both types of flowering lncRNAs. By mapping these targets to the known flowering pathways in model plants, we found some pivotal flowering homologs, including 2 florigen encoding genes (FT (FLOWERING LOCUS T) and TSF (TWIN SISTER of FT) homologs), 3 circadian clock related genes (EARLY FLOWERING 3 (ELF3), LATE ELONGATED HYPOCOTYL (LHY) and ELONGATED HYPOCOTYL 5 (HY5) homologs), 2 photoreceptor genes (PHYTOCHROME A (PHYA) and CRYPTOCHROME1 (CRY1) homologs), 1 B-BOX type CONSTANS (CO) homolog and 1 RELATED TO ABI3/VP1 (RAV1) homolog, were specifically affected by NB and competed by the positive and negative-flowering lncRNAs. We speculated that these potential flowering lncRNAs may mediate quinoa NB responses by modifying the expression of the floral homologous genes.

Conclusions

Together, the findings in this study will deepen our understanding of the roles of lncRNAs in NB responses, and provide valuable information for functional characterization in future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07605-2.

Genome-wide identification and characterization of long non-coding RNAs involved in flag leaf senescence of rice

KEY MESSAGE

This study showed the systematic identification of long non-coding RNAs (lncRNAs) involving in flag leaf senescence of rice, providing the possible lncRNA-mRNA regulatory relationships and lncRNA-miRNA-mRNA ceRNA networks during leaf senescence. LncRNAs have been reported to play crucial roles in diverse biological processes. However, no systematic identification of lncRNAs associated with leaf senescence in plants has been studied. In this study, a genome-wide high throughput sequencing analysis was performed using rice flag leaves developing from normal to senescence. A total of 3953 lncRNAs and 38757 mRNAs were identified, of which 343 lncRNAs and 9412 mRNAs were differentially expressed. Through weighted gene co-expression network analysis (WGCNA), 22 continuously down-expressed lncRNAs targeting 812 co-expressed mRNAs and 48 continuously up-expressed lncRNAs targeting 1209 co-expressed mRNAs were considered to be significantly associated with flag leaf senescence. Gene Ontology results suggested that the senescence-associated lncRNAs targeted mRNAs involving in many biological processes, including transcription, hormone response, oxidation-reduction process and substance metabolism. Additionally, 43 senescence-associated lncRNAs were predicted to target 111 co-expressed transcription factors. Interestingly, 8 down-expressed lncRNAs and 29 up-expressed lncRNAs were found to separately target 12 and 20 well-studied senescence-associated genes (SAGs). Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 6 down-expressed lncRNAs possibly regulated 51 co-expressed mRNAs through 15 miRNAs, and 14 up-expressed lncRNAs possibly regulated 117 co-expressed mRNAs through 21 miRNAs. Importantly, by expression validation, a conserved miR164-NAC regulatory pathway was found to be possibly involved in leaf senescence, where lncRNA MSTRG.62092.1 may serve as a ceRNA binding with miR164a and miR164e to regulate three transcription factors. And two key lncRNAs MSTRG.31014.21 and MSTRG.31014.36 also could regulate the abscisic-acid biosynthetic gene BGIOSGA025169 (OsNCED4) and BGIOSGA016313 (NAC family) through osa-miR5809. The possible regulation networks of lncRNAs involving in leaf senescence were discussed, and several candidate lncRNAs were recommended for prior transgenic analysis. These findings will extend the understanding on the regulatory roles of lncRNAs in leaf senescence, and lay a foundation for functional research on candidate lncRNAs.

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.

Roles of Non-Coding RNAs in Response to Nitrogen Availability in Plants

Nitrogen (N) is an essential nutrient for plant growth and development; therefore, N deficiency is a major limiting factor in crop production. Plants have evolved mechanisms to cope with N deficiency, and the role of protein-coding genes in these mechanisms has been well studied. In the last decades, regulatory non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs), and long ncRNAs (lncRNAs), have emerged as important regulators of gene expression in diverse biological processes. Recent advances in technologies for transcriptome analysis have enabled identification of N-responsive ncRNAs on a genome-wide scale. Characterization of these ncRNAs is expected to improve our understanding of the gene regulatory mechanisms of N response. In this review, we highlight recent progress in identification and characterization of N-responsive ncRNAs in Arabidopsis thaliana and several other plant species including maize, rice, and Populus.

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.

Transcriptional analysis reveals potential genes and regulatory networks involved in salicylic acid-induced flowering in duckweed (Lemna gibba)

Duckweed is a simple aquatic floating plant having great potential in sewage treatment and bioenergy production. Duckweed rarely flowers in nature, which greatly limits its germplasm collection, conservation, and heterosis usage. Salicylic acid (SA) can efficiently induce flowering of duckweed (e.g., Lemna gibba); however, the related genes and regulatory networks remain unclear. In this work, we demonstrated that L. gibba flowering induced by SA was photoperiod-dependent, stress-involved, and abscisic acid (ABA)-disrupted. Totally 202, 78, and 413 differentially expressed (DE) genes were up-regulated, while 429, 72, and 307 were down-regulated at flower induction, flower initiation, and flowering stages, respectively. At the flower induction stage, the down-regulated genes were mainly involved in cell wall, auxin and ABA, light reaction, and abiotic stress, while the up-regulated genes were involved in development, brassinosteroid, major CHO metabolism, and redox. At the flower initiation stage, the down-regulated genes were enriched in light reaction and lipid metabolism, whereas the up-regulated genes were enriched in starch degradation and Ca²⁺ signaling. At the flowering stage, the down-regulated genes were significantly enriched in photosynthesis, gibberellic acid, starch synthesis, nitrogen metabolism, and redox, while the up-regulated genes were enriched in cell wall, jasmonic acid, secondary metabolism, and Ca²⁺ signaling. Besides, 46 transcription factors and 13 flowering-related DE genes were identified. Finally, a possible floral pathway, where LgTEM1, LgSVP, and LgFT1 might play critical roles in SA-induced flowering in L. gibba, was discussed. These findings provide a useful foundation for further investigation of genes and regulatory networks of SA-induced flowering in duckweed.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

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.

Identification and characterization of mRNAs and lncRNAs of a barley shrunken endosperm mutant using RNA-seq

Barley shrunken endosperm mutants have been extensively reported. However, knowledge of the underlying molecular mechanisms of these mutants remains limited. Here, a pair of near isogenic lines (normal endosperm: Bowman and shrunken endosperm: sex1) was subjected to transcriptome analysis to identify mRNAs and lncRNAs related to endosperm development to further dissect its mechanism of molecular regulation. A total of 2123 (1140 up- and 983 down-regulated) unique differentially expressed genes (DEGs) were detected. Functional analyses showed that these DEGs were mainly involved in starch and sucrose metabolism, biosynthesis of secondary metabolites, and plant hormone signal transduction. A total of 343 unique target genes were identified for 57 differentially expressed lncRNAs (DE lncRNAs). These DE lncRNAs were mainly involved in glycerophospholipid metabolism, starch and sucrose metabolism, hormone signal transduction, and stress response. In addition, key lncRNAs were identified by constructing a co-expression network of the target genes of DE lncRNAs. Transcriptome results suggested that mRNA and lncRNA played a critical role in endosperm development. The shrunken endosperm in barley seems to be closely related to plant hormone signal transduction, starch and sucrose metabolism, and cell apoptosis. This study provides a foundation for fine mapping, elucidates the molecular mechanism of shrunken endosperm mutants, and also provides a reference for further studies of lncRNAs during the grain development of plants.

Plant Noncoding RNAs: Hidden Players in Development and Stress Responses

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