Whole-transcriptome sequencing reveals a vernalization-related ceRNA regulatory network in chinese cabbage (Brassica campestris L. ssp. pekinensis)

Source leaves are regulated by sink strengths through non-coding RNAs and alternative polyadenylation in cucumber (Cucumis sativus L.)

BACKGROUND

The yield of major crops is generally limited by sink capacity and source strength. Cucumber is a typical raffinose family oligosaccharides (RFOs)-transporting crop. Non-coding RNAs and alternative polyadenylation (APA) play important roles in the regulation of growth process in plants. However, their roles on the sink‒source regulation have not been demonstrated in RFOs-translocating species.

RESULTS

Here, whole-transcriptome sequencing was applied to compare the leaves of cucumber under different sink strength, that is, no fruit-carrying leaves (NFNLs) and fruit-carrying leaves (FNLs) at 12th node from the bottom. The results show that 1101 differentially expressed (DE) mRNAs, 79 DE long non-coding RNAs (lncRNAs) and 23 DE miRNAs were identified, which were enriched in photosynthesis, energy production and conversion, plant hormone signal transduction, starch and carbohydrate metabolism and protein synthesis pathways. Potential co-expression networks like, DE lncRNAs-DE mRNAs/ DE miRNAs-DE mRNAs, and competing endogenous RNA (ceRNA) regulation models (DE lncRNAs-DE miRNAs-DE mRNAs) associated with sink‒source allocation, were constructed. Furthermore, 37 and 48 DE genes, which enriched in MAPK signaling and plant hormone signal transduction pathway, exist differentially APA, and SPS (CsaV3_2G033300), GBSS1 (CsaV3_5G001560), ERS1 (CsaV3_7G029600), PNO1 (CsaV3_3G003950) and Myb (CsaV3_3G022290) may be regulated by both ncRNAs and APA between FNLs and NFNLs, speculating that ncRNAs and APA are involved in the regulation of gene expression of cucumber sink‒source carbon partitioning.

CONCLUSIONS

These results reveal a comprehensive network among mRNAs, ncRNAs, and APA in cucumber sink-source relationships. Our findings also provide valuable information for further research on the molecular mechanism of ncRNA and APA to enhance cucumber yield.

Regulatory Mechanism of Proanthocyanidins in Grape Peels Using vvi-miR828a and Its Target Gene VvMYBPA1

Anthocyanins and proanthocyanidins are considered to be essential secondary metabolites in grapes and are used to regulate metabolic processes, while miRNAs are involved in their synthesis of anthocyanins and proanthocyanidins to regulate metabolic processes. The present research work was carried out to investigate the underlying regulatory mechanism of target genes in the grape cultivars 'Italia' and 'Benitaka'. miRNA and transnscriptomic sequencing technology were employed to characterize both the profiles of miRNAs and the transcripts of grape peels at 10 and 11 weeks post flowering (10 wpf and 11 wpf). The results revealed that the expression level of vvi-miR828a in 'Italia' at 10 and 11 wpf was significantly higher than that in 'Benitaka'. miRNA-seq analysis predicted MYBPA1 to be the target gene of vvi-miR828a. In transcriptome analysis, the expression level of the VvMYBPA1 gene in 'Benitaka' was significantly higher than that in 'Italia'; in addition, the TPM values (expression levels) of VvMYBPA1 and miR828a also showed an evident negative correlation. The determination of the proanthocyanidin (PA) content in 'Italia' and 'Benitaka' peels at 11 wpf demonstrated that the PA content of 'Benitaka' was significantly higher than that of 'Italia'. The outcomes of RT-qRCR analysis exhibited that the expression levels of the VdPAL, VdCHS, VdCHI, VdDFR, VdMYB5b, VdANR, and VdMYBPA1 genes related anthocyanin and proanthocyanidin pathways were reduced, while the expression levels of all of the above genes were increased after the transient expression of the VvMYBPA1 vector into grape leaves. The results of the transient overexpression experiment of vvi-miR828a before the veraison period of strawberry fruits showed that vvi-miR828a can significantly slow down the coloration of strawberries. The vvi-miR828a negatively regulates the accumulation of proanthocyanidins in grape fruits by inhibiting the expression of VvMYBPA1.

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.

The Sink-Source Relationship in Cucumber (Cucumis sativus L.) Is Modulated by DNA Methylation

