TargetFinder: a software for antisense oligonucleotide target site selection based on MAST and secondary structures of target mRNA

Scenarios for the emergence of new microRNA genes in the plant Arabidopsis halleri

MicroRNAs (miRNAs) are central players in the regulation of gene expression in eukaryotes. The repertoires of miRNA genes vary drastically even among closely related species, indicating that they are evolutionarily labile. However, the processes by which they originate over the course of evolution and the nature of their progenitors across the genome remain poorly understood. Here, we analyzed miRNA genes in Arabidopsis halleri, a plant species where we recently documented a large number of species-specific miRNA genes, likely to represent recent events of emergence. Analysis of sequence homology across the genome indicates that a diversity of sources contributes to the emergence of new miRNA genes, including inverted duplications from protein-coding genes, rearrangements of transposable element (TE) sequences, and duplications of preexisting miRNA genes. Our observations indicate that the origin from protein-coding genes was less common than was previously considered. In contrast, we estimate that almost half of the new miRNA genes likely emerged from TEs. Miniature inverted-repeat TEs (MITEs) seem to be particularly important contributors to new miRNA genes, with the Harbinger and Mariner TE superfamilies representing disproportionate sources for their emergence. We further analyzed the recent expansion of a miRNA family derived from MuDR elements and the duplication of miRNA genes formed by two hAT transposons. Overall, our results illustrate the rapid pace at which new regulatory elements can arise from the modification of preexisting sequences in a genome and highlight the central role of certain categories of TEs in this process.

Transcriptomic analysis of two Chinese wheat landraces with contrasting Fusarium head blight resistance reveals miRNA-mediated defense mechanisms

Introduction

Fusarium head blight (FHB), caused primarily by Fusarium graminearum (Fg), poses a significant threat to wheat production. It is necessary to deeply understand the molecular mechanisms underlying FHB resistance in wheat breeding.

Methods

In this study, the transcriptomic responses of two Chinese wheat landraces-Wuyangmai (WY, resistant) and Chinese Spring (CS, susceptible)-to F. graminearum infection were examined using RNA sequencing (RNA-seq). Differential expression of mRNAs, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) was analyzed at 3 and 5 days post-Fg inoculation (dpi).

Results

The results showed that WY exhibited a targeted miRNA response, primarily modulating defense-related pathways such as glutathione metabolism and phenylpropanoid biosynthesis, which are crucial for oxidative stress regulation and pathogen defense response. In contrast, CS displayed a broader transcriptional response, largely linked to general metabolic processes rather than immune activation. Notably, the up-regulation of genes involved in oxidative stress and immune defense in WY confirmed its enhanced resistance to FHB. The integrated analysis of miRNA-mRNA interactions highlighted miRNAs as central regulators of defense mechanisms in WY, particularly at later stages of infection. These miRNAs targeted genes involved in immune responses, while lncRNAs and circRNAs played a more limited role in the regulation of defense responses. The GO and KEGG pathway enrichment analyses further revealed that WY enriched for plant-pathogen interaction and secondary metabolite biosynthesis pathways, which are crucial for pathogen resistance. In contrast, CS prioritized metabolic homeostasis, suggesting a less effective defense strategy.

Discussion

Overall, this study underscores the critical role of miRNA-mediated regulation in FHB resistance in WY. These insights into miRNA-mediated regulatory mechanisms provide a molecular basis for breeding FHB-resistant wheat varieties and highlight miRNA-mRNA interactions as promising targets for enhancing disease resilience.

Characterization of Exosome-like Nanoparticles from Saffron Tepals and Their Immunostimulatory Activity

Plant exosomes exhibit high stability and easy absorption, and have emerged as promising bioactive tools due to their potential health benefits and biomedical applications. Saffron tepals contain abundant metabolites with potential therapeutic properties and were used for exosome extraction by ultracentrifugation and gradient purification. The exosomes showed an average particle size of 151.5 ± 79.6 nm and exhibited a spherical morphology. Five well-conserved miRNAs-miR157, miR166, miR168, miR396, and miR398-were identified in the exosomes, which are involved in the coordination of growth and physiological plant responses with endogenous and environmental abiotic and biotic signals, and their potential targets in mammals are upregulated in specific cancer types and associated with inflammation. Proteome analysis revealed an enrichment of proteasome proteins, ribosomal proteins, and proteins involved in the cytoskeleton, transport across the membrane (ABC transporters), and vesicle trafficking (RAB GTPases, TM9SF and Coatomer subunits). Metabolite analyses showed mainly anthocyanins. The exosomes have selective stimulatory activity on macrophages, increasing the expression of surface molecules (CD80 and CD86), and cytokines (IL-1β, IL-6, and TNF-α), but not the levels of IL-10. Overall, these results indicated that saffron flowers are an effective and abundant source of exosomes as new nanomedicines for human health.

Deciphering the impact of microRNAs in plant biology: a review of computational insights and experimental validation

Exploring the complex world of microRNA (miRNA) biogenesis and functions in plants is essential for understanding their diverse regulatory mechanisms. This review highlights the processes involved in miRNA biogenesis and their crucial roles in growth and development of plant, stress responses, and nutrient homeostasis. miRNAs play a central role in various developmental processes, including the transition from the juvenile to adult stage, the growth of shoot apical meristem, leaf and floral morphogenesis, and the determination of flowering time. By presenting the current state of research, we focus on the vital role of computational tools and databases in deciphering the regulatory networks controlled by miRNAs, which helps us navigate the intricate world of plant biology. Furthermore, it stresses the importance of experimental validation techniques in confirming computational predictions, ensuring that miRNA research is reliable and robust. As the field continues to grow, this review emphasizes the urgent need for integrated approaches, to deepen our knowledge of plant miRNA biology and its implications. These insights will pave the way for advancements in crop improvement, stress resilience, and biotechnological innovations.

Transcriptome analysis and functional identification of transfer RNA-derived fragments in grape leaves exposed to UV-C radiation

Transfer RNA-derived fragments (tRFs) are noncoding small RNAs derived from transfer RNAs (tRNAs) in microorganisms, animals and plants. In plants, tRFs are known to respond to environmental stimuli, including heat, oxidative stress and UV radiation; however, their specific functions in horticultural plants, such as grapevine, remain poorly understood. In this study, we used RNA-seq to identify differentially expressed genes (DEGs) in grape leaves exposed to UV-C radiation. A total of 1329 and 8055 of genes were differentially expression after 1 and 6 h of UV-C treatment, respectively. We identified a large number of secondary metabolism-related genes in the DEGs, including genes involved in stilbene and flavonoid biosynthesis. Noticeably, the stilbene biosynthesis-related gene was induced earlier than the other genes in the phenylalanine metabolic pathway. We also conducted small RNA-seq and identified differentially expressed (DE) miRNAs and their targets. To explore whether the tRFs involved in UV-C response, further analysis of the small RNA-seq data revealed 23 down-regulated and 41 up-regulated DE tRFs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that the target genes of these tRFs are involved in multiple biological processing, including hormone signal transduction and metabolite synthesis. To validate the function of tRFs, tRF39 and tRF45 were selected and overexpressed in tobacco leaves, and the expression levels of their target genes were inhibited. Our study suggests that the tRFs may regulate multiple biological processes in response to UV-C exposure in grapevine. Our findings provide a foundation for further elucidating the regulatory mechanisms of tRFs in horticultural crops.

Circ_0001944 Targets the miR-1292-5p/FBLN2 Axis to Facilitate Sorafenib Resistance in Hepatocellular Carcinoma by Impeding Ferroptosis

BACKGROUND

Sorafenib, an orally active potent tyrosine kinase inhibitor (TKI), represented a primary treatment in patients with advanced hepatocellular carcinoma (HCC). Unfortunately, sorafenib resistance was regarded as a huge obstacle for HCC treatment.

METHODS

RNA-sequencing including circRNA Sequencing (circRNA-Seq) for circular RNAs (circRNAs), miRNA Sequencing (miRNA-Seq) for microRNAs (miRNAs), as well as mRNA Sequencing (mRNA-Seq) for mRNAs in sorafenib-resistant HCC cells vs sorafenib-sensitive HCC cells, were performed. Then, interaction correlation analysis between differentially expressed circRNAs and miRNAs and their target genes in Huh7/SOR and SMMC7721/SOR cells was exhibited. The "circRNA-miRNA-mRNA" network was constructed through the Cytoscape software application, Circular RNA Interactome and Targetscan prediction, RNA binding protein immunoprecipitation (RIP), RNA pull-down, and Dual luciferase reporter assay. Furthermore, mRNA-Seq, Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for the downstream genes involved in the "circRNA-miRNA-mRNA" network was implemented. Iron detection assay, Lipid peroxidation quantification assay, ROS measurement assay, CCK-8 assay, and tumor challenge in vivo were used to determine the mechanisms promoting sorafenib resistance in HCC, where the "circRNA-miRNA-mRNA" network is clearly involved in.

