scholarly journals Construction of Prognostic Risk Model of 5-Methylcytosine-Related Long Non-Coding RNAs and Evaluation of the Characteristics of Tumor-Infiltrating Immune Cells in Breast Cancer

2021 ◽  
Vol 12 ◽  
Author(s):  
Zhidong Huang ◽  
Junjing Li ◽  
Jialin Chen ◽  
Debo Chen

Purpose: The role of 5-methylcytosine-related long non-coding RNAs (m5C-lncRNAs) in breast cancer (BC) remains unclear. Here, we aimed to investigate the prognostic value, gene expression characteristics, and correlation between m5C-lncRNA risk model and tumor immune cell infiltration in BC.Methods: The expression matrix of m5C-lncRNAs in BC was obtained from The Cancer Genome Atlas database, and the lncRNAs were analyzed using differential expression analysis as well as univariate and multivariate Cox regression analysis to eventually obtain BC-specific m5C-lncRNAs. A risk model was developed based on three lncRNAs using multivariate Cox regression and the prognostic value, accuracy, as well as reliability were verified. Gene set enrichment analysis (GSEA) was used to analyze the Kyoto Encyclopedia of Genes and Genomes signaling pathway enrichment of the risk model. CIBERSORT algorithm and correlation analysis were used to explore the characteristics of the BC tumor-infiltrating immune cells. Finally, reverse transcription-quantitative polymerase chain reaction was performed to detect the expression level of three lncRNA in clinical samples.Results: A total of 334 differential m5C-lncRNAs were identified, and three BC-specific m5C-lncRNAs were selected, namely AP005131.2, AL121832.2, and LINC01152. Based on these three lncRNAs, a highly reliable and specific risk model was constructed, which was proven to be closely related to the prognosis of patients with BC. Therefore, a nomogram based on the risk score was built to assist clinical decisions. GSEA revealed that the risk model was significantly enriched in metabolism-related pathways and was associated with tumor immune cell infiltration based on the analysis with the CIBERSORT algorithm.Conclusion: The efficient risk model based on m5C-lncRNAs associated with cancer metabolism and tumor immune cell infiltration could predict the survival prognosis of patients, and AP005131.2, AL121832.2, and LINC01152 could be novel biomarkers and therapeutic targets for BC.

BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Shaokun Wang ◽  
Li Pang ◽  
Zuolong Liu ◽  
Xiangwei Meng

Abstract Background The change of immune cell infiltration essentially influences the process of colorectal cancer development. The infiltration of immune cells can be regulated by a variety of genes. Thus, modeling the immune microenvironment of colorectal cancer by analyzing the genes involved can be more conducive to the in-depth understanding of carcinogenesis and the progression thereof. Methods In this study, the number of stromal and immune cells in malignant tumor tissues were first estimated by using expression data (ESTIMATE) and cell-type identification with relative subsets of known RNA transcripts (CIBERSORT) to calculate the proportion of infiltrating immune cell and stromal components of colon cancer samples from the Cancer Genome Atlas database. Then the relationship between the TMN Classification and prognosis of malignant tumors was evaluated. Results By investigating differentially expressed genes using COX regression and protein-protein interaction network (PPI), the candidate hub gene serine protease inhibitor family E member 1 (SERPINE1) was found to be associated with immune cell infiltration. Gene Set Enrichment Analysis (GSEA) further projected the potential pathways with elevated SERPINE1 expression to carcinogenesis and immunity. CIBERSORT was subsequently utilized to investigate the relationship between the expression differences of SERPINE1 and immune cell infiltration and to identify eight immune cells associated with SERPINE1 expression. Conclusion We found that SERPINE1 plays a role in the remodeling of the colon cancer microenvironment and the infiltration of immune cells.


Author(s):  
Jinguo Zhang ◽  
Benjie Shan ◽  
Lin Lin ◽  
Jie Dong ◽  
Qingqing Sun ◽  
...  

