A Comprehensive Analysis of the Glutathione Peroxidase 8 (GPX8) in Human Cancer

Objective: Nowadays, cancer is still a leading public health problem all over the world. Several studies have reported the GPX8 could be correlated with the poor prognostic of Gastric Cancer and Breast Cancer. However, the prognostic potential of GPX8 in pan-cancer remains unclear. In this work, we aimed to explore the prognostic and immunological role of GPX8 in human cancer and confirm the oncogenic value in GBM.

Methods: The data of TCGA, CPTAC and GEO databases were adopted for the survival analysis. Based on the RNAseq and Methylation450 data of TCGA, the R language and package “ggplot2” were used to analyze the DNA methylation at the region of the promoter of GPX8 in tumors. The genetic alteration of GPX8 from TCGA cancers was investigated in cBioPortal. The R package “GSVA” and “ssGSEA” were employed to evaluate the correlation of GPX8 expression with the immune infiltration. The KEGG website was used for pathway analysis. The STRING website and GEPIA were performed to predict GPX8-binding proteins. The R package “ggplot2” and “clusterprofile” were used to analyze and visualize the GO and KEGG analysis. A normal human astrocyte cell line and three GBM cell lines were cultured under suitable conditions. The shRNA was transferred to cells by Lipofectamine 3000. The qRT-PCR and WB were adopted to detect the expression of GPX8. The wound-healing assay and transwell assay were taken to analyze the invasive and metastatic abilities. The tumor tissues and paracancerous ones were collected from patients with GBM. WB assay was employed to analyze the expression of GPX8 protein.

Results: GPX8 was a valuable diagnostic biomarker in multiple cancers, including GBM/LGG (glioblastoma multiforme/Brain lower grade glioma), KIRC (kidney renal clear cell carcinoma), KIRP (kidney renal papillary cell carcinoma) and STAD (stomach adenocarcinoma). Moreover, we observed a correlation between the expression of GPX8 and the reduced DNA methylation at the promoter region in several tumors, such as GBM/LGG. Our results indicated a positive correlation between the GPX8 expression and immune infiltration. In addition, the enrichment analysis demonstrated that antioxidant activity was mainly involved in the functional mechanism of GPX8. In particular, we first confirmed the up-regulated of GPX8 in GBM cells and observed the suppression of migrative and invasive phenotypes by knockdown of GPX8. Furthermore, we confirmed the expression of GPX8 was higher in GBM tumor tissues than paracancerous ones.

Conclusion: Our study showed a correlation of GPX8 expression with clinical prognosis, DNA methylation and immune infiltrates. Furthermore, we first confirmed GPX8 was highly expressed in GBM cells and contributed to migration and invasion. These results provided a predictive biomarker and an inclusive understanding of the GPX8 expression in multiple tumors types, especially in GBM.

Introduction

Nowadays, cancer is a main public health problem all over the world. Partly because of the benefit of significant development in the recognition, diagnosis, and treatment, the mortality ratio has improved in many common cancer types (1). However, compared to the UK and USA, China has a larger scale of cancer types with a poorer 5-year overall survival rate. The mortality and incidence ratio in patients with cancer is still extremely high (2, 3). Therefore, it is necessary to search for new biomarkers and treatment targets for cancers.

Since tumorigenesis is a complex and variable process, multiple genetic changes could impact the malignant transformation and clinical outcomes in cancer patients (4, 5). The pan-cancer expression analysis, which contains massive data from publicly funded databases, has been deemed an essential and effective means for researchers to analyze the correlation between gene expression and clinical prognosis. For example, previous studies have demonstrated that the SND1, HER2 and LAYN could act as predictive biomarkers for cancer prognosis (6–8).

It is generally accepted that the immune system activated by tumor cells could cause immune surveillance and defense to inhibit tumor progression (9). Nevertheless, as a double-edged sword, more and more evidence indicates that the multiple immune cells recruited by progressed tumors could stimulate an inflammatory environment and promote the secretion of cytokines and chemokines, which contribute to tumor growth, metastasis, invasion and pathological angiogenesis (10–12). Research of immune-related mechanisms is a promising direction and could prompt the development of cancer treatment (2). For instance, the programmed death-1 (PD-1), an immune checkpoint molecule, was been recognized as a target for cancer treatment (13). Multiple medicines targeted PD-1 have been approved by the FDA for the usage in malignant melanoma (13–15). Moreover, the tumor-infiltrating Immune cells have been verified to affect the prognosis (16).

Glutathione peroxidases (GPxs) (EC 1.11.1.9) are a kind of enzyme family that could protect organisms from oxidative damage by the reactive oxygen species (ROS) and has an important influence on substance metabolism (17–20). As a member of the GPX family, GPX8 is a transmembrane protein localized in the ER (21). It has been shown GPX8 can reduce oxidized PDI and prevent Ero1α-derived H2O2 leaking by interacting with Ero1α, a member of the sulfhydryl oxidase Ero1 family (21, 22). Increasing evidence has verified that GPX8 could be correlated with the poor prognostic of cancers and involved in diverse physiological processes and tumorigenesis, including gastric cancer, lung cancer and others (23–25). Furthermore, reports have indicated that ROS could modulate the function of immune cells (26–28). However, although the emerging genes are related to the prognosis, we have never retrieved any report with a pan-cancer analysis of GPX8. Thus, in this work, we aimed to explore the prognostic and immunological role of GPX8 in human cancer and confirm the oncogenic value in GBM.

Materials and Methods

GPX8 Expression in Human Cancers