The optimization of the sink-source relationship is of great importance for crop yield regulation. Cucumber is a typical raffinose family oligosaccharide (RFO)-transporting crop. DNA methylation is a common epigenetic modification in plants, but its role in sink-source regulation has not been demonstrated in RFO-translocating species. Here, whole-genome bisulfite sequencing (WGBS-seq) was conducted to compare the nonfruiting-node leaves (NFNLs) and leaves of fruit setting (FNLs) at the 12th node by removing all female flowers in other nodes of the two treatments. We found considerable differentially methylated genes enriched in photosynthesis and carbohydrate metabolic processes. Comparative transcriptome analysis between FNLs and NFNLs indicated that many differentially expressed genes (DEGs) with differentially methylated regions were involved in auxin, ethylene and brassinolide metabolism; sucrose metabolism; and RFO synthesis pathways related to sink-source regulation. Moreover, DNA methylation levels of six sink-source-related genes in the pathways mentioned above decreased in leaves after 5-aza-dC-2'-deoxycytidine (5-Aza-dC, a DNA methyltransferase inhibitor) treatment on FNLs, and stachyose synthase (CsSTS) gene expression, enzyme activity and stachyose content in RFO synthesis pathway were upregulated, thereby increasing fruit length and dry weight. Taken together, our findings proposed an up-to-date inference for the potential role of DNA methylation in the sink-source relationship, which will provide important references for further exploring the molecular mechanism of DNA methylation in improving the yield of RFO transport plants.

Comparative RNA-sequencing analysis of the prostate in a mouse model of benign prostatic hyperplasia with bladder outlet obstruction

In ageing men, benign prostatic hyperplasia (BPH) is a chronic disease that leads to progressive lower urinary tract symptoms (LUTS) caused by obstruction of the bladder outlet (BOO). Patients with LUTS (such as increased frequency and urgency of urination) and complications of BOO (such as hydronephrosis and bladder stones) are at risk of serious health problems. BPH causes a rapidly rising burden of LUTS far exceeding that of other urological conditions. Treatment outcomes are unsatisfactory for BPH largely due to the lacking of fully understanding of the pathogenesis. Hormonal imbalances related to androgen and oestrogen can cause BPH, but the exact mechanism is still unknown, even the animal model is not fully understood. Additionally, there are no large-scale data to explain this mechanism. A BPH mouse model was established using mixed slow-release pellets of testosterone (T) and estradiol (E2), and we measured gene expression in mouse prostate tissue using RNA-seq, verified the results using qRT‒PCR, and used bioinformatics methods to analyse the differentially expressed genes (DEGs).

Identification and Validation of Genes Exhibiting Dynamic Alterations in Response to Bleomycin-Induced Pulmonary Fibrosis

Idiopathic pulmonary fibrosis (IPF) carries a high mortality rate and has a poor prognosis. The pathogenesis of pulmonary fibrosis (PF) is highly related to dysregulation of multiple RNAs. This study aims to identify and validate dysregulated RNAs that exhibited dynamic alterations in response to bleomycin (BLM)-induced PF. The results will provide therapeutic targets for patients suffering from IPF. Whole transcriptomic profiles of BLM-induced PF were obtained through high-throughput RNA sequencing. miRNA profiling was downloaded from GSE45789 database in the Gene Expression Omnibus (GEO). We identified the differentially expressed RNAs (DERNAs) that exhibited dynamic alterations in response to BLM-induced PF. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway enrichment analysis were conducted to discovery regulatory processes of PF. Weighted gene co-expression network analysis (WGCNA), protein-protein interaction (PPI) analysis, and co-expression analysis were performed to identify key genes and pathogenic pattern during the progression of PF. MiRanda, miRcode, and TargetScan were utilized to predict target relationships in the potential competing endogenous RNA (ceRNA) network. The results were verified by qRT-PCR analysis. In the context of BLM-induced PF, this study identified a total of 167 differentially expressed messenger RNAs (DEmRNAs), 115 differentially expressed long non-coding RNAs (DElncRNAs), 45 differentially expressed circular RNAs (DEcircRNAs), and 87 differentially expressed microRNAs (DEmiRNAs). These RNA molecules showed dynamic alterations in response to BLM-induced PF. These DEmRNAs exhibited a predominant association with the biological processes pertaining to the organization of extracellular matrix. A regulatory network was built in PF, encompassing 31 DEmRNAs, 18 DE lncRNAs, 13 DEcircRNAs, and 13 DEmiRNAs. Several DERNA molecules were subjected to validate using additional BLM-induced PF model. The outcomes of this validation process shown a strong correlation with the results obtained from RNA sequencing analysis. The GSE213001 dataset was utilized to validate the expression levels and diagnostic efficacy of four specific hub mRNAs (CCDC80, CLU, COL5A1, and COL6A3) in individuals diagnosed with PF. In this study, we identified and validated several RNA molecules that exhibited dynamic alternations in response to BLM-induced PF. These dysregulated RNAs participated in the pathogenesis of PF and can be used as therapeutic targets for early-stage IPF. Although more work must be done to confirm the results, our study may provide directions for future studies.

