Dataset Information

Construction and Analysis of a circRNA-Mediated ceRNA Network in Lung Adenocarcinoma.

ABSTRACT:

Background

Circular RNAs (circRNAs), a new class of regulatory noncoding RNAs, are involved in gene regulation and may play a role in cancer development. The aim of this study was to identify circRNAs involved in lung adenocarcinoma (LUAD) using bioinformatics analysis.

Methods

CircRNA (GSE101684, GSE101586), miRNA (GSE135918), and mRNA (GSE130779) microarray datasets were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed circRNAs (DECs), miRNAs (DEMs), and mRNAs (DEGs) in LUAD. Circinteractome and StarBase were used to predict miRNAs and mRNAs, respectively. A circRNA-miRNA-mRNA-ceRNA network was constructed. Patient survival was analyzed using UALCAN, and a sub-network was established. Real-time quantitative PCR (qRT-PCR) was used to verify the expressed of DECs between LUAD tissues and paired adjacent normal tissues.

Results

Hsa_circ_0072088 was identified as a differentially expressed (upregulated) circRNA in the two datasets. Intersection analysis identified hsa-miR-532-3p and hsa-miR-942 as the two miRNAs with the highest potential for binding to hsa_circ_0072088. Differential expression analysis and target gene prediction were performed to build a ceRNA network of hsa_circ_0072088 using Circinteractome/StarBase 3.0. Intersection analysis showed that TMEM52, IL24, POF1B, KIF1A, NHS, LBH, HIST2H2BE, ABCC3, PYCR1, CD79A, IGF2BP3, ANKRD17, GTSE1, MKI67, CLSPN, PLAU, LUC7L, MAGIX, GPATCH4, and ABAT were potential downstream mRNAs. The association between the expression level of 20 DEGs and LUAD patient survival was analyzed using UALCAN and GEPIA, which showed that IGF2BP3, MKI67, CD79A, and ABAT were related to patient survival. Hsa_circ_0072088 was verified upregulated by qRT-PCR.

Conclusion

The circRNA hsa_circ_0072088, the DEMs (hsa-miR-532-3p and hsa-miR-942-5p), and the DEGs (IGF2BP3, MKI67, CD79A, and ABAT) may serve as prognostic markers in LUAD.

SUBMITTER: Wang Z

PROVIDER: S-EPMC8197624 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundA growing body of evidence indicates that long non-coding RNAs (lncRNAs) can act as competitive endogenous RNAs (ceRNAs) to bind to microRNAs (miRNAs), thereby affecting and regulating the expression of target genes. The lncRNA-miRNA-mRNA ceRNA network has been theorized to play an indispensable role in many types of tumors. However, the role of the lncRNA-related ceRNA regulatory network in lung adenocarcinoma (LUAD) remains unclear.MethodsWe downloaded the RNAseq and miRNAseq data of LUAD from The Cancer Genome Atlas (TCGA) data portal and identified differentially expressed lncRNAs (DElncRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) between LUAD and corresponding paracancerous tissues by using the edgeR package of R software. We constructed the lncRNA-miRNA-mRNA ceRNA network by using Cytoscape (version 3.7.2) on the basis of the interaction generated from the miRcode, miRTarBase, miRDB, and TargetScan databases. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed with DAVID 6.8 bioinformatics resources and plotted by using the ggplot2 package in R. The effect of genes on LUAD prognosis was assessed by applying the survival package in R in accordance with the Kaplan-Meier curve.ResultsIn total, 1645 DElncRNAs, 117 DEmiRNAs, and 2729 DEmRNAs were identified in LUAD. The LUAD-specific ceRNA network was composed of 157 nodes and 378 edges (329 DElncRNA-DEmiRNA interactions and 49 DEmiRNA-DEmRNA interactions). GO and KEGG pathway annotations suggested that the LUAD-specific ceRNA network was related to tumor-related molecular functions and pathways. Seven lncRNAs (DISC1-IT1, SYNPR-AS1, H19, LINC00460, LINC00518, DSCR10, and STEAP2-AS1), one miRNA (hsa-mir-31), and 16 mRNAs (ATAD2, OSCAR, KIF23, E2F7, PFKP, MCM4, CEP55, CBX2, CCNE1, CLSPN, CCNB1, CDC25A, EZH2, CHEK1, SLC7A11, and PBK) were revealed to be significantly correlated with overall survival.ConclusionIn this study, we described the potential regulatory mechanism of the progression of LUAD. We proposed a new lncRNA-miRNA-mRNA ceRNA network that could help further explore the molecular mechanisms of LUAD.

