scholarly journals Identification of potential therapeutic targets for atherosclerosis by analysing the gene signature related to different immune cells and immune regulators in atheromatous plaques

2021 ◽  
Vol 14 (1) ◽  
Author(s):  
Yang Shen ◽  
Li-rong Xu ◽  
Xiao Tang ◽  
Chang-po Lin ◽  
Dong Yan ◽  
...  

Abstract Background Atherosclerosis is a chronic inflammatory disease that affects multiple arteries. Numerous studies have shown the inherent immune diversity in atheromatous plaques and suggest that the dysfunction of different immune cells plays an important role in atherosclerosis. However, few comprehensive bioinformatics analyses have investigated the potential coordinators that might orchestrate different immune cells to exacerbate atherosclerosis. Methods Immune infiltration of 69 atheromatous plaques from different arterial beds in GSE100927 were explored by single-sample-gene-set enrichment analysis (presented as ssGSEA scores), ESTIMATE algorithm (presented as immune scores) and CIBERSORT algorithm (presented as relative fractions of 22 types of immune cells) to divide these plaques into ImmuneScoreL cluster (of low immune infiltration) and ImmuneScoreH cluster (of high immune infiltration). Subsequently, comprehensive bioinformatics analyses including differentially-expressed-genes (DEGs) analysis, protein–protein interaction networks analysis, hub genes analysis, Gene-Ontology-terms and KEGG pathway enrichment analysis, gene set enrichment analysis, analysis of expression profiles of immune-related genes, correlation analysis between DEGs and hub genes and immune cells were conducted. GSE28829 was analysed to cross-validate the results in GSE100927. Results Immune-related pathways, including interferon-related pathways and PD-1 signalling, were highly enriched in the ImmuneScoreH cluster. HLA-related (except for HLA-DRB6) and immune checkpoint genes (IDO1, PDCD-1, CD274(PD-L1), CD47), RORC, IFNGR1, STAT1 and JAK2 were upregulated in the ImmuneScoreH cluster, whereas FTO, CRY1, RORB, and PER1 were downregulated. Atheromatous plaques in the ImmuneScoreH cluster had higher proportions of M0 macrophages and gamma delta T cells but lower proportions of plasma cells and monocytes (p < 0.05). CAPG, CECR1, IL18, IGSF6, FBP1, HLA-DPA1 and MMP7 were commonly related to these immune cells. In addition, the advanced-stage carotid plaques in GSE28829 exhibited higher immune infiltration than early-stage carotid plaques. Conclusions Atheromatous plaques with higher immune scores were likely at a more clinically advanced stage. The progression of atherosclerosis might be related to CAPG, IGSF6, IL18, CECR1, FBP1, MMP7, FTO, CRY1, RORB, RORC, PER1, HLA-DPA1 and immune-related pathways (IFN-γ pathway and PD-1 signalling pathway). These genes and pathways might play important roles in regulating immune cells such as M0 macrophages, gamma delta T cells, plasma cells and monocytes and might serve as potential therapeutic targets for atherosclerosis.

2021 ◽  
Vol 12 ◽  
Author(s):  
Juexing Li ◽  
Lei Zhou ◽  
Zhenhua Li ◽  
Shangneng Yang ◽  
Liangyue Tang ◽  
...  

Sepsis-induced cardiomyopathy (SIC), with a possibly reversible cardiac dysfunction, is a potential complication of septic shock. Despite quite a few mechanisms including the inflammatory mediator, exosomes, and mitochondrial dysfunction, having been confirmed in the existing research studies we still find it obscure about the overall situation of gene co-expression that how they can affect the pathological process of SIC. Thus, we intended to find out the crucial hub genes, biological signaling pathways, and infiltration of immunocytes underlying SIC. It was weighted gene co-expression network analysis that worked as our major method on the ground of the gene expression profiles: hearts of those who died from sepsis were compared to hearts donated by non-failing humans which could not be transplanted for technical reasons (GSE79962). The top 25 percent of variant genes were abstracted to identify 10 co-expression modules. In these modules, brown and green modules showed the strongest negative and positive correlation with SIC, which were primarily enriched in the bioenergy metabolism, immunoreaction, and cell death. Next, nine genes (LRRC39, COQ10A, FSD2, PPP1R3A, TNFRSF11B, IL1RAP, DGKD, POR, and THBS1) including two downregulated and seven upregulated genes which were chosen as hub genes that meant the expressive level of which was higher than the counterparts in control groups. Then, the gene set enrichment analysis (GSEA) demonstrated a close relationship of hub genes to the cardiac metabolism and the necroptosis and apoptosis of cells in SIC. Concerning immune cells infiltration, a higher level of neutrophils and B cells native and a lower level of mast cells resting and plasma cells had been observed in patients with SIC. In general, nine candidate biomarkers were authenticated as a reliable signature for deeper exploration of basic and clinical research studies on SIC.