Whole transcriptome analysis and construction of a ceRNA regulatory network related to leaf and petiole development in Chinese cabbage (Brassica campestris L. ssp. pekinensis)

BACKGROUND

The growth and development of leaves and petioles have a significant effect on photosynthesis. Understanding the molecular mechanisms underlying leaf and petiole development is necessary for improving photosynthetic efficiency, cultivating varieties with high photosynthetic efficiency, and improving the yield of crops of which the leaves are foodstuffs. This study aimed to identify the mRNAs, long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs) related to leaf and petiole development in Chinese cabbage (Brassica campestris L. ssp. pekinensis). The data were used to construct a competitive endogenous RNA (ceRNA) network to obtain insights into the mechanisms underlying leaf and petiole development.

RESULTS

The leaves and petioles of the 'PHL' inbred line of Chinese cabbage were used as research materials for whole transcriptome sequencing. A total of 10,646 differentially expressed (DE) mRNAs, 303 DElncRNAs, 7 DEcircRNAs, and 195 DEmiRNAs were identified between leaves and petioles. Transcription factors and proteins that play important roles in leaf and petiole development were identified, including xyloglucan endotransglucosylase/hydrolase, expansion proteins and their precursors, transcription factors TCP15 and bHLH, lateral organ boundary domain protein, cellulose synthase, MOR1-like protein, and proteins related to plant hormone biosynthesis. A ceRNA regulatory network related to leaf and petiole development was constructed, and 85 pairs of ceRNA relationships were identified, including 71 DEmiRNA-DEmRNA, 12 DEmiRNA-DElncRNA, and 2 DEmiRNA-DEcircRNA pairs. Three LSH genes (BrLSH1, BrLSH2 and BrLSH3) with significant differential expression between leaves and petioles were screened from transcriptome data, and their functions were explored through subcellular localization analysis and transgenic overexpression verification. BrLSH1, BrLSH2 and BrLSH3 were nuclear proteins, and BrLSH2 inhibited the growth and development of Arabidopsis thaliana.

CONCLUSIONS

This study identifies mRNAs and non-coding RNAs that may be involved in the development of leaves and petioles in Chinese cabbage, and establishes a ceRNA regulatory network related to development of the leaves and petioles, providing valuable genomic resources for further research on the molecular mechanisms underlying leaf and petiole development in this crop species.

Whole-transcriptome analyses identify key differentially expressed mRNAs, lncRNAs, and miRNAs associated with male sterility in watermelon

Male sterility is a valuable trait for watermelon breeding, as watermelon hybrids exhibit obvious heterosis. However, the underlying regulatory mechanism is still largely unknown, especially regarding the related non-coding genes. In the present study, approximately 1035 differentially expressed genes (DEGs), as well as 80 DE-lncRNAs and 10 DE-miRNAs, were identified, with the overwhelming majority down-regulated in male-sterile floral buds. Enrichment analyses revealed that the general phenylpropanoid pathway as well as its related metabolisms was predicted to be altered in a mutant compared to its fertile progenitor. Meanwhile, the conserved genetic pathway DYT1-TDF1-AMS-MS188-MS1, as well as the causal gene ClAMT1 for the male-sterile mutant Se18, was substantially disrupted during male reproductive development. In addition, some targets of the key regulators AMS and MS188 in tapetum development were also down-regulated at a transcriptional level, such as ABCG26 (Cla004479), ACOS5 (Cla022956), CYP703A2 (Cla021151), PKSA (Cla021099), and TKPR1 (Cla002563). Considering lncRNAs may act as functional endogenous target mimics of miRNAs, competitive endogenous RNA networks were subsequently constructed, with the most complex one containing three DE-miRNAs, two DE-lncRNAs, and 21 DEGs. Collectively, these findings not only contribute to a better understanding of genetic regulatory networks underlying male sterility in watermelon, but also provide valuable candidates for future research.

