Integrated Bioinformatical Analysis Identified HTRA4 As An Immune-Related Prognostic Biomarker In The Gastric Cancer Microenvironment

Author(s):  
Ce Ji ◽  
Jingwei Dai ◽  
Zhen Wang ◽  
Xiangchen Hu ◽  
Youwei Kou

Abstract HtrA serine peptidase 4 (HTRA4), which belongs to the AB family of proteases/chapterones, participates in multiple signal pathways and plays an important regulatory role in some malignancies; however, its role in the prognosis and immune infiltrates of gastric cancer (GC) remains unclear. HTRA4 expression was thus investigated in tumor tissues, and its association with immune infiltrates was assessed, to determine its prognostic roles in GC patients. Data for patients with GC was collected from three GEO datasets. This study has investigated if high expression levels of HTRA4 are associated with a poor prognosis in GC patients, especially in T2,T3,N1,NI&N2&N3,M1,stage3 and HER2+ subgroups. HTRA4 was positively correlated with CD8+ T cells, CD4+ memory activated T cells, gamma delta T cells, and M2 macrophages in the GC microenvironment. Gene Ontology and Gene set Enrichment Analysis showed that HTRA4 was enriched in some immune-related pathways. High expression levels of HTRA4 were significantly correlated with tumors in the T stage, N stage, and clinical stage, as well as those of histological grade. In summary, HTRA4 was found to play a vital role in the progression of GC and may be a predictive biomarker for the survival of GC patients.

2021 ◽  
Author(s):  
beibei xu ◽  
Endian Zheng ◽  
Yi Huang ◽  
Liang Zheng ◽  
Qiaoli Lan ◽  
...  

Abstract BackgroundCircular RNA (circRNA) has been shown to be an important regulator in gastric cancer (GC). However, functions and regulatory mechanisms of circRNA-related competitive endogenous RNA (ceRNA) in GC have not been established.MethodsCircRNA data and clinical data were downloaded from the GEO and TCGA databases. The ceRNA and Protein-Protein Interaction (PPI) networks were constructed through bioinformatics analysis. Function enrichment analysis was performed. Additionally, correlations between expression levels of the top 10 hub genes and immune cell infiltration levels, histopathological grade and clinical stage were determined to establish their clinical values. The differentially expressed circRNA (DEcircRNA) was validated by quantitative real-time PCR (qRT-PCR).ResultsScreening of the GEO and TCGA databases revealed a total of 1627 DEcircRNAs, 6516 DEmRNAs, and 1451 DEmiRNAs. The ceRNA interaction network contained 2 circRNAs, 3 miRNAs and 55 mRNAs. Meanwhile, Gene Ontology (GO) analysis revealed a total of 323 biological processes (BP) terms, 53 cellular components (CC) terms, 51 molecular functions (MF) terms, while the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed 4 signaling pathways. Gene Set Enrichment Analysis (GSEA) analysis revealed that EPHA4, NCAM1 and NRXN1 were positively correlated with the axon guidan and adhesion molecules pathways. Most of top 10 hub genes were positively correlated with B cells, CD8+ T cells, CD4+ T cells, Neutrophils and Dendritic Cell infiltration. Correlation analysis between hub genes and clinical phenotypes revealed that elevated expressions of EPHA4 and KCNA1 indicated poor tissue differentiation and were associated with clinically advanced stages of GC. The qRT-PCR results revealed that the expression of has_circ_0002504 was significantly down-regulated in 3 GC cell lines which was consistent with the results of our bioinformatics analysis.ConclusionsHas_circ_0001998 and has_circ_0002504 are potential diagnostic biomarkers for GC, and the high expressions of both EPHA4 and KCNA1 may predict poor prognosis.


2021 ◽  
Author(s):  
Ting Liu ◽  
Zheng Gong ◽  
Hong Zhang ◽  
Yi Wan ◽  
Ming-Han Ren ◽  
...  

