The Landscape of the Tumor Microenvironment in Skin Cutaneous Melanoma Reveals a Prognostic and Immunotherapeutically Relevant Gene Signature

The tumorigenesis of skin cutaneous melanoma (SKCM) remains unclear. The tumor microenvironment (TME) is well known to play a vital role in the onset and progression of SKCM. However, the dynamic mechanisms of immune regulation are insufficient. We conducted a comprehensive analysis of immune cell infiltration in the TME. Based on the differentially expressed genes (DEGs) in clusters grouped by immune infiltration status, a set of hub genes related to the clinical prognosis of SKCM and tumor immune infiltration was explored.Methods: We analyzed immune cell infiltration in two independent cohorts and assessed the relationship between the internal pattern of immune cell infiltration and SKCM characteristics, including clinicopathological features, potential biological pathways, and gene mutations. Genes related to the infiltration pattern of TME immune cells were determined. Furthermore, the unsupervised clustering method (k-means) was used to divide samples into three different categories according to TME, which were defined as TME cluster-A, -B, and -C. DEGs among three groups of samples were analyzed as signature genes. We further distinguished common DEGs between three groups of samples according to whether differences were significant and divided DEGs into the Signature gene-A group with significant differences and the Signature gene-B group with insignificant differences. The Signature gene-A gene set mainly had exon skipping in SKCM, while the Signature gene-B gene set had no obvious alternative splicing form. Subsequently, we analyzed genetic variations of the two signatures and constructed a competing endogenous RNA (ceRNA) regulatory network. LASSO Cox regression was used to determine the immune infiltration signature and risk score of SKCM. Finally, we obtained 13 hub genes and calculated the risk score based on the coefficient of each gene to explore the impact of the high- and low-risk scores on biologically related functions and prognosis of SKCM patients further. The correlation between the risk score and clinicopathological characteristics of SKCM patients indicated that a low-risk score was associated with TME cluster-A classification (p < 0.001) and metastatic SKCM (p < 0.001). Thirteen hub genes also showed different prognostic effects in pan-cancer. The results of univariate and multivariate Cox analyses revealed that risk score could be used as an independent risk factor for predicting the prognosis of SKCM patients. The nomogram that integrated clinicopathological characteristics and immune characteristics to predict survival probability was based on multivariate Cox regression. Finally, 13 hub genes that showed different prognostic effects in pan-cancers were obtained. According to immunohistochemistry staining results, Ube2L6, SRPX2, and IFIT2 were expressed at higher levels, while CLEC4E, END3, and KIR2DL4 were expressed at lower levels in 25 melanoma specimens.Conclusion: We performed a comprehensive assessment of the immune-associated TME. To elucidate the potential development of immune-genomic features in SKCM, we constructed an unprecedented set of immune characteristic genes (EDN3, CLEC4E, SRPX2, KIR2DL4, UBE2L6, and IFIT2) related to the immune landscape of TME. These genes are related to different prognoses and drug responses of SKCM. The immune gene signature constructed can be used as a robust prognostic biomarker of SKCM and a predictor of an immunotherapy effect.

Download Full-text

Immuno-Genomic Landscape of Cutaneous Melanoma Identifies a Prognostic and Immunotherapeutically Relevant Gene Signature

10.21203/rs.3.rs-617208/v1 ◽

2021 ◽

Author(s):

Ronghua Yang ◽

Yidan Sun ◽

Tianqi Chen ◽

Jiehua Li ◽

Xiaobing Pi ◽

...

Keyword(s):

Risk Score ◽

Cutaneous Melanoma ◽

Immune Cell ◽

Gene Signature ◽

Low Risk ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Hub Genes ◽

Gene Set ◽

Immune Infiltration