2020 ◽  
Author(s):  
Bihui Han ◽  
Yanxiu Meng ◽  
Yumei Fan ◽  
Bing Liu ◽  
Jiajie Hou ◽  
...  

Abstract BackgroundHepatocellular carcinoma (HCC) is one of the most common malignancies and ranks as the second leading cause of cancer-related mortality worldwide. Heat shock factor 2 (HSF2) is a transcription factor that plays a critical role in development, particularly corticogenesis and spermatogenesis. However, studies on the expression and prognostic value of HSF2 and its association with tumor-infiltrating immune cells in HCC are still rare. MethodsThe TCGA, Oncomine, UALCAN, HCCDB and HPA databases were used to investigate HSF2 expression in HCC. Kaplan-Meier plotter, GEPIA and HCCDB databases were used to evaluate the association of HSF2 with the prognosis of HCC patients. Genetic alteration of HSF2 was examined by the cBioPortal database. The mechanism was investigated with Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GESA), and the relationship between HSF2 expression and immune infiltration was explored through the TIMER database and CIBERSORT algorithm.Results In the present study, we found that HSF2 expression was significantly upregulated in HCC compared with normal liver tissues. High HSF2 expression was associated with poor survival in HCC patients. GO, KEGG and GESA analyses demonstrated that HSF2 was associated with various signaling pathways, including the immune response. Notably, HSF2 expression was significantly correlated with the infiltration levels of different immune cells. HSF2 expression also displayed a significant correlation with multiple immune marker sets in HCC. ConclusionsIn summary, we explored the clinical significance of HSF2 and provided a therapeutic basis for the early diagnosis, prognostic judgment, and immunotherapy of HCC.


2021 ◽  
Vol 18 (6) ◽  
pp. 9336-9356
Author(s):  
Sidan Long ◽  
◽  
Shuangshuang Ji ◽  
Kunmin Xiao ◽  
Peng Xue ◽  
...  

<abstract> <sec><title>Background</title><p>LTB4 receptor 1 (LTB4R), as the high affinity leukotriene B4 receptor, is rapidly revealing its function in malignancies. However, it is still uncertain.</p> </sec> <sec><title>Methods</title><p>We investigated the expression pattern and prognostic significance of LTB4R in pan-cancer across different databases, including ONCOMINE, PrognoScan, GEPIA, and Kaplan-Meier Plotter, in this study. Meanwhile, we explored the significance of LTB4R in tumor metastasis by HCMDB. Then functional enrichment analysis of related genes was performed using GeneMANIA and DAVID. Lastly, utilizing the TIMER datasets, we looked into the links between LTB4R expression and immune infiltration in malignancies.</p> </sec> <sec><title>Results</title><p>In general, tumor tissue displayed higher levels of LTB4R expression than normal tissue. Although LTB4R had a negative influence on pan-cancer, a high expression level of LTB4R was protective of LIHC (liver hepatocellular carcinoma) patients' survival. There was no significant difference in the distribution of LTB4R between non-metastatic and metastatic tumors. Based on Gene Set Enrichment Analysis, LTB4R was implicated in pathways involved in inflammation, immunity, metabolism, and cancer diseases. The correlation between immune cells and LTB4R was found to be distinct across cancer types. Furthermore, markers of infiltrating immune cells, such as Treg, T cell exhaustion and T helper cells, exhibited different LTB4R-related immune infiltration patterns.</p> </sec> <sec><title>Conclusion</title><p>The LTB4R is associated with immune infiltrates and can be used as a prognostic biomarker in pan-cancer.</p> </sec> </abstract>


