Dataset Information

Comprehensive circRNA-microRNA-mRNA network analysis revealed the novel regulatory mechanism of Trichosporon asahii infection

ABSTRACT:

Background

Invasive Trichosporon asahii (T. asahii) infection frequently occurs with a high mortality in immunodeficient hosts, but the pathogenesis of T. asahii infection remains elusive. Circular RNAs (circRNAs) are a type of endogenous noncoding RNA that participate in various disease processes. However, the mechanism of circRNAs in T. asahii infection remains completely unknown.

Methods

RNA sequencing (RNA-seq) was performed to analyze the expression profiles of circRNAs, microRNAs (miRNAs), and mRNAs in THP-1 cells infected with T. asahii or uninfected samples. Some of the RNA-seq results were verified by RT-qPCR. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were used to analyze the differentially expressed mRNAs. A circRNA-miRNA-mRNA network was constructed and verified by dual-luciferase reporter assay and overexpression experiments.

Results

A total of 46 circRNAs, 412 mRNAs and 47 miRNAs were differentially expressed at 12?h after T. asahii infection. GO and KEGG analyses showed that the differentially expressed mRNAs were primarily linked to the leukocyte migration involved in the inflammatory response, the Toll-like receptor signaling pathway, and the TNF signaling pathway. A competing endogenous RNA (ceRNA) network was constructed with 5 differentially expressed circRNAs, 5 differentially expressed miRNAs and 42 differentially expressed mRNAs. Among them, hsa_circ_0065336 was found to indirectly regulate PTPN11 expression by sponging miR-505-3p.

Conclusions

These data revealed a comprehensive circRNA-associated ceRNA network during T. asahii infection, thus providing new insights into the pathogenesis of the T. asahii-host interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40779-021-00311-w.

SUBMITTER: Zhang M

PROVIDER: S-EPMC7941914 | biostudies-literature | 2021 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Trichosporon asahii is an emerging opportunistic pathogen that causes potentially fatal disseminated trichosporonosis. The global prevalence of coronavirus disease 2019 (COVID-19) poses an increasing fungal infection burden caused by T. asahii. Allicin is the main biologically active component with broad-spectrum antimicrobial activity in garlic. In this study, we performed an in-depth analysis of the antifungal characteristics of allicin against T. asahii based on physiological, cytological, and transcriptomic assessments. In vitro, allicin inhibited the growth of T. asahii planktonic cells and biofilm cells significantly. In vivo, allicin improved the mean survival time of mice with systemic trichosporonosis and reduced tissue fungal burden. Electron microscopy observations clearly demonstrated damage to T. asahii cell morphology and ultrastructure caused by allicin. Furthermore, allicin increased intracellular reactive oxygen species (ROS) accumulation, leading to oxidative stress damage in T. asahii cells. Transcriptome analysis showed that allicin treatment disturbed the biosynthesis of cell membrane and cell wall, glucose catabolism, and oxidative stress. The overexpression of multiple antioxidant enzymes and transporters may also place an additional burden on cells, causing them to collapse. Our findings shed new light on the potential of allicin as an alternative treatment strategy for trichosporonosis. IMPORTANCE Systemic infection caused by T. asahii has recently been recognized as an important cause of mortality in hospitalized COVID-19 patients. Invasive trichosporonosis remains a significant challenge for clinicians, due to the limited therapeutic options. The present work suggests that allicin holds great potential as a therapeutic candidate for T. asahii infection. Allicin demonstrated potent in vitro antifungal activity and potential in vivo protective effects. In addition, transcriptome sequencing provided valuable insights into the antifungal effects of allicin.

Project description:Viral noncoding RNAs have acquired increasing prominence as important regulators of infection and mediators of pathogenesis. Circular RNAs (circRNAs) generated by backsplicing events have been identified in several oncogenic human DNA viruses. Here, we show that Merkel cell polyomavirus (MCV), the etiologic cause of ?80% of Merkel cell carcinomas (MCCs), also expresses circular RNAs. By RNase R-resistant RNA sequencing, four putative circRNA backsplice junctions (BSJs) were identified from the MCV early region (ER). The most abundantly expressed MCV circRNA, designated circMCV-T, is generated through backsplicing of all of ER exon II to form a 762-nucleotide (nt) circular RNA molecule. Curiously, circMCV-T, as well as two other less abundantly expressed putative MCV circRNAs, overlaps in a complementary fashion with the MCV microRNA (miRNA) locus that encodes MCV-miR-M1. circMCV-T is consistently detected in concert with linear T antigen transcripts throughout infection, suggesting a crucial role for this RNA molecule in the regulatory functions of the early region, known to be vital for viral replication. Knocking out the hairpin structure of MCV-miR-M1 in genomic early region expression constructs and using a new high-efficiency, recombinase-mediated, recircularized MCV molecular clone demonstrates that circMCV-T levels decrease in the presence of MCV-miR-M1, underscoring the interplay between MCV circRNA and miRNA. Furthermore, circMCV-T partially reverses the known inhibitory effect of MCV-miR-M1 on early gene expression. RNase R-resistant RNA sequencing of lytic rat polyomavirus 2 (RatPyV2) identified an analogously located circRNA, stipulating a crucial, conserved regulatory function of this class of RNA molecules in the family of polyomaviruses.IMPORTANCE Covalently closed circular RNAs were recently described in the human DNA tumor viruses Epstein-Barr virus (EBV), Kaposi's sarcoma-associated herpesvirus (KSHV), and human papillomavirus (HPV). Here, we show that MCV, another DNA tumor virus, generates circRNAs from its early regulatory region in concert with T antigen linear transcripts. MCV circMCV-T interacts with another MCV noncoding RNA, miR-M1, to functionally modulate early region transcript expression important for viral replication and long-term episomal persistence. This work describes a dynamic regulatory network integrating circRNA/miRNA/mRNA biomolecules and underscores the intricate functional modulation between several classes of polyomavirus-encoded RNAs in the control of viral replication.