Project description:BackgroundAortic dissection (AD) is a rare and lethal disorder with its genetic basis remains largely unknown. Many studies have confirmed that circRNAs play important roles in various physiological and pathological processes. However, the roles of circRNAs in AD are still unclear and need further investigation. The present study aimed to elucidate the underlying molecular mechanisms of circRNAs regulation in AD based on the circRNA-associated competing endogenous RNA (ceRNA) network.MethodsExpression profiles of circRNAs (GSE97745), miRNAs (GSE92427), and mRNAs (GSE52093) were downloaded from Gene Expression Omnibus (GEO) databases, and the differentially expressed RNAs (DERNAs) were subsequently identified by bioinformatics analysis. CircRNA-miRNA-mRNA ceRNA network, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to predict the potential functions of circRNA-associated ceRNA network. RNA was isolated from human arterial blood samples after which qRT-PCR was performed to confirm the DERNAs.ResultsWe identified 14 (5 up-regulated and 9 down-regulated) differentially expressed circRNAs (DEcircRNAs), 17 (8 up-regulated and 9 down-regulated) differentially expressed miRNAs (DEmiRNAs) and 527 (297 up-regulated and 230 down-regulated) differentially expressed mRNAs (DEmRNAs) (adjusted P-value <0.05 and | log2FC | > 1.0). KEGG pathway analysis indicated that DEmRNAs were related to focal adhesion and extracellular matrix receptor interaction signaling pathways. Simultaneously, the present study constructed a ceRNA network based on 1 circRNAs (hsa_circRNA_082317), 1 miRNAs (hsa-miR-149-3p) and 10 mRNAs (MLEC, ENTPD7, SLC16A3, SLC7A8, TBC1D16, PAQR4, MAPK13, PIK3R2, ITGA5, SERPINA1). qRT-PCR demonstrated that hsa_circRNA_082317 and ITGA5 were significantly up-regulated, and hsa-miR-149-3p was dramatically down-regulated in AD (n = 3).ConclusionThis is the first study to demonstrate the circRNA-associated ceRNA network is altered in AD, implying that circRNAs may play important roles in regulating the onset and progression and thus may serve as potential biomarkers for the diagnosis and treatment of AD.

Project description:BackgroundCompetitive Endogenous RNA (ceRNA) may be closely associated with tumor progression. However, studies on ceRNAs and immune cells in LUAD are scarce.MethodThe profiles of gene expression and clinical data of LUAD patients were extracted from the TCGA database. Bioinformatics methods were used to evaluate differentially-expressed genes (DEGs) and to form a ceRNA network. Preliminary verification of clinical specimens was utilized to detect the expressions of key biomarkers at the tissues. Cox and Lasso regressions were used to identify key genes, and prognosis prediction nomograms were formed. The mRNA levels of 9 genes in the risk score model in independent clinical LUAD samples were detected by qRT-PCR. The interconnection between the risk of cancer and immune cells was evaluated using the CIBERSORT algorithm, while the conformation of notable tumor-infiltrating immune cells (TIICs) in the LUAD tissues of the high and low risk groups was assessed using the RNA transcript subgroup in order to identify tissue types. Finally, co-expression study was used to examine the interconnection between the key genes in the ceRNA networks and the immune cells.ResultA ceRNA network of 115 RNAs was established, and nine key genes were identified to construct a Cox proportional-hazard model and create a prognostic nomogram. This risk-assessment model might serve as an independent factor to forecast the prognosis of LUAD, and it was consistent with the preliminary verification of clinical specimens. Survival analysis of clinical samples further validated the potential value of high risk groups in predicting LUAD prognosis. Five immune cells were identified with significant differences in the LUAD tissues of the high and low risk groups. Besides, two pairs of biomarkers associated with the growth of LUAD were found, i.e., E2F7 and macrophage M1 (R = 0.419, p = 1.4e- 08) and DBF4 and macrophage M1 (R = 0.282, p < 2.2 e- 16).ConclusionThis study identified several important ceRNAs, i.e. (E2F7 and BNF4) and TIICs (macrophage M1), which might be related to the development and prognosis of LUAD. The established risk-assessment model might be a potential tool in predicting LUAD of prognosis.