Abstract BackgroundThe tumorigenesis of Skin cutaneous melanoma (SKCM) is still a mystery. Our study conducted a comprehensive analysis of the immune cell infiltration in the TME of SKCM. Based on the differential expression genes in the cluster grouped by the immune infiltration status, a set of hub genes related to the clinical prognosis of SKCM and tumor immune infiltration were explored.MethodsWe analyzed the immune cell infiltration in two independent cohorts, and then assessed the relationship between the internal pattern of immune cell infiltration and SKCM characteristics, including clinicopathological features, potential biological pathways and gene mutations. We further divided the three clusters of differential genes into two groups with different unique biological processes. The Signature gene-A gene set was mainly manifested as exon skipping (ES) in SKCM patients, while the Signature gene-B gene set has no obvious alternative splicing form. Subsequently, we not only analyzed the genetic variation of the two signatures, but also constructed a ceRNA regulatory network..LASSO Cox regression was utilized to find the immune infiltration signature and the risk score of SKCM. ResultWe finally obtained 13 Hub genes, and calculated the risk score based on the coefficient of each gene to further explore the impact of the high and low-risk score on the biologically related functions and prognosis of SKCM patients.The correlation between the risk score and the clinicopathological characteristics of SKCM patients indicated that the low risk score was associated with TMECluster-A classification (P <0.001) and metastatic SKCM (P <0.001). We finally obtained 13 Hub genes which showed different prognostic effects in pan-cancers. The IHC staining results showed that Ube2L6, SRPX2, IFIT2 were higher expression while CLEC4E, END3, KIR2DL4 were lower expression in 25 melanoma specimens.ConclusionWe performed a comprehensive assessment of SKCM's immune environment and constructed a set of unprecedented immune signatures related to the immune landscape (EDN3、CLEC4E、SRPX2、KIR2DL4、UBE2L6、IFIT2), which are correlated with the different prognosis and drug response of SKCM. The immune gene signature we constructed can be used as a robust prognostic biomarker of SKCM and a predictor of immunotherapy effect.

Download Full-text

A Ferroptosis-Related Signature Robustly Predicts Clinical Outcomes and Associates With Immune Microenvironment for Thyroid Cancer

Frontiers in Medicine ◽

10.3389/fmed.2021.637743 ◽

2021 ◽

Vol 8 ◽

Author(s):

Mingqin Ge ◽

Jie Niu ◽

Ping Hu ◽

Aihua Tong ◽

Yan Dai ◽

...

Keyword(s):

Thyroid Cancer ◽

High Risk ◽

Risk Score ◽

Immune Cell ◽

Cox Regression ◽

Low Risk ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Immune Microenvironment ◽

Cox Regression Analysis

Objective: This study aimed to construct a prognostic ferroptosis-related signature for thyroid cancer and probe into the association with tumor immune microenvironment.Methods: Based on the expression profiles of ferroptosis-related genes, a LASSO cox regression model was established for thyroid cancer. Kaplan-Meier survival analysis was presented between high and low risk groups. The predictive performance was assessed by ROC. The predictive independency was validated via multivariate cox regression analysis and stratified analysis. A nomogram was established and verified by calibration curves. The enriched signaling pathways were predicted via GSEA. The association between the signature and immune cell infiltration was analyzed by CIBERSORT. The ferroptosis-related genes were validated in thyroid cancer tissues by immunohistochemistry and RT-qPCR.Results: A ferroptosis-related eight gene model was established for predicting the prognosis of thyroid cancer. Patients with high risk score indicated a poorer prognosis than those with low risk score (p = 1.186e-03). The AUCs for 1-, 2-, and 3-year survival were 0.887, 0.890, and 0.840, respectively. Following adjusting other prognostic factors, the model could independently predict the prognosis (p = 0.015, HR: 1.870, 95%CI: 1.132–3.090). A nomogram combining the signature and age was constructed. The nomogram-predicted probability of 1-, 3-, and 5-year survival approached the actual survival time. Several ferroptosis-related pathways were enriched in the high-risk group. The signature was distinctly associated with the immune cell infiltration. After validation, the eight genes were abnormally expressed between thyroid cancer and control tissues.Conclusion: Our findings established a prognostic ferroptosis-related signature that was associated with the immune microenvironment for thyroid cancer.

