Transcriptomic analysis of ncRNAs and mRNAs interactions during drought stress in switchgrass

Switchgrass (Panicum virgatum L.) plays a pivotal role as a bioenergy feedstock in the production of cellulosic ethanol and contributes significantly to enhancing ecological grasslands and soil quality. The utilization of non-coding RNAs (ncRNAs) has gained momentum in deciphering the intricate genetic responses to abiotic stress in various plant species. Nevertheless, the current research landscape lacks a comprehensive exploration of the responses of diverse ncRNAs, including long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs), to drought stress in switchgrass. In this study, we employed whole transcriptome sequencing to comprehensively characterize the expression profiles of both mRNA and ncRNAs during episodes of drought stress in switchgrass. Our analysis identified a total of 12,511 mRNAs, 59 miRNAs, 38 circRNAs, and 368 lncRNAs that exhibited significant differential expression between normal and drought-treated switchgrass leaves. Notably, the majority of up-regulated mRNAs displayed pronounced enrichment within the starch and sucrose metabolism pathway, as validated through KEGG analysis. Co-expression analysis illuminated that differentially expressed (DE) lncRNAs conceivably regulated 1308 protein-coding genes in trans and 7110 protein-coding genes in cis. Furthermore, both cis- and trans-target mRNAs of DE lncRNAs exhibited enrichment in four common KEGG pathways. The intricate interplay between lncRNAs and circRNAs with miRNAs via miRNA response elements was explored within the competitive endogenous RNA (ceRNA) network framework. As a result, we constructed elaborate regulatory networks, including lncRNA-novel_miRNA480-mRNA, lncRNA-novel_miRNA304-mRNA, lncRNA/circRNA-novel_miRNA122-PvSS4, and lncRNA/circRNA-novel_miRNA14-PvSS4, and subsequently validated the functionality of the target gene, starch synthase 4 (PvSS4). Furthermore, through the overexpression of PvSS4, we ascertained its capacity to enhance drought tolerance in yeast. However, it is noteworthy that PvSS4 did not exhibit any discernible impact under salt stress conditions. These findings, as presented herein, not only contribute substantively to our understanding of ceRNA networks but also offer a basis for further investigations into their potential functions in response to drought stress in switchgrass.

Identification and characterization of CircRNA-associated CeRNA networks in moso bamboo under nitrogen stress

BACKGROUND

Nitrogen is a macronutrient element for plant growth and development. Circular RNAs (circRNAs) serve as pivotal regulators for the coordination between nutrient supply and plant demand. Moso bamboo (Phyllostachys edulis) is an excellent plant with fast growth, and the mechanism of the circRNA-target module in response to nitrogen remains unclear.

RESULTS

Deep small RNA sequencing results of moso bamboo seedlings under different concentrations of KNO₃ (N0 = 0 mM, N6 = 6 mM, N18 = 18 mM) were used to identify circRNAs. A total of 549 circRNAs were obtained, of which 309 were generated from corresponding parental coding genes including 66 new ones. A total of 536 circRNA-parent genes were unevenly distributed in 24 scaffolds and were associated with root growth and development. Furthermore, 52 differentially expressed circRNAs (DECs) were obtained, including 24, 33 and 15 DECs from three comparisons of N0 vs. N6, N0 vs. N18 and N6 vs. N18, respectively. Based on integrative analyses of the identified DECs, differentially expressed mRNAs (DEGs), and miRNAs (DEMs), a competitive endogenous RNA (ceRNA) network was constructed, including five DECs, eight DEMs and 32 DEGs. A regulatory module of PeSca_6:12,316,320|12,372,905-novel_miR156-PH02Gene35622 was further verified by qPCR and dual-luciferase reporter assays.

CONCLUSION

The results indicated that circRNAs could participate in multiple biological processes as miRNA sponges, including organ nitrogen compound biosynthesis and metabolic process regulation in moso bamboo. Our results provide valuable information for further study of circRNAs in moso bamboo under fluctuating nitrogen conditions.

Comprehensive Analysis of circRNAs, miRNAs, and mRNAs Expression Profiles and ceRNA Networks in Decidua of Unexplained Recurrent Spontaneous Abortion

The diagnosis and treatment of unexplained recurrent spontaneous abortion (URSA) are subject to debate, because the exact underlying mechanisms remain unclear. To address this issue, we elucidated the expression profiles of dysregulated circRNAs, miRNAs, and mRNAs and constructed circRNA-associated competitive endogenous RNA (ceRNA) networks by comparing the decidua of URSA with that of normal early pregnancy (NEP) using RNA-sequencing. In total, 550 mRNAs, 88 miRNAs, and 139 circRNAs were differentially expressed (DE) in decidua of URSA. Functional annotation revealed that DE mRNAs as well as potential target genes of DE miRNAs and DE circRNAs are mainly involved in immunologic function, such as antigen processing and presentation, allograft rejection, and T cell receptor signaling pathway. In addition, the top hub genes, including CCL4, DDX58, CXCL10, CXCL9, MX1, CD44, RPS2, SOCS3, RPS3A, and CXCL11, were identified. The mRNAs involved in ceRNA network were enriched in complement and coagulation cascades and protein processing in the endoplasmic reticulum. We found that circRNAs in the ceRNA network, which acted as decoys for hsa-miR-204-5p, were positively correlated with MFGE8 expression. Collectively, the results demonstrated that circRNAs, miRNAs, and mRNAs were aberrantly expressed in the decidua of patients with URSA and played a potential role in the development of URSA. Thus, the establishment of the ceRNA network may profoundly affect the diagnosis and therapy of URSA in the future.