2020 ◽  
Author(s):  
Peipei Gao ◽  
Ting Peng ◽  
Canhui Cao ◽  
Shitong Lin ◽  
Ping Wu ◽  
...  

Abstract Background: Claudin family is a group of membrane proteins related to tight junction. There are many studies about them in cancer, but few studies pay attention to the relationship between them and the tumor microenvironment. In our research, we mainly focused on the genes related to the prognosis of ovarian cancer, and explored the relationship between them and the tumor microenvironment of ovarian cancer.Methods: The cBioProtal provided the genetic variation pattern of claudin gene family in ovarian cancer. The ONCOMINE database and Gene Expression Profiling Interactive Analysis (GEPIA) were used to exploring the mRNA expression of claudins in cancers. The prognostic potential of these genes was examined via Kaplan-Meier plotter. Immunologic signatures were enriched by gene set enrichment analysis (GSEA). The correlations between claudins and the tumor microenvironment of ovarian cancer were investigated via Tumor Immune Estimation Resource (TIMER).Results: In our research, claudin genes were altered in 363 (62%) of queried patients/samples. Abnormal expression levels of claudins were observed in various cancers. Among them, we found that CLDN3, CLDN4, CLDN6, CLDN10, CLDN15 and CLDN16 were significantly correlated with overall survival of patients with ovarian cancer. GSEA revealed that CLDN6 and CLDN10 were significantly enriched in immunologic signatures about B cell, CD4 T cell and CD8 T cell. What makes more sense is that CLDN6 and CLDN10 were found related to the tumor microenvironment. CLDN6 expression was negatively correlated with immune infiltration level in ovarian cancer, and CLDN10 expression was positively correlated with immune infiltration level in ovarian cancer. Further study revealed the CLDN6 expression level was negatively correlated with gene markers of various immune cells in ovarian cancer. And, the expression of CLDN10 was positive correlated with gene markers of immune cells in ovarian cancer.Conclusions: CLDN6 and CLDN10 were prognostic biomarkers, and correlated with immune infiltration in ovarian cancer. Our results revealed new roles for CLDN6 and CLDN10, and they were potential therapeutic targets in the treatment of ovarian cancer.


2021 ◽  
Author(s):  
Sheng Fang ◽  
Xiao Fang ◽  
Xin Xu ◽  
Lin Zhong ◽  
An-quan Wang ◽  
...  

Abstract Relevance Rheumatoid arthritis (RA) is a systemic autoimmune disease with an aggressive, chronic synovial inflammation as the main pathological change. However, the specific etiology, pathogenesis, and related biomarkers in diagnosis and treatment are still not fully elucidated. This study attempts to provide new perspectives and insights into RA at the genetic, molecular, and cellular levels through the tenet of personalized medicine. Methods Gene expression profiles of four individual knee synovial tissues were downloaded from a comprehensive gene expression database, R language was used to screen for significantly differentially expressed genes (DEGs), Gene Ontology Enrichment Analysis, Kyoto Gene Encyclopedia, and Gene Set Enrichment Analysis were performed to analyze the biological functions and signaling pathways of these DEGs, STRING online database was used to establish protein-protein interaction networks, Cytoscape software to obtain ten hub genes, Goplot to get six inflammatory immune-related hub genes, and CIBERSORT algorithm to impute immune infiltration. Results Molecular pathways that play important roles in RA were obtained: Toll-like receptors, AMPK, MAPK, TNF, FoxO, TGF-beta, PI3K and NF-κB pathways, Ten hub genes: Ccr1, Ccr2, Ccr5, Ccr7, Cxcl5, Cxcl6, Cxcl13, Ccl13, Adcy2, and Pnoc. among which Adcy2 and Pnoc have not been reported in RA studies, suggesting that they may be worthy targets for further study. It was also found that among the synoviocytes in RA, the proportions of plasma cells, CD8 T cells, follicular helper T cells, monocytes, γ delta T cells, and M0 macrophages were higher, while the proportions of CD4 memory resting T cells, regulatory T cells (Tregs), activated NK cells, resting dendritic cells, M1 macrophages, eosinophils, activated mast cells, resting mast cells were lower in proportion, and each cell played an important role in RA. Conclusions This study may help understand the key genes, molecular pathways, the role of inflammatory immune infiltrating cells in RA’s pathogenesis and provide new targets and ideas for the diagnosis and personalized treatment of RA.