Download Full-text

Development of an Immune-Related LncRNA Prognostic Signature for Glioma

Frontiers in Genetics ◽

10.3389/fgene.2021.678436 ◽

2021 ◽

Vol 12 ◽

Author(s):

Yudong Cao ◽

Hecheng Zhu ◽

Jun Tan ◽

Wen Yin ◽

Quanwei Zhou ◽

...

Keyword(s):

Regression Analysis ◽

Clinical Outcome ◽

Immune Cell ◽

Cox Regression ◽

Gene Set Enrichment Analysis ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Prognostic Signature ◽

Cox Regression Analysis ◽

Immune Infiltration

IntroductionGlioma is the most common primary cancer of the central nervous system with dismal prognosis. Long noncoding RNAs (lncRNAs) have been discovered to play key roles in tumorigenesis in various cancers, including glioma. Because of the relevance between immune infiltrating and clinical outcome of glioma, identifying immune-related lncRNAs is urgent for better personalized management.Materials and methodsSingle-sample gene set enrichment analysis (ssGSEA) was applied to estimate immune infiltration, and glioma samples were divided into high immune cell infiltration group and low immune cell infiltration group. After screening differentially expressed lncRNAs in two immune groups, least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to construct an immune-related prognostic signature. Additionally, we explored the correlation between immune infiltration and the prognostic signature.ResultsA total of 653 samples were appropriate for further analyses, and 10 lncRNAs were identified as immune-related lncRNAs in glioma. After univariate Cox regression and LASSO Cox regression analysis, six lncRNAs were identified to construct a prognostic signature for glioma, which could be taken as independent prognostic factors in both univariate and multivariate Cox regression analyses. Moreover, risk score was significantly correlated with all the 29 immune-related checkpoint expression (p < 0.05) in ssGSEA except neutrophils (p = 0.43).ConclusionThe study constructed an immune-related prognostic signature for glioma, which contributed to improve clinical outcome prediction and guide immunotherapy.

Download Full-text

Identification of Crucial Methylation Genes Associated With Immune Cell Infiltration in Twin-to-twin Transfusion Syndrome Using Co-methylation Analysis

10.21203/rs.3.rs-210035/v1 ◽

2021 ◽

Author(s):

Chuang Li ◽

Yuan Wang ◽

Caixia Liu ◽

Shaowei Yin

Keyword(s):

Gene Expression ◽

Immune Cell ◽

Interaction Network ◽

Cell Infiltration ◽

Future Research ◽

Regulation Mechanism ◽

Immune Cell Infiltration ◽

Hub Genes ◽

Immune Infiltration ◽

Neurodevelopmental Abnormalities

Abstract Background: DNA methylation (DNAm), is an important transcriptional regulation mechanism, relevant to various diseases. Twin-to-twin transfusion syndrome (TTTS) is a complication in twin pregnancies resulting from disproportionate blood circulation. Survivors of TTTS show a high risk of neurodevelopmental abnormalities, particularly in the hippocampus, which is important in learning and memory. Here, we investigate gene expression and DNAm in hippocampus tissues of TTTS specimens. Methods: DNAm and gene expression levels were compared among the three groups: 10 recipients, 10 donors, and 10 matched control, using methylation microarray. We further explored the immune infiltration of six immune cell sub-populations using EpiDISH analysis. The methylated sites related to immune cell infiltration were identified using the WGCNA package. We explored the core methylation genes in the protein-protein interaction network using the MCODE plugin in Cytoscape software. Results: There were 188 differential methylation sites among three groups. Based on WGCNA, we found that the turquoise module containing 174 CpG sites is significantly related to the immune infiltration level. And four hub genes correlated with immune infiltration level, namely, PTPRJ, FYN, LYN, and AKT1, and were identified using gene sub-network analysis. Conclusions: We identify the four hub methylation genes related to immune infiltration in the TTTS. The molecular function of hub genes is still explored in the future research.

