scholarly journals Significance of hub genes and immune cell infiltration identified by bioinformatics analysis in pelvic organ prolapse

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e9773
Author(s):  
Ying Zhao ◽  
Zhijun Xia ◽  
Te Lin ◽  
Yitong Yin

Objective Pelvic organ prolapse (POP) refers to the decline of pelvic organ position and dysfunction caused by weak pelvic floor support. The aim of the present study was to screen the hub genes and immune cell infiltration related to POP disease. Methods Microarray data of 34 POP tissues in the GSE12852 gene expression dataset were used as research objects. Weighted gene co-expression network analysis (WGCNA) was performed to elucidate the hub module and hub genes related to POP occurrence. Gene function annotation was performed using the DAVID tool. Differential analysis based on the GSE12852 dataset was carried out to explore the expression of the selected hub genes in POP and non-POP tissues, and RT-qPCR was used to validate the results. The differential immune cell infiltration between POP and non-POP tissues was investigated using the CIBERSORT algorithm. Results WGCNA revealed the module that possessed the highest correlation with POP occurrence. Functional annotation indicated that the genes in this module were mainly involved in immunity. ZNF331, THBS1, IFRD1, FLJ20533, CXCR4, GEM, SOD2, and SAT were identified as the hub genes. Differential analysis and RT-qPCR demonstrated that the selected hub genes were overexpressed in POP tissues as compared with non-POP tissues. The CIBERSORT algorithm was employed to evaluate the infiltration of 22 immune cell types in POP tissues and non-POP tissues. We found greater infiltration of activated mast cells and neutrophils in POP tissues than non-POP tissues, while the infiltration of resting mast cells was lower in POP tissues. Moreover, we investigated the relationship between the type of immune cell infiltration and hub genes by Pearson correlation analysis. The results indicate that activated mast cells and neutrophils had a positive correlation with the hub genes, while resting mast cells had a negative correlation with the hub genes. Conclusions Our research identified eight hub genes and the infiltration of three types of immune cells related to POP occurrence. These hub genes may participate in the pathogenesis of POP through the immune system, giving them a certain diagnostic and therapeutic value.

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Xiaoqing Shen ◽  
Zhujuan Yang ◽  
Songwei Feng ◽  
Yi Li

Abstract Background While large-scale genomic analyses symbolize a precious attempt to decipher the molecular foundation of uterine leiomyosarcoma (ULMS), bioinformatics results associated with the occurrence of ULMS based totally on WGCNA and CIBERSORT have not yet been reported. This study aimed to screen the hub genes and the immune cell infiltration pattern in ULMS by bioinformatics methods. Methods Firstly, the GSE67463 dataset, including 25 ULMS tissues and 29 normal myometrium (NL) tissues, was downloaded from the public database. The differentially expressed genes (DEGs) were screened by the ‘limma’ package and hub modules were identified by weighted gene co-expression network analysis (WGCNA). Subsequently, gene function annotations were performed to investigate the biological role of the genes from the intersection of two groups (hub module and DEGs). The above genes were calculated in the protein–protein interaction (PPI) network to select the hub genes further. The hub genes were validated using external data (GSE764 and GSE68295). In addition, the differential immune cell infiltration between UL and ULMS tissues was investigated using the CIBERSORT algorithm. Finally, we used western blot to preliminarily detect the hub genes in cell lines. Results WGCNA analysis revealed a green-yellow module possessed the highest correlation with ULMS, including 1063 genes. A total of 172 DEGs were selected by thresholds set in the ‘limma’ package. The above two groups of genes were intersected to obtain 72 genes for functional annotation analysis. Interestingly, it indicated that 72 genes were mainly involved in immune processes and the Neddylation pathway. We found a higher infiltration of five types of cells (memory B cells, M0-type macrophages, mast cells activated, M1-type macrophages, and T cells follicular helper) in ULMS tissues than NL tissues, while the infiltration of two types of cells (NK cells activated and mast cells resting) was lower than in NL tissues. In addition, a total of five genes (KDR, CCL21, SELP, DPT, and DCN) were identified as the hub genes. Internal and external validation demonstrated that the five genes were over-expressed in NL tissues compared with USML tissues. Finally, the correlation analysis results indicate that NK cells activated and mast cells activated positively correlated with the hub genes. However, M1-type macrophages had a negative correlation with the hub genes. Moreover, only the DCN may be associated with the Neddylation pathway. Conclusion A series of evidence confirm that the five hub genes and the infiltration of seven types of immune cells are related to USML occurrence. These hub genes may affect the occurrence of USML through immune-related and Neddylation pathways, providing molecular evidence for the treatment of USML in the future.


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.