Abstract BackgroundGastric cancer (GC) is the fifth most common cancer worldwide. Previous studies have suggested that the tumor microenvironment (TME) plays an important role in the development and prognosis of GC. In this study, we aimed to identify genes in tumor-infiltrating immune cells (TICs) that influence the progression and prognosis of GC. MethodsWe used the ESTIMATE algorithm to calculate the scores of the stromal and immune components of the TME in 407 GC samples collected from The Cancer Genome Atlas (TCGA) database.The differentially expressed genes (DEGs) were intersected by a protein-protein interaction (PPI) network and analyzed by univariate Cox regression.Further analysis showed the correlation between MCEMP1 and the clinicopathological characteristics of GC patients (clinical stage, distant metastasis) and survival.Then we used Gene set enrichment analysis (GSEA) and CIBERSORT analysis to examine the relationship between MCEMP1 and the TME.ResultsThe analysis revealed that the expression of MCEMP1 was positively correlated with the clinicopathological characteristics of GC patients (clinical stage, distant metastasis) and negatively correlated with survival. Gene set enrichment analysis (GSEA) indicated that gene sets in the MCEMP1 high expression group were concentrated mainly in immune-related pathways. CIBERSORT analysis of the proportion of TICs revealed that neutrophils and M2 macrophages were positively correlated with MCEMP1 expression, suggesting that MCEMP1 is responsible for preservation of the immune-dominant status of the TME. ConclusionHigh MCEMP1 expression might be a biomarker of a poor prognosis in GC patients and provide a clue regarding the different statuses of the TME, offering additional insight into therapy for GC.


2021 ◽  
Author(s):  
Long-Jiang Chen ◽  
Lu-Lu Zhai ◽  
Wei Wang ◽  
Lun Wu ◽  
Li-Chao Yao ◽  
...  

Abstract While previous studies have suggested that B3GNT3 is associated with tumorigenesis and progression of several tumors, its expression level and clinical significance in pancreatic adenocarcinoma (PAAD) remains unclear. Our study aimed to investigate the role of B3GNT3 in PAAD. B3GNT3 RNA sequencing and clinicopathological data were collected from the TCGA, GEO and GTEx databases. We assessed the expression and prognostic value of B3GNT3 in PAAD using R program and attached packages. Additionally, we investigate the correlation between B3GNT3 expression and tumor-infiltrating immune cells using CIBERSORT and the “correlation” module of GEPIA. Finally, gene set enrichment analysis (GSEA) was used to elucidate B3GNT3 related signaling pathways in PAAD. Results showed that B3GNT3 expression was significantly higher in tumor tissues compared to normal tissues (P <0.05). Increased B3GNT3 expression was correlated with advanced histologic grade and stage (Ⅰ Vs Ⅱ). Patients with high B3GNT3 expression had a worse OS (HR = 1.713, P = 0.0005). Moreover, a negative correlation between increased B3GNT3 expression and immune infiltrating level of naive B cells, CD8 T cells, and CD4 memory activated T cells was revealed by CIBERSORT analysis. Then, further analysis verified the correlation using the “correlation” module of GEPIA. Finally, GSEA suggested that functional enrichment of B3GNT3 was mainly involved in pathways in cancer, p53 signaling pathway, TGF beta signaling pathway, catabolic and transport processes of proteins, etc. Collectively, these results suggested that overexpression of B3GNT3 might affect the infiltration of immune cells and could act as a potential prognostic biomarker of PAAD.


2020 ◽  
Vol 10 ◽  
Author(s):  
Ce Ji ◽  
Li-Sha Sun ◽  
Fei Xing ◽  
Nan Niu ◽  
Hong-Li Gao ◽  
...  