RESULTS 

circ_0001944 and circ_0078607 with upregulation and 2 downregulated expressed circRNAs (circ_0002874 and circ_0069981), as well as 11 upregulated miRNAs including miR-193a-5p, miR-197-3p, miR-27a-5p, miR-551b-5p, miR-335-3p, miR-767-3p, miR-767-5p, miR-92a-1-5p, miR-92a-3p, miR-3940-3p, and miR-664b-3p and 3 downregulated expressed miRNAs (miR-1292-5p, let-7c-5p, and miR-99a-5p) in sorafenib-resistant HCC cells were determined. Among these non-coding RNAs (ncRNAs), circ_0001944 and miR-1292-5p should not be drop out of sight; circ_0001944 has been proved to target miR-1292-5p to inhibit its expression in HCC. Subsequent findings also raise that miR-1292-5p directly targeted the 3'-noncoding region (3'-UTR) of Fibulin 2 (FBLN2) mRNA. Furthermore, circ_0001944 targets the miR-1292-5p/FBLN2 axis to inhibit cell ferroptosis in which the indicated regulators associated with iron overload and lipid peroxidation were "rearranged". Most importantly, circ_0001944 advanced sorafenib resistance in HCC through mitigating ferroptosis, where the miR-1292-5p/FBLN2 axis cannot be left unrecognized.

CONCLUSION 

Circ_0001944 is a putative target for reversing sorafenib resistance in HCC. Our findings are expected to provide new targets and new directions for sorafenib sensitization in the treatment of HCC.

Unveiling RNA structure-mediated regulations of RNA stability in wheat

Despite the critical role of mRNA stability in post-transcriptional gene regulation, research on this topic in wheat, a vital agricultural crop, remains unclear. Our study investigated the mRNA decay landscape of durum wheat (Triticum turgidum L. ssp. durum, BBAA), revealing subgenomic asymmetry in mRNA stability and its impact on steady-state mRNA abundance. Our findings indicate that the 3' UTR structure and homoeolog preference for RNA structural motifs can influence mRNA stability, leading to subgenomic RNA decay imbalance. Furthermore, single-nucleotide variations (SNVs) selected for RNA structural motifs during domestication can cause variations in subgenomic mRNA stability and subsequent changes in steady-state expression levels. Our research on the transcriptome stability of polyploid wheat highlights the regulatory role of non-coding region structures in mRNA stability, and how domestication shaped RNA structure, altering subgenomic mRNA stability. These results illustrate the importance of RNA structure-mediated post-transcriptional gene regulation in wheat and pave the way for its potential use in crop improvement.

Clinical Applications of Antisense Oligonucleotides in Cancer: A Focus on Glioblastoma

Antisense oligonucleotides (ASOs) are promising drugs capable of modulating the protein expression of virtually any target with high specificity and high affinity through complementary base pairing. However, this requires a deep understanding of the target sequence and significant effort in designing the correct complementary drug. In addition, ASOs have been demonstrated to be well tolerated during their clinical use. Indeed, they are already used in many diseases due to pathogenic RNAs of known sequences and in several neurodegenerative diseases and metabolic diseases, for which they were given marketing authorizations (MAs) in Europe and the United States. Their use in oncology is gaining momentum with several identified targets, promising preclinical and clinical results, and recent market authorizations in the US. However, many challenges remain for their clinical use in cancer. It seems necessary to take a step back and review our knowledge of ASOs and their therapeutic uses in oncology. The objectives of this review are (i) to summarize the current state of the art of ASOs; (ii) to discuss the therapeutic use of ASOs in cancer; and (iii) to focus on ASO usage in glioblastoma, the challenges, and the perspective ahead.

Comparative genome-wide analysis of circular RNAs in Brassica napus L.: target-site versus non-target-site resistance to herbicide stress

Circular RNAs (circRNAs), a class of non-coding RNA molecules, are recognized for their unique functions; however, their responses to herbicide stress in Brassica napus remain unclear. In this study, the role of circRNAs in response to herbicide treatment was investigated in two rapeseed cultivars: MH33, which confers non-target-site resistance (NTSR), and EM28, which exhibits target-site resistance (TSR). The genome-wide circRNA profiles of herbicide-stressed and non-stressed seedlings were analyzed. The findings indicate that NTSR seedlings exhibited a greater abundance of circRNAs, shorter lengths of circRNAs and their parent genes, and more diverse functions of parent genes compared with TSR seedlings. Compared to normal-growth plants, the herbicide-stressed group exhibited similar trends in the number of circRNAs, functions of parent genes, and differentially expressed circRNAs as observed in NTSR seedlings. In addition, a greater number of circRNAs that function as competing microRNA (miRNA) sponges were identified in the herbicide stress and NTSR groups compared to the normal-growth and TSR groups, respectively. The differentially expressed circRNAs were validated by qPCR. The differntially expressed circRNA-miRNA networks were predicted, and the mRNAs targeted by these miRNAs were annotated. Our results suggest that circRNAs play a crucial role in responding to herbicide stress, exhibiting distinct responses between NTSR and TSR in rapeseed. These findings offer valuable insights into the mechanisms underlying herbicide resistance in rapeseed.

Analysis of differential mRNA and miRNA expression induced by heterogeneous grafting in Gleditsia sinensis

Gleditsia sinensis Lam. is a multifaceted plant with medicinal, edible, chemical, timber, and ornamental applications. However, the effect of rootstocks on scions after grafting is still unclear. This study examined the mRNA and miRNA transcriptome among homografts, heterografts, and seedlings. GO enrichment analysis between seedlings and homograft/heterograft combinations revealed that biosynthesis, degradation, and transport were enriched. The KEGG enrichment results showed that plant hormone signal transduction and the plant MAPK signaling pathway were enriched in both seedlings and heterograft combinations. Through weighted correlation network analysis (WGCNA), the hub genes related to the content of plant hormones were obtained. Taking G. sinensis as the scion, there were 4594, 2887, 3429, and 5959 mRNAs that were specifically expressed in the grafted plants of G. sinensis/G. fera, G. sinensis/G. delavayi, G. sinensis/G. microphylla, and G. sinensis/G. japonica, respectively. The specifically expressed mRNA genes may participate in such processes and pathways as the rhythmic process, circadian rhythm, gibberellic-acid-mediated signaling pathway, and peptide-based amino acid modification. Additionally, 3, 16, 2, and 15 specifically expressed miRNAs were identified. This study examines the impact of grafting on gene expression in Gleditsia plants and establishes a foundation for the development of new resources and rootstock breeding.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

Deciphering the role of non-coding RNAs involved in sorafenib resistance

Sorafenib is an important treatment strategy for advanced hepatocellular carcinoma (HCC). Unfortunately, drug resistance has become a major obstacle in sorafenib application. In this study, whole transcriptome sequencing (WTS) was conducted to compare the paired differences between non-coding RNAs (ncRNAs), including long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), microRNAs (miRNAs), and mRNAs, in sorafenib-resistant and parental cells. The overlap of differentially expressed ncRNAs (DENs) between the SMMC7721/S and Huh7/S cells and their parental cells was determined. 2 upregulated and 3 downregulated lncRNAs, 2 upregulated and 1 downregulated circRNAs, as well as 10 upregulated and 2 downregulated miRNAs, in both SMMC7721/S and Huh7/S cells, attracted more attention. The target genes of these DENs were then identified as the overlaps between the differentially expressed mRNAs achieved using the WTS analysis and the predicted genes of DENs obtained using the "co-localization" or "co-expression," miRanda, and RNAhybrid analysis. Consequently, the potential regulatory network between overlapping DENs and their target genes in both SMMC7721/S and Huh7/S cells was explored. The "lncRNA-miRNA-mRNA" and "circRNA-miRNA-mRNA" networks were constructed based on the competitive endogenous RNA (ceRNA) theory using the Cytoscape software. In particular, lncRNA MED17-203-miRNA (miR-193a-5p, miR-197-3p, miR-27a-5p, miR-320b, miR-767-3p, miR-767-5p, miR-92a-3p, let-7c-5p)-mRNA," "circ_0002874-miR-27a-5p-mRNA" and "circ_0078607-miR-320b-mRNA" networks were first introduced in sorafenib-resistant HCC. Furthermore, these networks were most probably connected to the process of metabolic reprogramming, where the activation of the PPAR, HIF-1, Hippo, and TGF-β signaling pathways is governed. Alternatively, the network "circ_0002874-miR-27a-5p-mRNA" was also involved in the regulation of the activation of TGF-β signaling pathways, thus advancing Epithelial-mesenchymal transition (EMT). These findings provide a theoretical basis for exploring the mechanisms underlying sorafenib resistance mediated by metabolic reprogramming and EMT in HCC.