2021 ◽  
Vol 2021 ◽  
pp. 1-10
Author(s):  
Li Zhang ◽  
Yunlong Yang ◽  
Dechun Geng ◽  
Yonghua Wu

Background. Osteoporosis is characterized by low bone mass, deterioration of bone tissue structure, and susceptibility to fracture. New and more suitable therapeutic targets need to be discovered. Methods. We collected osteoporosis-related datasets (GSE56815, GSE99624, and GSE63446). The methylation markers were obtained by differential analysis. Degree, DMNC, MCC, and MNC plug-ins were used to screen the important methylation markers in PPI network, then enrichment analysis was performed. ROC curve was used to evaluate the diagnostic effect of osteoporosis. In addition, we evaluated the difference in immune cell infiltration between osteoporotic patients and control by ssGSEA. Finally, differential miRNAs in osteoporosis were used to predict the regulators of key methylation markers. Results. A total of 2351 differentially expressed genes and 5246 differentially methylated positions were obtained between osteoporotic patients and controls. We identified 19 methylation markers by PPI network. They were mainly involved in biological functions and signaling pathways such as apoptosis and immune inflammation. HIST1H3G, MAP3K5, NOP2, OXA1L, and ZFPM2 with higher AUC values were considered key methylation markers. There were significant differences in immune cell infiltration between osteoporotic patients and controls, especially dendritic cells and natural killer cells. The correlation between MAP3K5 and immune cells was high, and its differential expression was also validated by other two datasets. In addition, NOP2 was predicted to be regulated by differentially expressed hsa-miR-3130-5p. Conclusion. Our efforts aim to provide new methylation markers as therapeutic targets for osteoporosis to better treat osteoporosis in the future.


2020 ◽  
Author(s):  
Jukun Song ◽  
Song He ◽  
Wei Wang ◽  
Jiaming Su ◽  
Dongbo Yuan ◽  
...  

Abstract Background Immune infiltration of Prostate cancer (PCa) was highly related to clinical outcomes. However, previous works failed to elucidate the diversity of different immune cell types that make up the function of the immune response system. The aim of the study was to uncover the composition of TIICs in PCa utilizing the CIBERSORT algorithm and further reveal the molecular characteristics of PCa subtypes. Method In the present work, we employed the CIBERSORT method to evaluate the relative proportions of immune cell profiling in PCa and adjacent samples, normal samples. We analyzed the correlation between immune cell infiltration and clinical information. The tumor-infiltrating immune cells of the TCGA PCa cohort were analyzed for the first time. The fractions of 22 immune cell types were imputed to determine the correlation between each immune cell subpopulation and clinical feature. Three types of molecular classification were identified via R-package of “CancerSubtypes”. The functional enrichment was analyzed in each subtype. The submap and TIDE algorithm were used to predict the clinical response to immune checkpoint blockade, and GDSC was employed to screen chemotherapeutic targets for the potential treatment of PCa. Results In current work, we utilized the CIBERSORT algorithm to assess the relative proportions of immune cell profiling in PCa and adjacent samples, normal samples. We investigated the correlation between immune cell infiltration and clinical data. The tumor-infiltrating immune cells in the TCGA PCa cohort were analyzed. The 22 immune cells were also calculated to determine the correlation between each immune cell subpopulation and survival and response to chemotherapy. Three types of molecular classification were identified. Each subtype has specific molecular and clinical characteristics. Meanwhile, Cluster I is defined as advanced PCa, and is more likely to respond to immunotherapy. Conclusions Our results demonstrated that differences in immune response may be important drivers of PCa progression and response to treatment. The deconvolution algorithm of gene expression microarray data by CIBERSOFT provides useful information about the immune cell composition of PCa patients. In addition, we have found a subtype of immunopositive PCa subtype and will help to explore the reasons for the poor effect of PCa on immunotherapy, and it is expected that immunotherapy will be used to guide the individualized management and treatment of PCa patients.


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.


Sign in / Sign up

Export Citation Format

Share Document