HtrA serine peptidase 3 (HTRA3) participates in multiple signal pathways and plays an important regulatory role in various malignancies; however, its role on prognosis and immune infiltrates in gastric cancer (GC) remains unclear. The study investigated HTRA3 expression in tumor tissues and its association with immune infiltrates, and determined its prognostic roles in GC patients. Patients with GC were collected from the cancer genome atlas (TCGA). We compared the expression of HTRA3 in GC and normal gastric mucosa tissues with Wilcoxon rank sum test. And logistic regression was used to evaluate the relationship between HTRA3 and clinicopathological characters. Gene ontology (GO) term analysis, Gene set enrichment analysis (GSEA), and single-sample Gene Set Enrichment Analysis (ssGSEA) was conducted to explain the enrichmental pathways and functions and quantify the extent of immune cells infiltration for HTRA3. Kaplan-Meier analysis and Cox regression were performed to evaluate the correlation between HTRA3 and survival rates. A nomogram, based on Cox multivariate analysis, was used to predict the impact of HTRA3 on prognosis. High HTRA3 expression was significantly correlated with tumor histological type, histological grade, clinical stage, T stage, and TP53 status (P &lt; 0.05). HTRA3-high GC patients had a lower 10-year progression-free interval [PFI; hazard ratio (HR): 1.46; 95% confidence interval (CI): 1.02–2.08; P = 0.038], disease-specific survival (DSS; HR: 1.65; CI: 1.08–2.52; P = 0.021) and overall survival (OS; HR: 1.59; CI: 1.14–2.22; P = 0.006). Multivariate survival analysis showed that HTRA3 was an independent prognostic marker for PFI (HR: 1.456; CI: 1.021–2.078; P = 0.038), DSS (HR: 1.650; CI: 1.079–2.522; P = 0.021) and OS [hazard ratio (HR): 1.590; 95% confidence interval (CI):1.140–2.219; P = 0.006]. The C-indexes and calibration plots of the nomogram based on multivariate analysis indicated an effective predictive performance for GC patients. GSEA showed that High HTRA3 expression may activate NF-κB pathway, YAP1/WWTR1/TAZ pathway, and TGFβ pathway. There was a negative correlation between the HTRA3 expression and the abundances of adaptive immunocytes (T helper cell 17 cells) and a positive correlation with abundances of innate immunocytes (natural killer cells, macrophages etc.). HTRA3 plays a vital role in GC progression and prognosis and could be a moderate biomarker for prediction for survival after gastrectomy.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Tian Xia ◽  
Lian Meng ◽  
Zhijuan Zhao ◽  
Yujun Li ◽  
Hao Wen ◽  
...  

Abstract Background Rhabdomyosarcoma (RMS) is a malignant soft-tissue tumour. In recent years, the tumour microenvironment (TME) has been reported to be associated with the development of tumours. However, the relationship between the occurrence and development of RMS and TME is unclear. The purpose of this study is to identify potential tumor microenvironment-related biomarkers in rhabdomyosarcoma and analyze their molecular mechanisms, diagnostic and prognostic significance. Methods We first applied bioinformatics method to analyse the tumour samples of 125 patients with rhabdomyosarcoma (RMS) from the Gene Expression Omnibus database (GEO). Differential genes (DEGs) that significantly correlate with TME and the clinical staging of tumors were extracted. Immunohistochemistry (IHC) was applied to validate the expression of mitotic arrest deficient 2 like 1 (MAD2L1) and cyclin B2 (CCNB2) in RMS tissue. Then, we used cell function and molecular biology techniques to study the influence of MAD2L1 and CCNB2 expression levels on the progression of RMS. Results Bioinformatics results show that the RMS TME key genes were screened, and a TME-related tumour clinical staging model was constructed. The top 10 hub genes were screened through the establishment of a protein–protein interaction (PPI) network, and then Gene Expression Profiling Interactive Analysis (GEPIA) was conducted to measure the overall survival (OS) of the 10 hub genes in the sarcoma cases in The Cancer Genome Atlas (TCGA). Six DEGs of statistical significance were acquired. The relationship between these six differential genes and the clinical stage of RMS was analysed. Further analysis revealed that the OS of RMS patients with high expression of MAD2L1 and CCNB2 was worse and the expression of MAD2L1 and CCNB2 was related to the clinical stage of RMS patients. Gene set enrichment analysis (GSEA) revealed that the genes in MAD2L1 and CCNB2 groups with high expression were mainly related to the mechanism of tumour metastasis and recurrence. In the low-expression MAD2L1 and CCNB2 groups, the genes were enriched in the metabolic and immune pathways. Immunohistochemical results also confirmed that the expression levels of MAD2L1 (30/33, 87.5%) and CCNB2 (33/33, 100%) were remarkably higher in RMS group than in normal control group (0/11, 0%). Moreover, the expression of CCNB2 was related to tumour size. Downregulation of MAD2L1 and CCNB2 suppressed the growth, invasion, migration, and cell cycling of RMS cells and promoted their apoptosis. The CIBERSORT immune cell fraction analysis indicated that the expression levels of MAD2L1 and CCNB2 affected the immune status in the TME. Conclusions The expression levels of MAD2L1 and CCNB2 are potential indicators of TME status changes in RMS, which may help guide the prognosis of patients with RMS and the clinical staging of tumours.