Whole-transcriptome profiling and identification of cold tolerance-related ceRNA networks in japonica rice varieties

Introduction

Low-temperature stress negatively impacts rice yield, posing a significant risk to food security. While previous studies have explored the physiological and linear gene expression alterations in rice under low-temperature conditions, the changes in competing endogenous RNA (ceRNA) networks remain largely unexamined.

Methods

We conducted RNA sequencing on two japonica rice varieties with differing cold-tolerance capabilities to establish ceRNA networks. This enabled us to investigate the transcriptional regulatory network and molecular mechanisms that rice employs in response to low-temperature stress.

Results

We identified 364 differentially expressed circular RNAs (circRNAs), 224 differentially expressed microRNAs (miRNAs), and 12,183 differentially expressed messenger RNAs (mRNAs). WRKY family was the most prominent transcription factor family involved in cold tolerance. Based on the expression patterns and targeted relationships of these differentially expressed RNAs, we discerned five potential ceRNA networks related to low-temperature stress in rice: osa-miR166j-5p from the miR166 family was associated with cold tolerance; osa-miR528-3p and osa-miR156j-3p were linked to stress response; and osa-miR156j-3p was involved in the antioxidant system. In addition, Os03g0152000 in the antioxidant system, as well as Os12g0491800 and Os05g0381400, correlated with the corresponding stress response and circRNAs in the network. A gene sequence difference analysis and phenotypic validation of Os11g0685700 (OsWRKY61) within the WRKY family suggested its potential role in regulating cold tolerance in rice.

Discussion and conclusion

We identified Os11g0685700 (OsWRKY61) as a promising candidate gene for enhancing cold tolerance in japonica rice. The candidate miRNAs, mRNAs, and circRNAs uncovered in this study are valuable targets for researchers and breeders. Our findings will facilitate the development of cold-tolerant rice varieties from multiple angles and provide critical directions for future research into the functions of cold-tolerance-related miRNAs, mRNAs, and circRNAs in rice.

Identification and characterization of circular RNAs involved in the fertility stability of cotton CMS-D2 restorer line under heat stress

BACKGROUND

As a vital type of noncoding RNAs, circular RNAs (circRNAs) play important roles in plant growth and development and stress response. However, little is known about the biological roles of circRNAs in regulating the stability of male fertility restoration for cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) cotton under high-temperature (HT) stress.

RESULTS

In this study, RNA-sequencing and bioinformatics analysis were performed on pollen grains of isonuclear alloplasmic near-isogenic restorer lines NH [N(Rf1rf1)] and SH [S(Rf1rf1)] with obvious differences in fertility stability under HT stress at two environments. A total of 967 circRNAs were identified, with 250 differentially expressed under HT stress. We confirmed the back-splicing sites of eight selected circRNAs using divergent primers and Sanger sequencing. Tissue-specific expression patterns of five differentially expressed circRNAs (DECs) were also verified by RT-PCR and qRT-PCR. Functional enrichment and metabolic pathway analysis revealed that the parental genes of DECs were significantly enriched in fertility-related biological processes such as pollen tube guidance and cell wall organization, as well as the Pentose and glucuronate interconversions, Steroid biosynthesis, and N-Glycan biosynthesis pathways. Moreover, we also constructed a putative circRNA-mediated competing endogenous RNA (ceRNA) network consisting of 21 DECs, eight predicted circRNA-binding miRNAs, and their corresponding 22 mRNA targets, especially the two ceRNA modules circRNA346-miR159a-MYB33 and circRNA484-miR319e-MYB33, which might play important biological roles in regulating pollen fertility stability of cotton CMS-D2 restorer line under HT stress.

CONCLUSIONS

Through systematic analysis of the abundance, characteristics and expression patterns of circRNAs, as well as the potential functions of their parent genes, our findings suggested that circRNAs and their mediated ceRNA networks acted vital biological roles in cotton pollen development, and might be also essential regulators for fertility stability of CMS-D2 restorer line under heat stress. This study will open a new door for further unlocking complex regulatory mechanisms underpinning the fertility restoration stability for CMS-D2 in cotton.

Overexpression of Nta-miR6155 confers resistance to Phytophthora nicotianae and regulates growth in tobacco (Nicotiana tabacum L.)

Tobacco black shank induced by Phytophthora nicotianae causes significant yield losses in tobacco plants. MicroRNAs (miRNAs) play a pivotal role in plant biotic stress responses and have great potential in tobacco breeding for disease resistance. However, the roles of miRNAs in tobacco plants in response to P. nicotianae infection has not been well characterized. In this study, we found that Nta-miR6155, a miRNA specific to Solanaceae crops, was significantly induced in P. nicotianae infected tobacco. Some of predicted target genes of Nta-miR6155 were also observed to be involved in disease resistance. To further investigate the function of miR6155 in tobacco during P. nicotianae infection, Nta-miR6155 overexpression plants (miR6155-OE) were generated in the Honghua Dajinyuan tobacco variety (HD, the main cultivated tobacco variety in China). We found that the Nta-miR6155 overexpression enhanced the resistance in tobacco towards P. nicotianae infections. The level of reactive oxygen species (ROS) was significantly lower and antioxidant enzyme activities were significantly higher in miR6155-OE plants than those in control HD plants during P. nicotianae infection. In addition, we found that the accumulation of salicylic acid and the expression of salicylic acid biosynthesis and signal transduction-related genes is significantly higher in miR6155-OE plants in comparison to the control HD plants. Furthermore, we found that Nta-miR6155 cleaved target genes NtCIPK18 to modulate resistance towards P. nicotianae in tobacco plants. Additionally, phenotypic analysis of miR6155-OE plants showed that Nta-miR6155 could inhibit the growth of tobacco by suppressing nitrogen uptake and photosynthesis. In conclusion, our findings indicated that miR6155 plays a crucial role in the regulation of growth and resistance against P. nicotianae infections in tobacco plants.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

RNA-Sequencing Analysis Indicates That N-Cadherin Promotes Prostate Cancer Progression by the Epigenetic Modification of Key Genes

N-cadherin (cadherin-2 [CDH2]) is widely known as the promoter of prostate cancer (PCa) invasion and castration resistance. However, the biological mechanism of N-cadherin in PCa progression is unclear. In this study, we overexpressed N-cadherin in LNCaP cells and downregulated N-cadherin in PC3 cells by lentiviral transduction. Then, differentially expressed genes (DEGs) and dysregulated biological functions were investigated through RNA sequencing (RNA-seq) analyses. We found 13 long noncoding RNA (lncRNA) transcripts, 72 messenger RNA (mRNA) transcripts, and 3 integrated genes were dysregulated by N-cadherin. In the disease enrichment, bone cancer, and neurodegenerative and nervous system diseases were associated with N-cadherin in the circular RNA (circRNA; PC3 versus [vs.,/] LNCaP [PC3/LNCaP] comparison) and DEG analysis (LNCaP_oe_CDH2 vs. LNCaP_oe_NC [LNCaP_oe_CDH2/NC] comparison). Epigenetic reprogramming, such as nucleic acid binding, and chromatin and histone modifications, was enriched in Gene Ontology (GO) analysis (DEGs in LNCaP_oe_CDH2/NC and PC3_sh_NC/CDH2, and host genes of circRNA in PC3/LNCaP). Transcriptional misregulation in cancer, post-translational protein modification, gene expression, and generic transcription pathways were dysregulated in the pathway enrichment analysis (host genes of circRNA in PC3/LNCaP, and DEGs in LNCaP_oe_CDH2/NC and PC3_sh_NC/CDH2). Verifying DEGs through TCGA-PRAD dataset revealed six oncogenes (ARHGEF1, GRAMD1A, GTF2H4, MAPK8IP3, POLD1, and PTBP1) that were commonly upregulated by N-cadherin and in advanced PCa stages. In summary, we identified several oncogenes and biological functions associated with N-cadherin expression in PCa cells. N-cadherin may trigger epigenetic reprogramming in PCa cells to promote tumor progression.

Small RNA-Seq to Unveil the miRNA Expression Patterns and Identify the Target Genes in Panax ginseng

Panax ginseng, renowned for its medicinal properties, relies on adventitious roots and hairy roots as crucial sources for the production of ginsenosides. Despite the widespread utilization of ginseng, investigations into its miRNAs have remained scarce. To address this gap, two samples of ginseng adventitious roots and ginseng hairy roots were collected, and subsequent construction and sequencing of small RNA libraries of ginseng adventitious roots and hairy roots were performed using the Illumina HiSeq X Ten platform. The analysis of the sequencing data unveiled total miRNAs 2432. The miR166 and miR396 were the most highly expressed miRNA families in ginseng. The miRNA expression analysis results were used to validate the qRT-PCR. Target genes of miRNA were predicted and GO function annotation and KEGG pathway analysis were performed on target genes. It was found that miRNAs are mainly involved in synthetic pathways and biological processes in plants, which include metabolic and bioregulatory processes. The plant miRNAs enriched KEGG pathways are associated with some metabolism, especially amino acid metabolism and carbohydrate metabolism. These results provide valuable insights miRNAs and their roles in metabolic processes in ginseng.

