Dataset Information

A construction and comprehensive analysis of ceRNA networks and infiltrating immune cells in papillary renal cell carcinoma.

ABSTRACT:

Background

As 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.

Methods

CeRNA 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.

Results

COL1A1, 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).

Conclusion

The 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.

SUBMITTER: Fan Y

PROVIDER: S-EPMC8607257 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

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.

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.

Project description:Background: Kidney renal clear cell carcinoma (KIRC) has the highest invasion, mortality and metastasis of the renal cell carcinomas and seriously affects patient's quality of life. However, the composition of the immune microenvironment and regulatory mechanisms at transcriptomic level such as ceRNA of KIRC are still unclear. Methods: We constructed a ceRNA network associated with KIRC by analyzing the long non-coding RNA (lncRNA), miRNA and mRNA expression data of 506 tumor tissue samples and 71 normal adjacent tissue samples downloaded from The Cancer Genome Atlas (TCGA) database. In addition, we estimated the proportion of 22 immune cell types in these samples through "The Cell Type Identification by Estimating Relative Subsets of RNA Transcripts." Based on the ceRNA network and immune cells screened by univariate Cox analysis and Lasso regression, two nomograms were constructed to predict the prognosis of patients with KIRC. Receiver operating characteristic curves (ROC) and calibration curves were employed to assess the discrimination and accuracy of the nomograms. Consequently, co-expression analysis was carried out to explore the relationship between each prognostic gene in a Cox proportional hazards regression model of ceRNA and each survival-related immune cell in a Cox proportional hazards regression model of immune cell types to reveal the potential regulatory mechanism. Results: We established a ceRNA network consisting of 12 lncRNAs, 25 miRNAs and 136 mRNAs. Two nomograms containing seven prognostic genes and two immune cells, respectively, were successfully constructed. Both ROC [area under curves (AUCs) of 1, 3, and 5-year survival in the nomogram based on ceRNA network: 0.779, 0.747, and 0.772; AUCs of 1, 3, and 5-year survivals in nomogram based on immune cells: 0.603, 0.642, and 0.607] and calibration curves indicated good accuracy and clinical application value of both models. Through co-correlation analysis between ceRNA and immune cells, we found both LINC00894 and KIAA1324 were positively correlated with follicular helper T (Tfh) cells and negatively correlated with resting mast cells. Conclusion: Based on the ceRNA network and tumor-infiltrating immune cells, we constructed two nomograms to predict the survival of KIRC patients and demonstrated their value in improving the personalized management of KIRC.