Breast cancer (BC) represents a molecularly and clinically heterogeneous disease. Recent progress in immunotherapy has provided a glimmer of hope for several BC subtypes. The relationship between N6-methyladenosine (m6A) modification and long non-coding RNAs (LncRNAs) is still largely unexplored in BC. Here, with the intention to dissect the landscape of m6A-related lncRNAs and explore the immunotherapeutic value of the m6A-related lncRNA signature, we identified m6A-related lncRNAs by co-expression analysis from The Cancer Genome Atlas (TCGA) and stratified BC patients into different subgroups. Furthermore, we generated an m6A-related lncRNA prognostic signature. Four molecular subtypes were identified by consensus clustering. Cluster 3 preferentially had favorable prognosis, upregulated immune checkpoint expression, and high level of immune cell infiltration. Twenty-one m6A-related lncRNAs were applied to construct the m6A-related lncRNA model (m6A-LncRM). Survival analysis and receiver operating characteristic (ROC) curves further confirmed the prognostic value and prediction performance of m6A-LncRM. Finally, high- and low-risk BC subgroups displayed significantly different clinical features and immune cell infiltration status. Overall, our study systematically explored the prognostic value of the m6A-related LncRNAs and identified a high immunogenicity BC subtype. The proposed m6A-related LncRNA model might serve as a robust prognostic signature and attractive immunotherapeutic targets for BC treatment.


2021 ◽  
Vol 12 ◽  
Author(s):  
Qianhui Xu ◽  
Shaohuai Chen ◽  
Yuanbo Hu ◽  
Wen Huang

BackgroundIncreasing evdence supports the suggestion that the immune cell infiltration (ICI) patterns play a pivotal role in tumor progression in breast cancer (BRCA). Nonetheless, there has been no comprehensive analysis of the ICI patterns effects on the clinical outcomes and immunotherapy.MethodsMultiomic data for BRCA samples were downloaded from TCGA. ESTIMATE algorithm, ssGSEA method, and CIBERSORT analysis were used to uncover the landscape of the tumor immune microenvironment (TIME). BRCA subtypes based on the ICI pattern were identified by consensus clustering and principal-component analysis was performed to obtain the ICI scores to quantify the ICI patterns in individual tumors. Their prognostic value was validated by the Kaplan-Meier survival curves. Gene set enrichment analysis (GSEA) was applied for functional annotation. Immunophenoscore (IPS) was employed to explore the immunotherapeutic role of the ICI scores. Finally, the mutation data was analyzed by using the “maftools” R package.ResultsThree different immune infiltration patterns with a distinct prognosis and biological signature were recognized among 1,198 BRCA samples. The characteristics of TIME under these three patterns were highly consistent with three known immune profiles: immune- excluded, immune-desert, and immune-inflamed phenotypes, respectively. The identification of the ICI patterns within individual tumors based on the ICI score, developed under the ICI-related signature genes, contributed into dissecting biological processes, clinical outcome, immune cells infiltration, immunotherapeutic effect, and genetic variation. High ICI score subtype, characterized with a suppression of immunity, suggested an immune-exhausted phenotype. Abundant effective immune cells were discovered in the low ICI score patients, which corresponded to an immune-activated phenotype and might present an immunotherapeutic advantage. Immunophenoscore was implemented as a surrogate of immunotherapeutic outcome, low-ICI scores samples obtained a significantly higher immunophenoscore. Enrichment of the JAK/STAT and VEGF signal pathways were activated in the ICI low-score subgroup. Finally, the synergistic effect between the ICI score and the tumor mutation burden (TMB) was confirmed.ConclusionThis work comprehensively elucidated that the ICI patterns served as an indispensable player in complexity and diversity of TIME. Quantitative identification of the ICI patterns in individual tumor will contribute into mapping the landscape of TIME further optimizing precision immunotherapy.


2021 ◽  
Vol 12 ◽  
Author(s):  
Chuang Zhang ◽  
Danni Li ◽  
Ruoxi Yu ◽  
Ce Li ◽  
Yujia Song ◽  
...  

