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

2021 ◽  
Vol 12 ◽  
Author(s):  
Rui Huang ◽  
Jinying Liu ◽  
Hui Li ◽  
Lierui Zheng ◽  
Haojun Jin ◽  
...  

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.

2021 ◽  
Author(s):  
Chuang Li ◽  
Yuan Wang ◽  
Caixia Liu ◽  
Shaowei Yin

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.


Author(s):  
Sitong Zhou ◽  
Yidan Sun ◽  
Tianqi Chen ◽  
Jingru Wang ◽  
Jia He ◽  
...  

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.


2021 ◽  
Vol 11 ◽  
Author(s):  
Xingwei Xie ◽  
Shanshan Jiang ◽  
Xiang Li

Nuf2 participates in the regulation of cell apoptosis and proliferation by regulating the binding of centromere and spindle microtubules to achieve the correct separation of chromosomes. Previous reports have suggested that Nuf2 may play a role in various human cancers. However, the mechanism and function of Nuf2 in the development of Hepatocellular carcinoma (HCC) remains uncertain. This study investigated the prognostic potential of Nuf2 and its relation with immune cell infiltration in HCC. Nuf2 expression in tumor cells was examined using the TIMER and Oncomine databases, and its prognostic potential was assessed via the Kaplan-Meier plotter and GEPIA databases. The relationships between Nuf2 and tumor immune infiltration were analyzed using TIMER. The relationships between Nuf2 and biomarkers of tumor immune infiltration were analyzed using TIMER and GEPIA. Here we revealed that Nuf2 expression increased in tumor tissues containing HCC, and this correlated with poor relapse-free survival, disease-specific survival, progression-free survival, and overall survival in patients with HCC regardless of grades, genders, races, drinking behaviors and other clinical factors. Additionally, high expression of Nuf2 was positively correlated with differential immune cell infiltration and various immune biomarkers. Our work demonstrated that Nuf2 could be a potential prognostic biomarker and could be related to tumor immune cell infiltration in HCC.


Hereditas ◽  
2021 ◽  
Vol 158 (1) ◽  
Author(s):  
Yiqi Li ◽  
Jue Qi ◽  
Jiankang Yang

Abstract Objective Melanoma accounts for 80% of skin cancer deaths. The pathogenesis of melanoma is regulated by gene networks. Thus, we aimed here to identify gene networks and hub genes associated with melanoma and to further identify their underlying mechanisms. Methods GTEx (normal skin) and TCGA (melanoma tumor) RNA-seq datasets were employed for this purpose. We conducted weighted gene co-expression network analysis (WGCNA) to identify key modules and hub genes associated with melanoma. Log-rank analysis and multivariate Cox model analysis were performed to identify prognosis genes, which were validated using two independent melanoma datasets. We also evaluated the correlation between prognostic gene and immune cell infiltration. Results The blue module was the most relevant for melanoma and was thus considered the key module. Intersecting genes were identified between this module and differentially expressed genes (DEGs). Finally, 72 genes were identified and verified as hub genes using the Oncomine database. Log-rank analysis and multivariate Cox model analysis identified 13 genes that were associated with the prognosis of the metastatic melanoma group, and RTP4 was validated as a prognostic gene using two independent melanoma datasets. RTP4 was not previously associated with melanoma. When we evaluated the correlation between prognostic gene and immune cell infiltration, we discovered that RTP4 was associated with immune cell infiltration. Further, RTP4 was significantly associated with genes encoding components of immune checkpoints (PDCD1, TIM-3, and LAG3). Conclusions RTP4 is a novel prognosis-related hub gene in cutaneous melanoma. The novel gene RTP4 identified here will facilitate the exploration of the molecular mechanism of the pathogenesis and progression of melanoma and the discovery of potential new target for drug therapy.