Interactome of Long Non-Coding RNAs: Transcriptomic Expression Patterns and Shaping Cancer Cell Phenotypes

RNA biology has gained extensive recognition in the last two decades due to the identification of novel transcriptomic elements and molecular functions. Cancer arises, in part, due to the accumulation of mutations that greatly contribute to genomic instability. However, the identification of differential gene expression patterns of wild-type loci has exceeded the boundaries of mutational study and has significantly contributed to the identification of molecular mechanisms that drive carcinogenic transformation. Non-coding RNA molecules have provided a novel avenue of exploration, providing additional routes for evaluating genomic and epigenomic regulation. Of particular focus, long non-coding RNA molecule expression has been demonstrated to govern and direct cellular activity, thus evidencing a correlation between aberrant long non-coding RNA expression and the pathological transformation of cells. lncRNA classification, structure, function, and therapeutic utilization have expanded cancer studies and molecular targeting, and understanding the lncRNA interactome aids in defining the unique transcriptomic signatures of cancer cell phenotypes.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

Characterization of microRNAs and Target Genes in Musa acuminata subsp. burmannicoides, var. Calcutta 4 during Interaction with Pseudocercospora musae

Endogenous microRNAs (miRNAs) are small non-coding RNAs that perform post-transcriptional regulatory roles across diverse cellular processes, including defence responses to biotic stresses. Pseudocercospora musae, the causal agent of Sigatoka leaf spot disease in banana (Musa spp.), is an important fungal pathogen of the plant. Illumina HiSeq 2500 sequencing of small RNA libraries derived from leaf material in Musa acuminata subsp. burmannicoides, var. Calcutta 4 (resistant) after inoculation with fungal conidiospores and equivalent non-inoculated controls revealed 202 conserved miRNAs from 30 miR-families together with 24 predicted novel miRNAs. Conserved members included those from families miRNA156, miRNA166, miRNA171, miRNA396, miRNA167, miRNA172, miRNA160, miRNA164, miRNA168, miRNA159, miRNA169, miRNA393, miRNA535, miRNA482, miRNA2118, and miRNA397, all known to be involved in plant immune responses. Gene ontology (GO) analysis of gene targets indicated molecular activity terms related to defence responses that included nucleotide binding, oxidoreductase activity, and protein kinase activity. Biological process terms associated with defence included response to hormone and response to oxidative stress. DNA binding and transcription factor activity also indicated the involvement of miRNA target genes in the regulation of gene expression during defence responses. sRNA-seq expression data for miRNAs and RNAseq data for target genes were validated using stem-loop quantitative real-time PCR (qRT-PCR). For the 11 conserved miRNAs selected based on family abundance and known involvement in plant defence responses, the data revealed a frequent negative correlation of expression between miRNAs and target host genes. This examination provides novel information on miRNA-mediated host defence responses, applicable in genetic engineering for the control of Sigatoka leaf spot disease.

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.

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.

Non-coding RNAs in human health and disease: potential function as biomarkers and therapeutic targets

Human diseases have been a critical threat from the beginning of human history. Knowing the origin, course of action and treatment of any disease state is essential. A microscopic approach to the molecular field is a more coherent and accurate way to explore the mechanism, progression, and therapy with the introduction and evolution of technology than a macroscopic approach. Non-coding RNAs (ncRNAs) play increasingly important roles in detecting, developing, and treating all abnormalities related to physiology, pathology, genetics, epigenetics, cancer, and developmental diseases. Noncoding RNAs are becoming increasingly crucial as powerful, multipurpose regulators of all biological processes. Parallel to this, a rising amount of scientific information has revealed links between abnormal noncoding RNA expression and human disorders. Numerous non-coding transcripts with unknown functions have been found in addition to advancements in RNA-sequencing methods. Non-coding linear RNAs come in a variety of forms, including circular RNAs with a continuous closed loop (circRNA), long non-coding RNAs (lncRNA), and microRNAs (miRNA). This comprises specific information on their biogenesis, mode of action, physiological function, and significance concerning disease (such as cancer or cardiovascular diseases and others). This study review focuses on non-coding RNA as specific biomarkers and novel therapeutic targets.

RNA silencing proteins and small RNAs in oomycete plant pathogens and biocontrol agents

Introduction

Oomycetes cause several damaging diseases of plants and animals, and some species also act as biocontrol agents on insects, fungi, and other oomycetes. RNA silencing is increasingly being shown to play a role in the pathogenicity of Phytophthora species, either through trans-boundary movement of small RNAs (sRNAs) or through expression regulation of infection promoting effectors.

Methods

To gain a wider understanding of RNA silencing in oomycete species with more diverse hosts, we mined genome assemblies for Dicer-like (DCL), Argonaute (AGO), and RNA dependent RNA polymerase (RDRP) proteins from Phytophthora plurivora, Ph. cactorum, Ph. colocasiae, Pythium oligandrum, Py. periplocum, and Lagenidium giganteum. Moreover, we sequenced small RNAs from the mycelium stage in each of these species.

Results and discussion

Each of the species possessed a single DCL protein, but they differed in the number and sequence of AGOs and RDRPs. SRNAs of 21nt, 25nt, and 26nt were prevalent in all oomycetes analyzed, but the relative abundance and 5' base preference of these classes differed markedly between genera. Most sRNAs mapped to transposons and other repeats, signifying that the major role for RNA silencing in oomycetes is to limit the expansion of these elements. We also found that sRNAs may act to regulate the expression of duplicated genes. Other sRNAs mapped to several gene families, and this number was higher in Pythium spp., suggesting a role of RNA silencing in regulating gene expression. Genes for most effector classes were the source of sRNAs of variable size, but some gene families showed a preference for specific classes of sRNAs, such as 25/26 nt sRNAs targeting RxLR effector genes in Phytophthora species. Novel miRNA-like RNAs (milRNAs) were discovered in all species, and two were predicted to target transcripts for RxLR effectors in Ph. plurivora and Ph. cactorum, indicating a putative role in regulating infection. Moreover, milRNAs from the biocontrol Pythium species had matches in the predicted transcriptome of Phytophthora infestans and Botrytis cinerea, and L. giganteum milRNAs matched candidate genes in the mosquito Aedes aegypti. This suggests that trans-boundary RNA silencing may have a role in the biocontrol action of these oomycetes.

The pattern of DNA methylation alteration, and its association with the expression changes of non-coding RNAs and mRNAs in Moso bamboo under abiotic stress

Epigenetic changes play an important role in plant growth and development and in stress response. However, DNA methylation pattern and its relationship with the expression changes of non-coding RNAs and mRNAs of Moso bamboo in response to abiotic stress is still largely unknown. In this work, we used whole-genome bisulfite sequencing in combination with whole-transcriptome sequencing to analyze the DNA methylation and transcription patterns of mRNAs and non-coding RNAs in Moso bamboo under abiotic stresses such as cold, heat, ultraviolet (UV) and salinity. We found that CHH methylation in the promoter region was positively correlated with gene expression, while CHG and CHH methylations in the gene body regions were negatively associated with gene expression. Moreover, CG and CHG methylations in the promoter regions were negatively correlated with the transcript abundance of long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs). Similarly, the methylation levels of three contexts in the genic regions were negatively correlated with the transcript abundance of lncRNAs and miRNAs but positively correlated with that of circRNAs. In addition, we suggested that the reduction of 21-nt and 24-nt small interfering RNA (siRNA) expression tended to increase methylation levels in the genic regions. We found that stress-responsive genes such as CRPK1, HSFB2A and CIPK were differentially methylated and expressed. Our results also proposed that DNA methylation may regulate the expression of the transcription factors (TFs) and plant hormone signalling genes such as IAA9, MYC2 and ERF110 in response to abiotic stress. This study firstly reports the abiotic stress-responsive DNA methylation pattern and its involvement of expression of coding RNAs and non-coding RNAs in Moso bamboo. The results expand the knowledge of epigenetic mechanisms in Moso bamboo under abiotic stress and support in-depth deciphering of the function of specific non-coding RNAs in future studies.

Identification and differential analysis of noncoding RNAs in response to drought in Phyllostachys aureosulcata f. spectabilis