BackgroundGastric cancer (GC) still represents the third leading cause of cancer-related death worldwide. Peritoneal relapse (PR) is the most frequent metastasis occurring among patients with advanced gastric cancer. Increasingly more evidence have clarified the tumor immune microenvironment (TIME) may predict survival and have clinical significance in GC. However, tumor-transcriptomics based immune signatures derived from immune profiling have not been established for predicting the peritoneal recurrence of the advanced GC.MethodsIn this study, we depict the immune landscape of GC by using transcriptome profiling and clinical characteristics retrieved from GSE62254 of Gene Expression Omnibus (GEO). Immune cell infiltration score was evaluated via single-sample gene set enrichment (ssGSEA) analysis algorithm. The least absolute shrinkage and selection operator (LASSO) Cox regression algorithm was used to select the valuable immune cells and construct the final model for the prediction of PR. The receiver operating characteristic (ROC) curve and the Kaplan-Meier curve were used to check the accuracy of PRIs. Gene Set Enrichment Analysis (GSEA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to explore the molecular pathways associated with PRIs.ResultsA peritoneal recurrence related immune score (PRIs) with 10 immune cells was constructed. Compared to the low-PRIs group, the high-PRIs group had a greater risk. The upregulation of the focal adhesion signaling was observed in the high-PRIs subtype by GSEA and KEGG. Multivariate analysis found that both in the internal training cohort and the internal validation cohort, PRIs was a stable and independent predictor for PR. A nomogram that integrated clinicopathological features and PRIs to predict peritoneal relapse was constructed. Subgroup analysis indicated that the PRIs could obviously distinguish peritoneal recurrence in different molecular subtypes, pathological stages and Lauren subtypes, in which PRIs of Epithelial-Mesenchymal Transitions (EMT) subtype, III-IV stage and diffuse subtype are higher respectively.ConclusionOverall, we performed a comprehensive evaluation of the immune landscape of GC and constructed a predictive PR model based on the immune cell infiltration. The PRIs represents novel promising feature of predicting peritoneal recurrence of GC and sheds light on the improvement of the personalized management of GC patients after surgery.


2022 ◽  
Vol 12 ◽  
Author(s):  
Lan-Xin Mu ◽  
You-Cheng Shao ◽  
Lei Wei ◽  
Fang-Fang Chen ◽  
Jing-Wei Zhang

Purpose: This study aims to reveal the relationship between RNA N6-methyladenosine (m6A) regulators and tumor immune microenvironment (TME) in breast cancer, and to establish a risk model for predicting the occurrence and development of tumors.Patients and methods: In the present study, we respectively downloaded the transcriptome dataset of breast cancer from Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database to analyze the mutation characteristics of m6A regulators and their expression profile in different clinicopathological groups. Then we used the weighted correlation network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), and cox regression to construct a risk prediction model based on m6A-associated hub genes. In addition, Immune infiltration analysis and gene set enrichment analysis (GSEA) was used to evaluate the immune cell context and the enriched gene sets among the subgroups.Results: Compared with adjacent normal tissue, differentially expressed 24 m6A regulators were identified in breast cancer. According to the expression features of m6A regulators above, we established two subgroups of breast cancer, which were also surprisingly distinguished by the feature of the immune microenvironment. The Model based on modification patterns of m6A regulators could predict the patient’s T stage and evaluate their prognosis. Besides, the low m6aRiskscore group presents an immune-activated phenotype as well as a lower tumor mutation load, and its 5-years survival rate was 90.5%, while that of the high m6ariskscore group was only 74.1%. Finally, the cohort confirmed that age (p < 0.001) and m6aRiskscore (p < 0.001) are both risk factors for breast cancer in the multivariate regression.Conclusion: The m6A regulators play an important role in the regulation of breast tumor immune microenvironment and is helpful to provide guidance for clinical immunotherapy.


2021 ◽  
Vol 12 ◽  
Author(s):  
Zhuomao Mo ◽  
Daiyuan Liu ◽  
Dade Rong ◽  
Shijun Zhang