2021 ◽  
Author(s):  
Ronghua Yang ◽  
Yidan Sun ◽  
Tianqi Chen ◽  
Jiehua Li ◽  
Xiaobing Pi ◽  
...  

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.


2021 ◽  
Author(s):  
Wei ZHOU ◽  
Yaoyu LIU ◽  
Qinghong HU ◽  
Jiuyao ZHOU ◽  
Hua LIN

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.


2021 ◽  
Author(s):  
shenglan li ◽  
Zhuang Kang ◽  
jinyi Chen ◽  
Can Wang ◽  
Zehao Cai ◽  
...  

Abstract Background Medulloblastoma is a common intracranial tumor among children. In recent years, research on cancer genome has established four distinct subtypes of medulloblastoma: WNT, SHH, Group3, and Group4. Each subtype has its own transcriptional profile, methylation changes, and different clinical outcomes. Treatment and prognosis also vary depending on the subtype. Methods Based on the methylation data of medulloblastoma samples, methylCIBERSORT was used to evaluate the level of immune cell infiltration in medulloblastoma samples and identified 10 kinds of immune cells with different subtypes. Combined with the immune database, 293 Imm-DEGs were screened. Imm-DEGs were used to construct the co-expression network, and the key modules related to the level of differential immune cell infiltration were identified. Three immune hub genes (GAB1, ABL1, CXCR4) were identified according to the gene connectivity and the correlation with phenotype in the key modules, as well as the PPI network involved in the genes in the modules. Results The subtype marker was recognized according to the immune hub, and the subtype marker was verified in the external data set, the methylation level of immune hub gene among different subtypes was compared and analyzed, at the same time, tissue microarray was used for immunohistochemical verification, and a multi-factor regulatory network of hub gene was constructed. Conclusions Identifying subtype marker is helpful to accurately identify the subtypes of medulloblastoma patients, and can accurately evaluate the treatment and prognosis, so as to improve the overall survival of patients.


2021 ◽  
Author(s):  
Di Cao ◽  
Jun Wang ◽  
Yan Lin ◽  
Guangwei Li

Abstract Background: The therapeutic efficacy of immune checkpoint inhibitor therapy is highly influenced by tumor mutation burden (TMB). The relationship between TMB and prognosis in lower-grade glioma is still unclear. We aimed to explore the relationships and mechanisms between them in lower-grade glioma.Methods: We leveraged somatic mutation data from The Cancer Genome Atlas (TCGA). Clinical cases were divided into high- and low-TMB groups based on the median of TMB. Infiltrating immune cells were analyzed using CIBERSORT and differential expression analysis between the prognostic groups performed. The key genes were identified as intersecting between immune-related genes. Cox regression and survival analysis were performed on the intersecting genes. A total of 7 hub genes were identified. The effect of somatic copy number alterations (SCNA) of the hub genes on immune cell infiltration was analyzed using TIMER, which was used to determine the risk factors and immune infiltration status in LGG. Subsequently, based on hub genes, a TMB Prognosis Index (TMBPI) model was constructed to predict the risk in LGG patients. Besides, this model was validated using data from TCGA and Chinese Glioma Genome Atlas (CGGA).Results: High-TMB favored worse prognosis (P<0.001) and macrophage infiltration was an independent risk factor (P<0.001). In the high-TMB group (P=0.033, P=0.009), the proportion of macrophages M0 and M2 increased and monocytes decreased (P=0.006). Besides, the SCNA of the hub genes affected the level of immune cell infiltration by varying degrees among which IGF2BP3, NPNT, and PLA2G2A had a significant impact. The AUC of the ROC curve at 1-, 3- and 5-years were all above 0.74.Conclusions: This study implies that high-TMB correlated with unfavorable prognosis in lower-grade glioma. And high-TMB may have an impact on prognosis by changing tumor microenvironment, caused by the SCNAs of genes. The TMBPI model accurately predicted prognosis in LGG patients.


Sign in / Sign up

Export Citation Format

Share Document