The role of noncoding RNAs (ncRNAs) in plant resistance to abiotic stresses is increasingly being discovered. Drought stress is one of the most common stresses that affecting plant growth, and high intensity drought has a significant impact on the normal growth of plants. In this study, a high-throughput sequencing was performed on plant tissue samples of Phyllostachys aureosulcata f. spectabilis C. D. Chu et C. S. Chao by drought treatment for 0, 2, 4 and 6 days. The sequencing results were analysed bioinformatically. We detected 336,946 RNAs among all 12 samples, including 192,098 message RNAs (mRNAs), 142,761 long noncoding RNAs (lncRNAs), 1,670 circular RNAs (circRNAs), and 417 microRNAs (miRNAs). We detected 2,419 differentially expressed (DE) ncRNAs, including 213 DE circRNAs, 2,088 DE lncRNAs and 118 DE miRNAs. Then, we used Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to functionally predict DE ncRNAs. The results showed that most DE ncRNAs are involved in the response to drought stress, mainly in biochemical reactions involved in some metabolites, as well as in organelle activities. In addition, we validated two random circRNAs and demonstrated their circularity. We also found a stable internal reference gene available for Phyllostachys aureosulcata f. spectabilis and validated the accuracy of this experiment by quantitative real-time polymerase chain reaction (qRT-PCR).

Recent Progress and Challenges in the Development of Antisense Therapies for Myotonic Dystrophy Type 1

Myotonic dystrophy type 1 (DM1) is a dominant genetic disease in which the expansion of long CTG trinucleotides in the 3' UTR of the myotonic dystrophy protein kinase (DMPK) gene results in toxic RNA gain-of-function and gene mis-splicing affecting mainly the muscles, the heart, and the brain. The CUG-expanded transcripts are a suitable target for the development of antisense oligonucleotide (ASO) therapies. Various chemical modifications of the sugar-phosphate backbone have been reported to significantly enhance the affinity of ASOs for RNA and their resistance to nucleases, making it possible to reverse DM1-like symptoms following systemic administration in different transgenic mouse models. However, specific tissue delivery remains to be improved to achieve significant clinical outcomes in humans. Several strategies, including ASO conjugation to cell-penetrating peptides, fatty acids, or monoclonal antibodies, have recently been shown to improve potency in muscle and cardiac tissues in mice. Moreover, intrathecal administration of ASOs may be an advantageous complementary administration route to bypass the blood-brain barrier and correct defects of the central nervous system in DM1. This review describes the evolution of the chemical design of antisense oligonucleotides targeting CUG-expanded mRNAs and how recent advances in the field may be game-changing by forwarding laboratory findings into clinical research and treatments for DM1 and other microsatellite diseases.

Multi-omics reveals the key and specific miRNA-mRNA modules underlying salt tolerance in wild emmer wheat (Triticum dicoccoides L.)

Background

Salt stress is one of the most destructive environmental factors limiting crop growth and development. MicroRNAs (miRNAs) are a class of conserved endogenous small non-coding RNAs, playing the crucial role in regulating salt response and tolerance in plants. However, the miRNAs in wild emmer wheat, especially the key and specific salt-responsive miRNAs are not well studied.

Results

Here, we performed small RNA, transcriptome, and degradome sequencing of both of salt-tolerance (ST) and salt-sensitive (SS) wild emmer genotypes to identify the miRNA-mRNA modules associating with salt tolerance. Totally, 775 miRNAs, including 361 conserved known miRNAs and 414 novel miRNAs were detected. Differential expression analysis identified 93 salt-responsive miRNAs under salt stress. Combined with RNA-seq and degradome sequencing analysis, 224 miRNA-mRNA modules displayed the complete opposite expression trends between ST and SS genotypes, most of which functionally enriched into ROS homeostasis maintaining, osmotic pressure modulating, and root growth and development. Finally, the qRT-PCR and a large-scale yeast functional screening were also performed to initially validate the expression pattern and function of candidate genes.

Conclusions

This study reported the key and specific miRNA-mRNA modules associated with salt tolerance in wild emmer, which lay the foundation for improving the salt tolerance in cultivated emmer and bread wheat through miRNA engineering approach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08945-3.

Integrated analysis of competitive endogenous ribose nucleic acids (ceRNAs)-related regulatory networks in invasive and non-invasive non-functioning pituitary adenomas (NFPAs)

Background

This study aims to identify the differentially expressed (DE) non-coding ribose nucleic acids (ncRNAs), messenger RNA (mRNA) expression profiles, and competitive endogenous RNA (ceRNA)-related regulatory networks in invasive and non-invasive nonfunctioning pituitary adenomas (NFPAs).

Methods

A full-transcriptome sequencing of invasive and non-invasive NFPAs is carried out to evaluate the expression profiles of circular RNAs (circRNAs), long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and mRNA expression profiles.

Results

The screening criteria resulted in 118 DEcircRNAs (88 up-regulated and 30 down-regulated), 105 DElncRNAs (68 up-regulated and 37 down-regulated), 43 DEmiRNAs (22 up-regulated and 21 down-regulated), and 268 DEmRNAs (194 up-regulated and 74 down-regulated). Accordingly, a ceRNA regulatory network related to invasive NFPA is constructed. Further, the Gene Ontology and Kyoto Encylopedia of Genes and Genomes analyses showed that circRNAs and lncRNAs in the network are related to chromatin remodeling, participating in the Janus kinase/signal transducer and activator of transcription (JAK-STAT) and calcium signaling pathways. Hsa-miR-1248 showed exceptional connectivity in the ceRNA regulatory network, which could be closely related to the invasiveness of NFPAs.

Conclusions

Together, these findings clarified the regulatory mechanisms of invasive and non-invasive NFPAs, providing innovative research avenues and therapeutic targets for invasive NFPAs.

4D genetic networks reveal the genetic basis of metabolites and seed oil-related traits in 398 soybean RILs

Background

The yield and quality of soybean oil are determined by seed oil-related traits, and metabolites/lipids act as bridges between genes and traits. Although there are many studies on the mode of inheritance of metabolites or traits, studies on multi-dimensional genetic network (MDGN) are limited.

Results

In this study, six seed oil-related traits, 59 metabolites, and 107 lipids in 398 recombinant inbred lines, along with their candidate genes and miRNAs, were used to construct an MDGN in soybean. Around 175 quantitative trait loci (QTLs), 36 QTL-by-environment interactions, and 302 metabolic QTL clusters, 70 and 181 candidate genes, including 46 and 70 known homologs, were previously reported to be associated with the traits and metabolites, respectively. Gene regulatory networks were constructed using co-expression, protein–protein interaction, and transcription factor binding site and miRNA target predictions between candidate genes and 26 key miRNAs. Using modern statistical methods, 463 metabolite–lipid, 62 trait–metabolite, and 89 trait–lipid associations were found to be significant. Integrating these associations into the above networks, an MDGN was constructed, and 128 sub-networks were extracted. Among these sub-networks, the gene–trait or gene–metabolite relationships in 38 sub-networks were in agreement with previous studies, e.g., oleic acid (trait)–GmSEI–GmDGAT1a–triacylglycerol (16:0/18:2/18:3), gene and metabolite in each of 64 sub-networks were predicted to be in the same pathway, e.g., oleic acid (trait)–GmPHS–d-glucose, and others were new, e.g., triacylglycerol (16:0/18:1/18:2)–GmbZIP123–GmHD-ZIPIII-10–miR166s–oil content.

Conclusions

This study showed the advantages of MGDN in dissecting the genetic relationships between complex traits and metabolites. Using sub-networks in MGDN, 3D genetic sub-networks including pyruvate/threonine/citric acid revealed genetic relationships between carbohydrates, oil, and protein content, and 4D genetic sub-networks including PLDs revealed the relationships between oil-related traits and phospholipid metabolism likely influenced by the environment. This study will be helpful in soybean quality improvement and molecular biological research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02191-1.

Identification of miRNAs and Target Genes at Key Stages of Sexual Differentiation in Androdioecious Osmanthus fragrans

Androdioecy is the crucial transition state in the evolutionary direction of hermaphroditism to dioecy, however, the molecular mechanisms underlying the formation of this sex system remain unclear. While popular in China for its ornamental and cultural value, Osmanthus fragrans has an extremely rare androdioecy breeding system, meaning that there are both male and hermaphroditic plants in a population. To unravel the mechanisms underlying the formation of androdioecy, we performed small RNA sequencing studies on male and hermaphroditic O. fragrans. A total of 334 miRNAs were identified, of which 59 were differentially expressed. Functional categorization revealed that the target genes of differentially expressed miRNAs were mainly involved in the biological processes of reproductive development and the hormone signal transduction pathway. We speculated that the miRNA160, miRNA167, miRNA393 and miRNA396 families may influence the sex differentiation in O. fragrans. Overall, our study is the first exploration of miRNAs in the growth and development process of O. fragrans, and is also the first study of androdioecious plants from the miRNA sequencing perspective. The analysis of miRNAs and target genes that may be involved in the sex differentiation process lay a foundation for the ultimate discovery of the androdioecious molecular mechanism in O. fragrans.

A Combination of a Genome-Wide Association Study and a Transcriptome Analysis Reveals circRNAs as New Regulators Involved in the Response to Salt Stress in Maize

