- 1Department of Hepatobiliary Surgery, The Affiliated Hospital of Guizhou Medical University, Guiyang, China
- 2School of Clinical Medicine, Guizhou Medical University, Guiyang, China
- 3Hunan YoBon Biotechnology Limited Company, Changsha, China
- 4Department of Hepatobiliary Surgery, Renmin Hospital of Wuhan University, Wuhan, China
Hepatocellular carcinoma (HCC) is one of the most common gastrointestinal malignancies. It is not easy to be diagnosed in the early stage and is prone to relapse, with a very poor prognosis. And immune cell infiltration and tumor microenvironment play important roles in predicting therapeutic response and prognosis of HCC. Machado-Joseph domain-containing proteases (MJDs), as a gene family extensively involved in tumor progression, has pro-cancer and anti-cancer effects. However, the relationship between MJDs family members and immune cell infiltration and tumor microenvironment in HCC remains unclear. Therefore, cBio Cancer Genomics Portal (cBioPortal), The Cancer Genome Atlas (TCGA), UALCAN, Human Protein Atlas (HPA), MethSurv, and Tumor Immune Estimation Resource (TIMER) databases were performed to investigate the mRNA expression, DNA methylation, clinicopathologic features, immune cell infiltration and other related functions of MJDs family members in HCC. The results indicated that the expression of ATXN3, JOSD1, and JOSD2 was dramatically increased in HCC tissues and cell lines, and was correlated with histological grade, specimen type, TP53 mutation, lymph node metastatic, gender, and age of patients with HCC. Meanwhile, these genes also showed clinical value in improving the overall survival (OS), disease-specific survival (DSS), progression free survival (PFS), and relapse-free survival (RFS) in patients with HCC. The prognostic model indicated that the worse survival was associated with overall high expression of MJDs members. Next, the results suggested that promotor methylation levels of the MJDs family were closely related to these family mRNA expression levels, clinicopathologic features, and prognostic values in HCC. Moreover, the MJDs family were significantly correlated with CD4+ T cells, CD8+ T cells, B cells, neutrophils, macrophages, and DCs. And MJDs family members’ expression were substantially associated with the levels of several lymphocytes, immunomoinhibitors, immunomostimulators, chemokine ligands, and chemokine receptors. In addition, the expression levels of MJDs family were significantly correlated with cancer-related signaling pathways. Taken together, our results indicated that the aberrant expression of MJDs family in HCC played a critical role in clinical feature, prognosis, tumor microenvironment, immune-related molecules, mutation, gene copy number, and promoter methylation level. And MJDs family may be effective immunotherapeutic targets for patients with HCC and have the potential to be prognostic biomarkers.
Introduction
Hepatocellular carcinoma (HCC) is the fifth most common cancer, and is also the second causing reason for cancer-related deaths. According to National Cancer Center (NCC), HCC is the first cause of death for men aged 15–59 years in China (Zheng et al., 2022). In recent years, with the continuous improvement of diagnostic techniques and systemic treatment, the prognosis of HCC patients has been greatly improved. Unfortunately, the 5-year survival rate of HCC is still less than 19% (Feng et al., 2020). Thus, we need to further explore markers for early diagnosis and related therapeutic targets of HCC to improve the prognosis of patients.
The MJDs family, the sub-classes of deubiquitinases, consists of four members: ATXN3, ATXN3L, JOSD1, and JOSD2. The MJDs family members is of critical importance in tumorigenesis and tumor progression (Zeng et al., 2020). For example, downregulation of ATXN3 promotes sensitivity of neuroblastoma cells to Perifosine and MK-2206, but reduces sensitivity to etoposide and cisplatin in these cells (Gong et al., 2021); ATXN3 enhances breast cancer (BRCA) metastasis by deubiquitinating KLF4 (Zou et al., 2019); ATXN3 promotes mRNA expression of EIF5A2 by stabilizing EIF5A2 to decrease its ubiquitination and degradation in anaplastic thyroid carcinoma (Zhuang et al., 2021). Down-regulation of ATXN3L suppresses breast cancer cell proliferation by directly binding to KLF5 (Ge et al., 2015). Meanwhile, ATXN3L is also presented to promote the migration of NSCLC cells (Buus et al., 2009). JOSD1 suppresses mitochondrial apoptosis to induce the chemotherapy drug resistance in the gynaecological tumor by deubiquitinating and stabilizing MCL1 (Wu et al., 2020). In acute myeloid leukemia (AML) cells, JOSD1 interacts with JAK2-V617F and promotes its expression (Yang et al., 2022). JOSD2 deficiency inhibits tumor cell proliferation by reducing glycolysis, and its mRNA expression is related to the worst prognosis in NSCLC (Krassikova et al., 2021). JOSD2 directly binds to and reduces ubiquitination levels of CTNNB1, thus, augmenting Wnt pathway transduction in HCC. In addition, Over-expression of JOSD2 is positively associated with poor prognosis in HCC patients (Huang et al., 2022). However, the possible biological functions of MJDs family members in HCC and its prognostic value remain unclear.
In our study, we analyzed the mRNA expression alterations in MJDs family members, clinical characteristics, prognostic values, promoter methylation levels, MJDs members’ gene changes, their association with immune infiltration cells, and immunomodulators. The study shows the underlying biological functions and the prognostic values of MJDs family members, which will facilitate the early diagnosis, targeted therapy, and immunotherapy of HCC patients.
Materials and methods
UALCAN
UALCAN is a synthetically web-based database for analyzing TCGA data. In this study, HCC sample types, ages, genders, grades, tumor stages, lymph node metastasis TP53 mutations, and promotor methylation levels were from the UALCAN database.
GTEx
This database integrates multi-omics data based on various normal human tissues. According to the Genotype-Tissue Expression (GTEx) databases in normal tissues, the expression of MJDs family members were analyzed.
cBioPortal
cBioPortal is an integrated web tool to analyze the TCGA database. We used this tool to analyze the gene structure variation of MJDs family members, co-expression, correlation between methylation levels and mRNA expression levels.
MethSurv
MethSurv is a web tool to study the relationship between aggregation of methylation sites and patient survival. We used this tool to analyze the methylation levels of CpG island in members of the MJDs family and their relationship with patient survival.
Kaplan-Meier analysis
The association between the expressed level of MJDs family members and the outcome of HCC patients was analyzed with the Kaplan-Meier. In the study, we evaluated outcomes with HCC patients through means of disease-specific survival (DSS), overall survival (OS), and progression free survival (PFS) as well as relapse-free survival (RFS) curves.
Construction of prognostic signature model
Downloading the STAR-counts data and clinical details of patients with HCC from the TCGA database (https://portal.gdc.com), normalizing the data to log2 (TPM+1), and then keeping the data with samples of RNAseq data and clinical details allowed us to divide the patients into high-risk and low-risk groups for further analysis. The timeROC analysis was carried out to determine the accuracy of the prediction model, and the log rank was utilized to assess the KM survival analysis to compare the survival difference between the two groups mentioned above. The 10-fold cross-validation method and the least absolute shrinkage and selection operator (LASSO) regression algorithm were used to pick the features. R’s glmnet package was used for the analysis mentioned above. The log rank test and univariate Cox regression were used to provide p-values and hazard ratios (HR) with 95 percent confidence intervals (CI) for Kaplan-Meier curves. R software version 4.0.3 was used to carry out all of the aforementioned analysis techniques and R packages (R 4.0 Foundation for Statistical Computing, 2020). Statistics were deemed significant at p < 0.05.
GSCALite
The GSCALite is a versatile database analysis web tool that we use to analyze the signaling pathways activated or inhibited by MJDs family members.
TIMER immune score of MJDs
The TCGA database (https://portal.gdc.com), which contains RNAseq data, was used to gather the clinical data and related data for patients with liver cancer. Immunedeconv, a R software package that incorporates the TIMER algorithm, has been rigorously benchmarked and each approach has been proven to have its own distinct performance and advantages. Immunedeconv was used to reliably evaluate immunity scores. The ggplot2 and pheatmap R (v4.0.3) packages were used to produce the aforementioned findings.
TISIDB
TISIDB is an integrated database for the study of tumor-immune system interactions, allowing for the analysis and annotation of 10 types of data. The correlations between MJDs family members with lymphocytes and immunomodulators in HCC were analyzed by using the TISIDB database.
Immune checkpoints analysis
371 liver cancer patients’ RNAseq data and relative clinical data were gathered from the TCGA database (https://portal.gdc.com). Immune checkpoint genes include SIGLEC15, TIGIT, CD274, HAVCR2, PDCD1, CTLA4, LAG3, and PDCD1LG2. To track the expression of genes linked to immunological checkpoints, the expression values of these eight genes were identified. The R packages ggplot2 and pheatmap were used to produce the aforementioned findings.
Cell culture
Normal liver cell line HL-7702 and HCC cell lines (MHCC-97H, HepG2, PLC/PRF/5, Huh7, and HCCLM-3) were provided from the Department of Biliary-Pancreatic Surgery of Cell Line Resource, Tongji Hospital (Hubei, Wuhan, China). All HCC cell lines were cultured in DMEM (Meisen, Zhejiang, China) and HL-7702 was cultured in RPMI-1640 (Meisen, Zhejiang, China). Both DMEM and RPMI-1640 were added to 10% of decomplemented fetal bovine serum (Gibco,USA) and 1% of penicillin/streptomycin (BI, Israel). Moreover, the mentioned above cells were incubated at a suitable temperature (37°C) and 5% carbon dioxide (5% CO2).
RNA extraction and qRT-PCR assay
According to the protocol, total RNA was collected from HL-7702 and HCC cell lines using TRIzol reagent. The extracted total RNA was then subjected to RT-PCR using HiScript III RT SuperMix (Vazyme, Nanjing, China) and the corresponding primers to obtain complementary DNA (cDNA). Using the qPCR Master Mix (SYBR Green, Vazyme, Nanjing, China), qRT-PCR was performed in a CFX96 Real-Time PCR detection system (Bio-Rad, USA). To analyze the relative expression of genes was calculated using the 2−ΔΔCt method. The primer sequences of the genes were listed in Supplementary Table S1.
Human Protein Atlas (HPA)
The HPA database integrates transcriptomics and proteomics, which includes different profiles of normal and abnormal tissue types and various cell types. The immunohistochemical (IHC) staining images of HCC tissues and normal liver tissues of MJDs family members from the HPA database.
Results
Pan-cancer analysis of MJDs family members’ expression
To explore the expression levels of MJDs family members in HCC, we used the UALCAN database to find different transcriptional levels of MJDs members in pan-cancer. We identified three genes of the MJDs family with abnormal expression in HCC tissues. ATXN3, JOSD1, and JOSD2 were upregulated, while ATXN3L was only expressed in testicular tissue and testicular carcinoma tissue (Figure 1). The expressions of MJDs family members in normal tissues was then obtained from GTEx database, ATXN3L was only expressed in testicular tissue, and the other three genes were expressed in all tissues (Supplementary Figure S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g001.jpg)
FIGURE 1. (A–D) The expression levels of MJDs family members in different types of cancers. The mRNA expression levels of ATXN3, ATXN3L, JOSD1, and JOSD2 were compared with various TCGA cancers.
Correlation between MJDs family Members’ expression
Based on UALCAN database, the correlation analysis between MJDs family members’ expression and histological grades, sample types, TP53 mutations, lymph node metastatic status, genders, and years were further shown.
As shown in Figure 2A, mRNA expression levels of ATXN3 were significantly related to histological grades, sample types, TP53 mutations, lymph node metastatic status, genders, and years. Specifically, In HCC patients, the expression of ATXN3 showed an increasing tendency in stages 1–3. The expression of ATXN3 in HCC tissues was notably higher than that in normal liver tissues. HCC specimens with non-mutant TP53 and mutant TP53 had observably higher ATXN3 expression levels than the control group. In the lymph node metastasis group, the mRNA expression levels of ATXN3 in both N0 and N1 groups were observably higher than that in the normal control. The mRNA expression level of ATXN3 in N1 was higher than that in N0 without statistical difference. In gender groups, the expression level of ATXN3 in both male and female groups was higher than that in the control group, which was statistically significant. However, without statistical difference, the expression level of ATXN3 in females was higher than that in males. Furthermore, ATXN3 transcriptional levels were related to different age groups. As age increases, ATXN3 expression also enhanced differently in HCC patients. The mRNA expression of ATXN3 was significantly associated with different HCC subclasses, with hepatocellular carcinoma having the highest levels and fibrolamellar carcinoma showing the lowest levels.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g002.jpg)
FIGURE 2. (A–D) The transcription level of MJDs family members in the subgroup of HCC patients, and classified by according to histological grades, sample types, TP53 mutations, lymph node metastatic status, genders, years and histological subtypes.
Similarly, the UALCAN Database was used to explore the transcriptional levels of JOSD1 and JOSD2 in HCC. We also found that JOSD1 and JOSD2 expressions remarkably difference in varied sample types, histological grades, TP53 mutations, lymph node metastatic status, genders, and years of HCC (Figures 2B,C). Therefore, the expression of ATXN3, JOSD1, and JOSD2 can be used as potential diagnostic indicators of HCC.
Promoter methylation level of MJDs family members in HCC
To explore the mechanisms for the up-regulation of mRNA expression levels, we further used UALCAN database to analyze MJDs family members’ promoter methylation levels and their correlation with HCC patients of clinical characteristics. We found that the ATXN3 methylation β value of the HCC group and the normal group was less than 0.2, which was a completely unmethylated state. The results suggested the ATXN3 methylation levels in HCC tissues were lower than in corresponding normal groups, regardless of TP53 mutations, histological grades, lymph node metastatic, genders, and years (Figure 3A). This was the same as the up-regulated ATXN3 expression in HCC patients. In addition, we found that the ATXN3L and JOSD2 promoter methylation levels were significantly lower in HCC tissues than that in normal tissues, despite histological grades, TP53 mutations, lymph node metastatic, genders, and years (Figures 3B,D). Interestingly, we found that the methylation β value of ATXN3L in various subgroups of years, genders, histological grades, lymph node metastasis and TP53 mutation was greater than 0.6, presenting a complete methylation state. In age groups, the methylation level decreased with the increasing of age, but it was still in a complete methylation state. In tumor stages, β values of G1, G2, G3, and G4 decreased, indicating the down-regulation of methylation levels. However, JOSD1 methylation levels in HCC tissues was higher than that in corresponding normal tissues without statistic difference, in spite of histological grades, TP53 mutations, lymph node metastatic, genders, and years (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g003.jpg)
FIGURE 3. (A–D) Correlation analysis between methylation level of MJDs family members and HCC patients clinical characteristics. The MJDs methylation level in HCC patients classified by histological grades, sample types, TP53 mutations, lymph node metastatic status, genders, and years were shown.
Moreover, we evaluated the CpG sites of promoter methylation of MJDs family members in hepatocellular carcinoma. In the heat map of methylation levels at CpG sites, we found that ATXN3 had 12 CpG islands, among which cg12034871 and cg26081025 were completely methylated, and the rest were completely unmethylated (Supplementary Figure S2A). Meanwhile, we also found that ATXN3L had four CpG islands, with partial methylation and complete unmethylation in each sample (Supplementary Figure S2B). In addition, JOSD1 had 17 CpG islands among which cg03088955 was completely methylation (Supplementary Figure S2C). There were six CpG islands in JOSD2, and only cg18708810 and cg13521229 had complete and partial methylation cases, respectively (Supplementary Figure S2D).
Furthermore, we analyzed the correlations between the degree of methylation of ATXN3, JOSD1, and JOSD2 and the mRNA expression levels with the cBioPortal database. The results demonstrated that the Spearman coefficients of ATXN3, JOSD1, and JOSD2 were -0.31,-0.35, and -0.29, respectively, indicating negative association between the methylation levels and the mRNA transcriptional levels of these three genes in HCC patients (Figures 4A–C). Meanwhile, the correlations between the fraction genome alteration and the mRNA expression levels were analyzed, and we found that the Spearman coefficients of ATXN3 was -0.21, which was a negative correlation, while the Spearman coefficients of both JOSD1 and JOSD2 were 0.13 and 0.08, respectively, which were positive correlations (Figure 4 D, E, F). Thus, the results suggested that lower levels of ATXN3, JOSD1, and JOSD2 promoter methylation levels might lead to higher mRNA transcriptional levels of these genes in HCC. There may be other mechanisms of mRNA up-regulation of these three genes in HCC patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g004.jpg)
FIGURE 4. (A–C) The promoter methylation levels of MJDs family members were inversely correlated with the expression of MJDs family members (D) ATXN3 mRNA expression levels were positively associated to the fraction of the genome altered (E,F) The expressions of JOSD1 and JOSD2 were negatively related to the fraction of the genome altered.
Analysis of MJDs alterations in HCC
To explore the reasons for the upregulation of mRNA expression levels, we analyzed the structural variation of the MJDs family gene in HCC by the cBioPortal. MJDs family members have different trends in different HCC subclasses. The sub-types of with the highest proprotion of all alterations were Hepatocellular carcinoma plus Intrahepatic Cholangiocarcinoma, Hepatocellular Adenoma, and Hepatocellular carcinoma for ATXN3, ATXN3L, JOSD1, and JOSD2, respectively (Figures 5A–E). Figures 5F–I displayed the specific mutation sites in MJDs DNA sequences, in which the green dots represent missense mutations. Interestingly, Genetic alterations include missense mutation, amplification, deep deletion, and mRNA high or low (Figure 5J), among which 11% ATXN3, 2.5% ATXN3L, 6% JOSD1, and 4% JOSD2 showed structural variation. These results suggested that MJDs family members had considerable genetic stability as potential HCC patients’ diagnostic biomarkers.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g005.jpg)
FIGURE 5. The genetic alteration and mutation of MJDs family members in HCC (cBioPortal) (A–E) The proportions of alternations and mutations of MJDs family members were indicated in various subgroups of HCC (F–I). Schematic representation of gene mutation sites of MJDs on the coding strand (J)The genetic alternations and mutations of ATXN3, ATXN3L, JOSD1, and JOSD2 were showed in HCC patients (E).
Moreover, the frequency distribution of MJDs family members’ CNV patients in different stage and grade groups was shown. This suggested a high incidence and early event of the MJDs family members’ CNV alteration in HCC (Figures 6A–D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g006.jpg)
FIGURE 6. MJDs family members genomic alterations in HCC (cBioPortal) (A–D) CNV frequencies of MJDs family members were distributed in different of stages and grades of subgroups.
Survival analysis of MJDs family members in HCC
To evaluate the prognostic effects of MJDs and MJDs family members’ DNA methylation, we used Kaplan–Meier and MethSurv to explore the correlation between MJDs family members’ mRNA expressions and prognoses in HCC patients. The main parameters of survival analysis include DSS, OS, PFS, and RFS. The data from Kaplan–Meier Plotter was displayed in Figure 7, suggesting that highly expressed ATXN3 and ATXN3L levels were a remarkable correlation with better prognosis in HCC patients. Interestingly, we found that high-expression ATXN3 and ATXN3L had better DSS, OS, PFS, and RFS in HCC patients. On the contrary, high expression levels of JOSD1 and JOSD2 implicated a worse prognosis in HCC patients. Moreover, we found that highly expressed JOSD1 showed a significant association with worse DSS, OS, PFS, and RFS in overall HCC patients. Figure 7 and Table 1 revealed the more detailed prognostic information about MJDs members in HCC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g007.jpg)
FIGURE 7. (A‐G) The prognosis of MJDs expression level in HCC patients. The HCC specimens were divided into two groups according to the median expression levels of MJDs family members. The disease-specific survival (DSS, n = 362), overall survival (OS, n = 364), progression-free survival (PFS, n = 370), and relapse-free survival (RFS, n = 316) were compared between patients across high and low expression of MJDs family members. TCGA database screening results for 371 HCC patient samples, including risk score, survival time, and survival status. The middle shows the survival matching to the Risk score of various samples, while the top shows the Risk score scatter plot from low to high, with different colors denoting distinct risk categories. Time and survival status scatter plot distribution; the expression heat map of the signature’s genes is shown in the bottom figure (E) The KM survival curve distribution of the risk model in the data set, in which different groups are tested by log rank, HR (High risk) represents the risk coefficient between the high-risk group and the low-risk group samples, and HR > 1 indicates that the model is a risk model. 95% CI represents the HR confidence interval (F) The risk model’s ROC curve and AUC at various points; the greater the AUC value, the better the model’s capacity for forecasting (G).
We used dimensionality reduction based on Lasso regression to analyze the prognosis of MJDs members gene on patients with HCC and created a predictive model to assess OS. In Figure 7E, we found that when the three genes of ATXN3, JOSD1, and JOSD2 were co-expressed, the number of patients died and the risk score increased. The Kaplan–Meier Plotter showed that high risk score group of ATXN3, JOSD1, and JOSD2 overall high expression had worse OS, compared with low risk score group (Figure 7F). According to the ROC curve of the risk model established at different times, the AUC of 1-Years, 3-Years, 5-Years was 0.651, 0.571, 0.606, which indicated that the model has strong predictive ability (Figure 7G). This prognostic model indicated that worse survival was associated with overall high expression of MJDs members.
Subsequently, the MethSurv database was used to further reveal the relationship between methylation sites of MJDs and prognosis in HCC patients (Figure 8, Supplementary Figure S3). We found that the methylation level of cg26081025 in ATXN3 (HR = 3.161, p-value = 1.6e-06), cg07186939 in ATXN3L (HR = 1.621, p-value = 0.0062), cg25697769 in JOSD1(HR = 1.828, p-value = 0.004), and cg18708810 in JOSD2(HR = 1.735, p-value = 0.012) were a risk factor on survival. However, the methylation level of cg06711259 (HR = 0.455, p-value = 3.1e-05), cg19658332 (HR = 0.612, p-value = 0.014) cg27610821 (HR = 0.665, p-value = 0.023), and cg01138530 (HR = 0.67, p-value = 0.026) in JOSD1 were a protective factor on survival. Finally, the results suggested that eight CpG sites in ATXN3, the one CpG site in ATXN3L, 10 CpG sites in JOSD1 and the one CpG site in JOSD2 were related to the prognosis in HCC patients (Figure 8, Supplementary Figure S3, and Supplementary Table S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g008.jpg)
FIGURE 8. Kaplan-Meier survival curves comparing the high or low expression of MJDs methylation sites in HCC (MethSurv). Kaplan-Meier survival analysis of (A) ATXN3 and (B) ATXN3L methylation sites were shown.
MJDs family members may regulate the development, distribution, and maturation of immune cells, but the correlation between MJDs family members’ expression levels and immune-cell infiltration in HCC remains poorly unclear.
We used the TIMER score to assess this correlation. According to the expression levels of the ATXN3, JOSD1, and JOSD2 genes, the HCC patient samples in the TCGA database were split into 186 high-expression groups and 185 low-expression groups. In Figure 9A, the high-expression ATXN3 group had higher expression of CD4+T cells, CD8+T cells, macrophages, myeloid dendritic cells, B cells, and neutrophils, compared with the low-expression ATXN3 group. Then, JOSD1 expression was positively correlated with B cells, CD8+T cells, CD4+T cells, macrophages, neutrophils, and dendritic cells (Figure 9B). In addition, the change in JOSD2 expression was positively correlated with B cells, CD8+T cells, CD4+T cells, macrophages, neutrophils, and dendritic cells (Figure 9C). Therefore, these results implicated that MJDs family members expression may predict the infiltration of various kinds of immune cells in HCC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g009.jpg)
FIGURE 9. TIMER immune score of high expression group and low expression group of MJDs gene (A‐C) Different colors represent different immune cell types, the abscissa represents the group, and the ordinate represents the percentage of immune cells in a single group. This graph compares the percentage abundance of tumor-infiltrating immune cells in the high expression group (G1) and low expression group (G2) of MJDs genes. *p < 0.05, **p < 0.01, and ***p < 0.001.
In addition, correlations between the MJDs family members’ somatic copy number alterations (SCNA) and tumor immune infiltration levels in HCC patients were explored. The study showed that the SCNA of ATXN3 and JOSD1 had no significant correlations with the infiltration levels of CD4+ T cells, CD8+ T cells, B cells, neutrophils, macrophages, and DCs; the SCNA of ATXN3L had significant correlations with CD4+ T cells, CD8+ T cells, B cells, neutrophils, macrophages, and DCs; moreover, the SCNA of JOSD2 had significant correlations with CD4+ T cells (Supplementary Figure S4). These findings revealed that changes in the SCNA of the ATXN3L and JOSD2 may reflect tumor infiltration levels in HCC.
Relationship between MJDs members and immune molecules
The link between MJDs members and immune checkpoint molecules was then thoroughly investigated. The HCC patient samples in the TCGA database were divided into high and low expression groups of ATXN3, JOSD1, JOSD2, and the distribution of immune checkpoints in different gene expression samples was studied (Figure 10). The gene expression of ATXN3 was positively correlated with the expression of CD274 (p<0.001), PDCD1LG2 (p < 0.01), SIGLEC15 (p<0.01); the gene expression of JOSD1 had positive correlations with CD274 (p<0.001), CTLA4 (p < 0.001), HAVCR2(p < 0.001), LAG3 (p < 0.05), PDCD1 (p < 0.001), PDCD1LG2 (p < 0.01), TIGIT (p < 0.01); and the gene expression of JOSD2 had positive correlations with CTLA4 (p < 0.01), HAVCR2(p < 0.01), PDCD1 (p < 0.001),TIGIT (p < 0.05). These findings showed a positive correlation between the expression of immunological checkpoints and the elevation of MJDs member genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g010.jpg)
FIGURE 10. The relationship between MJDs gene expression and immune checkpoints (A‐C) The HCC patients in the TCGA database were classified into high expression group (G1) and low expression group (G2) based on the various levels of ATXN3, JOSD1 and JOSD2 gene expression, and the association between the two groups and the expression of immune checkpoints was presented, respectively. *p < 0.05, **p < 0.01, and ***p < 0.001.
We also explored relationships between MJDs expression and various immune signatures with the TISIDB database. Associations between MJDs expression and different immune markers were obtained from the TISIDB database. Figure 11 and Supplementary Figure S5 showed the correlations between MJDs and TILs which included MHC-TAP1, Act_CD4_abundance, Act_DC_abundance, and CD56bright. Immunomodulators can be divided into three groups: immunoinhibitors, immunostimulators, and MHC molecules. In Figure 11A and Table 2, HAVCR2, LGALS9, MICB, TNFRSF4, TNFRSF9, TNFSF9, MHC_TAP1, Act_CD4, and Act_DC were negatively associated with the ATXN3 expression, while KDR, IL6R, and ICOSLG positively correlated with the ATXN3 expression. In Figure 11B and Supplementary Table S3, JOSD1 negatively correlated with CD160, CD56bright, and CCL14, while CCL28 was positively related to this gene. In the Supplementary Figure S5 and Supplementary Table S4, JOSD2 significantly positively correlated with LGALS9, PVRL2, TNFRSF4, TNFRSF14, TNFRSF18, TNFRSF25, MHC_HLA-A, and MHC_HLA-DMA, while CD274, CD28, CCR4, IL6R, KDR, and TGFBR1 negatively correlated with JOSD2 expression. Thus, these results demonstrated that ATXN3, JOSD1, and JOSD2 were related to various immune molecules in HCC, which implicated a significant role in immune escape.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g011.jpg)
FIGURE 11. Correlation between MJDs level and lymphocytes, immunomodulators, and chemokines expression in HCC from TISIDB database. The correlation between alterations in expression levels of (A) ATXN3, and (B) JOSD1 with the expression of immune cells marker in HCC were performed.
Functions and pathways for MJDs family members in HCC
We used the LinkedOmics database to analyze the differentially expressed genes associated with ATXN3, JOSD1, and JOSD2 in HCC (Figure 12A, Supplementary Figures S6A,7A). LinkedOmics and the cBioPortal database were screened for the top 200 co-expressed genes to obtain crossover genes, respectively (Figure 12B, Supplementary Figure S6B,7B). The above intersecting genes were obtained and analyzed by GO and KEGG to predict functions and pathways of MJDs family members. The KEGG pathways showed that for ATXN3, ubiquitin mediated proteolysis, fanconi anemia pathway, aminoacyl-tRNA biosynthesis and RNA transport were significantly up-regulated, whereas valine, leucine and isoleucine degradation, staphylococcus aureus infection, glycine serine and threonine metabolism were significantly down-regulated (Figure 12C). For JOSD1, allograft rejection, Th17 cell differentiation, and tryptophan metabolism were significantly up-regulated, whereas DNA replication and fanconi anemia pathway were significantly down-regulated (Supplementary Figure S6C). For JOSD2, allograft rejection, Th17 cell differentiation, and ascorbate and aldarate metabolism were significantly up-regulated, whereas ribosome, oxidative phosphorylation, and protein export were significantly down-regulated (Supplementary Figure S7C). Interestingly, GO analysis showed that the co-expressed genes of MJDs family members were the same in Biological process (BP), Cellular component (CC) and Molecular Function (MF). Specifically, BP terms indicated regulation, metabolic process, and response to a stimulus. CC terms implicated the membrane, nucleus, and membrane-enclosed lumen. MF terms suggested protein binding, ion binding, and nucleic acid binding (Figure 12D, Supplementary Figure S6D,7D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g012.jpg)
FIGURE 12. Different genes and pathways connected to MJDs family members expression in HCC (A) Volcano plot suggested that the differential expression of genes correlated with ATXN3 in HCC (B) The top 200 co-expressed genes were extracted from the LinkedOmics and cBioPortal, and then the above genes were intersected (C,D) The prediction of GO analysis and KEGG pathway of MJDs family members were showed.
Furthermore, we analyzed the co-expression of ATXN3, JOSD1 and JOSD2 in hepatocellular carcinoma using the GSCALite database, and analyzed the signaling pathways involved by MJDs family member respectively (Supplementary Figure S8). We found that ATXN3 inhibited apoptosis, while JOSD1 and JOSD2 promoted apoptosis, JOSD2 inhibited DNA damage response, JOSD1 and JOSD2 promoted EMT, and JOSD2 suppressed RTK and androgen pathways.
Moreover, we used the cBioPortal database to explore the co-expression (Table 3) and mutual exclusion (Table 4) between individual members of the MJDs family. Our predicted results showed that there was no co-expression and no mutual exclusion among MJDs family members, but further experiments such as immunoprecipitation or immunofluorescence were needed to verify.
Validation of the expression of MJDs family members in HCC tissue and HCC cell lines
Next, we explored the MJDs family members’ expression levels between HCC samples and normal liver samples from the HPA database. The IHC staining of MJDs family members showed that the mRNA expression levels of ATXN3, ATXN3L, and JOSD1 were enhanced compared with the liver tissues (Figures 13A–C), which consistent with above results. The IHC staining of JOSD2 was pending HCC tumor tissue analysis in the HPA.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/965805/fgene-13-965805-HTML/image_m/fgene-13-965805-g013.jpg)
FIGURE 13. The expression levels of MJDs family members in HCC tissue and cell lines. The protein expression of (A) ATXN3, (B) ATXN3L, and (C) JOSD1 in HCC and normal liver tissues were obtained by immunohistochemical (IHC) data in the HPA database. The IHC image of JOSD2 was absent. The expression levels of (D) ATXN3, (E) JOSD1, and (F) JOSD2 in human immortalized HL-7702 and different HCC cell lines were conformed using RTqPCR.Ct values of ATXN3L were absent in HL-7702, MHCC-97 and HepG2 cell lines. *p < 0.05, **p < 0.01, and ***p < 0.001.
To verify above results with bioinformatics analysis, we used qRT-PCR to confirm the expression levels of MJDs in five HCC cell lines which included MHCC-97H, HepG2, PLC/PRF/5, Huh7, and HCCLM-3, respectively, compared with the HL-7702. The results suggested that, except ATXN3L, the expression levels of MJDs were elevated in multiple HCC cell lines (Figure 13 D-F), in consistent with the above results, and further experiments were needed to support it.
Discussion
Accumulating evidence indicated that the abnormal expression or mutation of deubiquitinases (DUBs), with the function of promoting the occurrence and development of malignant tumors, was emerging as therapeutic targets of cancer (Lv et al., 2020; Park and Baek, 2022). MJDs family, consisting of ATXN3, ATXN3L, JOSD1, and JOSD2, was the smallest family of the DUBs (Harrigan et al., 2018). MJDs family had a common Josephin domain, containing around 180 amino acids. According to recent studies, MJDs family members were aberrant expressions in some cancer and were related to biological functions and prognoses in cancers (Zeng et al., 2020). However, the aberrant mRNA expressions of MJDs family members and their biological functions and mechanisms in HCC were poorly understood. Thus, we further explored the role of MJDs family members in HCC, particularly, including immune microenvironment and promoter methylation levels.
Recent studies reported that oral squamous cell carcinoma (OSCC) (Song et al., 2021), testicular cancer (TC) (Shi et al., 2018), colon cancer (Li et al., 2019), breast cancer (BC) (Zou et al., 2019), anaplastic thyroid carcinoma (ATC) (Zhuang et al., 2021), neuroblastoma (NB) (Gong et al., 2021), non-small cell lung cancer (NSCLC) (Sacco et al., 2014), and breast cancer stem cells (Zhu et al., 2019)were closely related to ATXN3. For example, for oral squamous cell carcinoma, it had been studied that down-regulated miR-619-5p promoted the proliferation, migration, and invasion abilities by enhancing ATXN3 expression in OSCC cisplatin-resistant cells. The miR-619-5p/ATXN3 axis provided the potential molecular mechanism of OSCC, especially for promoting drug resistance of cisplatin of OSCC. And ATXN3 may be a putative therapeutic target for OSCC patients (Song et al., 2021). In addition, ATXN3 knockdown can inhibit cell growth. Meanwhile, down-regulated ATXN3 can promote PTEN expression and inactivate the AKT/mTOR pathway (Shi et al., 2018). Previous research reported that ATXN3L was expressed at higher levels in breast cancer. ATXN3L directly interacted with KLF5, thus, reducing KLF5’s ubiquitination and degradation (Ge et al., 2015). JOSD1 was found to be up-regulated in the head and neck squamous cell carcinoma (HNSCC) (Jing et al., 2021), gynaecological cancer (GC) cells (Wu et al., 2020), and acute myeloid leukemia (AML) cells (Yang et al., 2022). Over-expressed JOSD2 was positively correlated with the worse prognosis in hepatocellular carcinoma (Huang et al., 2022), cholangiocarcinoma (CCA) (Qian et al., 2021), and NSCLC (Krassikova et al., 2021).
Our results showed that ATXN3, JOSD1, and JOSD2 were markedly up-regulated in HCC samples, while ATXN3L was expressed only in testicular and testicular tumors. We found that the higher mRNA expression levels of JOSD1 and JOSD2 in HCC patients had a significantly worse prognosis. Interestingly, we analyzed that the prognostic value of MJDs family members’ expression levels in HCC patients revealed that higher mRNA expression of AXTN3 and ATXN3L have better OS, DSS, PFS, and RFS, which was not consistent with previous studies. The prognostic model indicated that worse survival was associated with overall high expression of MJDs members. Furthermore, expression levels of ATXN3, JOSD1, and JOSD2 were significantly related to different clinicopathological parameters in HCC. Thus, the above results strongly revealed that ATXN3, JOSD1, and JOSD2 deserved further investigation of potential targets in HCC.
Many studies indicated that methylation levels were closely correlated with tumor development, cancer progression, and chemotherapy and immunotherapy sensitivity of HCC (Al-Abdulla et al., 2019; Dong et al., 2021; Zhong et al., 2021; Liu et al., 2022) Thus, to explore the mechanism for the aberrant expression of MJDs family members in HCC, we further analyzed the correlations between these genes expression levels and promoter methylation levels with the UALCAN database. In our study, the promoter methylation levels of ATXN3 and JOSD2 were lower expressed in these genes’ high mRNA expression group in TCGA HCC samples, while JOSD1 was highly expressed in the JOSD1 high expression group. Moreover, we used the cBioPortal database to analyze the correlations between promoter methylation degrees of ATXN3, JOSD1, and JOSD2 and the transcriptional levels. We found that the Spearman coefficients of ATXN3, JOSD1, and JOSD2 were -0.31,-0.35, and -0.29, respectively, suggesting significantly negative correlations, which was consistent with the above results. In addition, it was worth noting that the promoter methylation levels of ATXN3L were significantly lower in HCC tissues than in normal tissues, and the methylation β value of ATXN3L was greater than 0.6, presenting a full methylation status. Furthermore, the previous research showed that the methylation mediated inactivation of genes was common in HCC compared with normal samples and connected to clinical features and prognose (He et al., 2020; Xu et al., 2021). Consistently, our results demonstrated that promoter methylation levels of MJDs family members were up-regulated in clinical features of HCC. Meanwhile, we further analyzed relationships of 39 loci within MJDs family members’ prognosis and confirmed these genes’ DNA methylation of twenty loci located on promoter CpG islands, which connected to worst survival. The above results indicated that abnormal methylation might be one of the crucial reasons for the aberrant expression of MJDs gene family. However, other epigenetic modifications, involving gene mutations, CNV, and SNV might also play important roles in the abnormal expression of these genes.
To further elucidate the potential mechanism of the MJDs family, co-expressed genes with ATXN3, JOSD1, and JOSD2 in HCC samples from the TCGA database were extracted for GO and KEGG enrichment analysis to study the functions of the MJDs family in HCC. As expected, the biological functions of MJDs family members primarily correlated with ubiquitin mediated proteolysis, Th17 cell differentiation, valine, leucine and isoleucine degradation, glycine serine and threonine metabolism, and oxidative phosphorylation. It was reported that MJDs family members involve multiple signaling pathways, including ubiquitin pathway (Zhuang et al., 2021), apoptosis (Li et al., 2019), Hedgehog signaling pathway (Xu et al., 2021), JAK2 (Yang et al., 2022), and Akt/PI3K signaling (Sacco et al., 2014). In addition, a previous study reported that the serine-biosynthesis pathway and tyrosine-relative kinase inhibitors played important regulatory roles in tumor progression and immunotherapy in liver cancer patients (Sangro et al., 2021). These findings suggested that ATXN3, JOSD1, and JOSD2 could participate in the progression of HCC through primary regulating genes post-translational modification pathways, and the immune environment.
It was critical to investigate the relationships between MJDs family members’ expression levels and the multiple immune-cell infiltration in HCC. Based on the TIMER database, obviously, correlations were shown between MJDs family members’ expression and the infiltrating levels of B cells, CD4+ T cells, CD8+ T cells, neutrophils, dendritic cells (DCs), and macrophages. Relevant evidence had shown that the presence of tumor-infiltrating lymphocytes (TILs) might relate to the prognosis of HCC patients (Li et al., 2020; Nie et al., 2021). Immune cell based immunotherapy (Liu et al., 2021), including T Cells (Woller et al., 2021), neutrophils (Geh et al., 2022), and dendritic cells (Matsui et al., 2021), played crucial roles in the immunological therapy of HCC. In general, previous studies and our research revealed that MJDs family members might be involved in the regulation of the infiltration cell recruitment of immune cells in the immune microenvironment, which could be targeted for immunotherapy in HCC patients.
The combination of immune therapy and the classic target therapies may be the emerging trend of clinical comprehensive treatment. For our study, we revealed that the transcriptional expression of MJDs family members was correlated with lymphocytes, immunomodulators (immunomoinhibitor and immunomostimulator), and chemokines in HCC from the TISIDB database. In this research, HAVCR2, LGALS9, MICB, Act_CD4, and Act_DC were negatively associated with the ATXN3 expression, while KDR, IL6R, and ICOSLG positively correlated with the ATXN3 expression. A recent study identified the binding of HAVCR2 to LGALS9 contributes to Th1 cell death via apoptosis (Zhu et al., 2005). Meanwhile, HAVCR2 was an inducible NK cell receptor that promotes IFN-γ production in response to LGALS9 (Gleason et al., 2012). Tocilizumab, an IL6R antagonist, was approved for the management of CAR T-cell-related Cytokine Release Syndrome (Patel et al., 2021). Resistance to checkpoint blockade is mainly mediated via down-regulated expression of MHC-I in tumor cells. The MHC -I chain-related polypeptide A/B (MICA/B), expressed in many human cancers, serves as ligands to activate the NKG2D receptor on the NK cells and T cells (Raulet et al., 2013). Inhibition of proteolytic shedding of MICA/B leads NK cells mediated immunity against tumors resistant to T cells (Ferrari de Andrade et al., 2020). Similarly, JOSD1 negatively correlated with CD160 and CCL14, while CCL28 was positively related to this gene. In addition, JOSD2 positively correlated with TNFRSF18, HLA-A, and HLA-DMA, while CD274, TGFBR1, CCR4, and CD28 negatively correlated with JOSD2 expression. Emerging studies have found that chemokine ligands (CCL) and chemokine receptors (CCR) are closely connected to the immune therapy of HCC. Notably, CCL14 expression in HCC negatively correlated with PD-1, HAVCR2, and CTLA-4, suggesting its role in regulating tumor immunity (Gu et al., 2020). CCR4 could enhance anti-tumor immunity, mainly by targeting and blocking the infiltration of regulatory T cells (Tregs)into the tumor microenvironment and inhibiting the stability of the TIL-Treg pool (Gao et al., 2022). Lenvatinib showed lower CD274(PD-L1) expression and Tregs infiltration in recurrent HCC compared with primary HCC(Yi et al., 2021). It has been excitingly reported that methylation at certain CpG sites was an indicator of immune infiltration of cancer and might predict patient response to checkpoint inhibitors (Bacolod et al., 2019). Combining previous studies with our results suggested that ATXN3, JOSD1, and JOSD2 might be a scheme to improve the outcome of immunotherapy. Therefore, it was necessary to further explore more detailed mechanism of transcriptional levels/DNA methylation of MJDs family members in regulating the tumor micro-environment.
In summary, our results suggested that the MJDs family members’ aberrant expression in HCC played an important role in clinical features, prognosis, tumor micro-environment, immune-related molecules, mechanisms, and functions. Therefore, our study will give support to elucidate the clinical significance of MJDs family members’ expression levels and DNA methylation in HCC, and further explore the mechanisms of the MJDs family member’s roles in HCC, and especially provide insight into potentially therapeutic targets to improve HCC patient prognosis. In addition, we found that the online tools were based on different algorithms, leading to inconsistent analysis results. Thus, our results need to be further verified in experiments.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
LZ conceived this study and drafted the manuscript. GC, TL, XL, and CY analyzed the data, collected the data, reviewed and revised the manuscript. JJ reviewed the manuscript. All authors reviewed and approved the submitted manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (No. 81871965).
Conflict of interest
Author GC was employed by the company Hunan YoBon Biotechnology Limited.
The remaining 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/fgene.2022.965805/full#supplementary-material
References
Al-Abdulla, R., Lozano, E., Macias, R. I. R., Monte, M. J., Briz, O., O'Rourke, C. J., et al. (2019). Epigenetic events involved in organic cation transporter 1-dependent impaired response of hepatocellular carcinoma to sorafenib. Br. J. Pharmacol. 176 (6), 787–800. doi:10.1111/bph.14563
Bacolod, M. D., Barany, F., and Fisher, P. B. (2019). Can CpG methylation serve as surrogate markers for immune infiltration in cancer? Adv. Cancer Res. 143, 351–384. doi:10.1016/bs.acr.2019.03.007
Buus, R., Faronato, M., Hammond, D. E., Urbé, S., and Clague, M. J. (2009). Deubiquitinase activities required for hepatocyte growth factor-induced scattering of epithelial cells. Curr. Biol. 19 (17), 1463–1466. doi:10.1016/j.cub.2009.07.040
Dong, Z.-R., Ke, A.-W., Li, T., Cai, J.-B., Yang, Y.-F., Zhou, W., et al. (2021). CircMEMO1 modulates the promoter methylation and expression of TCF21 to regulate hepatocellular carcinoma progression and sorafenib treatment sensitivity. Mol. Cancer 20 (1), 75. doi:10.1186/s12943-021-01361-3
Feng, D., Wang, M., Hu, J., Li, S., Zhao, S., Li, H., et al. (2020). Prognostic value of the albumin-bilirubin grade in patients with hepatocellular carcinoma and other liver diseases. Ann. Transl. Med. 8 (8), 553. doi:10.21037/atm.2020.02.116
Ferrari de Andrade, L., Kumar, S., Luoma, A. M., Ito, Y., Alves da Silva, P. H., Pan, D., et al. (2020). Inhibition of MICA and MICB shedding elicits NK-Cell-Mediated immunity against tumors resistant to cytotoxic T cells. Cancer Immunol. Res. 8 (6), 769–780. doi:10.1158/2326-6066.CIR-19-0483
Gao, Y., You, M., Fu, J., Tian, M., Zhong, X., Du, C., et al. (2022). Intratumoral stem-like CCR4+ regulatory T cells orchestrate the immunosuppressive microenvironment in HCC associated with hepatitis B. J. Hepatol. 76 (1), 148–159. doi:10.1016/j.jhep.2021.08.029
Ge, F., Chen, W., Qin, J., Zhou, Z., Liu, R., Liu, L., et al. (2015). Ataxin-3 like (ATXN3L), a member of the Josephin family of deubiquitinating enzymes, promotes breast cancer proliferation by deubiquitinating Krüppel-like factor 5 (KLF5). Oncotarget 6 (25), 21369–21378. doi:10.18632/oncotarget.4128
Geh, D., Leslie, J., Rumney, R., Reeves, H. L., Bird, T. G., and Mann, D. A. (2022). Neutrophils as potential therapeutic targets in hepatocellular carcinoma. Nat. Rev. Gastroenterol. Hepatol. 19 (4), 257–273. doi:10.1038/s41575-021-00568-5
Gleason, M. K., Lenvik, T. R., McCullar, V., Felices, M., O'Brien, M. S., Cooley, S. A., et al. (2012). Tim-3 is an inducible human natural killer cell receptor that enhances interferon gamma production in response to galectin-9. Blood 119 (13), 3064–3072. doi:10.1182/blood-2011-06-360321
Gong, B., Zhang, J., Hua, Z., Liu, Z., Thiele, C. J., and Li, Z. (2021). Downregulation of ATXN3 enhances the sensitivity to AKT inhibitors (perifosine or MK-2206), but decreases the sensitivity to chemotherapeutic drugs (etoposide or cisplatin) in neuroblastoma cells. Front. Oncol. 11, 686898. doi:10.3389/fonc.2021.686898
Gu, Y., Li, X., Bi, Y., Zheng, Y., Wang, J., Li, X., et al. (2020). CCL14 is a prognostic biomarker and correlates with immune infiltrates in hepatocellular carcinoma. Aging 12 (1), 784–807. doi:10.18632/aging.102656
Harrigan, J. A., Jacq, X., Martin, N. M., and Jackson, S. P. (2018). Deubiquitylating enzymes and drug discovery: Emerging opportunities. Nat. Rev. Drug Discov. 17 (1), 57–78. doi:10.1038/nrd.2017.152
He, H., Chen, D., Cui, S., Wu, G., Piao, H., Wang, X., et al. (2020). HDNA methylation data-based molecular subtype classification related to the prognosis of patients with hepatocellular carcinoma. BMC Med. Genomics 13 (1), 118. doi:10.1186/s12920-020-00770-5
Huang, Y., Zeng, J., Liu, T., Xu, Q., Song, X., and Zeng, J. (2022). Deubiquitinating enzyme JOSD2 promotes hepatocellular carcinoma progression through interacting with and inhibiting CTNNB1 degradation. Cell. Biol. Int. 46, 1089–1097. doi:10.1002/cbin.11812
Jing, C., Liu, D., Lai, Q., Li, L., Zhou, M., Ye, B., et al. (2021). JOSD1 promotes proliferation and chemoresistance of head and neck squamous cell carcinoma under the epigenetic regulation of BRD4. Cancer Cell. Int. 21 (1), 375. doi:10.1186/s12935-021-02060-1
Krassikova, L., Zhang, B., Nagarajan, D., Queiroz, A. L., Kacal, M., Samakidis, E., et al. (2021). The deubiquitinase JOSD2 is a positive regulator of glucose metabolism. Cell. Death Differ. 28 (3), 1091–1109. doi:10.1038/s41418-020-00639-1
Li, D., Zhang, T., Lai, J., Zhang, J., Wang, T., Ling, Y., et al. (2019). MicroRNA-25/ATXN3 interaction regulates human colon cancer cell growth and migration. Mol. Med. Rep. 19 (5), 4213–4221. doi:10.3892/mmr.2019.10090
Li, L., Shen, L., Ma, J., Zhou, Q., Li, M., Wu, H., et al. (2020). Evaluating distribution and prognostic value of new tumor-infiltrating lymphocytes in HCC based on a scRNA-seq study with CIBERSORTx. Front. Med. 7, 451. doi:10.3389/fmed.2020.00451
Liu, C., Shao, J., Dong, Y., Xu, Q., Zou, Z., Chen, F., et al. (2021). Advanced HCC patient benefit from neoantigen reactive T cells based immunotherapy: A case report. Front. Immunol. 12, 685126. doi:10.3389/fimmu.2021.685126
Liu, Y., Cheng, H., Cheng, C., Zheng, F., Zhao, Z., Chen, Q., et al. (2022). ZNF191 alters DNA methylation and activates the PI3K-AKT pathway in hepatoma cells via transcriptional regulation of DNMT1. Cancer Med. 11 (5), 1269–1280. doi:10.1002/cam4.4535
Lv, X.-Y., Duan, T., and Li, J. (2020). The multiple roles of deubiquitinases in liver cancer. Am. J. Cancer Res. 10 (6), 1647–1657.
Matsui, H. M., Hazama, S., Nakajima, M., Xu, M., Matsukuma, S., Tokumitsu, Y., et al. (2021). Novel adjuvant dendritic cell therapy with transfection of heat-shock protein 70 messenger RNA for patients with hepatocellular carcinoma: A phase I/II prospective randomized controlled clinical trial. Cancer Immunol. Immunother. 70 (4), 945–957. doi:10.1007/s00262-020-02737-y
Nie, H., He, T., Wang, L., and Zhang, L. (2021). Expression and prognostic value of tumor-infiltrating lymphocytes and PD-L1 in hepatocellular carcinoma. Onco. Targets. Ther. 14, 1377–1385. doi:10.2147/OTT.S289720
Park, H.-B., and Baek, K.-H. (2022). E3 ligases and deubiquitinating enzymes regulating the MAPK signaling pathway in cancers. Biochim. Biophys. Acta. Rev. Cancer 1877 (3), 188736. doi:10.1016/j.bbcan.2022.188736
Patel, K., Gooley, T. A., Bailey, N., Bailey, M., Hegerova, L., Batchelder, A., et al. (2021). Use of the IL-6R antagonist tocilizumab in hospitalized COVID-19 patients. J. Intern. Med. 289 (3), 430–433. doi:10.1111/joim.13163
Qian, M., Yan, F., Wang, W., Du, J., Yuan, T., Wu, R., et al. (2021). Deubiquitinase JOSD2 stabilizes YAP/TAZ to promote cholangiocarcinoma progression. Acta Pharm. Sin. B 11 (12), 4008–4019. doi:10.1016/j.apsb.2021.04.003
Raulet, D. H., Gasser, S., Gowen, B. G., Deng, W., and Jung, H. (2013). Regulation of ligands for the NKG2D activating receptor. Annu. Rev. Immunol. 31, 413–441. doi:10.1146/annurev-immunol-032712-095951
Sacco, J. J., Yau, T. Y., Darling, S., Patel, V., Liu, H., Urbé, S., et al. (2014). The deubiquitylase Ataxin-3 restricts PTEN transcription in lung cancer cells. Oncogene 33 (33), 4265–4272. doi:10.1038/onc.2013.512
Sangro, B., Sarobe, P., Hervás-Stubbs, S., and Melero, I. (2021). Advances in immunotherapy for hepatocellular carcinoma. Nat. Rev. Gastroenterol. Hepatol. 18 (8), 525–543. doi:10.1038/s41575-021-00438-0
Shi, Z., Chen, J., Zhang, X., Chu, J., Han, Z., Xu, D., et al. (2018). Ataxin-3 promotes testicular cancer cell proliferation by inhibiting anti-oncogene PTEN. Biochem. Biophys. Res. Commun. 503 (1), 391–396. doi:10.1016/j.bbrc.2018.06.047
Song, A., Wu, Y., Chu, W., Yang, X., Zhu, Z., Yan, E., et al. (2021). Involvement of miR-619-5p in resistance to cisplatin by regulating ATXN3 in oral squamous cell carcinoma. Int. J. Biol. Sci. 17 (2), 430–447. doi:10.7150/ijbs.54014
Woller, N., Engelskircher, S. A., Wirth, T., and Wedemeyer, H. (2021). Prospects and challenges for T cell-based therapies of HCC. Cells 10 (7), 1651. doi:10.3390/cells10071651
Wu, X., Luo, Q., Zhao, P., Chang, W., Wang, Y., Shu, T., et al. (2020). JOSD1 inhibits mitochondrial apoptotic signalling to drive acquired chemoresistance in gynaecological cancer by stabilizing MCL1. Cell. Death Differ. 27 (1), 55–70. doi:10.1038/s41418-019-0339-0
Xu, Q., Hu, Y., Chen, S., Zhu, Y., Li, S., Shen, F., et al. (2021). Immunological significance of prognostic DNA methylation sites in hepatocellular carcinoma. Front. Mol. Biosci. 8, 683240. doi:10.3389/fmolb.2021.683240
Yang, J., Weisberg, E. L., Liu, X., Magin, R. S., Chan, W. C., Hu, B., et al. (2022). Small molecule inhibition of deubiquitinating enzyme JOSD1 as a novel targeted therapy for leukemias with mutant JAK2. Leukemia 36 (1), 210–220. doi:10.1038/s41375-021-01336-9
Yi, C., Chen, L., Lin, Z., Liu, L., Shao, W., Zhang, R., et al. (2021). Lenvatinib targets FGF receptor 4 to enhance antitumor immune response of anti-programmed cell death-1 in HCC. Hepatology 74 (5), 2544–2560. doi:10.1002/hep.31921
Zeng, C., Zhao, C., Ge, F., Li, Y., Cao, J., Ying, M., et al. (2020). Machado-joseph deubiquitinases: From cellular functions to potential therapy targets. Front. Pharmacol. 26 (11), 1311. doi:10.3389/fphar.2020.01311
Zheng, R., Zhang, S., Zeng, H., Wang, S., Sun, K., Chen, R., et al. (2022). Reconceptualizing solidarity as power from below. Philos. Stud. 2 (1), 1–25. doi:10.1007/s11098-022-01845-y
Zhong, G.-C., Zhao, Z.-B., Cheng, Y., Wang, Y.-B., Qiu, C., Mao, L.-H., et al. (2021). Epigenetic silencing of GCH1promotes hepatocellular carcinoma growth by activating superoxide anion-mediated ASK1/p38 signaling via inhibiting tetrahydrobiopterin de novo biosynthesis. Free Radic. Biol. Med. 168, 81–94. doi:10.1016/j.freeradbiomed.2021.03.025
Zhu, C., Anderson, A. C., Schubart, A., Xiong, H., Imitola, J., Khoury, S. J., et al. (2005). The Tim-3 ligand galectin-9 negatively regulates T helper type 1 immunity. Nat. Immunol. 6 (12), 1245–1252. doi:10.1038/ni1271
Zhu, R., Gires, O., Zhu, L., Liu, J., Li, J., Yang, H., et al. (2019). TSPAN8 promotes cancer cell stemness via activation of sonic Hedgehog signaling. Nat. Commun. 10 (1), 2863. doi:10.1038/s41467-019-10739-3
Zhuang, S., Xie, J., Zhen, J., Guo, L., Hong, Z., Li, F., et al. (2021). The deubiquitinating enzyme ATXN3 promotes the progression of anaplastic thyroid carcinoma by stabilizing EIF5A2. Mol. Cell. Endocrinol. 537, 111440. doi:10.1016/j.mce.2021.111440
Keywords: MJDs family members, hepatocellular carcinoma, methylation, immune cell infiltration, prognostic biomarker
Citation: Zhou L, Chen G, Liu T, Liu X, Yang C and Jiang J (2022) MJDs family members: Potential prognostic targets and immune-associated biomarkers in hepatocellular carcinoma. Front. Genet. 13:965805. doi: 10.3389/fgene.2022.965805
Received: 10 June 2022; Accepted: 10 August 2022;
Published: 09 September 2022.
Edited by:
Wei Xu, Shanghai Changzheng Hospital, ChinaReviewed by:
Guoqiang Sun, Nanjing Medical University, ChinaHai Li Bao, Shanghai Jiao Tong University, China
Copyright © 2022 Zhou, Chen, Liu, Liu, Yang and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianxin Jiang, rm002979@whu.edu.cn
†These authors have contributed equally to this work