Download Full-text

Identification and Verification of Immune-Related Gene Prognostic Signature Based on ssGSEA for Osteosarcoma

Frontiers in Oncology ◽

10.3389/fonc.2020.607622 ◽

2020 ◽

Vol 10 ◽

Author(s):

Bo Xiao ◽

Liyan Liu ◽

Aoyu Li ◽

Cheng Xiang ◽

Pingxiao Wang ◽

...

Keyword(s):

Immune Cell ◽

Cox Regression ◽

Gene Signature ◽

Cell Infiltration ◽

Significant Prognostic Factor ◽

Immune Cell Infiltration ◽

Related Gene ◽

Prognostic Signature ◽

Immune Related Gene ◽

Immune Related Genes

Osteosarcoma is the most common malignant bone tumor in children and adolescence. Multiple immune-related genes have been reported in different cancers. The aim is to identify an immune-related gene signature for the prospective evaluation of prognosis for osteosarcoma patients. In this study, we evaluated the infiltration of immune cells in 101 osteosarcoma patients downloaded from TARGET using the ssGSEA to the RNA-sequencing of these patients, thus, high immune cell infiltration cluster, middle immune cell infiltration cluster and low immune cell infiltration cluster were generated. On the foundation of high immune cell infiltration cluster vs. low immune cell infiltration cluster and normal vs. osteosarcoma, we found 108 common differentially expressed genes which were sequentially submitted to univariate Cox and LASSO regression analysis. Furthermore, GSEA indicated some pathways with notable enrichment in the high- and low-immune cell infiltration cluster that may be helpful in understanding the potential mechanisms. Finally, we identified seven immune-related genes as prognostic signature for osteosarcoma. Kaplan-Meier analysis, ROC curve, univariate and multivariate Cox regression further confirmed that the seven immune-related genes signature was an innovative and significant prognostic factor independent of clinical features. These results of this study offer a means to predict the prognosis and survival of osteosarcoma patients with uncovered seven-gene signature as potential biomarkers.

Download Full-text

Identification of Hub Genes and Their Correlation With Immune Infiltration Cells in Hepatocellular Carcinoma Based on GEO and TCGA Databases

Frontiers in Genetics ◽

10.3389/fgene.2021.647353 ◽

2021 ◽

Vol 12 ◽

Author(s):

Rui Huang ◽

Jinying Liu ◽

Hui Li ◽

Lierui Zheng ◽

Haojun Jin ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Cell ◽

Primary Liver Cancer ◽

Cell Infiltration ◽

Regulation Mechanism ◽

Immune Cell Infiltration ◽

Hub Genes ◽

Rank Analysis ◽

Immune Infiltration ◽

Follicular Helper

Hepatocellular carcinoma (HCC) is a primary liver cancer with extremely high mortality in worldwide. HCC is hard to diagnose and has a poor prognosis due to the less understanding of the molecular pathological mechanisms and the regulation mechanism on immune cell infiltration during hepatocarcinogenesis. Herein, by performing multiple bioinformatics analysis methods, including the RobustRankAggreg (RRA) rank analysis, weighted gene co-expression network analysis (WGCNA), and a devolution algorithm (CIBERSORT), we first identified 14 hub genes (NDC80, DLGAP5, BUB1B, KIF20A, KIF2C, KIF11, NCAPG, NUSAP1, PBK, ASPM, FOXM1, TPX2, UBE2C, and PRC1) in HCC, whose expression levels were significantly up-regulated and negatively correlated with overall survival time. Moreover, we found that the expression of these hub genes was significantly positively correlated with immune infiltration cells, including regulatory T cells (Treg), T follicular helper (TFH) cells, macrophages M0, but negatively correlated with immune infiltration cells including monocytes. Among these hub genes, KIF2C and UBE2C showed the most significant correlation and were associated with immune cell infiltration in HCC, which was speculated as the potential prognostic biomarker for guiding immunotherapy.