Background: Generally, hepatocellular carcinoma (HCC) exists in an immunosuppressive microenvironment that promotes tumor evasion. Hypoxia can impact intercellular crosstalk in the tumor microenvironment. This study aimed to explore and elucidate the underlying relationship between hypoxia and immunotherapy in patients with HCC.Methods: HCC genomic and clinicopathological datasets were obtained from The Cancer Genome Atlas (TCGA-LIHC), Gene Expression Omnibus databases (GSE14520) and International Cancer Genome Consortium (ICGC-LIRI). The TCGA-LIHC cases were divided into clusters based on single sample gene set enrichment analysis and hierarchical clustering. After identifying patients with immunosuppressive microenvironment with different hypoxic conditions, correlations between immunological characteristics and hypoxia clusters were investigated. Subsequently, a hypoxia-associated score was established by differential expression, univariable Cox regression, and lasso regression analyses. The score was verified by survival and receiver operating characteristic curve analyses. The GSE14520 cohort was used to validate the findings of immune cell infiltration and immune checkpoints expression, while the ICGC-LIRI cohort was employed to verify the hypoxia-associated score.Results: We identified hypoxic patients with immunosuppressive HCC. This cluster exhibited higher immune cell infiltration and immune checkpoint expression in the TCGA cohort, while similar significant differences were observed in the GEO cohort. The hypoxia-associated score was composed of five genes (ephrin A3, dihydropyrimidinase like 4, solute carrier family 2 member 5, stanniocalcin 2, and lysyl oxidase). In both two cohorts, survival analysis revealed significant differences between the high-risk and low-risk groups. In addition, compared to other clinical parameters, the established score had the highest predictive performance at both 3 and 5 years in two cohorts.Conclusion: This study provides further evidence of the link between hypoxic signals in patients and immunosuppression in HCC. Defining hypoxia-associated HCC subtypes may help reveal potential regulatory mechanisms between hypoxia and the immunosuppressive microenvironment, and our hypoxia-associated score could exhibit potential implications for future predictive models.


2021 ◽  
Vol 14 (8) ◽  
pp. 1151-1159
Author(s):  
Chen-Lu Liao ◽  
◽  
Xing-Yu Sun ◽  
Qi Zhou ◽  
Min Tian ◽  
...  

AIM: To investigate the role of tumor microenvironment (TME)-related long non-coding RNA (lncRNA) in uveal melanoma (UM), probable prognostic signature and potential small molecule drugs using bioinformatics analysis. METHODS: UM expression profile data were downloaded from the Cancer Genome Atlas (TCGA) and bioinformatics methods were used to find prognostic lncRNAs related to UM immune cell infiltration. The gene expression profile data of 80 TCGA specimens were analyzed using the single sample Gene Set Enrichment Analysis (ssGSEA) method, and the immune cell infiltration of a single specimen was evaluated. Finally, the specimens were divided into high and low infiltration groups. The differential expression between the two groups was analyzed using the R package ‘edgeR’. Univariate, multivariate and Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analyses were performed to explore the prognostic value of TME-related lncRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses were also performed. The Connectivity Map (CMap) data set was used to screen molecular drugs that may treat UM. RESULTS: A total of 2393 differentially expressed genes were identified and met the criteria for the low and high immune cell infiltration groups. Univariate Cox analysis of lncRNA genes with differential expression identified 186 genes associated with prognosis. Eight prognostic markers of TME-included lncRNA genes were established as potentially independent prognostic elements. Among 269 differentially expressed lncRNAs, 69 were up-regulated and 200 were down-regulated. Univariate Cox regression analysis of the risk indicators and clinical characteristics of the 8 lncRNA gene constructs showed that age, TNM stage, tumor base diameter, and low and high risk indices had significant prognostic value. We screened the potential small-molecule drugs for UM, including W-13, AH-6809 and Imatinib. CONCLUSION: The prognostic markers identified in this study are reliable biomarkers of UM. This study expands our current understanding of the role of TME-related lncRNAs in UM genesis, which may lay the foundations for future treatment of this disease.


2021 ◽  
Vol 11 ◽  
Author(s):  
Young-Sil An ◽  
Se-Hyuk Kim ◽  
Tae Hoon Roh ◽  
So Hyun Park ◽  
Tae-Gyu Kim ◽  
...  

