Fukai Wen1,2†Fanshuai Meng1,2†Xuewen Li3†Qingyu Li1,2†Jiaming Liu1†
Rui Zhang1,2Yunzheng Zhao1,2Yu Zhang1,2Xin Wang1,2Shuai Ju4Yifeng Cui1,2*
Zhaoyang Lu1,2*- 1Department of Hepatic Surgery, The First Affiliated Hospital of Harbin Medical University, Harbin, China
- 2Key Laboratory of Hepatosplenic Surgery, Ministry of Education, The First Affiliated Hospital of Harbin Medical University, Harbin, China
- 3The Department of Inpatient Central Operating Room, The First Affiliated Hospital of Harbin Medical University, Harbin, China
- 4Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Harbin Medical University, Harbin, China
Background: The protein-coding gene RAB22A, a member of the RAS oncogene family, is amplified or overexpressed in certain cancers. However, its action mechanism in hepatocellular carcinoma (HCC) remains unclear. Here, we aimed to examine the connection between RAB22A and survival prognosis in HCC and explore the biological significance of RAB22A.
Methods: A database-based pan-cancer expression analysis of RAB22A was performed. Kaplan–Meier analysis and Cox regression were performed to evaluate the association between RAB22A expression and survival prognosis in HCC. Using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA), various potential biological functions and regulatory pathways of RAB22A in HCC were discovered. Tumor immune infiltration was studied using the single sample gene set enrichment analysis (ssGSEA) method. N6-methyladenosine modifications and the regulatory network of competitive endogenous RNA (ceRNA) were verified in the TCGA cohort.
Results: RAB22A was upregulated in HCC samples and cell lines. A high RAB22A expression in HCC was strongly correlated with sex, race, age, weight, TNM stage, pathological stage, tumor status, histologic grade, TP53 mutation status, and alpha fetal protein (AFP) levels. Overexpression of RAB22A indicated a poor prognosis was related to overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI). GO and KEGG analyses revealed that the differentially expressed genes related to RAB22A might be involved in the proteasomal protein catabolic process, ncRNA processing, ribosome ribosomal subunit, protein serine/threonine kinase activity, protein serine kinase activity, Endocytosis, and non-alcoholic fatty liver disease. GSEA analyses revealed that the differentially expressed genes related to RAB22A might be involved in the T cell receptor, a co-translational protein, that binds to the membrane, axon guidance, ribosome, phagocytosis, and Eukaryotic translation initiation. RAB22A was correlated with N6-methyladenosine expression in HCC and established RAB22A-related ceRNA regulatory networks. Finally,RAB22A expression was positively connected the levels of infiltrating with T helper cells, Tcm cells, and Th2 cells,In contrast, we observed negatively correlations with cytotoxic cells, DCs, and pDCs cells.Moreover,RAB22A expression showed a strong correlation with various immunomarkergroups in HCC.
Conclusions: RAB22A is a potential therapeutic target for improving HCC prognosis and is closely related to immune cell infiltration.
1 Introduction
Hepatocellular carcinoma (HCC) is the sixth most diagnosed cancer and the fourth leading cause of cancer death worldwide, with approximately 841,000 new cases and 782,000 deaths annually (1). Many key factors, including infection with hepatitis B or C and contact with foods contaminated with aflatoxin, contribute to HCC development (2). Surgery is the typical treatment for HCC; however, the disease is prone to relapse and metastasis, making it difficult to cure (3). Therefore, there is an urgent need to identify new relevant biomarkers to improve the early diagnosis, prognostic assessment, and treatment of HCC.
RAB22A is a small GTPase that belongs to the RAB protein family, specifically, the RAB5 subfamily (4). This protein is mainly located in early endosomes, Golgi bodies, and late endosomes. RAB proteins are involved in the regulation of vesicular traffic and exosome formation (5). Studies have found that the RAB5 subfamily (including RAB5, RAB21, RAB22A, and RAB22B) is primarily involved in the endocytosis, transport, and metabolism of growth factor receptors and may thus be associated with cancer progression (6–8). RAB22A expression is elevated in several malignancies, including breast, colorectal, and osteosarcoma cancer (9–11). It accelerates the progression of malignant tumors via various mechanisms, for instance, miRNA downregulation (11), recycling of extracellular matrix metalloproteinase inducer (EMMPRIN) (12), and hypoxia-inducible factor (13). Nevertheless, the function of RAB22A in HCC remains unclear.
Furthermore, RAB22A has multiple immune functions and is a novel immunomodulatory factor. Accurate intracellular transport of MHC-I molecules in dendritic cells (DCs) and T lymphocytes depends on RAB22A function (14). RAB22A is also part of the accommodative immune response and is absorbed by a process that separates it from the envelope proteins and spreads it throughout the body (15). Previous research has identified RAB22A as the main endosomal target in pathogen infection and a critical regulator of microbial infection and intracellular transport (16). In summary, RAB22A may have a significant prognostic and immunological significance in HCC.
In the current study, we analyzed the expression of RAB22A in HCC and paracancerous tissues using multiple datasets and in vitro experiments. Additionally, we examined the connection between RAB22A and survival prognosis in HCC and explored the biological significance of RAB22A by performing enrichment and protein-protein interaction (PPI) network analyses and determining the correlation with immune cell infiltration. Furthermore, we constructed ceRNA regulatory networks involving RAB22A in HCC. Our study proposes a possible connection between RAB22A expression and the presence of immune infiltrates in HCC.
2 Materials and methods
2.1 Database source and processing
Gene expression and clinical data were extracted from multiple databases (Supplementary Table 1) and RAB22A expression levels from RNA-seq data (TPM) of patients with HCC were analyzed. The Supplementary Materials and Methods (17–19) presents detailed information on the included data.
2.2 Patients and clinical samples
The First Affiliated Hospital of Harbin Medical University provided 30 matched sets of HCC and nearby non-tumor liver samples from patients undergoing hepatectomy between February 2020 and June 2022. This project was approved by the First Affiliated Hospital of Harbin Medical University’s Ethics Committee.
2.3 Western blotting and quantitative real-time PCR
Total proteins and total RNA were extracted from HCC samples. Details of the experimental procedures are provided in the Supplementary Materials and Methods.
2.4 Tumor immune infiltration analysis
We used the single sample gene set enrichment analysis (ssGSEA) method (20, 21) and TIMER database (22, 23) to investigate the relationships between RAB22A expression and immune cell infiltration, as detailed in the Supplementary Materials and Methods.
2.5 Gene set enrichment analysis
Enrichment analyses of relevant functional pathways were performed using the GO and KEGG databases (Supplementary Tables 3, 4) and GSEA (Supplementary Tables 5, 6), as detailed in the Supplementary Materials and Methods (24, 25).
2.6 Prediction and construction of ceRNA networks
Multiple databases were used to predict and screen the lncRNA-miRNA-mRNA (RAB22A) ceRNA network online. Details are provided in the Supplementary Materials and Methods.
2.7 Protein interaction network and module analysis
We created the protein–protein interaction (PPI) network using the Search Tool for the Retrieval of Interacting Genes (STRING) database (Supplementary Table 6) (26, 27), as detailed in the Supplementary Materials and Methods.
2.8 Statistical analysis
The R package (version 3.6.3) was used for statistical analyses and plotting. RAB22A expression in unpaired and paired samples was analyzed using the Wilcoxon rank sum test, and Wilcoxon signed rank test, respectively, with the pROC (1.17.0.1) package for ROC analysis. The RAB22A expression level was analyzed by querying the GEO, TIMER, and UALCAN databases (18). Using the KM method and log-rank test, we compared the differences in 10-year OS, DSS, and PFI between patients with high RAB22A expression and those with low RAB22A expression in TCGA. Cox analysis was used to determine the correlation between RAB22A expression and clinical features. p < 0.05 was considered to indicate significance.
3 Results
3.1 RAB22A is upregulated in HCC
First, we examined the RAB22A expression levels in different malignancies by assessing TCGA databases. RAB22A was highly expressed in 33 malignant tumors, including HCC (Figure 1A). In addition, RAB22A was highly expressed in the GEO datasets GSE121248, GSE87630, GSE76427, GSE84005, GSE57957, and GSE39791 HCC samples (p < 0.001) (Figures 1B–G). Western blot analysis of human normal liver cells (L02) and HCC cells (Hep G2, SK-Hep1, Huh7, HCCLM3, and MHCC97-H) validated the high expression of RAB22A in HCC cell lines (Figure 1H). The same results were obtained through qRT-PCR (p < 0.001) (Figure 1K). Next, we extracted 30 pairs of proteins from HCC and adjacent tissues and analyzed them using western blotting, which revealed that RAB22A was highly expressed in the former (Figure 1I). Results of western blot analysis of the liver and adjacent tissues are shown in Supplementary Figure 1. The high RAB22A mRNA expression levels in HCC tissues were further substantiated using qRT-PCR (p < 0.001) (Figure 1L). Immunohistochemistry (IHC) results also verified that RAB22A was upregulated in HCC tissues (Figure 1J). Finally, a receiver operating characteristic (ROC) curve was created. The ROC curve enclosed by the axes is the area below the curve (AUC). The AUC for RAB22A was 0.891, suggesting its remarkable diagnostic value for HCC (Figure 1M).
Figure 1 Expression level of RAB22A in HCC was verified in TCGA and GEO databases and in vitro experiments. (A) Comparison of the expression levels of RAB22A in different cancerous and normal tissues. (B–G) CEO database analysis of RAB22A expression in HCC tissues. (H) Western blotting assay of RAB22A protein expression levels in L02, Hep G2, SK-Hep1, Huh7, HCCLM3, and MHCC97-H cell lines. (I) Western blotting assay of RAB22A protein expression levels in HCC and adjacent tissues. (J) RAB22A protein levels in normal liver and HCC were measured using IHC. (K) qRT-PCR assay of RAB22A mRNA expression levels in L02, Hep G2, SK-Hep1, Huh7, HCCLM3, and HCCH97-H cell lines. (L) qRT-PCR assay of RAB22A mRNA expression levels in 30 pairs of HCC and adjacent tissues. (M) ROC curves were created to investigate the value of RAB22A in identifying HCC tissues. *p < 0.05, **p < 0.01, ***p < 0.001, NS, no significance.
3.2 Association of RAB22A expression with clinical characteristics
Using the UALCAN database to perform subgroup analysis of numerous pathological characteristics, we found that RAB22A transcript levels were elevated in patients with HCC. (Figure 2A). The sub-group analysis of cancer stage, ethnicity, sex, age, weight, tumor grade, and TP53 mutation showed that the expression of RAB22A in HCC patients was significantly higher than that in the normal group (Figures 2B–H).
Figure 2 Box plot showing the relative transcription of RAB22A in individual cancer stages, race, gender, age, weight, tumor grade, and TP53 mutation status in a subgroup of patients with HCC. (A) RAB22A in normal and HCC tissues. (B) RAB22A in normal individuals or in patients with stages 1–4 liver cancer. (C) RAB22A in normal and LIHC samples based on patient ethnicity. (D) RAB22A in normal individuals and males and females with HCC. (E) RAB22A in healthy subjects of any age and patients aged 21–40, 41–60, 61–80, and 81–100 years with HCC. (F) RAB22A in healthy subjects of any weight and normal weight patients, extreme weight patients, obese patients, and extremely obese patients. (G) RAB22A in normal subjects and patients with different liver cancer tumor grades. (H) RAB22A in normal and TP53-mutant or TP53-non mutant patients. *p < 0.05, **p < 0.01, ***p < 0.001, NS, no significance.
Logistic regression analysis showed that the increased expression of RAB22A in HCC was significantly correlated with sex (OR = 0.627 for male vs. female, p = 0.036), weight (OR = 0.567 for weight > 70kg vs. ≤ 70kg, p = 0.009), histological grades (OR=1.611 for G3 and G4 vs G1 and G2, p = 0.028), and tumor status (OR = 1.619 for with tumors vs. tumor free, p = 0.026). Conversely, RAB22A expression was not associated with age, M stage, T stage, N stage, height, BMI, AFP, or vascular invasion (Table 1).
Table 1 Association between RAB22A expression and clinicopathologic parameters by Logistic regression.
Next, we collected data from TCGA database to determine the clinicopathological parameters of RAB22A in different patients with HCC. Detailed information on the clinical data is provided in Table 2. After excluding cases without the necessary clinical data, 374 cases with a median age of 61.5 (range: 49.25−70.00) years and male preponderance of 67% were included. High expression of RAB22A in HCC was positively associated with tumor status (tumor-free vs. with tumor, p = 0.033), sex (female vs. male, p = 0.047), weight (≤ 70 vs. > 70, p = 0.012), and histological grade (grades 3 and 4 vs. grades 1 and 2, p = 0.031). These results indicate that the overexpression of RAB22A in HCC is closely related to the clinicopathological characteristics.
Table 2 Correlation between clinicopathological variables and RAB22A expression.
3.3 Prognostic value of RAB22A in HCC
Kaplan–Meier survival curves were analyzed to determine the connection between RAB22A expression and overall survival (OS), disease-free survival (DSS), and progression-free interval (PFI) in the prognosis of patients with HCC. Increased levels of RAB22A expression were inversely related to prognosis (Figures 3A–C). Additionally, subgroup analysis was performed on patients with low RAB22A expression and AFP < 400, and these patients had better OS, DSS, and PFI prognosis (Figures 3D–F). However, the groups with AFP (ng/mL) > 400 showed no significant differences (Supplementary Figures 2A–C). The high expression of RAB22A in stage M0 liver cancer was associated with poor OS, DSS, and PFI in a subgroup of patients (Figures 3G–I). The subgroups of T3 versus T4, stages III vs. IV, and tumor versus tumor-free status had significantly worse OS (Supplementary Figures 1D–F). Finally, we compared predictive variables in patients with HCC obtained by univariate regression analysis to those obtained via multivariate survival analysis (OS) (Supplementary Table 7). Pathologic stage (stages I and II compared with stages III and IV; p <0.001), tumor size (T stages 1 and 2 versus T stages 3 and 4; p < 0.001), metastatic spread (M stages 0 and 1; p = 0.017), and tumor status (without or with tumor; p < 0.001) were highly significant in the univariate analysis. The multivariate analysis showed that with tumor (p = 0.014) was significant, suggesting that it is an independent risk factor.
Figure 3 Kaplan–Meier survival plots comparing the relationship between RAB22A and prognosis in HCC. (A–C) Survival curves of OS, DSS, and PFI between RAB22A-high and -low patients with HCC. (D–F) OS, DSS, and PFI survival curves of patients with HCC with high and low RAB22A expression of AFP (ng/mL) ≤ 400. (G–I) Survival curves comparing OS, DSS, and PFI in patients with HCC at the M0 stage with high and low expression of RAB22A.
3.4 GSEA and GO/KEGG enrichment analyses
GO and KEGG pathway co-expression analyses of RAB22A-related genes in liver cancer mRNA sequencing data with 371 patients from the TCGA were performed using the functional module of Linkedomics. The top 50 marker genes and their connections with RAB22A expression are displayed on the heat map (Figures 4A, B; Supplementary Table 8). These findings revealed a widespread effect of RAB22A on the transcriptome.
Figure 4 Enrichment of biofunction and associated gene analysis of RAB22A in HCC. (A, B) Heat map showing genes positively and negatively associated with RAB22A in liver cancer (top 50). Positively associated genes are indicated in red, while negatively associated genes are in green. (C) The enriched terms in GO categories in HCC. (D) KEGG pathway analysis based on RAB22A-associated DEGs. (E) The five most positively correlated pathways were revealed by GO term analysis. (F) KEGG pathway analysis revealed the five most positively correlated pathways. (G) The five most positively correlated pathways were identified via REACTOME pathway analysis. (H) The five most negatively correlated pathways were identified via GO term analysis. (I) KEGG pathway analysis identified the five most negatively correlated pathways. (J) The five most negatively correlated pathways were identified via REACTOME pathway analysis.
Next, we conducted an enrichment analysis using the GO and KEGG databases to support the concept that RAB22A-related DEGs play a biological role in HCC (Figures 4C, D). The results of GO analysis showed that these DEGs were related to biological processes (BP), cellular components (CC), and molecular functions (MF). In the GO analysis, DEGs were enriched in diverse biological pathways, including proteasomal protein catabolic process, ncRNA processing, ribosome ribosomal subunit, protein serine/threonine kinase activity, and protein serine kinase activity. In the KEGG analysis, DEGs were highly concentrated in endocytosis and non-alcoholic fatty liver disease. GSEA was used to analyze the biological functions related to RAB22A expression.
The later criteria were enrichment score | NSE | > 1 (p < 0.05), according to which the five most positively relevant signal pathways were selected. GO analysis revealed that RAB22A expression was strongly positively correlated with the processes of homophilic cell adhesion via plasma membrane adhesion, immunoglobulin, T cell receptor, and plasma membrane signaling receptor complex (Figure 4E). The expression of RAB22A was inversely linked to that of co-translational proteins that bind to the membrane, the cytosolic ribosome, the structural components of the ribosome, the ribosomal subunit, and the nonsense-mediated decay of nuclear-transcribed mRNA catabolic processes (Figure 4H). KEGG analysis revealed that RAB22A expression was most strongly negatively connected with axon guidance, extracellular matrix receptor interaction, focal adhesion, FCγR-mediated phagocytosis, and the interaction with neuroactive ligand receptors (Figure 4F). The ribosome, Parkinson’s disease, retinol metabolism, oxidative phosphorylation, and complement and coagulation cascades were the top five most negatively correlated pathways (Figure 4I). REACTOME pathway analysis determined that phospholipids play a role in phagocytosis, Fc gamma receptor FCGR-dependent phagocytosis, CD22-mediated B cell receptor (BCR) regulation, and FCGR activation, and that second messenger is activated by BCR antigen that all positively correlated with RAB22A expression (Figure 4G). Eukaryotic translation initiation, eukaryotic translation elongation, translocation, the response of eukaryotic initiation factor 2 alpha subunit kappa B cyclin N2 to amino acid deprivation, co-translational protein of SRP-dependent targeting to the membrane, and nonsense-mediated decay were all negatively correlated with RAB22A expression (Figure 4J).
3.5 PPI network analysis
The PPI network of co-expressed genes conforming to the STRING conditions was assembled and visualized using Cytoscape, and analysis of the interactions among 108 DEGs in the HCC group was conducted. A total of 51 proteins and 534 edges were screened (Figure 5A; Supplementary Table 9).
Figure 5 PPI network enrichment analysis. (A) The PPI network was built based on PPI pairs identified by the STRING dataset. (B) Hub gene clusters were selected from the PPI network (criteria of total scores ≥ 14,000). (C) Top 10 hub genes in the PPI network.
After screening 12 nodes and 212 edges, a primary gene cluster with a total score ≥ 14,000 was discovered (Figure 5B). Finally, we screened the top 10 central genes, namely RAB22A, RABGEF1, VPS45, VPS18, VPS11, MON1A, VPS39, VPS16, ZFYV20, and VPS8 (Figure 5C).
3.6 Role of RAB22A and m6A methylation regulators in HCC
M6A methylation affects the development of HCC (24–27). The expression of RAB22A was compared with that of the 23 M6A methylation genes reported in the literature to verify this conclusion (Figure 6A). RAB22A expression was closely connected with that of the 23 m6A-related genes in HCC (Figures 6B–X). Moreover, groups were formed according to RAB22A median expression. By analyzing the differences in the 23 m6A methylation genes in RAB22A between the high- and low-expression groups of patients with HCC, we observed that the expression levels of all genes in the RAB22A high-expression group were upregulated (Figure 6Y). Overall, we observed an obvious relationship between m6A methylation and RAB22A expression levels in HCC.
Figure 6 Correlation analysis of RAB22A expression levels with m6A-related gene expression in HCC tissues. (A) Correlation of RAB22A expression levels with m6A gene expression in HCC. (B–X) Scatter plot showing the relationship between RAB22A and the m6A gene. Differences in 23 M6A-related genes between the RAB22A high-expression group and RAB22A low-expression group in liver cancer patients (Y).**p < 0.01, ***p < 0.001, NS, no significance.
3.7 Construction of a triple regulatory network for RAB22A-associated ceRNA
Increasing evidence has demonstrated the regulatory effect on the lncRNA-miRNA-mRNA ceRNA network in HCC. The Venn diagram showed 41 overlapping miRNAs in the Targerscan, starBase, and MiRDB databases (Figure 7A). Five human-derived miRNAs (miR-328-3p, miR-3163, miR-2114-5p, miR-664b-3p, and miR-204-5p) were verified to negatively correlate with RAB22A expression (Figure 7B). The expression of RAB22A and target microRNAs is displayed as a scatter plot (Figures 7C–G). We consulted the Rnalnter and starBase databases to predict lncRNAs that can have a mutual effect on target miRNAs (miR-204-5p and miR-328-3p) (Figures 7H–I). The expression levels of lncRNAs and miRNAs were inversely correlated, which accounted for the mutual influence between the two. We used the starBase database to filter and identify lncRNAs that were adversely associated with the two target miRNAs in HCC. Nine HCC-related ceRNA regulatory networks were constructed (Figure 7J).
Figure 7 Prediction of ceRNA networks in HCC. (A) Venn diagram results showing 41 overlapping miRNAs in Targerscan, starBase, and MiRDB databases. (B) Five miRNAs screened for negative correlation with RAB22A expression. (C–G) Scatter plots showed that miRNAs were significantly correlated with mRNAs. (H, I) Prediction of lncRNAs bound to target miRNAs using miRNet and starBase online databases and displayed as a Venn diagram, including hsa-miR-204-5p and hsa-miR-328-3p. (J) Sankey diagram showing the RAB22A-related ceRNA regulatory network.
3.8 Association of RAB22A expression with immune cell infiltration
Using the ssGSEA method, we verified the strong connection between RAB22A and immune cells (Figure 8A). The expression of RAB22A was positively connected with T helper cells, Tcm cells, and Th2 cells (p < 0.001) but negatively with cytotoxic cells, DCs, and pDCs (p < 0.001) (Figures 8B–G). RAB22A may be heavily involved in the T-cell immune response to HCC. Moreover, RAB22A expression in HCC correlated with various immune cell markers (Table 3). In the M2 macrophages in HCC, we found that RAB22A expression was substantially relevant to the expression of the immunological markers CD163, VSIG4, and MS4A4A. These results indicate that RAB22A caused the macrophages in HCC to adopt an M2 phenotype. The expression of RAB22A was substantially linked to 66 immunological markers, including CD8A, CD3D, and T-bet, in an analysis of functional T-cell immunity indicators. Furthermore, RAB22A expression was linked to immunological markers for B cells, T cells, TAMs, and neutrophils (Table 3). The TIMER database was utilized to determine whether RAB22A expression in HCC was connected with immune cell invasion levels. The results indicated that the CNV of RAB22A was related to the level of neutrophil infiltration. (Figure 8H). Subsequently, the infiltration of macrophages, T helper cells, Tcm, and Th2 cells increased (p < 0.001) in the RAB22A high-expression group; however, cytotoxic cells, DCs, and pDCs decreased (p < 0.001) (Figure 8I). These results verified that the increased expression of RAB22A in HCC is tightly linked with the infiltration of immune cells.
Figure 8 Relationship between the expression of RAB22A and microenvironment of immune infiltrating cells in HCC. (A) Forest plot depicting the relationship between RAB22A expression levels and the relative abundance of the 24 immune cells. (B–G) Scatter plots showing the degree of differentiation of pDCs, T helper cells, DCs, Th2 cells, cytotoxic cells, and Tcm cells between the high and low RAB22A expression groups. (H) SCNA showed that the expression of RAB22A correlated with the degree of immune cell infiltration. (I) Scatter plot showing the correlation of 24 immune cells with RAB22A expression levels. *p < 0.05, **p < 0.01, ***p < 0.001, NS, no significance.
Table 3 Correlation analysis of RAB22A expression with immune cell biomarkers.
4 Discussion
RAB22A is a member of the RAS oncogene family that controls membrane properties and vesicle budding, delamination, movement, and fusion and is central to ensuring that cargo is transported to its correct destination. RAB22A is referred to in the early formation of endosomes and regulates vesicle transport (28, 29).
Furthermore, RAB22A is a critical oncogene that has a crucial impact on the course of many different forms of cancer (12, 30). RAB22A promotes the epithelial–mesenchymal transition of papillary thyroid cancer cells, thereby promoting their proliferation, migration, and invasion (31). CD147 is recycled by RAB22A to control lung carcinoma cell motility and invasion (13). In metastatic breast cancer, hypoxia facilitates MV production and HIF-dependent RAB22A gene expression (14). In addition, RAB22A is involved in a miRNA downregulation mechanism in which the overexpression of small GTPases promotes tumor growth and carcinogenesis. Several tumor models, including kidney, colorectal, glioma, and bile duct cancer, have utilized RAB22A as a target gene for miRNAs (32–34). Changes in RAB22A in HCC may be significant as hepatocytes always maintain high metabolic levels and active vesicular transport; nevertheless, the potential effect on RAB22A in HCC is unclear.
In the present study, we first found that RAB22A was upregulated in HCC and various malignant tumors by analyzing multiple databases. Subsequently, we verified the elevation of RAB22A expression in HCC cell lines and HCC samples using western blotting, qRT-PCR, and IHC in vitro. Overexpression of RAB22A in HCC tissues was closely associated with clinicopathologic features. The ROC curve analysis suggested RAB22A as a promising diagnostic biomarker for differentiating HCC from normal tissues. Moreover, the overexpression of RAB22A was interrelated with a poor prognosis of HCC, as indicated by OS, DSS, and PFI.
To elucidate the potential biological functions and regulatory pathways of RAB22A, we investigated genes encoding RAB22A-related proteins and co-expression genes in HCC tissues. mRNA sequencing data with HCC were evaluated in the TCGA database, while the DEGs associated with RAB22A in HCC were shown in a heat map. Insights gained from pathway enrichment analyses using GO and KEGG indicated that RAB22A has far-reaching effects on the transcriptome. Through enrichment pathway analysis, we verified that these DEGs were involved in proteasomal protein catabolic process, ncRNA processing, ribosomes, and ribosomal subunits, protein serine/threonine kinase activity, GTPase combining, herpes simplex virus type 1 infection, multiple neurodegenerative illnesses, and Alzheimer’s disease pathways. Next, we analyzed 30 signaling pathways positively and negatively correlated with RAB22A expression using GSEA. Overexpression of RAB22A was linked to processes such as cell adhesion via the plasma membrane (35), nonsense-mediated decay of nuclear-transcribed mRNA (36), axon guidance (37), oxidative phosphorylation (38), FCGR activation (39), and eukaryotic translation elongation (40) in a GSEA of HCC. Overall, we suggest that RAB22A may participate in various cellular immune functions and intracellular transport and may facilitate the advance of HCC by adjusting these signaling pathways.
Subsequently, we built a PPI network using Cytoscape. One central gene cluster (total score ≥ 14,000) and the top 10 central genes were filtered, namely RAB22A, RABGEF1, VPS45, VPS18, VPS11, MON1A, VPS39, VPS16, ZFYV20, and VPS. These findings provide important insights for subsequent study designs and experimental validations.
M6A methylation has been examined to elucidate the mechanisms of HCC since it has been proven to affect cancer via numerous mechanisms (41). m6A is a critical player in HCC (42, 43). Methyltransferases (the “Writers”), demethylases (the “Erasers”), and methylated reading proteins have access to the same m6A methylation (Readers). Methylation transferases, such as METTL3/14, WTAP, and KIAA1429, are primarily responsible for catalyzing the m6A alteration of adenosine on mRNA. Demethylases, such as FTO and ALKHB5, facilitate the demethylation of m6A. Methylation reading proteins, such as YTHDF 1-3 and YTHDC 1-3, recognize RNA methylation and play a role in regulatory processes, such as RNA translation, degradation, and miRNA processing (44). Further analysis of the connection between RAB22A expression and m6A methylation proteins revealed a positive and significant association between RAB22A expression and the expression of methylation transferases, demethylases, and methylated reading proteins. Patients with HCC have a poor prognosis because m6A-modified proteins are highly elevated in the disease, and their overexpression increases the disease progression. Several reports have verified that IGF2BP1, YTHDF1, and RBM15 are all highly elevated in HCC and contribute to its development and progression. These findings indicate that m6A may alter the RAB22A gene to enhance the consistency of its mRNA, hence boosting the occurrence and development of HCC. Evidence for lncRNA-miRNA-mRNA ceRNA networks’ regulatory role in cancers is mounting (45). Based on these predictions, we constructed a ceRNA regulatory network that predicted that RAB22A might affect several critical pathways of HCC regulatory mechanisms. We intend to conduct further experiments to validate this network.
Cancer cells that invade Immune cells, known as tumor-infiltrating immune cells (TIICs), play a key regulatory role in tumorigenesis and development (46). The HCC prognosis may be affected by the presence of TIICs, which are essential for HCC development (28–30). TIICs facilitate a tangled web of cellular interactions that boost the immunosuppressive milieu, facilitate immune escape, and ultimately aid in tumor progression. Changes in the immune environment of the liver can cause liver lesions, such as chronic inflammation and fibrosis/cirrhosis (22, 47). RAB22A is a regulator of immune functions. Independent studies have also shown that Th2 cells contribute to cancer development and progression (48, 49). Effector T helper cell subgroups are essential for coordinating immune responses to diverse infections and participate in the nosogenesis of numerous inflammatory disorders, including autoimmunity and allergies (50). pDCs are a sentinel cell type that can test pathogen-derived nucleic acids and reactions via the rapid and significant production of type I interferons, primarily in autoimmune diseases, immune deficiencies, and cancer (51). Cytotoxic T cells and DCs are also essential effectors of antitumor immunity (40, 52). These findings suggest that RAB22A plays an indispensable role in regulating immune cell infiltration in HCC.
We also demonstrated that RAB22A expression was significantly correlated with 66 immune markers. These results indicate that the upregulation of RAB22A expression in HCC is linked to immune cell infiltration. Understanding the function of RAB22A in immune activation will help to facilitate future research using various immune cell types and animal models.
Although our study identified the molecular mechanism of RAB22A in HCC through bioinformatics analysis, there remain limitations. Firstly, to elucidate the effect of RAB22A on HCC, several subjective factors, such as the treatment details received by patients and follow-up, should be considered simultaneously. However, some experiments were conducted in different centers, thereby limiting the information or causing inconsistency in the public database, which led to some errors. Secondly, the number of patients with cancer in the experimental control group was different from that in the current study; hence accessional studies are needed to eliminate the error caused by sample offset.
Thirdly, multicenter investigations based on communal databases seek to compensate for the paucity of single-center studies. However, retrospective studies have drawbacks, including inconsistent interventions and a lack of data. Since this study is retrospective, prospective investigations should be undertaken to eliminate analytical bias. Based on previous verifications, the results are robust; advancements in single-cell and spatial transcriptomics technologies allowed for the increased use of single-cell multi-omics technologies to gain insights into complex cellular ecosystems and biological processes. Currently, there is a gap in the rapidly growing single-cell multi-omics data, while effective methods for comprehensive analysis of these inherently sparse and heterogeneous data are limited. Therefore, new algorithms, such as SMGR (53) and spaCI (54), have been derived to address this gap. Single-cell multi-omics gene co-regulation algorithms provide multiple regulatory stages to study the control of cellular heterogeneity and complex biological mechanisms, which provide great clinical value for identifying mechanisms, targets, and predictors to enhance translational therapy. The spaCI algorithm can detect upstream transcription factors (TFS) mediating the L-R signaling axis, which provides insights into the underlying molecular mechanisms of the intercellular crosstalk. These emerging algorithms can be used to verify the biological mechanism of RAB22A in HCC.
In summary, we demonstrated, to the best of our knowledge, for the first time that RAB22A promotes carcinogenesis via m6A methylation and ceRNA network processes and is strongly linked with HCC development, poor survival, and immune infiltration.
Data availability statement
The article/Supplementary material contains the original contributions presented in the study. Any additional questions can be forwarded to the corresponding authors. The datasets presented in this study can be found in online repositories. The names of the repository can be found below: https://portal.gdc.cancer.gov.
Ethics statement
The studies involving human participants were reviewed and approved by First Affiliated Hospital of Harbin Medical University’s ethics committee. The patients/participants provided their written informed consent to participate in this study.
Author contributions
FW, FM, XL, QL, and LJ contributed equally to this work. ZL and YC designed this research. FW and FM drafted this manuscript. XL , XW and QL performed the data collection and analysis. JL, RZ, YunZ, YuZ and SJ participated in the data interpretation and study design. All authors approved the final manuscript.
Funding
This work was supported by The National Natural Scientific Foundation of China (Grant Nos. 81972230), The Heilongjiang Postdoctoral Science Foundation (Grant No. LBH-Z20178), The Scientific Foundation of the First Affiliated Hospital of Harbin Medical University (Grant No. 2021B03), the Excellent Youth Science Fund of the First Affiliated Hospital of Harbin Medical University (Grant No. 2021Y01), and the Chen Xiaoping Foundation for the Development of Science and Technology of Hubei Province (CXPJJH122002-092).
Acknowledgments
We thank the patients who provided the experimental samples and the teachers for providing our funding.
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/fimmu.2023.1086342/full#supplementary-material
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J Clin (2021) 71:209–49. doi: 10.3322/caac.21660
2. Konyn P, Ahmed A, Kim D. Current epidemiology in hepatocellular carcinoma. Expert Rev Gastroenterol Hepatol (2021) 15:1295–307. doi: 10.1080/17474124.2021.1991792
3. Cui Y, Sun D, Liu L. Upregulation of cystatin SN promotes hepatocellular carcinoma progression and predicts a poor prognosis. J Cell Physiol (2019) 234:22623–34. doi: 10.1002/jcp.28828
4. Stenmark H. Rab GTPases as coordinators of vesicle traffic. Nat Rev Mol Cell Biol (2009) 10:513–25. doi: 10.1038/nrm2728
5. Antonyak MA, Wilson KF, Cerione RAR. R(h)oads to microvesicles. Small GTPases (2012) 3:219–24. doi: 10.4161/sgtp.20755
6. Mesa R, Magadán J, Barbieri A, López C, Stahl PD, Mayorga LS. Overexpression of Rab22a hampers the transport between endosomes and the golgi apparatus. Exp Cell Res (2005) 304:339–53. doi: 10.1016/j.yexcr.2004.11.017
7. Weigert R, Yeung AC, Li J, Donaldson JG. Rab22a regulates the recycling of membrane proteins internalized independently of clathrin. Mol Biol Cell (2004) 15:3758–70. doi: 10.1091/mbc.e04-04-0342
8. Kauppi M, Simonsen A, Bremnes B, Vieira A, Callaghan J, Stenmark H, et al. The small GTPase Rab22 interacts with EEA1 and controls endosomal membrane trafficking. J Cell Sci (2002) 115:899–911. doi: 10.1242/jcs.115.5.899
9. Sun L, He M, Xu N, Xu D-H, Ben-David Y, Yang Z-Y, et al. Regulation of RAB22A by mir-193b inhibits breast cancer growth and metastasis mediated by exosomes. Int J Oncol (2018) 53:2705–14. doi: 10.3892/ijo.2018.4571
10. Yin Y, Zhang B, Wang W, Fei B, Quan C, Zhang J, et al. miR-204-5p inhibits proliferation and invasion and enhances chemotherapeutic sensitivity of colorectal cancer cells by downregulating RAB22A. Clin Cancer Res (2014) 20:6187–99. doi: 10.1158/1078-0432.CCR-14-1030
11. Zheng S, Jiang F, Ge D, Tang J, Chen H, Yang J, et al. LncRNA SNHG3/miRNA-151a-3p/RAB22A axis regulates invasion and migration of osteosarcoma. BioMed Pharmacother (2019) 112:108695. doi: 10.1016/j.biopha.2019.108695
12. Zhou Y, Wu B, Li J-H, Nan G, Jiang J-L, Chen Z-N. RAB22A enhances CD147 recycling and is required for lung cancer cell migration and invasion. Exp Cell Res (2017) 357:9–16. doi: 10.1016/j.yexcr.2017.04.020
13. Wang T, Gilkes DM, Takano N, Xiang L, Luo W, Bishop CJ, et al. Hypoxia-inducible factors and RAB22A mediate formation of microvesicles that stimulate breast cancer invasion and metastasis. Proc Natl Acad Sci USA (2014) 111:E3234–42. doi: 10.1073/pnas.1410041111
14. Mayorga LS, Cebrian I. RAB22A: A novel regulator of immune functions. Mol Immunol (2019) 113:87–92. doi: 10.1016/j.molimm.2018.03.028
15. Johnson DC, Baines JD. Herpesviruses remodel host membranes for virus egress. Nat Rev Microbiol (2011) 9:382–94. doi: 10.1038/nrmicro2559
16. Isberg RR, O’Connor TJ, Heidtman M. The legionella pneumophila replication vacuole: making a cosy niche inside host cells. Nat Rev Microbiol (2009) 7:13–24. doi: 10.1038/nrmicro1967
17. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: A portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (2017) 19:649–58. doi: 10.1016/j.neo.2017.05.002
18. Liu J, Lichtenberg T, Hoadley KA, Poisson LM, Lazar AJ, Cherniack AD, et al. An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell (2018) 173:400–16.e11. doi: 10.1016/j.cell.2018.02.052
19. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res (2012) 41:D991–5. doi: 10.1093/nar/gks1193
20. Hänzelmann 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
21. Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity (2013) 39:782–95. doi: 10.1016/j.immuni.2013.10.003
22. Li B, Severson E, Pignon J-C, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol (2016) 17:174. doi: 10.1186/s13059-016-1028-7
23. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res (2017) 77:e108–e10. doi: 10.1158/0008-5472.CAN-17-0307
24. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
25. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an r package for comparing biological themes among gene clusters. OMICS (2012) 16:284–7. doi: 10.1089/omi.2011.0118
26. von Mering C. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res (2003) 31:258–61. doi: 10.1093/nar/gkg034
27. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal (2013) 6:pl1. doi: 10.1126/scisignal.2004088
28. He H, Dai F, Yu L, She X, Zhao Y, Jiang J, et al. Identification and characterization of nine novel human small GTPases showing variable expressions in liver cancer tissues. Gene Expr (2002) 10:231–42. doi: 10.3727/000000002783992406
29. Magadán JG, Barbieri MA, Mesa R, Stahl PD, Mayorga LS. RAB22A regulates the sorting of transferrin to recycling endosomes. Mol Cell Biol (2006) 26:2595–614. doi: 10.1128/MCB.26.7.2595-2614.2006
30. Wang J, Luo X, Lu J, Wang X, Miao Y, Li Q, et al. RAB22A promotes the proliferation, migration, and invasion of lung adenocarcinoma via up-regulating PI3K/Akt/mTOR signaling pathway. Exp Cell Res (2022) 416:113179. doi: 10.1016/j.yexcr.2022.113179
31. Luo X, Wang J, Lu J, Wang X, Miao Y, Li Q, et al. RAB22A promotes epithelial-mesenchymal transition in papillary thyroid carcinoma by activating PI3K/AKT/mTOR signaling pathway. BioMed Res Int (2022) 2022:1–17. doi: 10.1155/2022/1874550
32. Liu Q, Lei C. LINC01232 serves as a novel biomarker and promotes tumour progression by sponging miR-204-5p and upregulating RAB22A in clear cell renal cell carcinoma. Ann Med (2021) 53:2153–64. doi: 10.1080/07853890.2021.2001563
33. Kong W, Li H, Xie L, Cui G, Gu W, Zhang H, et al. LncRNA MCF2L-AS1 aggravates the malignant development of colorectal cancer via targeting miR-105-5p/RAB22A axis. BMC Cancer (2021) 21:1069. doi: 10.1186/s12885-021-08668-w
34. Cui Y, Liang S, Liu L. ABCA8 is regulated by miR-374b-5p and inhibits proliferation and metastasis of hepatocellular carcinoma through the ERK/ZEB1 pathway. J Exp Clin Cancer Res (2020) 39:90. doi: 10.1186/s13046-020-01591-1
35. Piao G, Yuchao H, Lu C, Lisha Q, Dongming L, Ziye C, et al. Cytosolic phospholipase A2α modulates cell-matrix adhesion via the FAK/paxillin pathway in hepatocellular carcinoma. Cancer Biol Med (2019) 16:377. doi: 10.20892/j.issn.2095-3941.2018.0386
36. Zhang H, You Y, Zhu Z. The human RNA surveillance factor up-frameshift 1 inhibits hepatic cancer progression by targeting MRP2/ABCC2. BioMed Pharmacother (2017) 92:365–72. doi: 10.1016/j.biopha.2017.05.090
37. Liang JQ, Teoh N, Xu L, Pok S, Li X, Chu ESH, et al. Dietary cholesterol promotes steatohepatitis related hepatocellular carcinoma through dysregulated metabolism and calcium signaling. Nat Commun (2018) 9:4490. doi: 10.1038/s4146-018-06931-6
38. Yu H, Wang C, Ke S, Bai M, Xu Y, Lu S, et al. Identification of CFHR4 as a potential prognosis biomarker associated with lmmune infiltrates in hepatocellular carcinoma. Front Immunol (2022) 13:892750. doi: 10.3389/fimmu.2022.892750
39. Kudo Y, Sugimoto M, Arias E, Kasashima H, Cordes T, Linares JF, et al. PKCλ/ι loss induces autophagy, oxidative phosphorylation, and NRF2 to promote liver cancer progression. Cancer Cell (2020) 38:247–62.e11. doi: 10.1016/j.ccell.2020.05.018
40. Oh DY, Fong L. Cytotoxic CD4+ T cells in cancer: Expanding the immune effector toolbox. Immunity (2021) 54:2701–11. doi: 10.1016/j.immuni.2021.11.015
41. Sun T, Wu R, Ming L. The role of m6A RNA methylation in cancer. BioMed Pharmacother (2019) 112:108613. doi: 10.1016/j.biopha.2019.108613
42. Wang P, Wang X, Zheng L, Zhuang C. Gene signatures and prognostic values of m6A regulators in hepatocellular carcinoma. Front Genet (2020) 11:540186. doi: 10.3389/fgene.2020.540186
43. Zhou T, Li S, Xiang D, Liu J, Sun W, Cui X, et al. m6A RNA methylation-mediated HNF3γ reduction renders hepatocellular carcinoma dedifferentiation and sorafenib resistance. Sig Transduct Target Ther (2020) 5:296. doi: 10.1038/s41392-020-00299-0
44. Chen X-Y, Zhang J, Zhu J-S. The role of m6A RNA methylation in human cancer. Mol Cancer (2019) 18:103. doi: 10.1186/s12943-019-1033-z
45. Qi X, Zhang D-H, Wu N, Xiao J-H, Wang X, Ma W. ceRNA in cancer: possible functions and clinical implications. J Med Genet (2015) 52:710–8. doi: 10.1136/jmedgenet-2015-103334
46. Garnelo M, Tan A, Her Z, Yeong J, Lim CJ, Chen J, et al. Interaction between tumour-infiltrating b cells and T cells controls the progression of hepatocellular carcinoma. Gut (2017) 66:342–51. doi: 10.1136/gutjnl-2015-310814
47. Ruf B, Heinrich B, Greten TF. Immunobiology and immunotherapy of HCC: spotlight on innate and innate-like immune cells. Cell Mol Immunol (2021) 18:112–27. doi: 10.1038/s41423-020-00572-w
48. Orecchioni M, Bedognetti D, Newman L, Fuoco C, Spada F, Hendrickx W, et al. Single-cell mass cytometry and transcriptome profiling reveal the impact of graphene on human immune cells. Nat Commun (2017) 8:1109. doi: 10.1038/s41467-017-01015-3
49. Mantovani A, Romero P, Palucka AK, Marincola FM. Tumour immunity: effector response to tumour and role of the microenvironment. Lancet (2008) 371:771–83. doi: 10.1016/S0140-6736(08)60241-X
50. Zhu J. T Helper cell differentiation, heterogeneity, and plasticity. Cold Spring Harb Perspect Biol (2018) 10:a030338. doi: 10.1101/cshperspect.a030338
51. Reizis B. Plasmacytoid dendritic cells: Development, regulation, and function. Immunity (2019) 50:37–50. doi: 10.1016/j.immuni.2018.12.027
52. Martinek J, Wu T-C, Cadena D, Banchereau J, Palucka K. Interplay between dendritic cells and cancer cells. Int Rev Cell Mol Biol (2019) 348:179–215. doi: 10.1016/bs.ircmb.2019.07.008
53. Song Q, Zhu X, Jin L, Chen M, Zhang W, Su J. SMGR: a joint statistical method for integrative analysis of single-cell multi-omics data. NAR Genom Bioinform (2022) 4:lqac056. doi: 10.1093/nargab/lqac056
Keywords: RAB22A, hepatocellular carcinoma, cancer immune infiltrates, prognosis, biomarker, bioinformatics analysis
Citation: Wen F, Meng F, Li X, Li Q, Liu J, Zhang R, Zhao Y, Zhang Y, Wang X, Ju S, Cui Y and Lu Z (2023) Characterization of prognostic value and immunological roles of RAB22A in hepatocellular carcinoma. Front. Immunol. 14:1086342. doi: 10.3389/fimmu.2023.1086342
Received: 01 November 2022; Accepted: 20 February 2023;
Published: 03 March 2023.
Edited by:
Hao Zhang, Chongqing Medical University, ChinaCopyright © 2023 Wen, Meng, Li, Li, Liu, Zhang, Zhao, Zhang, Wang, Ju, Cui and Lu. 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: Yifeng Cui, cui88963342@hrbmu.edu.cn; Zhaoyang Lu, lzy76772005@hrbmu.edu.cn
†These authors have contributed equally to this work