2022 ◽  
Vol 12 ◽  
Author(s):  
Meng-Di Xia ◽  
Rui-Ran Yu ◽  
Dong-Ming Chen

BackgroundAntineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a systemic autoimmune disease that generally induces the progression of rapidly progressive glomerulonephritis (GN). The purpose of this study was to identify key biomarkers and immune-related pathways involved in the progression of ANCA-associated GN (ANCA-GN) and their relationship with immune cell infiltration.MethodsGene microarray data were downloaded from the Gene Expression Omnibus (GEO). Hub markers for ANCA-GN were mined based on differential expression analysis, weighted gene co-expression network analysis (WGCNA) and lasso regression, followed by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) of the differential genes. The infiltration levels of 28 immune cells in the expression profile and their relationship to hub gene markers were analysed using single-sample GSEA (ssGSEA). In addition, the accuracy of the hub markers in diagnosing ANCA-GN was subsequently evaluated using the receiver operating characteristic curve (ROC).ResultsA total of 651 differential genes were screened. Twelve co-expression modules were obtained via WGCNA; of which, one hub module (black module) had the highest correlation with ANCA-GN. A total of 66 intersecting genes were acquired by combining differential genes. Five hub genes were subsequently obtained by lasso analysis as potential biomarkers for ANCA-GN. The immune infiltration results revealed the most significant relationship among monocytes, CD4+ T cells and CD8+ T cells. ROC curve analysis demonstrated a prime diagnostic value of the five hub genes. According to the functional enrichment analysis of the differential genes, hub genes were mainly enhanced in immune- and inflammation-related pathways.ConclusionB cells and monocytes were closely associated with the pathogenesis of ANCA-GN. Hub genes (CYP3A5, SLC12A3, BGN, TAPBP and TMEM184B) may be involved in the progression of ANCA-GN through immune-related signal pathways.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yangming Hou ◽  
Yingjuan Xu ◽  
Dequan Wu

AbstractThe infiltration degree of immune and stromal cells has been shown clinically significant in tumor microenvironment (TME). However, the utility of stromal and immune components in Gastric cancer (GC) has not been investigated in detail. In the present study, ESTIMATE and CIBERSORT algorithms were applied to calculate the immune/stromal scores and the proportion of tumor-infiltrating immune cell (TIC) in GC cohort, including 415 cases from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were screened by Cox proportional hazard regression analysis and protein–protein interaction (PPI) network construction. Then ADAMTS12 was regarded as one of the most predictive factors. Further analysis showed that ADAMTS12 expression was significantly higher in tumor samples and correlated with poor prognosis. Gene Set Enrichment Analysis (GSEA) indicated that in high ADAMTS12 expression group gene sets were mainly enriched in cancer and immune-related activities. In the low ADAMTS12 expression group, the genes were enriched in the oxidative phosphorylation pathway. CIBERSORT analysis for the proportion of TICs revealed that ADAMTS12 expression was positively correlated with Macrophages M0/M1/M2 and negatively correlated with T cells follicular helper. Therefore, ADAMTS12 might be a tumor promoter and responsible for TME status and tumor energy metabolic conversion.