2021 ◽  
Vol 12 ◽  
Author(s):  
Jian Zhou ◽  
Menghui Zhang ◽  
Yan Zhang ◽  
Xi Shi ◽  
Linlin Liu ◽  
...  

Multiple myeloma (MM) is a malignant disease of plasma cells, which remains incurable because of its unclear mechanism and drug resistance. Herein, we aimed to explore new biomarkers and therapeutic targets in MM. After screening differentially expressed genes (DEGs) in GSE6477 and GSE13591 dataset, we performed Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses of DEGs using DAVID online database. The results indicated that the downregulated DEGs were mainly enriched in the immune-associated biological process. The protein–protein interaction network was constructed by STRING database, on which we performed module analysis and identified key genes. Gene set enrichment analysis (GSEA) and Kaplan–Meier analysis showed that RRM2 could be a novel biomarker in MM diagnosis. We further confirmed that novel RRM2 inhibitor osalmid inhibited MM cell proliferation and triggered cell cycle S phase arrest. Targeting RRM2 was expected to develop new therapeutic strategies for malignant MM.


Cells ◽  
2021 ◽  
Vol 10 (3) ◽  
pp. 488
Author(s):  
Yoshihisa Tokumaru ◽  
Masanori Oshi ◽  
Ankit Patel ◽  
Wanqing Tian ◽  
Li Yan ◽  
...  

The colon adenoma–carcinoma sequence is a multistep genomic-altering process that occurs during colorectal cancer (CRC) carcinogenesis. Organoids are now commonly used to model both non-cancerous and cancerous tissue. This study aims to investigate how well organoids mimic tissues in the adenoma–carcinoma sequence by comparing their transcriptomes. A total of 234 tissue samples (48 adenomas and 186 CRC) and 60 organoid samples (15 adenomas and 45 CRC) were analyzed. We found that cell-proliferation-related gene sets were consistently enriched in both CRC tissues and organoids compared to adenoma tissues and organoids by gene set enrichment analysis (GSEA). None of the known pathways in the colon adenoma–carcinoma sequence were consistently enriched in CRC organoids. There was no enrichment of the tumor microenvironment-related gene sets in CRC organoids. CRC tissues enriched immune-response-related gene sets, whereas CRC organoids did not. The proportions of infiltrating immune cells were different between tissues and organoids, whereas there was no difference between cancer and adenoma organoids. The amounts of cancer stem cells and progenitor cells were not different between CRC and adenoma organoids, whereas a difference was noted between CRC and adenoma tissues. In conclusion, we demonstrated that organoids model only part of the adenoma–carcinoma sequence and should be used with caution after considering their limitations.


Author(s):  
Cong Luo ◽  
Zhixiong Liu ◽  
Wenrui Ye ◽  
Fangkun Liu

Background: Tumor microenvironment, especially infiltrating immune cell, is crucial for solid tumors including glioma. However, the hub genes as well as their effects on patient prognosis and immunotherapy efficacy remain obscure.Methods: We employed a total of 952 lower grade glioma (LGG) patients from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases, and 24 samples in our hospital for subsequent analyses. Abundances of immune infiltrates were evaluated using CIBERSORT and ImmuCellAI. Their correlations with prognosis were assessed by log-rank test. Immune infiltration-related hub genes were obtained from overlapped differential expressed genes (DEGs) in various subsets of survival-related immune cell types. The risk signature was constructed by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis. The functional analyses were estimated by GVSA and Gene Set Enrichment Analysis (GSEA) algorithms. And protein–protein interaction enrichment analysis was carried out with the Metascape database integrating STRING, BioGrid, OmniPath, and InWeb_IM.Results: Among the 21 infiltrates, the abundances of five immune infiltrates were correlated with overall survival (OS) in LGG patients. Higher abundances of naïve CD4+ T cells (p = 0.002), activated mast cells (p = 0.015), and monocytes (p = 0.014) were correlated with better prognosis, while higher abundances of resting memory CD4+ T cells (p = 0.015) and M1 macrophages (p = 0.020) correlated with poorer OS. We finally obtained 44 hub genes and constructed an immune infiltration-related signature (IIRS). The IIRS correlates with clinicopathological characteristics and exhibited potential power in predicting the immunotherapy efficacy. The IRRS correlates with cancer related pathways, especially “epithelial-mesenchymal transition (EMT),” and cytotoxic T lymphocytes.Conclusion: Our study constructed and validated a novel signature for risk stratification and prediction of immunotherapy response in grade II and III gliomas, which was closely associated with glioma immune microenvironment and could serve as a promising prognostic biomarker for glioma patients.