Project description:BackgroundCircular RNA (circRNA) plays an important role in the regulation of gene expression and the occurrence of human diseases. However, studies on the role of circRNA in acute myocardial infarction (AMI) are limited. This study was performed to explore novel circRNA-related regulatory networks in AMI, aiming to better understand the molecular mechanism of circRNAs involvement in AMI and provide basis for further scientific research and clinical decision-making.MethodsThe AMI-related microarray datasets GSE160717 (circRNA), GSE31568 (miRNA), GSE61741 (miRNA), and GSE24519 (mRNA) were obtained from the Gene Expression Omnibus (GEO) database. After differential expression analysis, the regulatory relationships between these DERNAs were identified by online databases circBank, circInteractome, miRDB, miRWalk, Targetscan, and then two circRNA-miRNA-mRNA regulatory networks were constructed. Differentially expressed genes (DEGs) in this network were selected followed by enrichment analysis and protein-protein interaction (PPI) analysis. Hub genes were identified using Cytohubba plug-in of Cytoscape software. Hub genes and hub gene-related miRNAs were used for receiver operating characteristic curve (ROC) analysis to identify potential biomarkers. The relative expression levels of these biomarkers were further assessed by GSE31568 (miRNA) and GSE66360 (mRNA). Finally, on the basis of the above analysis, myocardial hypoxia model was constructed to verify the expression of Hub genes and related circRNAs.ResultsA total of 83 DEcircRNAs, 109 CoDEmiRNAs and 1204 DEGs were significantly differentially expressed in these datasets. The up-regulated circRNAs and down-regulated circRNAs were used to construct a circRNA-miRNA-mRNA regulatory network respectively. These circRNA-related DEGs were mainly enriched in the terms of "FOXO signaling pathway," "T cell receptor signaling pathway," "MAPK signaling pathway," "Insulin resistance," "cAMP signaling pathway," and "mTOR signaling pathway." The top 10 hub genes ATP2B2, KCNA1, GRIN2A, SCN2B, GPM6A, CACNA1E, HDAC2, SRSF1, ANK2, and HNRNPA2B1 were identified from the PPI network. Hub genes GPM6A, SRSF1, ANK2 and hub gene-related circRNAs hsa_circ_0023461, hsa_circ_0004561, hsa_circ_0001147, hsa_circ_0004771, hsa_circ_0061276, and hsa_circ_0045519 were identified as potential biomarkers in AMI.ConclusionIn this study, the potential circRNAs associated with AMI were identified and two circRNA-miRNA-mRNA regulatory networks were constructed. This study explored the mechanism of circRNA involvement in AMI and provided new clues for the selection of new diagnostic markers and therapeutic targets for AMI.

Project description:BackgroundOsteosarcoma is a common malignant primary bone tumor in adolescents and children. Numerous studies have shown that circRNAs were involved in the proliferation and invasion of various tumors. However, the role of circRNAs in osteosarcoma remains unclear. Here, we aimed to explore the regulatory network among circRNA-miRNA-mRNA in osteosarcoma.MethodsThe circRNA (GSE140256), microRNA (GSE28423), and mRNA (GSE99671) expression profiles of osteosarcoma were collected from the Gene Expression Omnibus (GEO) database. Differentially expressed circRNAs, miRNAs and mRNAs were identified. CircRNA-miRNA interactions and miRNA-mRNA interactions were determined by Circular RNA Interactome (CircInteractome) database and microRNA Data Integration Portal (mirDIP) database, respectively. Then, we constructed a regulatory network. Function enrichment analysis of miRNA and mRNA was performed by DIANA-miRPath v3.0 and Metascape database, respectively. mRNAs with significant prognostic value were identified based on expression profiles from The Cancer Genome Atlas (TCGA) database, and we constructed a subnetwork for them. To make the most of the network, we used the CLUE database to predict potential drugs for the treatment of osteosarcoma based on mRNA expression in the network. And we used the STITCH database to analyze and validate the interactions among these drugs and mRNAs, and to further screen for potential drugs.ResultsA total of 9 circRNAs, 19 miRNAs, 67 mRNAs, 54 pairs of circRNA-miRNA interactions and 110 pairs of miRNA-mRNA interactions were identified. A circRNA-miRNA-mRNA network was constructed. Function enrichment analysis indicated that these miRNAs and mRNAs in the network were involved in the process of tumorigenesis and immune response. Among these mRNAs, STC2 and RASGRP2 with significantly prognostic value were identified, and we constructed a subnetwork for them. Based on mRNA expression in the network, three potential drugs, quinacridine, thalidomide and zonisamide, were screened for the treatment of osteosarcoma. Among them, quinacridine and thalidomide have been proved to have anti-tumor effects in previous studies, while zonisamide has not been reported. And a corresponding drug-protein interaction network was constructed.ConclusionOverall, we constructed a circRNA-miRNA-mRNA regulatory network to investigate the possible mechanism in osteosarcoma, and predicted that quinacridine, thalidomide and zonisamide could be potential drugs for the treatment of osteosarcoma.