Download Full-text

A Comprehensive Prognostic and Immunological Analysis of a New Three-Gene Signature in Hepatocellular Carcinoma

Stem Cells International ◽

10.1155/2021/5546032 ◽

2021 ◽

Vol 2021 ◽

pp. 1-25

Author(s):

Jun Liu ◽

Jianjun Lu ◽

Wenli Li

Keyword(s):

Hepatocellular Carcinoma ◽

Risk Score ◽

Immune Checkpoint ◽

Immune Cell ◽

Prognostic Significance ◽

Gene Signature ◽

Cancer Genome ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Immune Checkpoint Molecules

There are few reports on the role of genes associated with the mRNA expression-based stemness index (mRNAsi) in the prognosis and immune regulation of hepatocellular carcinoma (HCC). This study is aimed at analyzing the expression profile and prognostic significance of a new mRNAsi-based three-gene signature in HCC. This three-gene signature was identified by analyzing mRNAsi data from the Cancer Genome Atlas (TCGA) HCC dataset. The prognostic value of the risk score based on the three-gene signature was evaluated by Cox regression and Kaplan-Meier analysis and then verified in the International Cancer Genome Consortium (ICGC) database. Meanwhile, the correlations between the risk score and immune cell infiltration patterns, microsatellite instability (MSI), tumor mutation burden (TMB), immune checkpoint molecules, hypoxia-related genes, immunotherapy response, and compounds targeting the gene signature were explored, respectively. The results showed that compared with normal liver tissues, the mRNAsi score of HCC tissues was significantly increased. PTDSS2, MRPL9, and SOCS were the genes most related to mRNAsi in HCC tissues. Survival analysis results suggested the risk score based on the three-gene signature was an independent predictor of the prognosis for patients with HCC. The nomogram combining the risk score and pathological stage showed a good predictive ability for the overall survival of patients with HCC patients. Meanwhile, the risk score was significantly related to immune cell infiltration patterns, MSI, TMB, several immune checkpoint molecules, and hypoxia-related genes. In addition, the risk score was associated with the immunotherapy response, and fifteen potential therapeutic drugs targeting the three-gene signature were identified. Therefore, we propose to use this three-gene signature including PTDSS2, MRPL9, and SOCS as a potential prognostic biomarker for HCC.

Download Full-text

The Landscape of Immune Cell Infiltration in the Glomerulus of Diabetic Nephropathy: Evidence Based on Bioinformatics

10.21203/rs.3.rs-554679/v1 ◽

2021 ◽

Author(s):

Wei ZHOU ◽

Yaoyu LIU ◽

Qinghong HU ◽

Jiuyao ZHOU ◽

Hua LIN

Keyword(s):

T Cells ◽

Diabetic Nephropathy ◽

Mast Cells ◽

Protein Interactions ◽

Immune Cell ◽

Gene Set Enrichment Analysis ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Hub Genes ◽

Immune Infiltration

