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

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.

2021 ◽  
Author(s):  
Qi Zhou ◽  
Xin Xiong ◽  
Min Tang ◽  
Yingqing Lei ◽  
Hongbin Lv

Abstract BackgroundDiabetic retinopathy (DR), a severe complication of diabetes mellitus (DM), is a global social and economic burden. However, the pathological mechanisms mediating DR are not well-understood. This study aimed to identify differentially methylated and differentially expressed hub genes (DMGs and DEGs, respectively) and associated signaling pathways, and to evaluate immune cell infiltration involved in DR. MethodsTwo publicly available datasets were downloaded from the Gene Expression Omnibus database. Transcriptome and epigenome microarray data and multi-component weighted gene coexpression network analysis (WGCNA) were utilized to determine hub genes within DR. One dataset was utilized to screen DEGs and to further explore their potential biological functions using functional annotation analysis. A protein-protein interaction network was constructed. Gene set enrichment and variation analyses (GSVA and GSEA, respectively) were utilized to identify the potential mechanisms mediating the function of hub genes in DR. Infiltrating immune cells were evaluated in one dataset using CIBERSORT. The Connectivity Map (CMap) database was used to predict potential therapeutic agents. ResultsIn total, 673 DEGs (151 upregulated and 522 downregulated genes) were detected. Gene expression was significantly enriched in the extracellular matrix and sensory organ development, extracellular matrix organization, and glial cell differentiation pathways. Through WGCNA, one module was found to be significantly related with DR (r=0.34, P =0.002), and 979 hub genes were identified. By comparing DMGs, DEGs, and genes in WGCNA, we identified eight hub genes in DR ( AKAP13, BOC, ACSS1, ARNT2, TGFB2, LHFPL2, GFPT2, TNFRSF1A ), which were significantly enriched in critical pathways involving coagulation, angiogenesis, TGF-β, and TNF-α-NF-κB signaling via GSVA and GSEA. Immune cell infiltration analysis revealed that activated natural killer cells, M0 macrophages, resting mast cells, and CD8 + T cells may be involved in DR. ARNT2, TGFB2, LHFPL2 , and AKAP13 expression were correlated with immune cell processes, and ZG-10, JNK-9L, chromomycin-a3, and calyculin were identified as potential drugs against DR. Finally, TNFRSF1A , GFPT2 , and LHFPL2 expression levels were consistent with the bioinformatic analysis. ConclusionsOur results are informative with respect to correlations between differentially methylated and expressed hub genes and immune cell infiltration in DR, providing new insight towards DR drug development and treatment.


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.


2020 ◽  
Author(s):  
Biao Huang ◽  
Wei Han ◽  
Zu-Feng Sheng ◽  
Guo-Liang Shen

Abstract Background Skin cutaneous melanoma (SKCM) is known as the most malignancy and treatment-resistant in human tumor, causing about 72% of deaths in skin carcinoma. However, the potential mechanism and new effective targets remain to be further elucidated. Available datasets such as Gene Expression Omnibus (GEO) can be utilized to search for novel therapeutic targets and prognostic biomarkers. Methods Three data sets were downloaded from GEO database . The differentially expressed genes (DEGs) were identified via Venn software. Protein‐protein interaction network of DEGs was developed and the module hub genes analysis was constructed by Cytoscape. Subsequently, multiple online tools and Kaplan-Meier survival curves were analyzed to detect underlying signaling pathways, gene expression, drug-gene interaction and prognostic value of hub genes. In addition, we explored the correlation between hub genes and immune cell infiltration. At last, the related miRNA, lncRNA networks were constructed by R software. Results A total of 308 DEGs and 12 hub genes were identified. Function and pathway enrichment results demonstrated a correlation between DEGs and the tumor microenvironment, immune response and melanoma tumorigenesis. Subsequently, we focused on assessing potential value of 12 hub genes. Seven hub genes ( CCL4, CCL5, NMU, GAL, CXCL9, CXCL10, CXCL13 ) were identified with significant overall survival for prognosis. What’s more, five of these seven hub genes were found to be related to clinical stages (P values<0.05). In addition, the most important pathways of hub genes include interleukin-10 signaling, peptide ligand-binding receptors, which play important roles in tumor microenvironment for immune activation or immunosuppressive by regulating the infiltration of immune cells. Our results revealed a strong positive correlation between gene expression (CCL4, CCL5, CXCL9, CXCL10 and CXCL13) and immune cell infiltration (B-cell, CD8+ T cells, CD4+ T cells, macrophages, Neutrophils, Dendritic cells). Interestingly, 8 of 12 hub genes (CXCL10, CCL4, CCL5, IL6, CXCL2, PTGER3, GAL, NPY1R) were also found in the predicted drug-gene interaction. The related miRNA, lncRNA for diagnosis and prognosis were found in networks. Conclusion In conclusion, CCL4, CCL5, NMU, GAL, CXCL9, CXCL10, CXCL13 were of high prognostic value and may be potential targets for the diagnosis and therapy of patients with melanoma.


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 &lt; 0.001) and metastatic SKCM (p &lt; 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 8 ◽  
Author(s):  
Yujia Yang ◽  
Yue Cai ◽  
Yuan Zhang ◽  
Xu Yi ◽  
Zhiqiang Xu

Atherosclerotic cardiovascular disease (ASCVD) caused by atherosclerosis (AS) is one of the highest causes of mortality worldwide. Although there have been many studies on AS, its etiology remains unclear. In order to carry out molecular characterization of different types of AS, we retrieved two datasets composed of 151 AS samples and 32 normal samples from the Gene Expression Omnibus database. Using the non-negative matrix factorization (NMF) algorithm, we successfully divided the 151 AS samples into two subgroups. We then compared the molecular characteristics between the two groups using weighted gene co-expression analysis (WGCNA) and identified six key modules associated with the two subgroups. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analysis were used to identify the potential functions and pathways associated with the modules. In addition, we used the cytoscape software to construct and visualize protein–protein networks so as to identify key genes in the modules of interest. Three hub genes including PTGER3, GNAI1, and IGFBP5 were further screened using the least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE) algorithms. Since the modules were associated with immune pathways, we performed immune cell infiltration analysis. We discovered a significant difference in the level of immune cell infiltration by naïve B cells, CD8 T cells, T regulatory cells (Tregs), resting NK cells, Monocytes, Macrophages M0, Macrophages M1, and Macrophages M2 between the two subgroups. In addition, we observed the three hub genes were positively correlated with Tregs but negatively correlated with Macrophages M0. We also found that the three key genes are differentially expressed between normal and diseased tissue, as well as in the different subgroups. Receiver operating characteristic (ROC) results showed a good performance in the validation dataset. These results may provide novel insight into cellular and molecular characteristics of AS and potential markers for diagnosis and targeted 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 ◽  
Vol 11 (1) ◽  
Author(s):  
Zitong Feng ◽  
Jingge Qu ◽  
Xiao Liu ◽  
Jinghui Liang ◽  
Yongmeng Li ◽  
...  