Salinization seriously threatens the normal growth of maize, especially at the seedling stage. Recent studies have demonstrated that circular RNAs (circRNAs) play vital roles in the regulation of plant stress resistance. Here, we performed a genome-wide association study (GWAS) on the survival rate of 300 maize accessions under a salt stress treatment. A total of 5 trait-associated SNPs and 86 candidate genes were obtained by the GWAS. We performed RNA sequencing for 28 transcriptome libraries derived from 2 maize lines with contrasting salt tolerance under normal and salt treatment conditions. A total of 1217 highly expressed circRNAs were identified, of which 371 were responsive to a salt treatment. Using PCR and Sanger sequencing, we verified the reliability of these differentially expressed circRNAs. An integration of the GWAS and RNA-Seq analyses uncovered two differentially expressed hub genes (Zm00001eb013650 and Zm00001eb198930), which were regulated by four circRNAs. Based on these results, we constructed a regulation model of circRNA/miRNA/mRNA that mediated salt stress tolerance in maize. By conducting hub gene-based association analyses, we detected a favorable haplotype in Zm00001eb198930, which was responsible for high salt tolerance. These results help to clarify the regulatory relationship between circRNAs and their target genes as well as to develop salt-tolerant lines for maize breeding.

Transcriptome analysis revealed the expression levels of genes related to abscisic acid and auxin biosynthesis in grapevine (Vitis vinifera L.) under root restriction

The root system is essential for the stable growth of plants. Roots help anchor plants in the soil and play a crucial role in water uptake, mineral nutrient absorption and endogenous phytohormone formation. Root-restriction (RR) cultivation, a powerful technique, confines plant roots to a specific soil space. In the present study, roots of one-year-old "Muscat Hamburg" grapevine under RR and control (nR) treatments harvested at 70 and 125 days after planting were used for transcriptome sequencing, and in total, 2031 (nR7 vs. nR12), 1445 (RR7 vs. RR12), 1532 (nR7 vs. RR7), and 2799 (nR12 vs. RR12) differentially expressed genes (DEGs) were identified. Gene Ontology (GO) enrichment analysis demonstrated that there were several genes involved in the response to different phytohormones, including abscisic acid (ABA), auxin (IAA), ethylene (ETH), gibberellins (GAs), and cytokinins (CTKs). Among them, multiple genes, such as PIN2 and ERF113, are involved in regulating vital plant movements by various phytohormone pathways. Moreover, following RR cultivation, DEGs were enriched in the biological processes of plant-type secondary cell wall biosynthesis, the defense response, programmed cell death involved in cell development, and the oxalate metabolic process. Furthermore, through a combined analysis of the transcriptome and previously published microRNA (miRNA) sequencing results, we found that multiple differentially expressed miRNAs (DEMs) and DEG combinations in different comparison groups exhibited opposite trends, indicating that the expression levels of miRNAs and their target genes were negatively correlated. Furthermore, RR treatment indeed significantly increased the ABA content at 125 days after planting and significantly decreased the IAA content at 70 days after planting. Under RR cultivation, most ABA biosynthesis-related genes were upregulated, while most IAA biosynthesis-related genes were downregulated. These findings lay a solid foundation for further establishing the network through which miRNAs regulate grapevine root development through target genes and for further exploring the molecular mechanism through which endogenous ABA and IAA regulate root architecture development in grapevine.

Integrative analysis of transcriptome and miRNAome reveals molecular mechanisms regulating pericarp thickness in sweet corn during kernel development

Pericarp thickness affects the edible quality of sweet corn (Zea mays L. saccharata Sturt.). Therefore, breeding varieties with a thin pericarp is important for the quality breeding of sweet corn. However, the molecular mechanisms underlying the pericarp development remain largely unclear. We performed an integrative analysis of mRNA and miRNA sequencing to elucidate the genetic mechanism regulating pericarp thickness during kernel development (at 15 days, 19 days, and 23 days after pollination) of two sweet corn inbred lines with different pericarp thicknesses (M03, with a thinner pericarp and M08, with a thicker pericarp). A total of 2,443 and 1,409 differentially expressed genes (DEGs) were identified in M03 and M08, respectively. Our results indicate that phytohormone-mediated programmed cell death (PCD) may play a critical role in determining pericarp thickness in sweet corn. Auxin (AUX), gibberellin (GA), and brassinosteroid (BR) signal transduction may indirectly mediate PCD to regulate pericarp thickness in M03 (the thin pericarp variety). In contrast, abscisic acid (ABA), cytokinin (CK), and ethylene (ETH) signaling may be the key regulators of pericarp PCD in M08 (the thick pericarp variety). Furthermore, 110 differentially expressed microRNAs (DEMIs) and 478 differentially expressed target genes were identified. miRNA164-, miRNA167-, and miRNA156-mediated miRNA-mRNA pairs may participate in regulating pericarp thickness. The expression results of DEGs were validated by quantitative real-time PCR. These findings provide insights into the molecular mechanisms regulating pericarp thickness and propose the objective of breeding sweet corn varieties with a thin pericarp.

Key regulatory pathways, microRNAs, and target genes participate in adventitious root formation of Acer rubrum L

Red maple (Acer rubrum L.) is a type of colorful ornamental tree with great economic value. Because this tree is difficult to root under natural conditions and the seedling survival rate is low, vegetative propagation methods are often used. Because the formation of adventitious roots (ARs) is essential for the asexual propagation of A. rubrum, it is necessary to investigate the molecular regulatory mechanisms of AR formation in A. rubrum. To address this knowledge gap, we sequenced the transcriptome and small RNAs (sRNAs) of the A. rubrum variety ‘Autumn Fantasy’ using high-throughput sequencing and explored changes in gene and microRNA (miRNA) expression in response to exogenous auxin treatment. We identified 82,468 differentially expressed genes (DEGs) between the treated and untreated ARs, as well as 48 known and 95 novel miRNAs. We also identified 172 target genes of the known miRNAs using degradome sequencing. Two key regulatory pathways (ubiquitin mediated proteolysis and plant hormone signal transduction), Ar-miR160a and the target gene auxin response factor 10 (ArARF10) were selected based on KEGG pathway and cluster analyses. We further investigated the expression patterns and regulatory roles of ArARF10 through subcellular localization, transcriptional activation, plant transformation, qRT-PCR analysis, and GUS staining. Experiments overexpressing ArARF10 and Ar-miR160a, indicated that ArARF10 promoted AR formation, while Ar-miR160a inhibited AR formation. Transcription factors (TFs) and miRNAs related to auxin regulation that promote AR formation in A. rubrum were identified. Differential expression patterns indicated the Ar-miR160a-ArARF10 interaction might play a significant role in the regulation of AR formation in A. rubrum. Our study provided new insights into mechanisms underlying the regulation of AR formation in A. rubrum.

Integrated small RNA profiling and degradome analysis of Anthurium andraeanum cultivars with different-colored spathes

MicroRNAs (miRNAs) are known to play vital roles in coloration of leaves, flowers, and fruits in plants. However, their functions in spathe coloration are poorly known. Anthurium andraeanum is a popular ornamental plant with various spathe colors. In this study, small RNA and degradome libraries from three A. andraeanum cultivars with different-colored spathes were constructed and sequenced. Illumina sequencing resulted in 94 conserved miRNAs, and 34 novel miRNAs in total were then identified based on precursor sequences and hairpin structures. Differential expression analysis showed that 52, 51, and 49 miRNAs were differentially expressed in comparisons of orange- versus white-colored spathe, purple- versus white-colored spathe, and purple- versus orange-colored spathe, respectively. The expression patterns of miRNAs and their corresponding targets involved in spathe coloration were further analyzed, and displayed that miR156b and miR529 were highly abundant in the spathes with higher anthocyanin content. These two miRNAs co-targeted a gene encoding SPL17, which may function as a negative regulator in anthocyanin accumulation. In addition, miR408 was also abundantly expressed in purple- and orange-colored spathes, and its typical targets were also identified. This comprehensive integrated analysis provides insight into the miRNA-mediated genetic regulation in spathe coloration of A. andraeanum.

Comparative Small RNA and Degradome Sequencing Provides Insights into Antagonistic Interactions in the Biocontrol Fungus Clonostachys rosea