Abstract BackgroundIncreasing evidence suggests that immune cell infiltration contributes to the pathogenesis and progression of diabetic nephropathy (DN). We aim to unveil the immune infiltration pattern in the glomerulus of DN and provide potential targets for immunotherapy. MethodsInfiltrating percentage of 22 types of immune cell in the glomerulus tissues were estimated by the CIBERSORT algorithm based on three transcriptome datasets mined from the GEO database. Differentially expressed genes (DEGs) were identified by the “limma” package. Then immune-related DEGs were identified by intersecting DEGs with immune-related genes (downloaded from Immport database). The protein-protein interactions of Immune-related DEGs were explored using the STRING database and visualized by Cytoscape. The enrichment analyses for KEGG pathways and GO terms were carried out by the gene set enrichment analysis (GSEA) method. Results9 types of immune cell were revealed to be significantly altered in the glomerulus tissues of DN (Up: B cells memory, T cells CD4 naive, Macrophages M2, Dendritic cells resting, Mast cells resting, Mast cells activated; Down: NK cells resting, Monocytes, Neutrophils). Correlation analysis revealed that immune infiltration act as a complicated and tightly regulated network, among which T cells gamma delta and T cells CD4 naive show the most synergistic effect (r = 0.58, p ＜ 0.001); meanwhile, T cells CD8 and T cells CD4 memory resting show the most competitive effect (r = - 0.67, p ＜ 0.001). Several pathways related to immune were significantly activated. Moreover, 6 hub genes with a medium to strong correlation with renal function (eGFR) were identified (ALB, EGF, FOS, CXCR1, CXCR2, CCL2). ConclusionIn the glomerulus of DN, the immune infiltration pattern changed significantly. A complicated and tightly regulated network of immune cells exists in the pathological of DN. The hub genes identified here will facilitate the development of immunotherapy.

Download Full-text

PD-L1 and Immune Infiltration of m6A RNA Methylation Regulators and Its miRNA Regulators in Hepatocellular Carcinoma

BioMed Research International ◽

10.1155/2021/5516100 ◽

2021 ◽

Vol 2021 ◽

pp. 1-16

Author(s):

Yingxue Lin ◽

Yinhui Yao ◽

Ying Wang ◽

Lingdi Wang ◽

Haipeng Cui

Keyword(s):

Risk Score ◽

Immune Cell ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Consensus Clustering ◽

Lasso Regression ◽

Rna Methylation ◽

Immune Infiltration ◽

Significant Difference ◽

M6a Rna Methylation

Background. The aim of this study was to systematically evaluate the relationship between the expression of m6A RNA methylation regulators and prognosis in HCC. Methods. We compared the expression of m6A methylation modulators and PD-L1 between HCC and normal in TCGA database. HCC samples were divided into two subtypes by consensus clustering of data from m6A RNA methylation regulators. The differences in PD-L1, immune infiltration, and prognosis between the two subtypes were further compared. The LASSO regression was used to build a risk score for m6A modulators. In addition, we identified miRNAs that regulate m6A regulators. Results. We found that fourteen m6A regulatory genes were significantly differentially expressed between HCC and normal. HCC samples were divided into two clusters. Of these, there are higher PD-L1 expression and poorer overall survival (OS) in cluster 1. There was a significant difference in immune cell infiltration between cluster 1 and cluster 2. Through the LASSO model, we obtained 12 m6A methylation regulators to construct a prognostic risk score. Compared with patients with a high-risk score, patients with a low-risk score had upregulated PD-L1 expression and worse prognosis. There was a significant correlation between risk score and tumor-infiltrating immune cells. Finally, we found that miR-142 may be the important regulator for m6A RNA methylation in HCC. Conclusion. Our results suggest that m6A RNA methylation modulators may affect the prognosis through PD-L1 and immune cell infiltration in HCC patients. In addition, the two clusters may be beneficial for prognostic stratification and improving immunotherapeutic efficacy.

Download Full-text

Development of a miR-3648-Targeted Immune-Related 4-Gene Signature as a Prognostic Biomarker in Esophageal Adenocarcinoma Based on ehe WGCNA Algorithm

10.21203/rs.3.rs-549122/v1 ◽

2021 ◽

Author(s):

Donglei Zhang ◽

Hang Yin ◽

Ping Xu ◽

Xiaozhe Qian

Keyword(s):

Esophageal Adenocarcinoma ◽

Immune Cell ◽

Cox Regression ◽

Gene Signature ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Cell Infiltration ◽

Prognostic Signature ◽

Cox Regression Analysis ◽

Normal Tissues

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.

Download Full-text