The Construction and Analysis of Tumor-Infiltrating Immune Cells and ceRNA Networks in Bladder Cancer
Ontology highlight
ABSTRACT: Background
Bladder cancer (BLCA) is the 11th most common malignancy worldwide. Although significant improvements have been made in screening, diagnosis, and precise management in recent years, the prognosis of BLCA remains bleak. Objectives
This study aimed to investigate the prognostic significance of tumor-infiltrating immune cells and construct ceRNA networks in BLCA patients. Methods
The expression data of BLCA patients were obtained from The Cancer Genome Atlas (TCGA) database. A competing endogenous RNA (ceRNA) network was constructed to identify the hub genes involved in the prognosis of BLCA. The CIBERSORT algorithm was utilized to investigate the infiltration levels of 22 subsets of immune cells. Ultimately, the nomogram was generated to visualize the survival probability of each patient, with the calibration curve being performed to assess its performance. Furthermore, the Pearson correlation test was used to explore the correlation between the identified hub genes in the ceRNA network and the prognostic-related immune cells. Results
A total of eight elements in the ceRNA network were considered as key members and correlated with the prognosis of BLCA, including ELN, SREBF1, DSC2, TTLL7, DIP2C, SATB1, hsa-miR-20a-5p, and hsa-miR-29c-3p. T cells CD8, T cells follicular helper (Tfh), and neutrophils were identified as independent prognostic factors in BLCA. The co-expression analysis showed that there was a significant correlation between the identified hub genes and immune cells. Conclusion
Our results suggest that the mechanism of hsa-miR-29c-3p regulates the expression of ELN and DSC2, and the infiltration of Tfh and neutrophils might play pivotal roles in the progression of BLCA.
SUBMITTER: Jiang A
PROVIDER: S-EPMC7775311 | biostudies-literature | 2020 Jan
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:PURPOSE:To construct and analyze tumor-infiltrating immune cell and ceRNA (competitive endogenous RNA) networks in metastatic adrenal cortical carcinoma (ACC). METHODS:A ceRNA network was established to identify the ceRNAs involved in metastasis of ACC based on 92 samples from TCGA, including 18 cases of metastasis and 74 cases of non-metastatic primary tumors. And the algorithm "cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT)" was used to quantify the proportion of immune cells in ACC. In addition, predictive nomograms based on the types of important immune cells or ceRNAs were constructed to predict ACC prognosis. Moreover, we evaluated the relationships between metastatic ACC-specific immune cells and ceRNA networks to identify the potential immune gene characteristics. RESULTS:Ten prognostic biomarkers were identified as key members of the ceRNA network and three tumor-infiltrating immune cells were identified by CIBERSORT algorithm. Some important co-expression patterns between immune cells and ceRNAs network indicate significant correlation between Macrophages M0 and hsa-miR-130b-3p (P < 0.001), Macrophages M0 and H2AFX (P = 0.003). CONCLUSIONS:The present study inferred that the metastasis-related ceRNAs of H2AFX, hsa-miR-130b-3p and Macrophages M0 might play important roles in ACC metastasis.
| S-EPMC7103591 | biostudies-literature
Project description:Background/Aims: As a malignant and melanocytic tumor, cutaneous melanoma is the devastating skin tumor with high rates of recurrence and metastasis. Bone is the common metastatic location, and bone metastasis may result in pathologic fracture, neurologic damage, and severe bone pain. Although metastatic melanoma was reported to get benefits from immunotherapy, molecular mechanisms and immune microenviroment underlying the melanoma bone metastasis and prognostic factors are still unknown. Methods: Gene expression profiling of 112 samples, including 104 primary melanomas and 8 bone metastatic melanomas from The Cancer Genome Atlas database, was assayed to construct a ceRNA network associated with bone metastases. Besides, we detected the fraction of 22 immune cell types in melanoma via the algorithm of "cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT)." Based on the significant ceRNAs or immune cells, we constructed nomograms to predict the prognosis of patients with melanoma. Ultimately, correlation analysis was implemented to discover the relationship between the significant ceRNA and immune cells to reveal the potential signaling pathways. Results: We constructed a ceRNA network based on the interaction among 8 pairs of long noncoding RNA-microRNA and 15 pairs of microRNA-mRNA. CIBERSORT and ceRNA integration analysis discovered that AL118506.1 has both significant prognostic value (P = 0.002) and high correlation with T follicular helper cells (P = 0.033). Meanwhile, T cells CD8 and macrophages M2 were negatively correlated (P < 0.001). Moreover, we constructed two satisfactory nomograms (area under curve of 3-year survival: 0.899; 5-year survival: 0.885; and concordance index: 0.780) with significant ceRNAs or immune cells, to predict the prognosis of patients. Conclusions: In this study, we suggest that bone metastasis in melanoma might be related to AL118506.1 and its role in regulating thrombospondin 2 and T follicular helper cells. Two nomograms were constructed to predict the prognosis of patients with melanoma and demonstrated their value in improving the personalized management.
| S-EPMC6774271 | biostudies-literature
Project description:Soft tissue sarcoma (STS) is one of the most challenging tumors for medical oncologists, with a high rate of recurrence after initial resection. In this study, a recurrent STS-specific competitive endogenous RNA (ceRNA) network including seven recurrence and overall survival (OS)-associated genes (LPP-AS2, MUC1, GAB2, hsa-let-7i-5p, hsa-let-7f-5p, hsa-miR-101-3p and hsa-miR-1226-3p) was established based on the gene expression profiling of 259 primary sarcomas and 3 local recurrence samples from the TCGA database. The algorithm "cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT)" was applied to estimate the fraction of immune cells in sarcomas. Based on 5 recurrence and OS-associated immune cells (NK cells activated, dendritic cells resting, mast cells resting, mast cells activated and macrophages M1), we constructed a recurrent STS-specific immune cells network. Both nomograms were identified to have good reliabilities (Area Under Curve (AUC) of 5-year survival is 0.724 and 0.773, respectively). Then the co-expression analysis was performed to identify the potential regulation network among recurrent STS-specific immune cells and ceRNAs. Hsa-miR-1226-3p and MUC1 were significantly correlated and dendritic cells resting was related to hsa-miR-1226-3p. Additionally, the expression of MUC1 and dendritic cell marker CD11c were also verified by immunohistochemistry (IHC) assay and multidimensional databases. In conclusion, this study illustrated the potential mechanism of hsa-miR-1226-3p regulating MUC1 and dendritic cells resting might play an important role in STS recurrence. These findings might provide potential prognostic biomarkers and therapeutic targets for recurrent STS.
| S-EPMC6914407 | biostudies-literature
Project description:BackgroundAs the second most common malignancy in adults, papillary renal cell carcinoma (PRCC) has shown an increasing trend in both incidence and mortality. Effective treatment for advanced metastatic PRCC is still lacking. In this study, we aimed to establish competitive endogenous RNA (ceRNA) networks related to PRCC tumorigenesis, and analyze the specific role of differentially expressed ceRNA components and infiltrating immune cells in tumorigenesis.MethodsCeRNA networks were established to identify the key ceRNAs related to PRCC tumorigenesis based on the 318 samples from The Cancer Genome Atlas database (TCGA), including 285 PRCC and 33 normal control samples. The R package, "CIBERSORT," was used to evaluate the infiltration of 22 types of immune cells. Then we identified the significant ceRNAs and immune cells, based on which two nomograms were obtained for predicting the prognosis in PRCC patients. Finally, we investigated the co-expression of PRCC-specific immune cells and core ceRNAs via Pearson correlation test.ResultsCOL1A1, H19, ITPKB, LDLR, TCF4, and WNK3 were identified as hub genes in ceRNA networks. Four prognostic-related tumor-infiltrating immune cells, including T cells CD4 memory resting, Macrophages M1, and Macrophages M2 were revealed. Pearson correlation test indicated that Macrophage M1 was negatively related with COL1A1 (p < 0.01) and LDLR (p < 0.01), while Macrophage M2 was positively related with COL1A1 (p < 0.01), TCF4 (p < 0.01), and H19 (p = 0.032). Two nomograms were conducted with favorable accuracies (area under curve of 1-year survival: 0.935 and 0.877; 3-year survival: 0.849 and 0.841; and 5-year survival: 0.818 and 0.775, respectively).ConclusionThe study constructed two nomograms suited for PRCC prognosis predicting. Moreover, we concluded that H19-miR-29c-3p-COL1A1 axis might promote the polarization of M2 macrophages and inhibit M1 macrophage activation through Wnt signaling pathway, collaborating to promote PRCC tumorigenesis and lead to poor overall survival of PRCC patients.
| S-EPMC8607257 | biostudies-literature
Project description:Ulcerative colitis (UC) is a common disease with great variability in severity, with a high recurrence rate and heavy disease burden. In recent years, the different biological functions of competing endogenous RNA (ceRNA) networks of long noncoding RNAs (lncRNAs) and microRNAs (miRs) have aroused wide concerns, the ceRNA network of ulcerative colitis (UC) may have potential research value, and these expressed noncoding RNAs may be involved in the molecular basis of inflammation recurrence and progression. This study analyzed 490 colon samples associated with UC from 4 gene expression microarrays from the GEO database and identified gene modules by weighted correlation network analysis (WGCNA). CIBERSORT detected tissue-infiltrating leukocyte profiling by deconvolution of microarray data. LncBase and multiMIR were used to identify lncRNA-miRNA-mRNA interaction. We constructed a ceRNA network which includes 4 lncRNAs (SH3BP5-AS1, MIR4435-2HG, ENTPD1-AS1, and AC007750.1), 5 miRNAs (miR-141-3p, miR-191-5p, miR-192-5p, miR-194-5p, and miR196-5p), and 52 mRNAs. Those genes are involved in interleukin family signals, neutrophil degranulation, adaptive immunity, and cell adhesion pathways. lncRNA MIR4435-2HG is a variable in the decision tree for moderate-to-severe UC diagnostic prediction. Our work identifies potential regulated inflammation-related lncRNA-miRNA-mRNA regulatory axes. The regulatory axes are dysregulated during the deterioration of UC, suggesting that it is a risk factor for UC progression.
| S-EPMC8277522 | biostudies-literature
Project description:Esophageal cancer (ESCA) is one of the most commonly diagnosed cancers in the world. Tumor immune microenvironment is closely related to tumor prognosis. The present study aimed at analyzing the competing endogenous RNA (ceRNA) network and tumor-infiltrating immune cells in ESCA. The expression profiles of mRNAs, lncRNAs, and miRNAs were downloaded from the Cancer Genome Atlas database. A ceRNA network was established based on the differentially expressed RNAs by Cytoscape. CIBERSORT was applied to estimate the proportion of immune cells in ESCA. Prognosis-associated genes and immune cells were applied to establish prognostic models basing on Lasso and multivariate Cox analyses. The survival curves were constructed with Kaplan-Meier method. The predictive efficacy of the prognostic models was evaluated by the receiver operating characteristic (ROC) curves. The differentially expressed mRNAs, lncRNAs, and miRNAs were identified. We constructed the ceRNA network including 23 lncRNAs, 19 miRNAs, and 147 mRNAs. Five key molecules (HMGB3, HOXC8, HSPA1B, KLHL15, and RUNX3) were identified from the ceRNA network and five significant immune cells (plasma cells, T cells follicular helper, monocytes, dendritic cells activated, and neutrophils) were selected via CIBERSORT. The ROC curves based on key genes and significant immune cells all showed good sensitivity (AUC of 3-year survival: 0.739, AUC of 5-year survival: 0.899, AUC of 3-year survival: 0.824, AUC of 5-year survival: 0.876). There was certain correlation between five immune cells and five key molecules. The present study provides an effective bioinformatics basis for exploring the potential biomarkers of ESCA and predicting its prognosis.
| S-EPMC8164107 | biostudies-literature
Project description:Background Diabetic foot ulcer (DFU) is a severe complication characterized by low-grade infectious inflammation and probably associated with specific competitive endogenous RNAs (ceRNAs) and infiltrating immune cells. Nonetheless, no reliable biomarkers are used for detecting infectious inflammation in DFU. Therefore, it is essential to explore potential biomarkers for the accurate diagnosis and treatment of DFU. Methods The gene expression profile was retrieved from Gene Expression Omnibus (GEO) database and divided into two groups, namely, standard samples and DFU samples. To establish the ceRNA networks, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were utilized to analyze differential expression genes (DEGs). The cell type identification was achieved by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm to screen-specific immune-infiltrating cells associated with DFU. Results A ceRNA network was constructed with 20 differential expression circRNA (DEcircRNAs), 11 differential expression microRNAs (DEmiRNAs), and 9 differential expression mRNAs (DEmRNAs). Functional enrichment analysis demonstrated that DFU was mainly enriched in vascular endothelial growth factor (VEGF) and T-cell receptor signaling. In addition, CIBERSORT estimation indicated that CD8+ T cells and Monocytes were significantly related to the expression of IL-6, a DFU-specific infectious inflammation factor. Conclusion This study identified that some significant ceRNAs (JUNB, GATA3, hsa-circ-0049271 and hsa-circ-0074559) and infiltrating immune cells (CD8+ T cells and monocytes) might be related to DFU infectious inflammation.
| S-EPMC9329527 | biostudies-literature
Project description:BackgroundBladder cancer (BC) is the ninth most common malignant tumor, accounting for an estimate of 549,000 new BC cases and 200,000 BC-related deaths worldwide in 2018. The prognosis of BC has not substantially improved despite significant advances in the diagnosis and treatment of the disease.MethodsThe RNA sequencing (RNA-seq) data and clinical information of BC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm was used to assess immune infiltration. The survival analyses were performed using the selected components of a ceRNA network and selected immune cell types by least absolute shrinkage and selection operator (LASSO) Cox regression to calculate the risk score. The accuracy of prognosis prediction was determined by receiver operating characteristic (ROC) curves, survival curves, and nomograms. Finally, the correlation analysis was performed to investigate the relationships between the signature components of the ceRNA network and the immune cell signature.ResultsTwo completed survival analyses included selected components of the ceRNA network (ELN, SREBF1, DSC2, TTLL7, DIP2C, SATB1, hsa-miR-20a-5p, and hsa-miR-29c-3p) and selected immune cell types (M0 macrophages, M2 macrophages, resting mast cells, and neutrophils). ROC curves, survival curves (all P values <0.05), nomograms, and calibration curves indicated that the accuracy of the two survival analyses was acceptable. Moreover, the correlations between TTLL7 and resting mast cells (R=0.24, P<0.001), DSC2 and resting mast cells (R=-0.23, P<0.001), ELN and resting mast cells (R=0.44, P<0.001), and hsa-miR-29c-3p and M0 macrophages (R=-0.29, P<0.001) were significant, indicating that interactions of these factors may play significant roles in the prognosis of BC.ConclusionsTTLL7, DSC2, ELN, hsa-miR-29c-3p, resting mast cells, and M0 macrophages may play an important role in the development of BC. However, additional studies are needed to confirm this hypothesis.
| S-EPMC8185653 | biostudies-literature
Project description:BackgroundMuscle-invasive bladder cancer (MIBC) is a lethal disease with poor treatment response and a high death rate. Immune cells infiltrating the tumor tissues have been shown to play a vital role in tumorigenesis and tumor progression, but their prognostic significance in MIBC remains unclear.ObjectivesTo explore the landscape and prognostic significance of tumor-infiltrating immune cells (TIICs) in MIBC, and to develop a model to improve the prognostic predictions of MIBC.Methods and materialsThe gene expression profile and clinical data of MIBC patients were downloaded from the Gene Expression Omnibus and The Cancer Genome Atlas portal. The fractions of 22 TIIC subtypes were calculated using the Cell Type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm. A TIICs-based model was constructed using least absolute shrinkage and selection operator (LASSO) Cox regression in a training cohort and validated in the validation cohort.ResultsTen types of TIICs demonstrated different infiltration abundance between MIBC and normal tissue. We also found 11 types of TIICs that were significantly associated with overall survival (OS). A TIICs-based model was established, consisting of 15 types of immune cells, and an immunoscore was calculated. Significant differences in OS were found between the high and low immunoscore groups, in both training (n = 343) and validation (n = 146) cohorts. The model could identify patients who would have worse OS despite having similar clinical characteristics. Furthermore, multivariate analysis identified the immunoscore as an independent risk factor (hazard ratio, 3.23; 95% confidence interval; 2.22-4.70) for OS in MIBC patients.ConclusionThe landscape of immune infiltration is different between MIBC and normal tissue. The TIICs-based model could provide promising predictive value to complement the existing staging system for predicting the OS of MIBC patients.
| S-EPMC7653560 | biostudies-literature
Project description:BackgroundAt present, the mechanisms behind atrial fibrillation (AF) pathogenesis are still unclear. We construct a ceRNA immunoregulatory network to further understand the mechanism of AF.MethodsFour AF mRNA datasets from the Gene Expression Omnibus (GEO) database were integrated by SVA method. AF-related immune genes (AF-IRGs) were selected via combining ImmPort database with the genes in the module most associated with AF obtained by a weighted gene coexpression network analysis (WGCNA). Then, circRNA and miRNA expressions from the GEO database were extracted and mapped with related databases. Next, an immune-related circRNA-miRNA-mRNA ceRNA network was constructed and hub genes were filtered from a protein-protein interaction (PPI) network, and the differentially expressed (DE) hub genes in AF were further screened. Additionally, immune infiltration was investigated in AF by using CIBERSORT. Subsequently, the relationships between DE hub genes and AF-related infiltrating immune cells were performed by using Pearson correlation coefficients. Ulteriorly, the immune-cells-related ceRNA subnetwork in AF was built.ResultsA total of 95 AF-IRGs were detected, and an immune-related ceRNA network in AF was constructed with 12 circRNAs, 7 miRNAs and 50 mRNAs. The immune infiltration analysis indicated that a higher level of neutrophils, as well as a lower level of T cells regulatory (Tregs) and NK cells activated in AF. Four DE hub genes (CXCL12, IL7R, TNFSF13B, CD8A) were associated with Tregs or NK cells activated immune cells (P < 0.05). Tregs or NK cells activated immune cells-related ceRNA subnetwork including 5 circRNAs (has_circ_0001190, has_circ_0006725, has_circ_0079284, has_circ_0005299, and has_circ_0002103), 4 miRNAs (has-miR-198, has-miR-623, has-miR-1246, and has-miR-339-3p) and 4 DE hub genes was eventually constructed in AF.ConclusionOur results provide new insights into the molecular mechanisms governing AF progression from the perspective of immune-related ceRNA network.
| S-EPMC8643171 | biostudies-literature