Necrotrophic mycoparasitism is an intricate process involving recognition, physical mycelial contact, and killing of host fungi (mycohosts). During such interactions, mycoparasites undergo a complex developmental process involving massive regulatory changes of gene expression to produce a range of chemical compounds and proteins that contribute to the parasitism of the mycohosts. Small RNAs (sRNAs) are vital components of posttranscriptional gene regulation, although their role in gene expression regulation during mycoparasitisms remain understudied. Here, we investigated the role of sRNA-mediated gene regulation in mycoparasitism by performing sRNA and degradome tag sequencing of the mycoparasitic fungus Clonostachys rosea interacting with the plant-pathogenic mycohosts Botrytis cinerea and Fusarium graminearum at two time points. The majority of differentially expressed sRNAs were downregulated during the interactions with the mycohosts compared to a C. rosea self-interaction control, thus allowing desuppression (upregulation) of mycohost-responsive genes. Degradome analysis showed a positive correlation between high degradome counts and antisense sRNA mapping and led to the identification of 201 sRNA-mediated potential gene targets for 282 differentially expressed sRNAs. Analysis of sRNA potential gene targets revealed that the regulation of genes coding for membrane proteins was a common response against both mycohosts. The regulation of genes involved in oxidative stress tolerance and cellular metabolic and biosynthetic processes was exclusive against F. graminearum, highlighting common and mycohost-specific gene regulation of C. rosea. By combining these results with transcriptome data collected during a previous study, we expand the understanding of the role of sRNA in regulating interspecific fungal interactions and mycoparasitism.

IMPORTANCE Small RNAs (sRNAs) are emerging as key players in pathogenic and mutualistic fungus-plant interactions; however, their role in fungus-fungus interactions remains elusive. In this study, we employed the necrotrophic mycoparasite Clonostachys rosea and the plant-pathogenic mycohosts Botrytis cinerea and Fusarium graminearum and investigated the sRNA-mediated gene regulation in mycoparasitic interactions. The combined approach of sRNA and degradome tag sequencing identified 201 sRNA-mediated putative gene targets for 282 differentially expressed sRNAs, highlighting the role of sRNA-mediated regulation of mycoparasitism in C. rosea. We also identified 36 known and 13 novel microRNAs (miRNAs) and their potential gene targets at the endogenous level and at a cross-species level in B. cinerea and F. graminearum, indicating a role of cross-species RNA interference (RNAi) in mycoparasitism, representing a novel mechanism in biocontrol interactions. Furthermore, we showed that C. rosea adapts its transcriptional response, and thereby its interaction mechanisms, based on the interaction stages and identity of the mycohost.

A Multi-Level Iterative Bi-Clustering Method for Discovering miRNA Co-regulation Network of Abiotic Stress Tolerance in Soybeans

Although growing evidence shows that microRNA (miRNA) regulates plant growth and development, miRNA regulatory networks in plants are not well understood. Current experimental studies cannot characterize miRNA regulatory networks on a large scale. This information gap provides an excellent opportunity to employ computational methods for global analysis and generate valuable models and hypotheses. To address this opportunity, we collected miRNA-target interactions (MTIs) and used MTIs from Arabidopsis thaliana and Medicago truncatula to predict homologous MTIs in soybeans, resulting in 80,235 soybean MTIs in total. A multi-level iterative bi-clustering method was developed to identify 483 soybean miRNA-target regulatory modules (MTRMs). Furthermore, we collected soybean miRNA expression data and corresponding gene expression data in response to abiotic stresses. By clustering these data, 37 MTRMs related to abiotic stresses were identified, including stress-specific MTRMs and shared MTRMs. These MTRMs have gene ontology (GO) enrichment in resistance response, iron transport, positive growth regulation, etc. Our study predicts soybean MTRMs and miRNA-GO networks under different stresses, and provides miRNA targeting hypotheses for experimental analyses. The method can be applied to other biological processes and other plants to elucidate miRNA co-regulation mechanisms.

A fungal microRNA-like RNA subverts host immunity and facilitates pathogen infection by silencing two host receptor-like kinase genes

Small RNAs (sRNAs) play important roles in various biological processes by silencing their corresponding target genes in most eukaryotes. However, cross-kingdom regulation mediated by fungal microRNA-like RNAs (milRNAs) in plant-pathogen interactions is still largely unknown. Using molecular, genetic, histological, and biochemical approaches, we found that the apple tree Valsa canker pathogen Valsa mali milRNA Vm-milR1 could suppress the host immunity by silencing two host receptor-like kinase genes, MdRLKT1 and MdRLKT2. Vm-milR1 was highly induced during V. mali infection. Deletion of Vm-milR1 precursor abolished the generation of Vm-milR1 and reduced the virulence of V. mali. Inoculation of Vm-milR1 deletion mutants induced the host defence responses, including reactive oxygen species (ROS) accumulation, callose deposition, and high expression of defence-related genes. Furthermore, Vm-milR1 was confirmed to be able to suppress the expression of MdRLKT1 and MdRLKT2 in a sequence-specific manner. Moreover, overexpression of either MdRLKT1 or MdRLKT2 enhanced apple resistance to V. mali by activating the host defence responses. Furthermore, knockdown of MdRLKT1 or MdRLKT2 compromised the host resistance to V. mali. Our study revealed that V. mali was equipped with Vm-milR1 as an sRNA effector to silence host receptor-like kinase genes, suppress the host defence responses, and facilitate pathogen infection.

Identification and Comparative Analysis of Genes and MicroRNAs Involved in the Floral Transition of the Xinjiang Early-Flowering Walnut (Juglans regia L.)

For tree crops, shortening the juvenile phase is a vital strategy to advance fruit bearing and enhance the breeding efficiency. Walnut (Juglans regia L.) seedlings usually take at least eight to 10 years to flower, but early-flowering (EF) types can flower one or two years after planting. In this study, RNA sequencing (RNA-Seq) and microRNA sequencing (miRNA-Seq) were used for a transcriptome-wide analysis of gene and miRNA expression in hybrids of the Xinjiang EF walnut variety ‘Xinwen 81’ and later-flowering (LF) walnut. Based on a high-quality chromosome-scale reference genome, a total of 3009 differentially expressed genes (DEGs) were identified, of which 933 were upregulated (accounting for 31%) and 2076 were downregulated (accounting for 69%). DEGs were functionally annotated, and the flowering-related genes FLOWERING LOCUS T (FT), SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), and LEAFY (LFY) showed remarkable upregulation in EF compared with in the LF walnut. In addition, miRNAs associated with floral transition were screened as candidates for flowering time regulation in the walnut. This work provides new insights into walnut floral transition, which may ultimately contribute to genetic improvement of the walnut.

Small Drugs, Huge Impact: The Extraordinary Impact of Antisense Oligonucleotides in Research and Drug Development

Antisense oligonucleotides (ASOs) are an increasingly represented class of drugs. These small sequences of nucleotides are designed to precisely target other oligonucleotides, usually RNA species, and are modified to protect them from degradation by nucleases. Their specificity is due to their sequence, so it is possible to target any RNA sequence that is already known. These molecules are very versatile and adaptable given that their sequence and chemistry can be custom manufactured. Based on the chemistry being used, their activity may significantly change and their effects on cell function and phenotypes can differ dramatically. While some will cause the target RNA to decay, others will only bind to the target and act as a steric blocker. Their incredible versatility is the key to manipulating several aspects of nucleic acid function as well as their process, and alter the transcriptome profile of a specific cell type or tissue. For example, they can be used to modify splicing or mask specific sites on a target. The entire design rather than just the sequence is essential to ensuring the specificity of the ASO to its target. Thus, it is vitally important to ensure that the complete process of drug design and testing is taken into account. ASOs' adaptability is a considerable advantage, and over the past decades has allowed multiple new drugs to be approved. This, in turn, has had a significant and positive impact on patient lives. Given current challenges presented by the COVID-19 pandemic, it is necessary to find new therapeutic strategies that would complement the vaccination efforts being used across the globe. ASOs may be a very powerful tool that can be used to target the virus RNA and provide a therapeutic paradigm. The proof of the efficacy of ASOs as an anti-viral agent is long-standing, yet no molecule currently has FDA approval. The emergence and widespread use of RNA vaccines during this health crisis might provide an ideal opportunity to develop the first anti-viral ASOs on the market. In this review, we describe the story of ASOs, the different characteristics of their chemistry, and how their characteristics translate into research and as a clinical tool.

Genome-Wide Identification and Comparative Profiling of MicroRNAs Reveal Flavonoid Biosynthesis in Two Contrasting Flower Color Cultivars of Tree Peony

MicroRNA (miRNA)-mediated gene regulation is involved in various physiological processes in plants. Flower color is one of the vital ornamental traits of tree peony (Paeonia suffruticosa Andr.). However, the yellow-flowered tree peony cultivars are particularly rare. To elucidate the miRNA-mediated gene regulatory mechanism underlying yellow pigmentation in tree peony, we combined pigment assessment, miRNA identification, expression analysis, and gene functional verification in two contrasting flower color cultivars "High Noon" and "Roufurong." Flavones/flavonols and anthocyanins were found to be the main contributors to the coloration of "High Noon" and "Roufurong" petals, respectively. Subsequently, miRNA analysis based on available genome data identified 9 differentially expressed miRNAs and 12 relevant target genes implicated in flavonoid biosynthesis. Their dynamic expression patterns determined the key role of mdm-miR156b-PsSPL2 module in yellow pigmentation of tree peony flowers. The sequence analysis and subcellular localization validated that PsSPL2 might function as a nuclear-localized transcription factor. Overexpression of PsSPL2 in tobacco resulted in a decrease of anthocyanin content and down-regulation of NtF3'H and NtDFR transcripts. PsSPL2-silenced petals exhibited lighter yellow color, and the contents of THC, Ap, and Ch decreased significantly. Meanwhile, expression levels of PsCHS, PsCHI, and PsF3H were significantly decreased in the petals with PsSPL2 silencing, while those of PsF3'H and PsDFR were remarkably increased. This study offers a novel insight into yellow pigmentation-related miRNA regulation network in tree peony, and further provides the valuable information on physiological changes during yellow coloring process of tree peony.