The expression of GPX8 gene in different types of cancers was explored in the TIMER database (http://timer.cistrome.org/) and ONCOMINE database (www.oncomine.org) with the conditions: fold change > 2 and P < 0.001. The R language (version 3.6.3) and R package “ggplot2” were performed to analyze and visualize the RNAseq data of TCGA processed uniformly by the toiling process (29). The CPTAC analysis of the UALCAN portal (http://ualcan.path.uab.edu/analysis-prot.html) was used to identify the protein expression level of GPX8 in six cancers, including breast cancer, colon cancer, clear cell RCC (Renal cell carcinoma), LUAD (Lung adenocarcinoma), UCEC (Uterine corpus endometrial carcinoma) and ovarian cancer. The Human Protein Atlas (HPA) (http://www.proteinatlas.org/) were employed to confirm the expression of GPX8 in normal or tumor tissue of colon, kidney, lung, endometrium.

Survival and ROC Curve Analysis

The TCGA datasets were employed to investigate the correlation of GPX8 expression with the prognosis of different cancers. The Kaplan-Meier Plotter (https://kmplot.com/analysis/) was performed to analyze the correlation between GPX8 expression and the Overall Survival and Disease Specific Survival. The PrognoScan (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html) was adopted to analyze the correlation between GPX8 expression and survival in GSE4412-GPL97, GSE4271-GPL97, GSE17536 and GSE31210. The R package “pROC” package and “ggplot2” were used to analyze and visualize the RNAseq data of TCGA and GTEx. The area under the ROC curves (AUC) was calculated for the determination of diagnosis and prognosis.

DNA Methylation and Genetic Alteration Analysis

To analyze and visualize the RNAseq and Methylation450 data of TCGA, the R language and package “ggplot2” were used. Moreover, the genetic alteration of GPX8 from TCGA cancers was explored in cBioPortal (https://www.cbioportal.org/) and displayed as alteration frequency.

Immune Infiltration Analysis

The “Gene_Corr” module of TIMER2 was used to explore the correlation between GPX8 and cancer-associated fibroblasts. The single sample gene set enrichment analysis (ssGSEA) of R package “GSVA” (30) was used to analyze the correlation between GPX8 and Immune Cells, including aDC [activated dendritic cells], B cells, CD8 T cells, Cytotoxic cells, DC [dendritic cells]; Eosinophils, iDC [immature dendritic cells], Macrophages, Mast cells, Neutrophils, NK CD56bright cells, NK CD56dim cells, NK cells, pDC [Plasmacytoid dendritic cells], T cells, T helper cells, Tcm [T central memory], Tem [T effector memory], Tfh [T follicular helper], Tgd [T gamma delta], Th1 cells, Th17 cells, Th2 cells and Treg.

GPX8-Related Gene Enrichment Analysis

The KEGG website (https://www.kegg.jp/) was used for pathway analysis. The STRING website (https://string-db.org/) and GEPIA (http://gepia.cancer-pku.cn/index.html) were performed to predict GPX8-binding proteins. The R package “ggplot2” and “clusterProfiler” (31) were used to analyze and visualize the GO and KEGG analysis. The GBM/LGG, KIRC, KIRC and STAD project RNAseq data from TCGA (https://portal.gdc.cancer.gov/) were applied to analyze the co-expression among GPX8, IKBIP, SERPINH1, PPIC, OSTC, TNFAIP8 and CRTAP.

Patients and Tissue Samples

GBM tissues and paracancerous ones were surgically resected from 3 patients undergoing surgery or biopsy, which were pathologically confirmed and stored at −80°C. This study got approval from the Ethics Committee of Guizhou Provincial People’s Hospital and was conducted after obtaining written informed consent of each subject.

Cell Culture and Transfection

A normal human astrocyte cell line (HA) was received from ScienCell (Carlsbad, USA) and maintained in Astrocyte Medium. A GBM cell line LN-299 was received from ATCC (Rockville, USA). The other two GBM cell lines A172 and U251 were received from the cell bank of the Chinese Academy of Science (Shanghai, China). These cells were maintained in DMEM medium with 10% FBS and cultured in suitable conditions (37°C; 5% CO2). All shRNAs (sh-NC,sh-GPX8-1 and sh-GPX8-2) were purchased from General Biosystems (Anhui, China) and transfected to cells by Lipofectamine 3000 (Invitrogen, USA). The Sequences of shRNAs were displayed in Table S1.

Quantitative Reverse Transcription-PCR (qRT-PCR)

Cells were treated with Trizol (Takara, Japan). Total RNA was reverse transcribed by the RT kit (Takara, Japan). SYBR Green PCR kit (Takara, Japan) was performed for qRT-PCR. The relative expression of GPX8 was analyzed by the 2^-△△Ct method and normalized with GAPDH. The primer sequences were exhibited in Table S2.

Western Blot Assay (WB)

Cells and tissue samples were collected for protein extraction. The BCA Protein assay kit (Beyotime, China) was used for protein concentration evaluation. Proteins were separated with a sodium dodecyl sulfate–polyacrylamide gel electrophoresis gel. The polyvinylidene fluoride membrane containing proteins was blocked with 5% milk. Then, specific primary antibodies (GPX8, Abcam; GAPDH, Abcam) were applied to the membrane at 4°C overnight. After incubating with secondary antibodies, an ECL Western Blotting Substrate (Solarbio, China) was used to detect the protein blots.

Wound-Healing Assay

LN-299 and U251 cells were cultured in plates for 24h. The pipette tip was employed to draw the surface of the cell layer. The microscope was used to capture the images at 0 and 48h after injury. The distance of the injury area at 48h was measured. A relative migration rate was calculated by normalizing the distance of the injury area at 0h.

Transwell Assay

Transwell chamber with matrigel was performed to transwell assay. Cells were mixed with serum-free media and injected into the upper layer. Outside of the transwell chamber was full with complete medium. After cells were cultured in suitable conditions for 48h, the transwell chambers were taken for fixing and staining by 4% paraformaldehyde and 0.1% crystal violet (Solarbio, China). The cells on the bottom of the chambers were counted.

Statistical Analysis

SPSS 22.0 (IBM, USA) was used for statistical analyses. The data were analyzed with the Student’s t-test or one-factor analysis. All results were displayed as mean ± SD. P <0.05 was regarded as statistically significant.

Results

GPX8 Expression in Pan-Cancer

To ascertain differences in GPX8 expression between the tumor and normal tissues, several databases were used, including Oncomine, TCGA, GEO and others. By searching in the Oncomine, the expression level of the GPX8 gene was found significantly up-regulated in many types of cancers, including BRCA (Breast invasive carcinoma), CHOL(Cholangiocarcinoma), COAD (Colon adenocarcinoma), GBM(Glioblastoma multiforme), HNSC(Head and Neck squamous cell carcinoma), KIRC(Kidney renal clear cell carcinoma), KIRP(Kidney renal papillary cell carcinoma), LIHC(Liver hepatocellular carcinoma), LUAD(Lung adenocarcinoma), LUSC(Lung squamous cell carcinoma), STAD(Stomach adenocarcinoma) (Figures 1A, B). The tumor and adjacent normal tissues for GPX8 across all TCGA tumors were shown in Figure 1C. After including the normal tissue of the GTEx dataset as controls, we further confirmed the GPX8 expression was significantly higher in BRCA, COAD, GBM/LGG, HNSC, KIRC, KIRP, LUAD and STAD compared with normal tissues (Figure 1D). Moreover, the result of CPTAC analysis demonstrated the protein expression of GPX8 was up-regulated and correlated with the pathological stages of colon cancer, clear cell RCC, LUAD and UCEC (Figures 2A, B). The data of the HPA database showed the different expression levels of GPX8 in normal or tumor tissue of the colon, kidney, lung, endometrium (Figure 2C). In a word, GPX8 was highly expressed in most cancers.

FIGURE 1

Figure 1 The expression of GPX8 gene in different cancers. (A) The expression of GPX8 in pan-cancer was explored in the Oncomine database. (www.oncomine.org). (B) The expression of GPX8 in pan-cancer was explored in TCGA through TIMER. (http://timer.cistrome.org/). (C) The expression of GPX8 in pan-cancer was analyzed using paired tumor/normal samples from TCGA and GTEx databases. (D) The expression of the GPX8 was analyzed in BRCA, COAD, GBM/LGG, HNSC, KIRC, KIRP, LUAD and STAD. *P < 0.05, **P < 0.005, ***P < 0.001. “ns” was regarded as no statistically significant.

FIGURE 2

Figure 2 The expression of GPX8 protein in different cancers. (A) Based on the CPTAC analysis (http://ualcan.path.uab.edu/analysis-prot.html), the expression of GPX8 protein was analyzed in Breast cancer, colon cancer, clear cell RCC, Lung adenocarcinoma, UCEC, Ovarian cancer. (B) The expression of GPX8 protein was analyzed by the main pathological stages of colon cancer, clear cell RCC, Lung adenocarcinoma, UCEC. (C) Based on the HPA database (http://www.proteinatlas.org/), the expression of GPX8 protein in normal or tumor tissue of colon, kidney, lung, the endometrium was displayed. ***P < 0.001.

Prognostic Potential of GPX8 in Cancers

Based on the TCGA databases, the results of survival analysis demonstrated that the GPX8 expression significantly affects prognosis in GBM/LGG, KIRC, KIRP, LUAD and STAD (Figures 3A, B). Moreover, the GEO datasets, including GSE4412-GPL97, GSE4271-GPL97, GSE17536 and GSE31210, showed that high GPX8 expression was correlated with poorer prognosis of patients with Brain cancer, Colorectal cancer and Lung cancer (Figure 3C). Furthermore, the receiver operating characteristic (ROC) curve was emplored to verify the diagnostic value of GPX8 in different cancers. As shown in Figure 3D, the GPX8 had a moderate diagnostic accuracy of BRCA, GBM/LGG, HNSC, KIRC, KIRP and STAD (AUCs were above 0.7 and even 0.8). To summarize, the Survival and ROC curve analyses indicated that GPX8 was a valuable diagnostic biomarker in multiple types of cancers, including GBM/LGG, KIRC, KIRP and STAD.

FIGURE 3

Figure 3 Survival prognosis analysis of cancers. (A) Overall Survival analysis of GPX8 genes in the TCGA dataset. (https://kmplot.com/analysis/). (B) Disease Specific Survival analysis of GPX8 genes in the TCGA dataset. (https://kmplot.com/analysis/). (C) Survival analysis of GPX8 genes in GSE4412-GPL97, GSE4271-GPL97, GSE17536 and GSE31210 (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html). (D) ROC analysis of GPX8 genes in TCGA database.

Correlation Between GPX8 Expression and Critical Clinical/Pathological Features

Table 1 provided the relations between clinical data and GPX8 expression in GBM/LGG, KIRC, KIRP and STAD. High GPX8 expression was correlated with the WHO grade and age of GBM/LGG, and the T stage and Pathologic stage of KIRC and STAD. However, the high expression of GPX8 was not significantly correlated with other clinical features of KIRP. Univariate analysis demonstrated that the OS was correlated with the WHO grade and T stage of GBM/LGG, the N stage, M stage and the Pathologic T stage of KIRC and KIRP and the N stage and M stage of STAD. In multivariate analysis, the OS was associated with the WHO grade of GBM/LGG, the M stage of KIRC, KIRP and STAD and the N stage of STAD (Table 2).

TABLE 1

Table 1 Correlation of the GPX8 expression and the clinical characteristic of patients with cancers.

TABLE 2

Table 2 Univariate and multivariate analyses of overall survival.

The Expression of GPX8 Was Correlated to DNA Methylation and Genetic Alteration

Through exploring in the RNAseq and Methylation450 data of TCGA, a negative correlation between GPX8 expression and GPX8 DNA methylation of the promoter region was observed in GBM/LGG, KIRC and STAD (Figure 4A). However, the different cancer types showed a different correlation coefficient. For instance, the correlation coefficient of two probes (cg08212230 and cg08508814) were -0.480 and -0.5 in GBM/LGG and were 0.15 and 0.11 in KIRC. Furthermore, a genetic alteration analysis showed the alteration frequency of GPX8 was >4%, and the primary type was amplification (Figure 4B). Moreover, based on the cbioportal database, the missense mutation of GPX8 has been found as the main type of genetic alteration, and S61Y/F alteration was the important site, which was detected in 2 cases of COAD, 1 case of BLCA and 1 case of UCEC (Figure 4C). The S61Y/F site in the 3D structure of GPX8 protein was shown in Figure 4D.

FIGURE 4

Figure 4 DNA methylation analysis and Mutation feature of GPX8 in cancers. (A) DNA methylation of GPX8 in BM/LGG, KIRC, KIRP and STAD. (B) The alteration frequency and mutation type of GPX8. (https://www.cbioportal.org/). (C) The alteration frequency with mutation type and mutation site are displayed. (https://www.cbioportal.org/). (D) The mutation site with the highest alteration frequency (S61Y/F) in the 3D structure of GPX8. (https://www.cbioportal.org/).

Immune Infiltration Analysis

Subsequently, the results of TIMER, TIDE, XCEL, MCPCOUNTER, EPIC algorithms indicated a positive correlation between the GPX8 expression and cancer-associated fibroblast in many human cancers, including GBM, LGG, KIRC, KIRP and STAD (Figure 5A). Moreover, based on the TCGA database, the result of R package “GSVA” was demonstrated the positive correlation of GPX8 and the different immune cells was observed in these cancers (Figures 5B–E), and the detail was displayed in Table 3. Furthermore, the scatterplot data from the TIMER database indicated that the GPX8 expression level was correlated with immune cells (Figure 6). To sum up, the expression of GPX8 was related to Immune infiltration.

FIGURE 5

Figure 5 Correlation analysis between GPX8 expression and tumor-infiltrating immune cells. (A) Correlation analysis between GPX8 expression and CAFs in pan-cancer (http://timer.comp-genomics.org/). (B–E) Correlation analysis between GPX8 expression and immune cells in GBM/LGG, KIRC, KIRP, STAD.

TABLE 3

Table 3 Correlation analysis between GPX8 and immune cells.

FIGURE 6

Figure 6 Correlation analysis between GPX8 expression and immune cells in TIMER (http://timer.cistrome.org/).

Enrichment Analysis of GPX8-Related Partners

The GPX8-binding PPI networks were created by the STRING online database (Figure 7A). The GO and KEGG analysis demonstrated that GPX8 and GPX8-binding proteins were mainly involved in Glutathione metabolism (Figure 7B). The relationship of these pathways was displayed (Figure 7C). The GBMLGG, KIRC, KIRC and STAD project RNAseq data from TCGA demonstrated that GPX8 was co-expression among IKBIP, SERPINH1, PPIC, OSTC, TNFAIP8 and CRTAP (Figure 7D).

FIGURE 7

Figure 7 The analysis of GPX8 and GPX8-related partners. (A) The GPX8-binding PPI analysis of GPX8. (https://string-db.org/, http://gepia.cancer-pku.cn/index.html). (B) The GO and KEGG analysis of GPX8 and GPX8-related partners. (https://www.kegg.jp/). (C) The relationship of pathways. (D) The co-expression analysis of GPX8 and GPX8-related genes in GBM/LGG, KIRC, KIRP and STAD. (https://portal.gdc.cancer.gov/). ***P < 0.001.

GPX8 Was Highly Expressed and Influenced the Migration and Invasion of GBM Cells

To confirm the function of GPX8, we selected GBM as a typical type of human cancer for further investigation. By comparing with the normal human astrocyte cell line, the result of qRT-PCR and WB demonstrates that GPX8 was highly expressed in GBM cells, including LN-299, A172 and U251 (Figure 8A). Then, we designed two shRNAs that targeted different sites of GPX8. The result of qRT-PCR and WB indicated both two sh-RNAs have excellent efficiency in knockdown of the GPX8 expression (Figure 8B). Subsequently, we observed the knockdown of GPX8 could inhibit the migration and invasion of GBM cells (Figures 8C, D). To confirm the expression of GPX8 in primary tumors, we collected the tumor tissues and paracancerous ones from patients with GBM. By comparing with the normal tissue of the same patient, the result of WB implied that GPX8 was highly expressed in GBM tissue (Figure 8E). The clinical parameters of the patients were displayed in Figure 8F.

FIGURE 8

Figure 8 GPX8 is up-regulated in GBM cells and related to migration and invasion in vitro. (A) qRT-PCR and WB confirmed the expression of GPX8 was higher in GBM cells. (B) qRT-PCR and WB verified the efficiency of sh-RNAs. (C) Wound-healing assays indicated the knockdown of GPX8 could restrain the migration of GBM cells. (D) Transwell assays showed the knockdown of GPX8 could inhibit the invasion of GBM cells. (E). WB confirmed the expression of GPX8 was higher in GBM tumor tissues than paracancerous ones. (F) The clinical parameters of the patients. Data were presented as mean ± SD. ***P < 0.001.

Discussion

In this work, we comprehensively investigated the molecular characteristics of the GPX8 gene in 33 different cancers from diverse databases, including TCGA, CPTAC and GEO, and aimed to explore the value of GPX8 in cancer prognosis, progression, and treatment.

Based on our results, compared to normal tissues, the mRNA of GPX8 was highly expressed in multiple types of cancer, especially in BRCA, COAD, GBM/LGG, HNSC, KIRC, KIRP, LUAD and STAD. Furthermore, the UALCAN database data confirmed that the protein expression level of GPX8 in tumor tissues of COAD, KIRC and LUAD were higher than normal, and the higher expression of GPX8 was related to the late clinical stage. Nevertheless, we failed to find significant differences between tumor and normal tissues in Breast cancer. Additionally, although the GPX8 protein was highly expressed in UCEC, the mRNA of GPX8 was lower than normal. As proteins could be modulated at multiple levels, the mRNA/protein expression is not always consistent (32). And the different data sources of mRNA and protein might lead to these differences.

Several studies have reported the GPX8 could be correlated with the poor prognostic of Gastric Cancer and Breast Cancer (24, 33). Herein, based on GEO, TCGA and GTEx data, we conducted a series of survival analyses which found the prognostic value of GPX8 seems to differ for different tumors. For GBM/LGG, KIRC, KIRP and STAD, the GPX8 expression was related to the prognosis of OS and DSS and had good diagnostic accuracy (AUC>0.7). For BRCA and HNSC, there was no significant correlation between GPX8 and OS or DSS. Thus, although the prognostic value of GPX8 in BRCA and HNSC is still debated, the current evidence could suggest the role of GPX8 in the clinical prognosis of GBM/LGG, KIRC, KIRP and STAD.

For cancers, DNA methylation is a crucial mechanism for epigenetic changes. Reports have demonstrated that abnormal DNA methylation has been involved in multiple types of tumorigenesis (34, 35). For example, hypermethylated FBXL7 is often observed in aggressive tumors (36). In our work, evidence suggested that the expression of GPX8 was correlated to the reduced DNA methylation at the promoter region in different types of tumors, especially for patients with GBM/LGG.

Reports have uncovered that the tumor microenvironment (TME), which consists of tumor cells and many other cells, plays a crucial role in cancer progression and could significantly influence clinical outcomes (37). As a considerable component of the TME, tumor-infiltrating immune cells are associated with the growth, invasion, and metastasis of nearly all tumors (2, 6). For example, the cancer-associated fibroblasts (CAFs), fibroblasts activated by tumors, play a crucial role in generating the interstitial matrix and contributing to advanced carcinomas’ desmoplastic stroma. Accumulating evidence has demonstrated that the CAFs were associated with poor prognosis in multiple types of cancer and could promote tumor aggressiveness and immune evasion by providing a cytokine and extracellular matrix (ECM) milieu of tumors (38–41). Furthermore, the multiple immune cells recruited by progressed tumors could stimulate an inflammatory environment and promote the secretion of cytokines and chemokines, which contribute to tumor growth, metastasis, invasion and pathological angiogenesis (10–12). Herein, our result demonstrated that the GPX8 expression might be related to cancer-associated fibroblasts and immune infiltration levels in certain tumors. Our study first implied the GPX8 expression might associate with the TME, especially the immune infiltration.

Additionally, through the enrichment analyses of GPX8-related genes, we recognized the “Glutathione Metabolism”, “Antioxidant Activity” and “Antioxidant Activity” may be the potential mechanism for GPX8 to influence the etiology(ies) and pathogenesis of cancers. Reports have demonstrated that oxidative stress is a critical metabolic feature in the TME due to stromal inflammatory cell recruitment, can promote the function of CAFs (42–44). As an important gene for preventing ROS accumulation, the expression of GPX8 may contribute to this progress and then influence the clinical outcomes.

Furthermore, we confirmed GPX8 was highly expressed in GBM cells and contributed to migration and invasion. Chen et al. had reported that the over-expression of GPX8 was observed in gastric cancer and could promote the malignant behavior of gastric cancer cells (23). Another report had demonstrated the GPX8 was up-regulated in breast cancer and associated with epithelial-mesenchymal transition (24). However, our study first indicated that GPX8 plays a cancer-promoting role in GBM cells, especially in migrative and invasive phenotypes. Reports revealed that fibrosis is driven by ROS and cytokines induced by immune cells (45, 46). The complex TME may be critical for invasion and metastasis (47, 48). Interestingly, all of these characteristics above are related to GPX8. It implied that GPX8 might be a potential target for the diagnosis and treatment of cancers.

Overall, our pan-cancer analyses of GPX8 have clarified the landscape of GPX8 expression with the clinical prognosis and immune cell infiltration, which provided a predictive biomarker in pan-cancer, especially in the GBM/LGG, KIRC, KIRP and STAD. In particular, we first confirmed the up-regulated of GPX8 in GBM cells and observed the suppression of migrative and invasive phenotypes by knockdown of GPX8. Guided by the bioinformatics analysis, our study further revealed the mechanisms of the diagnosis and prognosis of cancers. However, although our findings have indicated the correlation between the expression of GPX8 and TME, more experiments are required to explore, including the role of GPX8 in cell proliferation, apoptosis and the molecular mechanisms underlying. Furthermore, the validation of GPX8 expression in more tumor samples needs to be further studied, as well as the relationship between GPX8 expression and the survival curve of patients with GBM. Moreover, as GPX8, which does not have GSH binding domain, is reported to contain a weak GPX activity, the precise biological function of GPX8 remains to be explored (21, 49).

Conclusion

Our study showed a correlation of GPX8 expression with clinical prognosis, DNA methylation and immune infiltrates. Furthermore, we first confirmed GPX8 was highly expressed in GBM cells and contributes to migration and invasion. These results provided a predictive biomarker and an inclusive understanding of the GPX8 expression in multiple tumor types, especially in GBM.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: TIMER database (http://timer.cistrome.org/), ONCOMINE database (www.oncomine.org), The CPTAC analysis of the UALCAN portal (http://ualcan.path.uab.edu/analysis-prot.html), The Human Protein Atlas (HPA) (http://www.proteinatlas.org/), The Kaplan-Meier Plotter (https://kmplot.com/analysis/), The PrognoScan (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html), GSE4412-GPL97, GSE4271-GPL97, GSE17536, GSE17536 and GSE31210. The ONCOMINE database is no longer available, any data request should be directed to the corresponding authors.

Author Contributions

XT, ZY, DW, and ZR conceived the project. ZR, YH, QY, and JG performed the data analysis. BL and ZY participated in writing the manuscript. ZR and HH edited the language. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by National Natural Science Foundation of China (31660265, 31960208, 81060145, and 81560356), Youth Fund of Guizhou Provincial People’s Hospital (GZSYQN[2015]06 and GZSYQN[2019]15), Subsidy Foundation of National Natural Science Foundation of Guizhou Provincial People’s Hospital (Guizhou Science and Technology Platform (2017)5724), Science and Technology Foundation of Guizhou Province (Guizhou Science and Technology J Word (2015)2096), Subsidy Foundation of National Natural Science Foundation of Guizhou Provincial People’s Hospital(GPPH-NSFC-2019-1).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.812811/full#supplementary-material

Supplementary Table 1 | Primer sequences for qRT-PCR.

Supplementary Table 2 | Sequences of shRNAs.

References

1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2020. CA Cancer J Clin (2020) 70(1):7–30. doi: 10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Steven A, Seliger B. The Role of Immune Escape and Immune Cell Infiltration in Breast Cancer. Breast Care (Basel) (2018) 13(1):16–21. doi: 10.1159/000486585

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Feng RM, Zong YN, Cao SM, Xu RH. Current Cancer Situation in China: Good or Bad News From the 2018 Global Cancer Statistics? Cancer Commun (Lond) (2019) 39(1):22. doi: 10.1186/s40880-019-0368-6

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Ramazzotti D, Lal A, Wang B, Batzoglou S, Sidow A. Multi-Omic Tumor Data Reveal Diversity of Molecular Mechanisms That Correlate With Survival. Nat Commun (2018) 9(1):4453. doi: 10.1038/s41467-018-06921-8

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Repana D, Nulsen J, Dressler L, Bortolomeazzi M, Venkata SK, Tourna A, et al. The Network of Cancer Genes (NCG): A Comprehensive Catalogue of Known and Candidate Cancer Genes From Cancer Sequencing Screens. Genome Biol (2019) 20(1):1. doi: 10.1186/s13059-018-1612-0

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Cui X, Zhang X, Liu M, Zhao C, Zhang N, Ren Y, et al. A Pan-Cancer Analysis of the Oncogenic Role of Staphylococcal Nuclease Domain-Containing Protein 1 (SND1) in Human Tumors. Genomics (2020) 112(6):3958–67. doi: 10.1016/j.ygeno.2020.06.044

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Li Z, Chen S, Feng W, Luo Y, Lai H, Li Q, et al. A Pan-Cancer Analysis of HER2 Index Revealed Transcriptional Pattern for Precise Selection of HER2-Targeted Therapy. EBioMedicine (2020) 62:103074. doi: 10.1016/j.ebiom.2020.103074

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Pan JH, Zhou H, Cooper L, Huang JL, Zhu SB, Zhao XX, et al. LAYN Is a Prognostic Biomarker and Correlated With Immune Infiltrates in Gastric and Colon Cancers. Front Immunol (2019) 10:6. doi: 10.3389/fimmu.2019.00006

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Castro A, Pyke RM, Zhang X, Thompson WK, Day CP, Alexandrov LB, et al. Strength of Immune Selection in Tumors Varies With Sex and Age. Nat Commun (2020) 11(1):4128. doi: 10.1038/s41467-020-17981-0

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hamarsheh S, Gross O, Brummer T, Zeiser R. Immune Modulatory Effects of Oncogenic KRAS in Cancer. Nat Commun (2020) 11(1):5439. doi: 10.1038/s41467-020-19288-6

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Benner B, Scarberry L, Suarez-Kelly LP, Duggan MC, Campbell AR, Smith E, et al. Generation of Monocyte-Derived Tumor-Associated Macrophages Using Tumor-Conditioned Media Provides a Novel Method to Study Tumor-Associated Macrophages In Vitro. J Immunother Cancer (2019) 7(1):140. doi: 10.1186/s40425-019-0622-0

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Poh AR, Love CG, Masson F, Preaudet A, Tsui C, Whitehead L, et al. Inhibition of Hematopoietic Cell Kinase Activity Suppresses Myeloid Cell-Mediated Colon Cancer Progression. Cancer Cell (2017) 31(4):563–575 e565. doi: 10.1016/j.ccell.2017.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sade-Feldman M, Yizhak K, Bjorgaard SL, Ray JP, de Boer CG, Jenkins RW, et al. Defining T Cell States Associated With Response to Checkpoint Immunotherapy in Melanoma. Cell (2018) 175(4):998–1013.e1020. doi: 10.1016/j.cell.2018.10.038

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Callahan MK, Postow MA, Wolchok JD. Targeting T Cell Co-Receptors for Cancer Therapy. Immunity (2016) 44(5):1069–78. doi: 10.1016/j.immuni.2016.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

15. McDermott DF, Huseni MA, Atkins MB, Motzer RJ, Rini BI, Escudier B, et al. Clinical Activity and Molecular Correlates of Response to Atezolizumab Alone or in Combination With Bevacizumab Versus Sunitinib in Renal Cell Carcinoma. Nat Med (2018) 24(6):749–57. doi: 10.1038/s41591-018-0053-3

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhang Z, Li J, He T, Ding J. Bioinformatics Identified 17 Immune Genes as Prognostic Biomarkers for Breast Cancer: Application Study Based on Artificial Intelligence Algorithms. Front Oncol (2020) 10:330. doi: 10.3389/fonc.2020.00330

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Bansal A, Simon MC. Glutathione Metabolism in Cancer Progression and Treatment Resistance. J Cell Biol (2018) 217(7):2291–8. doi: 10.1083/jcb.201804161

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Brigelius-Flohe R, Flohe L. Regulatory Phenomena in the Glutathione Peroxidase Superfamily. Antioxid Redox Signal (2020) 33(7):498–516. doi: 10.1089/ars.2019.7905

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Traverso N, Ricciarelli R, Nitti M, Marengo B, Furfaro AL, Pronzato MA, et al. Role of Glutathione in Cancer Progression and Chemoresistance. Oxid Med Cell Longev (2013) 2013:972913. doi: 10.1155/2013/972913

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Jiao Y, Wang Y, Guo S, Wang G. Glutathione Peroxidases as Oncotargets. Oncotarget (2017) 8(45):80093–102. doi: 10.18632/oncotarget.20278

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Nguyen VD, Saaranen MJ, Karala AR, Lappi AK, Wang L, Raykhel IB, et al. Two Endoplasmic Reticulum PDI Peroxidases Increase the Efficiency of the Use of Peroxide During Disulfide Bond Formation. J Mol Biol (2011) 406(3):503–15. doi: 10.1016/j.jmb.2010.12.039

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ramming T, Hansen HG, Nagata K, Ellgaard L, Appenzeller-Herzog C. GPx8 Peroxidase Prevents Leakage of H2O2 From the Endoplasmic Reticulum. Free Radical Biol Med (2014) 70:106–16. doi: 10.1016/j.freeradbiomed.2014.01.018

CrossRef Full Text | Google Scholar

23. Chen H, Xu L, Shan ZL, Chen S, Hu H. GPX8 Is Transcriptionally Regulated by FOXC1 and Promotes the Growth of Gastric Cancer Cells Through Activating the Wnt Signaling Pathway. Cancer Cell Int (2020) 20(1):596. doi: 10.1186/s12935-020-01692-z

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Khatib A, Solaimuthu B, Ben Yosef M, Abu Rmaileh A, Tanna M, Oren G, et al. The Glutathione Peroxidase 8 (GPX8)/IL-6/STAT3 Axis Is Essential in Maintaining an Aggressive Breast Cancer Phenotype. Proc Natl Acad Sci USA (2020) 117(35):21420–31. doi: 10.1073/pnas.2010275117

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhang J, Liu Y, Guo Y, Zhao Q. GPX8 Promotes Migration and Invasion by Regulating Epithelial Characteristics in Non-Small Cell Lung Cancer. Thorac Cancer (2020) 11(11):3299–308. doi: 10.1111/1759-7714.13671

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Mimura K, Kua LF, Shimasaki N, Shiraishi K, Nakajima S, Siang LK, et al. Upregulation of Thioredoxin-1 in Activated Human NK Cells Confers Increased Tolerance to Oxidative Stress. Cancer Immunol Immunother (2017) 66(5):605–13. doi: 10.1007/s00262-017-1969-z

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Deng H, Yang W, Zhou Z, Tian R, Lin L, Ma Y, et al. Targeted Scavenging of Extracellular ROS Relieves Suppressive Immunogenic Cell Death. Nat Commun (2020) 11(1):4951. doi: 10.1038/s41467-020-18745-6

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ou W, Jiang L, Thapa RK, Soe ZC, Poudel K, Chang JH, et al. Combination of NIR Therapy and Regulatory T Cell Modulation Using Layer-by-Layer Hybrid Nanoparticles for Effective Cancer Photoimmunotherapy. Theranostics (2018) 8(17):4574–90. doi: 10.7150/thno.26758

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Vivian J, Rao AA, Nothaft FA, Ketchum C, Armstrong J, Novak A, et al. Toil Enables Reproducible, Open Source, Big Biomedical Data Analyses. Nat Biotechnol (2017) 35(4):314–6. doi: 10.1038/nbt.3772

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Hanzelmann S, Castelo R, Guinney J. GSVA: Gene Set Variation Analysis for Microarray and RNA-Seq Data. BMC Bioinform (2013) 14:7. doi: 10.1186/1471-2105-14-7

CrossRef Full Text | Google Scholar

31. Yu G, Wang LG, Han Y, He QY. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Kustatscher G, Grabowski P, Schrader TA, Passmore JB, Schrader M, Rappsilber J. Co-Regulation Map of the Human Proteome Enables Identification of Protein Functions. Nat Biotechnol (2019) 37(11):1361–71. doi: 10.1038/s41587-019-0298-5

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Zhang X, Zhan D, Li Y, Wang H, Cheng C, Yao Y, et al. Glutathione Peroxidase 8 as a Prognostic Biomarker of Gastric Cancer: An Analysis of The Cancer Genome Atlas (TCGA) Data. Med Sci Monit (2020) 26:e921775. doi: 10.12659/MSM.921775

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Gkountela S, Castro-Giner F, Szczerba BM, Vetter M, Landin J, Scherrer R, et al. Circulating Tumor Cell Clustering Shapes DNA Methylation to Enable Metastasis Seeding. Cell (2019) 176(1-2):98–112 e114. doi: 10.1016/j.cell.2018.11.046

PubMed Abstract | CrossRef Full Text | Google Scholar

35. So JY, Skrypek N, Yang HH, Merchant AS, Nelson GW, Chen WD, et al. Induction of DNMT3B by PGE2 and IL6 at Distant Metastatic Sites Promotes Epigenetic Modification and Breast Cancer Colonization. Cancer Res (2020) 80(12):2612–27. doi: 10.1158/0008-5472.CAN-19-3339

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Moro L, Simoneschi D, Kurz E, Arbini AA, Jang S, Guaragnella N, et al. Epigenetic Silencing of the Ubiquitin Ligase Subunit FBXL7 Impairs C-SRC Degradation and Promotes Epithelial-to-Mesenchymal Transition and Metastasis. Nat Cell Biol (2020) 22(9):1130–42. doi: 10.1038/s41556-020-0560-6

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Jabeen S, Zucknick M, Nome M, Dannenfelser R, Fleischer T, Kumar S, et al. Serum Cytokine Levels in Breast Cancer Patients During Neoadjuvant Treatment With Bevacizumab. Oncoimmunology (2018) 7(11):e1457598. doi: 10.1080/2162402X.2018.1457598

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Miao L, Liu Q, Lin CM, Luo C, Wang Y, Liu L, et al. Targeting Tumor-Associated Fibroblasts for Therapeutic Delivery in Desmoplastic Tumors. Cancer Res (2017) 77(3):719–31. doi: 10.1158/0008-5472.CAN-16-0866

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Hu G, Zhong K, Chen W, Wang S, Huang L. Podoplanin-Positive Cancer-Associated Fibroblasts Predict Poor Prognosis in Lung Cancer Patients. Onco Targets Ther (2018) 11:5607–19. doi: 10.2147/OTT.S175566

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Duperret EK, Trautz A, Ammons D, Perales-Puchalt A, Wise MC, Yan J, et al. Alteration of the Tumor Stroma Using a Consensus DNA Vaccine Targeting Fibroblast Activation Protein (FAP) Synergizes With Antitumor Vaccine Therapy in Mice. Clin Cancer Res (2018) 24(5):1190–201. doi: 10.1158/1078-0432.CCR-17-2033

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Curtis M, Kenny HA, Ashcroft B, Mukherjee A, Johnson A, Zhang Y, et al. Fibroblasts Mobilize Tumor Cell Glycogen to Promote Proliferation and Metastasis. Cell Metab (2019) 29(1):141–155 e149. doi: 10.1016/j.cmet.2018.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Maj T, Wang W, Crespo J, Zhang H, Wang W, Wei S, et al. Oxidative Stress Controls Regulatory T Cell Apoptosis and Suppressor Activity and PD-L1-Blockade Resistance in Tumor. Nat Immunol (2017) 18(12):1332–41. doi: 10.1038/ni.3868

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Parri M, Chiarugi P. Redox Molecular Machines Involved in Tumor Progression. Antioxid Redox Signal (2013) 19(15):1828–45. doi: 10.1089/ars.2012.5040

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Cirri P, Chiarugi P. Cancer-Associated-Fibroblasts and Tumour Cells: A Diabolic Liaison Driving Cancer Progression. Cancer Metastasis Rev (2012) 31(1-2):195–208. doi: 10.1007/s10555-011-9340-x

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Morry J, Ngamcherdtrakul W, Yantasee W. Oxidative Stress in Cancer and Fibrosis: Opportunity for Therapeutic Intervention With Antioxidant Compounds, Enzymes, and Nanoparticles. Redox Biol (2017) 11:240–53. doi: 10.1016/j.redox.2016.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Fielding CA, Jones GW, McLoughlin RM, McLeod L, Hammond VJ, Uceda J, et al. Interleukin-6 Signaling Drives Fibrosis in Unresolved Inflammation. Immunity (2014) 40(1):40–50. doi: 10.1016/j.immuni.2013.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Wu HW, Sheard MA, Malvar J, Fernandez GE, DeClerck YA, Blavier L, et al. Anti-CD105 Antibody Eliminates Tumor Microenvironment Cells and Enhances Anti-GD2 Antibody Immunotherapy of Neuroblastoma With Activated Natural Killer Cells. Clin Cancer Res (2019) 25(15):4761–74. doi: 10.1158/1078-0432.CCR-18-3358

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Zanotelli MR, Rahman-Zaman A, VanderBurgh JA, Taufalele PV, Jain A, Erickson D, et al. Energetic Costs Regulated by Cell Mechanics and Confinement Are Predictive of Migration Path During Decision-Making. Nat Commun (2019) 10(1):4185. doi: 10.1038/s41467-019-12155-z

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Herbette S, Roeckel-Drevet P, Drevet JR. Seleno-Independent Glutathione Peroxidases. More Than Simple Antioxidant Scavengers. FEBS J (2007) 274(9):2163–80. doi: 10.1111/j.1742-4658.2007.05774.x

PubMed Abstract | CrossRef Full Text | Google Scholar