AbstractEsophageal squamous cell carcinoma (ESCC) is a life-threatening thoracic tumor with a poor prognosis. The role of molecular alterations and the immune microenvironment in ESCC development has not been fully elucidated. The present study aimed to elucidate key candidate genes and immune cell infiltration characteristics in ESCC by integrated bioinformatics analysis. Nine gene expression datasets from the Gene Expression Omnibus (GEO) database were analysed to identify robust differentially expressed genes (DEGs) using the robust rank aggregation (RRA) algorithm. Functional enrichment analyses showed that the 152 robust DEGs are involved in multiple processes in the tumor microenvironment (TME). Immune cell infiltration analysis based on the 9 normalized GEO microarray datasets was conducted with the CIBERSORT algorithm. The changes in macrophages between ESCC and normal tissues were particularly obvious. In ESCC tissues, M0 and M1 macrophages were increased dramatically, while M2 macrophages were decreased. A robust DEG-based protein–protein interaction (PPI) network was used for hub gene selection with the CytoHubba plugin in Cytoscape. Nine hub genes (CDA, CXCL1, IGFBP3, MMP3, MMP11, PLAU, SERPINE1, SPP1 and VCAN) had high diagnostic efficiency for ESCC according to receiver operating characteristic (ROC) curve analysis. The expression of all hub genes except MMP3 and PLAU was significantly related to macrophage infiltration. Univariate and multivariate regression analyses showed that a 7-gene signature constructed from the robust DEGs was useful for predicting ESCC prognosis. Our results might facilitate the exploration of potential targeted TME therapies and prognostic evaluation in ESCC.


2018 ◽  
Vol 115 (50) ◽  
pp. E11701-E11710 ◽  
Author(s):  
Yoong Wearn Lim ◽  
Haiyin Chen-Harris ◽  
Oleg Mayba ◽  
Steve Lianoglou ◽  
Arthur Wuster ◽  
...  

Cancer immunotherapy has emerged as an effective therapy in a variety of cancers. However, a key challenge in the field is that only a subset of patients who receive immunotherapy exhibit durable response. It has been hypothesized that host genetics influences the inherent immune profiles of patients and may underlie their differential response to immunotherapy. Herein, we systematically determined the association of common germline genetic variants with gene expression and immune cell infiltration of the tumor. We identified 64,094 expression quantitative trait loci (eQTLs) that associated with 18,210 genes (eGenes) across 24 human cancers. Overall, eGenes were enriched for their being involved in immune processes, suggesting that expression of immune genes can be shaped by hereditary genetic variants. We identified the endoplasmic reticulum aminopeptidase 2 (ERAP2) gene as a pan-cancer type eGene whose expression levels stratified overall survival in a subset of patients with bladder cancer receiving anti–PD-L1 (atezolizumab) therapy. Finally, we identified 103 gene signature QTLs (gsQTLs) that were associated with predicted immune cell abundance within the tumor microenvironment. Our findings highlight the impact of germline SNPs on cancer-immune phenotypes and response to therapy; and these analyses provide a resource for integration of germline genetics as a component of personalized cancer 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.


Sign in / Sign up

Export Citation Format

Share Document