Project description:BackgroundMore and more studies have proven that circular RNAs (circRNAs) play vital roles in cancer development via sponging miRNAs. However, the expression pattern of competing endogenous RNA (ceRNA) in lung adenocarcinoma (LUAD) remains largely unclear. The current study explored functional roles and the regulatory mechanisms of circRNA as ceRNAs in LUAD and their potential impact on LUAD patient prognosis.MethodsIn this study, we systematically screened differential expression circRNAs (DEcircRNAs), miRNAs (DEmiRNAs) and mRNAs (DEGs) associated with LUAD. Then, DEcircRNAs, DEmiRNAs and DEGs were selected to construct a circRNA-miRNA-mRNA prognosis-related regulatory network based on interaction information from the ENCORI database. Subsequently, the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed on the genes in the network to predict the potential underlying mechanisms and functions of circRNAs in LUAD. In addition, Kaplan-Meier survival analysis was performed to evaluate clinical outcomes of LUAD patients, and drug sensitivity analysis was used to screen potential biomarkers for drug treatment of patients with LUAD.ResultsAs a result, 10 circRNAs were aberrantly expressed in LUAD tissues. The ceRNA network was built, which included 3 DEcircRNAs, 6 DEmiRNAs and 157 DEGs. The DEGs in the ceRNA network of hsa_circ_0049271 enriched in biological processes of cell proliferation and the Jak-STAT signaling pathway. We also detected 7 mRNAs in the ceRNA network of hsa_circ_0049271 that were significantly associated with the overall survival of LUAD patients (P < 0.05). Importantly, four genes (PDGFB, CCND2, CTF1, IL7R) identified were strongly associated with STAT3 activation and drugs sensitivity in GDSC.ConclusionsIn summary, a ceRNA network of hsa_circ_0049271 was successfully constructed, which including one circRNA, two miRNAs, and seven mRNAs. Seven mRNAs (PDGFB, TNFRSF19, CCND2, CTF1, IL11RA, IL7R and MAOA) were remarkably associated with the prognosis of LUAD patients. Among seven mRNA species, four genes (PDGFB, CCND2, CTF1, and IL7R) could be considered as drug targets in LUAD. Our research will provide new insights into the prognosis-related ceRNA network in LUAD.

Project description:BackgroundA competitive endogenous RNA (ceRNA) network was constructed to examine the potential mechanisms of circular RNAs (circRNAs) in lung adenocarcinoma (LUAD).MethodsLUAD-related data sets were obtained from the Gene Expression Omnibus (GEO) database and screened for differentially expressed circRNAs (DECs) and differentially expressed microRNAs (DEMs). We identified the target microRNAs (miRNAs) regulated by the DECs and the potential target genes of the miRNA. The basic structure of the DECs were analyzed and enrichment analyses were conducted to determine the function of the circRNA. The Kaplan-Meier method for survival analysis was used to examine the clinical data from The Cancer Genome Atlas (TCGA) database. A protein-protein interaction (PPI) network was constructed to determine the hub genes. The relative expression of the RNA molecules on the ceRNA axis was verified by quantitative real-time polymerase chain reaction (qRT-PCR) and Western blot analysis.ResultsA total of 17 DECs and 237 DEMs were selected for analysis. After reviewing the cancer-specific circRNA database (CSCD), 10 circRNAs were identified. The 432 target miRNAs were screened by circRNA interactome (CRI) and cross-referenced with the DEMs to obtain 126 miRNAs of interest. The expression of miR-3611, which is regulated by hsa_circ_0031968, was found to significantly affect the survival and prognosis of patients with LUAD (P≤0.05). The target gene function of hsa_circ_0031968 was determined to be mainly enriched in SMAD binding, and the signaling pathway was primarily enriched in miRNAs related to cancer. The TCGA database screened out 2,484 differentially expressed mRNAs (DEmRNAs) and intersection analysis with the target gene of miR-3611 revealed 1 gene, namely the proglucagon gene (GCG). Therefore, we chose the hsa_circ_0031968/miR-3611/GCG axis for further research. The expression of GCG was determined to be associated with a poorer survival rate and higher T stage in LUAD patients. Finally, 17 hub genes that interact with GCG were identified.ConclusionsThe ceRNA regulatory network hsa_circ_0031968/miR-3611/GCG was successfully constructed and this provided novel insights into the identification of biomarkers and the pathogenesis of LUAD. This knowledge will contribute to the early diagnosis and development of potential treatment for patients with LUAD.

Project description:Purpose:Accumulating evidence has indicated that circRNAs are closely involved in tumorigenesis and progression of human cancers. However, the molecular mechanism underlying function of circRNAs in breast cancer has not been thoroughly elucidated. Currently, we aimed to characterize the circRNA-related competing endogenous RNA (ceRNA) regulatory network in breast cancer and to construct prognostic model. Materials and Methods:First, we constructed circRNA expression profiles for paired breast cancer in a Chinese population using a human circRNA microarray. Expression profiles of circRNAs, miRNAs, and mRNAs were retrieved from our circRNA dataset, the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. We applied the limma and edgeR packages to identify differentially expressed RNAs. Weighted gene correlation network analysis (WGCNA) was used to identify critical modules of mRNAs. Next, a ceRNA network was established based on circRNA-miRNA and miRNA-mRNA intersections. Both Cox regression analysis and ROC curve analysis were performed to generate prognostic model. Additionally, we performed Gene Set Enrichment Analysis (GSEA) on prognostic signatures. Results:Total of 59 circRNAs, 98 miRNAs and 3966 mRNAs were identified as differentially expressed RNAs. We first identified 38 miRNA-mRNA pairs and 38 circRNA-miRNA pairs to construct the circRNA-miRNA-mRNA regulatory network and then generated a prognostic model based on 7 signatures (MMD, SLC29A4, CREB5, FOS, ANKRD29, MYOCD, and PIGR), and patients with high-risk scores presented poor prognosis. Several cancer-related pathways were enriched, including the TGF-? pathway, the focal adhesion pathway, and the JAK-STAT signaling pathway, and 20 prognostic ceRNA regulatory networks were subsequently identified. Conclusion:In all, we screened a series of dysregulated circRNAs, miRNAs, and mRNAs, and constructed circRNA-related ceRNA network in breast cancer. Our findings may help to deepen the understanding of circRNA-related regulatory mechanisms. Moreover, we generated a prognostic model that provided new insight into postoperative management for breast cancer.