2021 ◽  
Vol 104 (2) ◽  
pp. 003685042110113
Author(s):  
Yusha Chen ◽  
Xiaoqian Lin ◽  
Jinwen Zheng ◽  
Jiancui Chen ◽  
Huifeng Xue ◽  
...  

Apelin (APLN) is recently demonstrated a direct association with many malignant diseases. However, its effects on cervical cancer remain unclear. This study therefore aims to evaluate the association between APLN expression and cervical cancer using publicly available data from The Cancer Genome Atlas (TCGA). The Pearson χ2 test and Fish exact test, as well as logistic regression, were used to evaluate the relationship between clinicopathological factors in cervical cancer and the expression of APLN. Additionally, the Cox regression and Kaplan-Meier methods were conducted to analyze the Overall Survival (OS) of cervical cancer patients in TCGA. Finally, gene set enrichment analysis (GSEA) was performed to establish its biological functions. High expression of APLN in cervical cancer was significantly associated with a more advanced clinical stage (OR = 1.91 (1.21–3.05) for Stage II, Stage III, and Stage IV vs Stage I, p = 0.006). Additionally, it was associated with poor outcome after primary therapy (OR = 2.14 (1.03–4.59) for Progressive Disease (PD), Stable Disease (SD), and Partial Response (PR) vs Complete Remission (CR), p = 0.045) and high histologic grade (OR = 1.67 (1.03–2.72) for G3 and G4 vs G1 and G2, p = 0.037). Moreover, multivariate analysis showed that high expression of APLN was associated with a shorter OS. GSEA demonstrated that six KEGG pathways, including PPAR signaling, ECM-receptor interaction, focal adhesion, MAPK signaling, TGF-beta signaling, and Gap junction pathways were differentially enriched in the high expression APLN phenotype. The recent study suggests that APLN plays an important role in the progression of cervical cancer and might be a promising prognostic biomarker of the disease.


2021 ◽  
Vol 8 ◽  
Author(s):  
Jinfeng Zhu ◽  
Chen Luo ◽  
Jiefeng Zhao ◽  
Xiaojian Zhu ◽  
Kang Lin ◽  
...  

Background: Lysyl oxidase (LOX) is a key enzyme for the cross-linking of collagen and elastin in the extracellular matrix. This study evaluated the prognostic role of LOX in gastric cancer (GC) by analyzing the data of The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) dataset.Methods: The Wilcoxon rank-sum test was used to calculate the expression difference of LOX gene in gastric cancer and normal tissues. Western blot and immunohistochemical staining were used to evaluate the expression level of LOX protein in gastric cancer. Kaplan-Meier analysis was used to calculate the survival difference between the high expression group and the low expression group in gastric cancer. The relationship between statistical clinicopathological characteristics and LOX gene expression was analyzed by Wilcoxon or Kruskal-Wallis test and logistic regression. Univariate and multivariate Cox regression analysis was used to find independent risk factors affecting the prognosis of GC patients. Gene set enrichment analysis (GSEA) was used to screen the possible mechanisms of LOX and GC. The CIBERSORT calculation method was used to evaluate the distribution of tumor-infiltrating immune cell (TIC) abundance.Results: LOX is highly expressed in gastric cancer tissues and is significantly related to poor overall survival. Wilcoxon or Kruskal-Wallis test and Logistic regression analysis showed, LOX overexpression is significantly correlated with T-stage progression in gastric cancer. Multivariate Cox regression analysis on TCGA and GEO data found that LOX (all p &lt; 0.05) is an independent factor for poor GC prognosis. GSEA showed that high LOX expression is related to ECM receptor interaction, cancer, Hedgehog, TGF-beta, JAK-STAT, MAPK, Wnt, and mTOR signaling pathways. The expression level of LOX affects the immune activity of the tumor microenvironment in gastric cancer.Conclusion: High expression of LOX is a potential molecular indicator for poor prognosis of gastric cancer.


Sign in / Sign up

Export Citation Format

Share Document