Comparative Transcriptome Sequencing of Taro Corm Development With a Focus on the Starch and Sucrose Metabolism Pathway

Taro (Colocasia esculenta) is an important tuber crop and staple food. Taro corms have higher nutritional value and starch contents as compared to most of the other root/tuber crops. However, the growth and development of the taro rhizome have not been critically examined in terms of transcriptomic signatures in general or specific to carbohydrates (starch and sucrose) accumulation. In current study, we have conducted a comprehensive survey of transcripts in taro corms aged 1, 2, 3, 4, 5, and 8 months. In this context, we have employed a whole transcriptome sequencing approach for identification of mRNAs, CircRNAs, and miRNAs in corms and performed functional enrichment analysis of the screened differentially expressed RNAs. A total of 11,203 mRNAs, 245 CircRNAs, and 299 miRNAs were obtained from six developmental stages. The mRNAs included 139 DEGs associated with 24 important enzymes of starch and sucrose metabolism. The expression of genes encoding key enzymes of starch and sucrose metabolism pathway (GBSS, AGPase, UGPase, SP, SSS, βFRUCT and SuSy) demonstrated significant variations at the stage of 4 months (S4). A total of 191 CircRNAs were differentially expressed between the studied comparisons of growth stages and 99 of these were associated with those miRNA (or target genes) that were enriched in starch and sucrose metabolism pathway. We also identified 205 miRNAs including 46 miRNAs targeting DEGs enriched in starch and sucrose biosynthesis pathway. The results of current study provide valuable resources for future exploration of the molecular mechanisms involved in the starch properties of Taro.

Tr-milRNA1 Contributes to Lignocellulase Secretion under Heat Stress by Regulating the Lectin-Type Cargo Receptor Gene Trvip36 in Trichoderma guizhouence NJAU 4742

Background: MicroRNA plays an important role in multifarious biological processes by regulating their corresponding target genes. However, the biological function and regulatory mechanism of fungal microRNA-like RNAs (milRNAs) remain poorly understood. Methods: In this study, combined with deep sequencing and bioinformatics analysis, milRNAs and their targets from Trichoderma guizhouence NJAU 4742 were isolated and identified under solid-state fermentation (SSF) by using rice straw as the sole carbon source at 28 °C and 37 °C, respectively. Results: A critical milRNA, TGA1_S04_31828 (Tr-milRNA1), was highly expressed under heat stress (37 °C) and adaptively regulated lignocellulase secretion. Overexpression of Tr-milRNA1 (OE-Tr-milRNA1) did not affect vegetative growth, but significantly increased lignocellulose utilization under heat stress. Based on the bioinformatics analysis and qPCR validation, a target of Tr-milRNA1 was identified as Trvip36, a lectin-type cargo receptor. The expression of Tr-milRNA1 and Trvip36 showed a divergent trend under SSF when the temperature was increased from 28 °C to 37 °C. In addition, the expression of Trvip36 was suppressed significantly in Tr-milRNA1 overexpression strain (OE-Tr-milRNA1). Compared with the wild type, deletion of Trvip36 (ΔTrvip36) significantly improved the secretion of lignocellulases by reducing the retention of lignocellulases in the ER under heat stress. Conclusions: Tr-milRNA1 from NJAU 4742 improved lignocellulose utilization under heat stress by regulating the expression of the corresponding target gene Trvip36. These findings might open avenues for exploring the mechanism of lignocellulase secretion in filamentous fungi.

Small RNA sequencing reveals various microRNAs involved in piperine biosynthesis in black pepper (Piper nigrum L.)

Background

Black pepper (Piper nigrum L.), an important and long-cultivated spice crop, is native to South India and grown in the tropics. Piperine is the main pungent and bioactive alkaloid in the berries of black pepper, but the molecular mechanism for piperine biosynthesis has not been determined. MicroRNAs (miRNAs), which are classical endogenous noncoding small RNAs, play important roles in regulating secondary metabolism in many species, but less is known regarding black pepper or piperine biosynthesis.

Results

To dissect the functions of miRNAs in secondary metabolism especially in piperine biosynthesis, 110 known miRNAs, 18 novel miRNAs and 1007 individual targets were identified from different tissues of black pepper by small RNA sequencing. qRT-PCR and 5′-RLM-RACE experiments were conducted to validate the reliability of the sequencing data and predicted targets. We found 3 miRNAs along with their targets including miR166-4CL, miR396-PER and miR397-CCR modules that are involved in piperine biosynthesis.

Conclusion

MiRNA regulation of secondary metabolism is a common phenomenon in plants. Our study revealed new miRNAs that regulate piperine biosynthesis, which are special alkaloids in the piper genus, and they might be useful for future piperine genetic improvement of black pepper.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08154-4.

The Effect of Low Temperature Stress on the Leaves and MicroRNA Expression of Potato Seedlings

In recent years, with the wanton destruction of the ecological environment by humans and the frequent occurrence of extreme bad weather, many places that should have been warm and blooming in spring have instead experienced the phenomenon of the “April blizzard,” which has seriously affected China's crops, especially spring potato production in most areas. Potato cultivars, especially potato seedlings, are sensitive to frost, and low temperature frost has become one of the most important abiotic stresses affecting potato production. Potato cold tolerance is regulated by a complex gene network. Although some low temperature resistant microRNAs have been identified, little is known about the role of miRNAs in response to low temperature stress in potato. Therefore, the objective of this study is to clarify the influence of low temperature stress on the miRNA expression of potato by comparing the expression differences of miRNA in potato which was treated with different low temperatures. For the study, 307 known miRNAs belonging to 73 small RNA families and 211 novel miRNAs were obtained. When the temperature decreased, the number of both known and novel miRNA decreased, and the minimum temperature was −2°C. Most of the miRNAs respond to low temperature, drought, and disease stress; some conserved miRNAs were first found to respond to low temperature stress in potato, such as stu-miR530, stu-miR156d, and stu-miR167b. The Gene Ontology, Kyoto Encyclopedia of Genes, and Genomes pathway enrichment analysis of 442 different expression miRNAs target genes indicated that there existed diversified low temperature responsive pathways, but Abscisic Acid was found likely to play a central coordinating role in response to low temperature stress in many metabolism pathways. Quantitative real-time PCR assays indicated that the related targets were negatively regulated by the tested different expression miRNAs during low temperature stress. The results indicated that miRNAs may play an important coordination role in response to low temperature stress in many metabolic pathways by regulating abscisic acid and gibberellin, which provided insight into the roles of miRNAs during low temperature stress and would be helpful for alleviating low temperature stress and promoting low temperature resistant breeding in potatoes.

Identification of microRNA-like RNAs from Trichoderma asperellum DQ-1 during its interaction with tomato roots using bioinformatic analysis and high-throughput sequencing

MicroRNA-like small RNAs (milRNAs) and their regulatory roles in the interaction between plant and fungus have recently aroused keen interest of plant pathologists. Trichoderma spp., one of the widespread biocontrol fungi, can promote plant growth and induce plant disease resistance. To investigate milRNAs potentially involved in the interaction between Trichoderma and tomato roots, a small RNA (sRNA) library expressed during the interaction of T. asperellum DQ-1 and tomato roots was constructed and sequenced using the Illumina HiSeqTM 2500 sequencing platform. From 13,464,142 sRNA reads, we identified 21 milRNA candidates that were similar to other known microRNAs in the miRBase database and 22 novel milRNA candidates that possessed a stable microRNA precursor hairpin structure. Among them, three milRNA candidates showed different expression level in the interaction according to the result of stem-loop RT-PCR indicating that these milRNAs may play a distinct regulatory role in the interaction between Trichoderma and tomato roots. The potential transboundary milRNAs from T. asperellum and their target genes in tomato were predicted by bioinformatics analysis. The results revealed that several interesting proteins involved in plant growth and development, disease resistance, seed maturation, and osmotic stress signal transduction might be regulated by the transboundary milRNAs. To our knowledge, this is the first report of milRNAs taking part in the process of interaction of T. asperellum and tomato roots and associated with plant promotion and disease resistance. The results might be useful to unravel the mechanism of interaction between Trichoderma and tomato.