BackgroundThe purpose of this study was to investigate the correlation between 18F-fluorodeoxyglucose (FDG) uptake and infiltrating immune cells in metastatic brain lesions.MethodsThis retrospective study included 34 patients with metastatic brain lesions who underwent brain 18F-FDG positron emission tomography (PET)/computed tomography (CT) followed by surgery. 18F-FDG uptake ratio was calculated by dividing the standardized uptake value (SUV) of the metastatic brain lesion by the contralateral normal white matter uptake value. We investigated the clinicopathological characteristics of the patients and analyzed the correlation between 18F-FDG uptake and infiltration of various immune cells. In addition, we evaluated immune-expression levels of glucose transporter 1 (GLUT1), hexokinase 2 (HK2), and Ki-67 in metastatic brain lesions.ResultsThe degree of 18F-FDG uptake of metastatic brain lesions was not significantly correlated with clinical parameters. There was no significant relationship between the 18F-FDG uptake and degree of immune cell infiltration in brain metastasis. Furthermore, other markers, such as GLUT1, HK2, and Ki-67, were not correlated with degree of 18F-FDG uptake. In metastatic brain lesions that originated from breast cancer, a higher degree of 18F-FDG uptake was observed in those with high expression of CD68.ConclusionsIn metastatic brain lesions, the degree of 18F-FDG uptake was not significantly associated with infiltration of immune cells. The 18F-FDG uptake of metastatic brain lesions from breast cancer, however, might be associated with macrophage activity.


2021 ◽  
Author(s):  
Xiaofen Pan ◽  
Xingkui Tang ◽  
Minling Liu ◽  
Xijun Luo ◽  
Mengyuan Zhu ◽  
...  

Abstract BackgroundTumor microenvironment consists of tumor cells, immune cells and other matric components. Tumor infiltration immune cells are associated with prognosis. But all the current prognosis evaluation system dose not take tumor immune cells other matrix component into consideration. In the current study, we aimed to construct a prognosis predictive model based on tumor microenvironment.MethodCIBERSORT and ESTIMATE algorithms were used to reveal the immune cell infiltration landscape of colon cancer. Patients were classified into three clusters by ConsensusClusterPlus algorithm. Immune cell infiltration (ICI) scores of each patient were determine by principal-component analysis. Patients were divided to high and low ICI score groups. Survival, gene expression and somatic mutation of the two groups were compared.ResultsPatients with no lymph node invasion, no metastasis, T1-2 disease and stage I-II had higher ICI scores. Calcium signaling pathway, leukocyte transendothelial migration pathway, MAPK signaling pathway, TGF β pathway, and WNT signaling pathway were enriched in high ICI score group. Immune-checkpoint genes and immune-activity associated genes were significantly decreased in high ICI score. Patients in high ICI score group had better survival than low ICI score group. Prognostic value of ICI score was independent of TMB.ConclusionICI score might serve as an independent prognostic biomarker in colon cancer.


2022 ◽  
Vol 2022 ◽  
pp. 1-24
Author(s):  
Bin-Bin Da ◽  
Shuai Luo ◽  
Ming Huang ◽  
Fei Song ◽  
Rong Ding ◽  
...  

It has been demonstrated that the inflammatory response influences cancer development and can be used as a prognostic biomarker in various tumors. However, the relevance of genes associated with inflammatory responses in hepatocellular carcinoma (HCC) remains unknown. The Cancer Genome Atlas (TCGA) database was analyzed using weighted gene coexpression network analysis (WGCNA) and differential analysis to discover essential inflammatory response-related genes (IFRGs). Cox regression studies, both univariate and multivariate, were employed to develop a prognostic IFRGs signature. Additionally, Gene Set Enrichment Analysis (GSEA) was used to deduce the biological function of the IFRGs signature. Finally, we estimated immune cell infiltration using a single sample GSEA (ssGSEA) and x-cell. Our results revealed that, among the major HCC IFRGs, two (DNASE1L3 and KLKB1) were employed to create a predictive IFRG signature. The IFRG signature could correctly predict overall survival (O.S) as per Kaplan-Meier time-dependent roc curves analysis. It was also linked to pathological tumor stage and T stage and might be used as a prognostic predictor in HCC. GSEA analysis concluded that the IFRG signature might influence the immune response in HCC. Immunological cell infiltration and immune checkpoint molecule expression differed in the high-risk and low-risk groups. As a result of our findings, DNASILE may play a role in the tumor microenvironment. However, more research is necessary to confirm the role of DNASE1L3 and KLKB1.


Sign in / Sign up

Export Citation Format

Share Document