scholarly journals Dissecting the Role of N6-Methylandenosine-Related Long Non-coding RNAs Signature in Prognosis and Immune Microenvironment of Breast Cancer

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):  
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.


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):  
Guanghui Zhao ◽  
Junhua An ◽  
Qian Pu ◽  
Wenwen Geng ◽  
Haiyun Song ◽  
...  

The N6-methyladenosine (m6A) has been considered as a new layer of epitranscriptomic regulation on mRNA processing, stability, and translation. However, potential roles of m6A RNA methylation modification in tumor immune microenvironment (TIME) of breast cancer are yet fully understood. In this study, we comprehensively evaluated the genetic variations and transcript expressions of 15 m6A regulators in 1,079 breast cancer samples from the Cancer Genome Atlas (TCGA) database. We validated major regulators had significantly differential mRNA and protein expression in tumor tissue compared to normal tissues from 39 pairs of clinical breast cancer samples with different molecular subtypes, and especially high expression of m6A readers YTHDF1 and YTHDF3 predicted poor survival. Two clusters of breast cancer patients identified by the 15 m6A regulators’ pattern showed distinct overall survival, immune activation status, and immune cell infiltration, and clinical samples confirmed the diversity of lymphocytic infiltration. The profiles of these two clusters accorded with that of two classical cancer-immune phenotypes, immune-excluded and immune-inflamed phenotypes, it suggested that m6A regulators-based patterns might serve as crucial mediators of TIME in breast cancer. Moreover, the m6A phenotype-related gene signatures could also be survival predictor in breast cancer. Therefore, comprehensive evaluation of tumor m6A modification pattern will contribute to enhance our understanding of the characterization of immune cell infiltration in the tumor microenvironment and promote the responsiveness of breast cancer to immunotherapy.


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 21 (1) ◽  
Author(s):  
Hao Xu ◽  
Qianhui Xu ◽  
Lu Yin

Abstract Background Although immunotherapy for colon cancer has made promising progress, only a few patients currently benefit from it. A recent study revealed that infiltrating immune cells are highly relevant to tumor prognosis and influence the expression of immune-related genes. However, the characterization of immune cell infiltration (ICI) has not yet been comprehensively analyzed and quantified in colon adenocarcinoma (COAD). Methods The multiomic data of COAD samples were downloaded from TCGA. ESTIMATE algorithm, ssGSEA method and CIBERSORT analysis were conducted to estimate the subpopulations of infiltrating immune cells. COAD subtypes based on ICI pattern were identified by consensus clustering then principal-component analysis was performed to obtain ICI scores to quantify the ICI patterns in individual tumors. Kaplan–Meier analysis was employed to validate prognostic value. Gene set enrichment analysis (GSEA) was applied for functional annotation. Finally, the mutation data was analyzed by employing “maftools” package. Results Three bioinformatics algorithms were used to evaluate the ICI patterns from 538 patients with COAD. Two ICI subtypes were determined using consensus clustering, and the ICI score was constructed by performing principal component analysis. Our findings showed that a higher ICI score often indicated a more advanced tumor and worse prognosis. The high-ICI score subgroup had a higher stromal score and more M0 macrophages but fewer plasma cells and decreased CD8 T cell infiltration. In addition, patients with high ICI scores had significantly higher expression levels of HAVCR2 and PCDC1LG2. Real-time polymerase chain reaction (PCR) was conducted to determine the prognostic significances of ICI-related genes. Conclusions In conclusion, ICI score may be considered as an original and useful indicator for independent prognostic prediction and individual immune-related therapy.


Author(s):  
Qi Zhao ◽  
Junfeng Liu

Objective: Prolyl 4-hydroxylase, alpha polypeptide I (P4HA1), a key enzyme in collagen synthesis, comprises two identical alpha subunits and two beta subunits. However, the immunomodulatory role of P4HA1 in tumor immune microenvironment (TIME) remains unclear. This study aimed to evaluate the prognostic value of P4HA1 in pan-cancer and explore the relationship between P4HA1 expression and TIME.Methods: P4HA1 expression, clinical features, mutations, DNA methylation, copy number alteration, and prognostic value in pan-cancer were investigated using the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression data. Pathway enrichment analysis of P4HA1 was performed using R package “clusterProfiler.” The correlation between immune cell infiltration level and P4HA1 expression was analyzed using three sources of immune cell infiltration data, including ImmuCellAI database, TIMER2 database, and a published work.Results: P4HA1 was substantially overexpressed in most cancer types. P4HA1 overexpression was associated with poor survival in patients. Additionally, we discovered that P4HA1 expression was positively associated with infiltration levels of immunosuppressive cells, such as tumor-associated macrophages, cancer-associated fibroblasts, nTregs, and iTregs, and negatively correlated with CD8+ T and NK cells in pan-cancer.Conclusions: Our results highlighted that P4HA1 might serve as a potential prognostic biomarker in pan-cancer. P4HA1 overexpression is indicative of an immunosuppressive microenvironment. P4HA1 may be a potential target of immunotherapy.


2021 ◽  
Author(s):  
Donglei Zhang ◽  
Hang Yin ◽  
Ping Xu ◽  
Xiaozhe Qian

Abstract Background: Esophageal adenocarcinoma (EA) has a poor prognosis and is a typical immunogenic malignant tumor. Abnormal expression of miR-3648 has been reported in EA, but its value in prognosis prediction and immune cell infiltration imbalance mediation is still unknown. We aimed to mine immune-related genes (IRGs) targeted by miR-3648 and construct a multigene signature to improve the prognostic prediction of EA.Methods: The gene expression data of EA tumor or normal tissues from The Cancer Genome Atlas (TCGA) database and GTEx database were downloaded. Weighted gene coexpression network analysis (WGCNA), the CIBERSORT algorithm and Cox regression analysis were applied to identify IRGs and to construct a prognostic signature and nomogram.Results: miR-3648 was obviously highly expressed in EA tumor tissues and was correlated with patient survival time [hazard ratio (HR) = 1.28, 95% confidence interval (CI): 1.09-1.49, p = 0.002]. A total of 70 miR-3648-targeted genes related to immune cell infiltration were identified, and a novel 4-gene signature (C10orf55, DLL4, PANX2 and NKAIN1) was established. The prognostic signature-based risk score has superior capability to predict overall survival (AUC = 0.740 for 1 year; AUC = 0.717 for 3 years; AUC = 0.622 for 5 years). A higher score was indicative of a poorer prognosis than a lower score (HR = 1.69, 95% CI: 1.08-2.64, p = 0.20, adjusted by TNM stage).Conclusion: miR-3648 might play a crucial role in the progression of EA. The novel miR-3648-targeted immune-related 4-gene signature is expected to become a potential prognostic marker in EA.


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


Sign in / Sign up

Export Citation Format

Share Document