Identification, biogenesis, function, and mechanism of action of circular RNAs in plants

Circular RNAs (circRNAs) are a class of single-stranded, closed RNA molecules with unique functions that are ubiquitously expressed in all eukaryotes. The biogenesis of circRNAs is regulated by specific cis-acting elements and trans-acting factors in humans and animals. circRNAs mainly exert their biological functions by acting as microRNA sponges, forming R-loops, interacting with RNA-binding proteins, or being translated into polypeptides or proteins in human and animal cells. Genome-wide identification of circRNAs has been performed in multiple plant species, and the results suggest that circRNAs are abundant and ubiquitously expressed in plants. There is emerging compelling evidence to suggest that circRNAs play essential roles during plant growth and development as well as in the responses to biotic and abiotic stress. However, compared with recent advances in human and animal systems, the roles of most circRNAs in plants are unclear at present. Here we review the identification, biogenesis, function, and mechanism of action of plant circRNAs, which will provide a fundamental understanding of the characteristics and complexity of circRNAs in plants.

Whole-transcriptome analyses of Sorghum leaves identify key mRNAs and ncRNAs associated with GA3-mediated alleviation of salt stress

Sorghum has recently attracted much attention for its tolerance in high salt environment. However, the effect and regulatory mechanism of the gibberellic acid (GA₃)-mediated alleviation of salt stress in sorghum remains unclear. Herein, we reported that a GA₃ concentration of 50 mg/L is optimal for sorghum ("Jitian 3") development under salt stress. We conducted a whole-transcriptome analysis between GA₃-treated and control sorghum leaves under salt stress, and we identified 1002 differentially expressed (DE)-messenger RNAs (mRNAs), 81 DE-long non-coding RNAs (lncRNAs), 7 DE-circular RNAs (circRNAs), and 26 DE-microRNA (miRNAs) in sorghum following GA₃ treatment. We also identified a majority of DE-mRNAs and non-coding RNAs (ncRNAs) targets that serve essential roles in phenylpropanoid biosynthesis and plant hormone networks. In addition, we generated a competitive endogenous RNA (ceRNA)-miRNA-target gene network, and 3 circRNAs (circRNA_2746, circRNA_6515, circRNA_5622), 4 lncRNAs (XR_002450182.1, XR_002452422.1, XR_002448510.1, XR_002448296.1) and 4 genes (LOC8056546, LOC8062245, LOC8061469, LOC8071960) probably act as valuable candidates for the regulation of the GA₃-mediated alleviation of salt stress in sorghum. Our findings uncovered potential mRNA and non-coding RNAs that contribute to GA₃ regulation, thus offering a basis for the future investigation of underlying mechanisms of salt stress in sorghum.

Increase Crop Resilience to Heat Stress Using Omic Strategies

Varieties of various crops with high resilience are urgently needed to feed the increased population in climate change conditions. Human activities and climate change have led to frequent and strong weather fluctuation, which cause various abiotic stresses to crops. The understanding of crops' responses to abiotic stresses in different aspects including genes, RNAs, proteins, metabolites, and phenotypes can facilitate crop breeding. Using multi-omics methods, mainly genomics, transcriptomics, proteomics, metabolomics, and phenomics, to study crops' responses to abiotic stresses will generate a better, deeper, and more comprehensive understanding. More importantly, multi-omics can provide multiple layers of information on biological data to understand plant biology, which will open windows for new opportunities to improve crop resilience and tolerance. However, the opportunities and challenges coexist. Interpretation of the multidimensional data from multi-omics and translation of the data into biological meaningful context remained a challenge. More reasonable experimental designs starting from sowing seed, cultivating the plant, and collecting and extracting samples were necessary for a multi-omics study as the first step. The normalization, transformation, and scaling of single-omics data should consider the integration of multi-omics. This review reports the current study of crops at abiotic stresses in particular heat stress using omics, which will help to accelerate crop improvement to better tolerate and adapt to climate change.