2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Siqing Sun ◽  
Yutao Wang ◽  
Jianfeng Wang ◽  
Jianbin Bi

Abstract Background The Wnt signaling pathway is core to the growth of bladder tumors. Epithelial-to-mesenchymal transition (EMT) is significant for bladder tumor metastasis. Nevertheless, the relationship between the Wnt signaling pathway, outcomes of bladder cancer (BLCA), and the specific mechanisms driving immune infiltration have not been studied. Methods We obtained Wnt pathway-related gene mRNA and clinicopathological data from the Cancer Genome Atlas (TCGA). We obtained 34 genes that were greatly correlated with outcome using univariate Cox regression analysis and conducted a completely randomized data t-test to perform clinical staging. According to the single-sample gene set enrichment analysis (ssGSEA), the weighted correlation network analysis (WGCNA) was applied to identify relevant biological functions. Various subtypes were identified using consensus cluster analysis. Univariate Cox regression and least absolute shrinkage sum selection operator–Cox regression algorithm analysis were conducted on TCGA and Gene Expression Omnibus data to identify risk characteristics. The Kaplan–Meier method and receiver running feature curves were adopted to calculate overall survival. Single-sample gene set enrichment analysis (ssGSEA) was adopted for the assessment of the degree of immune infiltration. Then, we demonstrated the relationship between PPP2CB and EMT function in two cell lines. Results Thirty-four Wnt signaling pathway-related genes were risk factors for BLCA outcome, and their expression levels differed by clinical stage. The co-expression of WGCNA showed the relationship between the Wnt signaling pathway and biological functions and was closely associated with EMT. We divided BLCA patients into two subtypes using consensus clustering. Survival curves and clinical analysis showed that the Wnt pathway enriched group had worse outcomes. The Wnt signature showed the significance of the outcome for MAPK10, PPP2CB, and RAC3. Based on these genes, the degree of immune infiltration was evaluated. Cell function experiments suggested that PPP2CB drives the proliferation and migration of BLCA cells. Conclusion We found that Wnt signaling pathway-related genes can be used as prognostic risk factors for BLCA, and the Wnt signaling pathway is a cancer-promoting signaling pathway associated with EMT. We identified three critical genes: MAPK10, RAC3, and PPP2CB. The genes in these three Wnt signaling pathways are associated with tumor cell EMT and immune cell infiltration. The most important finding was that these three genes were independent prognostic factors for BLCA.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yumei Fan ◽  
Bing Liu ◽  
Fei Chen ◽  
Zhiyuan Song ◽  
Bihui Han ◽  
...  

Lung cancer has the highest death rate among cancers globally. Hepcidin is a fascinating regulator of iron metabolism; however, the prognostic value of hepcidin and its correlation with immune cell infiltration in lung cancer remain unclear. Here, we comprehensively clarified the prognostic value and potential function of hepcidin in lung cancer. Hepcidin expression was significantly increased in lung cancer. High hepcidin expression was associated with sex, age, metastasis, and pathological stage and significantly predicted an unfavorable prognosis in lung cancer patients. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) results suggested that hepcidin is involved in the immune response. Furthermore, hepcidin expression was positively correlated with the infiltration levels of immune cells and the expression of diverse immune cell marker sets. Importantly, hepcidin may affect prognosis partially by regulating immune infiltration in lung cancer patients. Hepcidin may serve as a candidate prognostic biomarker for determining prognosis associated with immune infiltration in lung cancer.


Sign in / Sign up

Export Citation Format

Share Document