Project description:BackgroundMilk fat is one of the main reference elements for evaluating milk quality and is a primary objective trait in dairy cattle breeding. In recent years, circular RNAs (circRNAs) have been found to play crucial roles in many biological processes. However, the function and expression profiles of circRNAs in milk fat synthesis in cows are not completely understood. We performed RNA sequencing to analyze the genome-wide expression of circRNA transcripts in bovine mammary epithelial cells (BMECs) from cows with extreme differences in milk fat percentage. We identified candidate differential circRNAs associated with milk fat metabolism using functional enrichment analysis and constructed a lipid metabolism-related competing endogenous RNA (ceRNA) interactive regulatory network.ResultsA total of 290 circRNAs were significantly differentially expressed (DE-circRNAs) in high milk fat percentage (HMF) cows compared to that in low milk fat percentage (LMF) cows. Of the 290 circRNAs, 142 were significantly upregulated and 148 were significantly downregulated. Enrichment analysis (Gene Ontology and Kyoto Encyclopedia of Genes and Genomes) identified four DE-circRNAs (circ_0001122, circ_0007367, circ_0018269, and circ_0015179) that potentially regulate milk fat metabolism. Among them, circ_0001122, circ_0007367, and circ_0015179 had relatively high expression levels in cow mammary gland tissue compared to other tissues (heart, liver, kidney, uterus, ovaries, and small intestine) of cows. The regulatory networks circ_0001122:miR-12043:LIPG, circ_0007367:miR-331-3p:CIDEA/PML, and circ_0018269:miR-11989:RORC/HPX are potential networks to explore the mechanism of milk fat regulation.ConclusionsThese results reveal the possible role of circRNAs in milk fat metabolism in dairy cows. Several important circRNAs and ceRNAs affecting milk fat synthesis were identified, providing insights into the complex biology of milk fat synthesis as well as a novel theoretical perspective for future research on lactation, milk quality, and breed improvement in dairy cows.

Project description:Background and purposeAccumulating evidence suggests that circular RNAs (circRNAs) are involved in immune and inflammatory processes after acute ischemic stroke (AIS). However, the roles of circRNA-mediated competing endogenous RNA (ceRNA) in modulating immune inflammation of AIS have not yet been determined. This study aimed to construct a circRNA-mediated immune-related ceRNA network and identify novel circRNAs in AIS.MethodsMicroarray data were downloaded from the GEO database and further analysed by R software. Then, we constructed a circRNA-mediated ceRNA network based on interaction information from the bioinformatics database. A topological property analysis of the ceRNA network was conducted to screen novel circRNAs. Finally, we further applied quantitative real-time polymerase chain reaction (qRT-PCR) to two independent sets.ResultsWe constructed an AIS immune-related ceRNA (AISIRC) network containing immune-related genes (IRGs), miRNAs, and circRNAs. Additionally, we extracted the subnetwork from the AISIRC network and screened six immune-related circRNAs. After identification and validation, we finally confirmed that plasma levels of circPTP4A2 and circTLK2 were significantly increased in AIS patients compared with both healthy control subjects (HCs) and transient ischemic attack (TIA) patients. Logistic regression and receiver-operating characteristic (ROC) curve analyses demonstrated that these two circRNAs may function as predictive and discriminative biomarkers for AIS. We also confirmed that plasma levels of circPTP4A2 were elevated in TIA patients compared with HCs and might be an independent risk factor for predicting TIA. Longitudinal analysis of circRNA expression up to 90 days after AIS indicated that the ability of circPTP4A2 and circTLK2 to monitor AIS dynamics was highly desirable.ConclusionIn summary, the circRNA-mediated immune-related ceRNA network was successfully constructed, and two circulating circRNAs (circPTP4A2 and circTLK2) improved sensitivity for the diagnosis of AIS and could be considered diagnostic biomarkers.