Molecular Associations and Clinical Significance of RAPs in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is an aggressive gastrointestinal malignancy with a high rate of mortality. Multiple studies have individually recognized members of RAP gene family as critical regulators of tumor progression in several cancers, including hepatocellular carcinoma. These studies suffer numerous limitations including a small sample size and lack of analysis of various clinicopathological and molecular features. In the current study, we utilized authoritative multi-omics databases to determine the association of RAP gene family expression and detailed molecular and clinicopathological features in hepatocellular carcinoma (HCC). All five RAP genes were observed to harbor dysregulated expression in HCC compared to normal liver tissues. RAP2A exhibited strongest ability to differentiate tumors from the normal tissues. RAP2A expression was associated with progressive tumor grade, TP53 and CTNNB1 mutation status. Additionally, RAP2A expression was associated with the alteration of its copy numbers and DNA methylation. RAP2A also emerged as an independent marker for patient prognosis. Further, pathway analysis revealed that RAP2A expression is correlated with tumor-infiltrating immune cell composition and oncogenic molecular pathways, such as cell cycle and cellular metabolism.

Introduction

Hepatocellular carcinoma (HCC) is the sixth leading cancer in incidence and the fourth most common cause of cancer mortality in the world (Bray et al., 2018). It is the most common type of primary liver cancer that usually arises on the background of chronic liver disease, hepatitis B or C virus infection, or nonalcoholic steatohepatitis (Bruix et al., 2011; Villanueva, 2019). For locally advanced cancers without cirrhosis, the 5-years survival rate of patients is only 36–70% and 60–70% with successful surgical resection or liver transplantation, respectively. Further, postoperative recurrence and metastasis are common in HCC, which pose a challenge in the management of this disease. Therefore, biomarkers to predict prognosis in HCC are highly needed. The common indicators of prognosis of HCC include tumor size, degree of cirrhosis, tumor differentiation and microvascular invasion (Villanueva, 2019). The recent emergence of high throughput sequencing data by multiple studies has enabled researchers to describe molecular features of HCC in detail and has provided several potential biomarkers for the prediction of patient prognosis (Wheeler and Roberts, 2017).

RAP proteins (Ras proximate proteins) are members of the Ras GTP binding family sharing 50–60% sequence homology with the Ras family. The diversity and specificity of Ras and RAP proteins are contributed by different sets of GEFs (guanine nucleotide exchange factors) and GAPs (GTPase-activating proteins). Five different genes of the RAP family, RAP1A, RAP1B, RAP2A, RAP2B, and RAP2C have been identified in the human genome (Bokoch, 1993). RAP proteins primarily function in cell adhesion, migration, and polarity (Bokoch, 1993; Ehrhardt et al., 2002; Di et al., 2015b; Qu et al., 2016; Meng et al., 2018). The effect of RAP activation depends on the context-specific interaction of RAP with their regulators and downstream effectors.

Oncogenic functions of RAP proteins have been well established in multiple cancer types, such as breast (Di et al., 2015a), lung (Fu et al., 2009; Wu et al., 2014; Xie et al., 2015; Peng et al., 2016), ovary (Che et al., 2015; Lu et al., 2016), stomach (Zhang J. et al., 2020), cervix (Li et al., 2018), prostate (Bigler et al., 2007) and brain (Wang et al., 2017). Accumulating evidence suggests that RAP proteins also play critical roles in hepatocellular carcinogenesis and tumor progression. Single nucleotide polymorphism (SNPs) in RAP1A gene rs494453 has been shown to associate with higher incidence and recurrence after liver transplantation (Mo et al., 2018; Zhang R. et al., 2020). Further, higher activity of the NF-κB/RAP1 signaling pathway is associated with tumorigenicity in HCC cells (Mo et al., 2018). Some studies have also provided strong links between RAP1A expression and liver inflammation, a risk factor for liver carcinogenesis. RAPGEF1, the GEF for RAP1A has also been shown to be overexpressed in HCC (Sequera et al., 2018). A previous study reported that HBV replication promotes liver carcinogenesis through upregulation of RAP1B (Sheng et al., 2014). Further, overexpression of RAP1B enhances the proliferation and migration of HCC cells by regulating Twist-1 gene expression (Tang et al., 2018). Overexpression of RAB2B has also been reported in HCC and its inhibition reduces cell proliferation and invasion (Zhang et al., 2017). Recently, Zheng et al. reported that HCC tissues exhibit significantly higher mRNA and protein expression of RAP2A, which is associated with tumor size, metastasis, pathological differentiation, and vascular invasion (Zheng et al., 2017). Furthermore, they also demonstrated that higher protein levels of RAP2A are independently associated with poor overall survival in HCC.

While the current literature suggests that RAP genes might play critical roles in the pathophysiology of HCC, these studies are limited by determining individual genes of the RAP signaling pathway, limited number of clinical samples used in different studies. Further, studies focused on determining the association of RAP genes with genetic alteration and molecular alterations remain limited. In the current study, we utilized authoritative multi-omics databases to determine the association of RAP gene family expression and detailed molecular and clinicopathological features. Furthermore, we also determined their association with multiple survival parameters to determine their prognostic value.

Materials and Methods

Data Retrieval

For mRNA expression analysis, RNA seq data of TCGA-LIHC dataset, which was originally sourced from Broad GDAC Firehose (http://gdac.broadinstitute.org/) (Wheeler and Roberts, 2017) was extracted using UCSC XENA webserver (Goldman et al., 2020). Clinicopathological and molecular characterstics of the TCGA-LIHC dataset has been given in Table 1. Microarray gene expression data from multiple studies was accessed through the TNMplot webserver (Bartha and Győrffy, 2020). This web server hosts data from multiple HCC studies, where gene expression data has been normalized for all available studies and can be used for comparison between the collective groups of all normal samples with tumor samples. Multi-Omics dataset of hepatocellular carcinoma released by Clinical Proteomic Tumor Analysis Consortium (CPTAC) (https://cptac-data-portal.georgetown.edu/cptacPublic/) was utilized to analyze both mRNA and protein levels of RAP genes.

TABLE 1

TABLE 1. Patient charcterstics in TCGA-LIHC dataset.

DNA Methylation Analysis

DNA methylation of RAP genes in TCGA cancer dataset was estimated and visualized using MEXPRESS web server (https://mexpress.be) (Koch et al., 2015; Koch et al., 2019) and TCGA Wanderer (Díez-Villanueva et al., 2015). The MEXPRESS web server uses DNA methylation data of cancer and normal tissues from TCGA datasets, which were originally developed on the Illumina Human Methylation 450 BeadChip platform. The predesignated methylation probes for each gene were taken into consideration.

Survival Analysis

Kaplan-Meier survival analysis was performed using the tool available in the KM-plotter (Nagy et al., 2018). For Kaplan-Meier analysis, patients were distributed in high and low expression groups based on median expression value as a cut-off point for each gene. For survival analysis using univariate and multivariate Cox proportionate hazard model, RAP2A gene expression was taken as a continuous variable with multiple survival parameters for the TCGA-LIHC dataset, as recommended (Liu et al., 2018).

Correlation and Pathway Enrichment Analysis

Similarly, whole transcriptome correlations of RAP2A in the TCGA-LIHC study were downloaded from the cBioPortal website (https://www.cbioportal.org/) (Cerami et al., 2012; Gao et al., 2013). After applying a filter for a cutoff of FDR corrected p-value of 0.05 for Spearman’s r-value, 10,980 genes with Spearman’s correlation q value <0.05 were filtered and used for gene set enrichment analysis in GSEA software (Broad Institute, http://www.broad.mit.edu/gsea/). Hallmark gene set (version 7.1) (Subramanian et al., 2005) from predefined molecular signature database was used as a reference gene set for pathway enrichment (Liberzon et al., 2015).

Tumor Immunity Associations

Tumor immune estimation score (TIMER) webserver (https://cistrome.shinyapps.io/timer/), which utilizes the RNA sequencing data from TCGA for estimation of correlation between gene expression and level of immune cells, present in the tumor samples (Li et al., 2017). We utilized TIMER to calculate the association between RAP2A gene expression with infiltration of six immune cells including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells in TCGA-LIHC datasets. Default parameters were used in the TIMER database for the gene-specific analysis module.

Further, CIBERSORT (Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts) analysis data of TCGA-LIHC were extracted from a previously published study (Chen et al., 2018; Thorsson et al., 2018). This provided relative fractions of 22 different immune cells from a mixture of gene expression profiles (TCGA-LIHC study) and was used to correlate with RAP2A expression using Spearman’s correlation test. A total of 360 HCC samples were available with both gene expression data and CIBERSORT analysis estimated fractions of immune cells. Heatmap of the immune cell profiling data was generated along with hierarchical clustering using HemI (Deng et al., 2014). The default parameters of hierarchical clustering using the average linkage method and Pearson distance were used.

Statistical Analysis

Data analysis was performed using Graphpad (version 6) and Stata software (version 11). Mann-Whitney U-test was used for comparison among histological subtypes, molecular subtype and grades (***p < 0.001; **p < 0.01; *p < 0.05; ns, p > 0.05). Pearson correlation analysis was used to determine the association of DNA methylation of RAP2A to its expression in the TCGA-LIHC dataset. Kaplan-Meier survival analysis was performed using the log-rank test. A p-value < 0.05 was considered statistically significant.

Results

Expression Pattern of RAPs in HCC

Further, RNA sequencing data from TCGA-LIHC study was utilized to compare RAP gene expression in tumor tissues with both tumor adjacent normal tissues from the same dataset and with non-tumor associated normal hepatic tissue. RAP1A, RAP1B, RAP2A and RAP2B exhibited significant higher expression in tumors compared to other two groups (Figures 1A,C,E,G). Although RAP2C expression was higher in tumors compared to adjacent normal tissues, but both these groups exhibited lower expression of RAP2C compared to normal tissues from GTEx (Figure 1I). Furthermore, tumor adjacent normal tissues also exhibited higher expression of RAP1A and RAP2B, while no difference was observed for RAP1B and RAP2A. Comparison between 50 paired normal and tumor tissues from TCGA-LIHC also revealed that all RAP genes exhibit higher expression in tumor tissues compared to tumor adjacent normal tissues (Figures 1B,D,F,H,J). Among all RAPs, RAP2A displayed most robust upregulation in tumor tissues in TCGA dataset (Figure 1F). Further, we utilized multiomics data of hepatocellular carcinoma developed by CPTAC study, where both mRNA and proteomic data was available. The expression analysis in CPTAC data also suggested that expression of RAP genes differ between normal and tumor tissues both at the mRNA and protein level. In CPTAC data also, RAP2A exhibited most robust upregulation of mRNA and protein levels in tumor tissues compared to normal tissues, while expression of RAP2C was found to be reduced in tumor tissues compared to normal tissues (Supplementary Figure S1).

FIGURE 1

FIGURE 1. Expression of RAP genes in tumor tissues compared with adjacent normal tissues and normal tissues from GTEx study. GTEx, Genotype-Tissue Expression project; TANCT, tumor adjacent non-cancerous tissue. ****p < 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05; ns, p > 0.05.

We further performed receiver operating characteristic (ROC) curve analysis to determine potential of RAP gene expression in differentiating tumor tissues from normal liver tissues. Interestingly, among five RAP genes, RAP2A exhibited highest area under curve (AUC) of 0.8676 in TCGA-LIHC mRNA data (Figure 2A). Similarly, analysis of CPTAC mRNA data also suggested highest AUC of RAP2A (AUC: 0.9173, Figure 2B) compared to other RAP genes. Interestingly, analysis of AUC in CPTAC protein expression data revealed that RAP2C exhibit highest AUC of 0.8445 followed by RAP2A with AUC of 0.8172 (Figure 2C).

FIGURE 2

FIGURE 2. (A) ROC curve for the utility of RAP gene expression to differentiate between liver tumor tissues and normal tissue group in (A) TCGA mRNA data (B) CPTAC mRNA data, and (C) CPTAC protein expression data. For (A), normal tissue group consisted of tumor-adjacent normal tissues from TCGA study and normal tissues from non-disease controls from the GTEx study. (D) Comparison of RAP2A gene expression among normal tissue, tumor tissue, and metastatic tissues assessed through TNM webtool.

Furthermore, expression data of RAP gene family in hepatocellular carcinoma tissues and normal tissues from multiple other datasets was assessed through TNMplot web server. The analysis revealed that RAP1A, RAP1B, RAP2A, and RAP2B genes exhibit significantly higher expression in HCC tissues compared to normal tissues in comparison of both available paired (Supplementary Figure S2, left panel) and unpaired tissues (Supplementary Figure S2, right panel). However, RAP2C did not exhibit significant difference in expression in paired tissue analysis (Supplementary Figures S2J). Considering robust upregulation of RAP2A in tumors, and its established involvement in cell migration, we compared expression of RAP2A in metastatic tissues with both normal and primary tissues, which revealed highest expression of RAP2A in metastatic tissues compared to other two groups (Figure 2D).

Association of RAP Family Expression and Clinicopathological Features in HCC

We further assessed the association of RAP genes with clinicopathological including pathological age, gender, stage, tumor grade, blood AFP levels. Among all RAP genes, higher expression of RAP1B was associated with advanced-stage (Figure 3A). Higher expression of RAP2A and RAP2B, and low expression of RAP1A was associated with advanced grade (Figure 3B). High RAP2A expression was associated with younger age (<50 years, Supplementary Figure S3A) and female gender (Supplementary Figure S3B). Higher expression of RAP2A was also associated with increased AFP levels (Figure 4A). A history of alcohol consumption was associated with lower levels of RAP2A and RAP2C expression (Figure 4B).

FIGURE 3

FIGURE 3. Expression of RAP genes in tumor tissues compared between different stage (A) and grade (B). ****p < 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05; ns, p > 0.05.

FIGURE 4

FIGURE 4. Association of RAP expression in TCGA-KIRC dataset with (A) AFP levels, and (B) alcohol history. ****p < 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05; ns, p > 0.05.

Association of RAP Family Expression and Genetic and Epigenetic Alterations in HCC

To further determine whether the expression of RAP genes is associated with genetic alterations in HCC, we compared their expression in tumors with mutated or wild type TP53, CTNNB1, ALB, PCLO, and LRP1B. TP53 mutation was observed to be associated with higher expression of RAP1A, RAP1B, RAP2A, and RAP2B (Figure 5A). Further, CTNNB1 mutation was significantly associated with reduced levels of RAP1B, RAP2A, and RAP2B expression (Figure 5B). No RAP gene exhibited association with PCLO and ALB mutation status (Supplementary Figures S4A,B respectively), while higher expression of RAP1A was associated with LRPB1 mutant tumors (Supplementary Figure S5).

FIGURE 5

FIGURE 5. Association of RAP expression in TCGA-KIRC dataset with (A) TP53 mutation, and (B) CTNNB1 mutation. ****p < 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05; ns, p > 0.05.

To determine the potential role of copy number alterations and DNA methylation in the regulation of RAP2A expression in HCC, we utilized a TCGA-LIHC study where copy number variation, DNA methylation, and gene expression data were available. The DNA methylation data in TCGA study was developed on “Illumina HumanMethylation450 Beadchip” platform, where representative CpG sites from different regions of most genes are captured. Interestingly, RAP2A gene expression was reduced in tumor tissues and exhibited a negative correlation with DNA methylation at several sites within the RAP2A promoter regions and gene body (Figures 6A,B). A similar association was also observed for normal liver tissues (Supplementary Figure S6). We observed that in both normal and tumor tissues, DNA methylation at an intragenic region represented by cg03608515 was most negatively correlated with gene expression, suggesting this region, but not promoter region is the major regulatory site for the expression (Figure 6C). Furthermore, a comparison of 47 paired normal and tumor tissue also revealed significantly reduced methylation levels of cg03608515 in tumor tissues., these results strongly suggest the role of DNA methylation in aberrant expression of RAP2A in HCC. Additionally, the expression of RAP2A was also positively correlated with its copy number (r = 0.450, p < 0.001). Further, analysis of CNV data revealed frequent alterations in RAP2A copy number in HCC tissues was associated with its higher expression with copy number gain (Figure 6D, Kruskal-Wallis test, p < 0.0001).

FIGURE 6

FIGURE 6. Association of mRNA expression of RAP2A with its copy number variation and DNA methylation in TCGA-LIHC dataset. (A) Comparison of DNA methylation level of RAP2A between tumor tissues and normal tissues. (B) Correlation of RAP2A mRNA expression of RAP2A with its copy number variation and DNA methylation in tumor tissues. (C) Comparison of DNA methylation level of RAP2A at an intragenic site associated probe cg03608515. (D) Comparison of RAP2A gene expression among different copy number based groups in TCGA-LIHC dataset. ***p < 0.001; **p < 0.01; *p < 0.05. Insignificant associations (p > 0.05), are faded.

Prognostic Significance of RAP Genes in Hepatocellular Carcinoma

To determine the association of RAP gene family expression with patient prognosis, we utilized the TCGA-LIHC dataset where information for overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI) was available. We performed survival analysis by constructing a Kaplan-Meier plot for all RAP genes using median expression levels for allotting patients into high and low groups. We observed that higher expression of RAP2A was significantly associated with poor OS (HR = 1.72, CI = 1.21–2.45, p = 0.0023, Figure 7A) and DSS (HR = 1.9, CI = 1.2–2.99, p = 0.005, Figure 7B), while no significant association was observed with DFI and PFI (Figures 7C,D, respectively). In light of the high positive correlation of RAP2A with other RAP genes, we also assessed their association with patient survival (Supplementary Figure S7). Among other RAPs, higher expression of RAP1A and RAP1B was also associated with poor overall survival (Supplementary Figures S7A,B). We further performed univariate and multivariate survival analysis for RAP2A expression and other clinicopathological features, such as age, gender, stage, grade, alcohol intake history, radiotherapy status, and embolization status using Cox proportionate hazard model. Interestingly, higher RAP2A expression was also associated with poor OS, DSS, and PFI in both univariate and multivariate survival analysis (Tables 2, 3, respectively). This suggested that RAP2A expression is independently associated with poor outcome in HCC patients.

FIGURE 7

FIGURE 7. Kaplan-Meier survival analysis of RAP2A in TCGA-LIHC dataset, including (A) OS, overall survival (B) DSS, disease specific survival (C) DFI, disease free interval and (D) PFI, progression free interval. HR, hazard ratio; ****p < 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05; ns, p > 0.05.

TABLE 2

TABLE 2. Univariate analysis for association of RAP2A expression with patient prognosis in HCC.

TABLE 3

TABLE 3. Multivariate analysis for association of RAP2A expression with patient prognosis in HCC.

RAP2A Associated Cellular Pathways in Hepatocellular Carcinoma

To determine RAP2A associated cancer-related pathways, gene expression data of the TCGA-LIHC study was used. GSEA analysis revealed that RAP2A expression is positively correlated with cell cycle associated pathways such as mitotic spindle (Figure 8A), G2M checkpoint (Figure 8B), and E2F targets (Figure 8C) besides protein secretion (Figure 8D). Further, negatively correlated genes were enriched in metabolism associated pathways, such as oxidative phosphorylation (Figure 8E), xenobiotic metabolism (Figure 8F), fatty acid metabolism (Figure 8G), bile acid metabolism (Figure 8H), adipogenesis (Figure 8I), reactive oxygen species (Figure 8J) and others such as coagulation (Figure 8K), peroxisome (Figure 8L), interferon-alpha response (Figure 8M), DNA repair (Figure 8N) and Myc target genes (Figure 8O).

FIGURE 8

FIGURE 8. Gene set enrichment analysis of RAP2A correlated genes in TCGA-LIHC dataset. (A–D) depicts positively enriched pathways (E–L) depicts negatively enriched pathways with normalized enrichment score (NES), false discovery rate (FDR), and p-value depicted inside the respective pathway.

Association of RAP2A Expression With Tumor Immunity

Considering the previously described role of RAP genes in immune cell functions (Carvalho et al., 2019a), we analyzed the association of RAP2A expression with the level of immune cell infiltration. Using the TIMER tool, we determine tumor purity normalized spearman correlation of RAP2A expression with infiltration level of six different immune cells. This analysis revealed a positive correlation between RAP2A expression with B cells (r = 0.3, p = 1.37e-08), CD8+ T cells (r = 0.237, p = 9.06e-06), CD4+ T cells (r = 0.474, p = 1.16e-20), macrophages (r = 0.469, p = 4.56e-20), neutrophils (r = 0.374, p = 7.19e-13), and dendritic cells (r = 0.401, p = 1.36e-14) in HCC (Figure 9A). Furthermore, we utilized CIBERSORT analysis to determine the association of RAP2A gene expression with the relative abundance of 22 different types of immune cells in the TCGA-LIHC dataset (Figure 9B, Supplementary Table S2). Among immune cells, RAP2A expression was positively correlated to CD4 Memory Resting T cells, resting dendritic cells, neutrophils, M0 type macrophages, and naïve B cells, while it exhibited negative correlations to monocytes, activated NK cells, CD4 naïve T cells, CD8 T cells.

FIGURE 9

FIGURE 9. Association of RAP2A gene with tumor immunity in TCGA-LIHC dataset. (A) TIMER analysis showing the correlation of RAP2A expression with an abundance of six different immune cell types in TCGA-LIHC dataset. (B) CIBERSORT analysis showing relative fractions of 22 different immune cell types in HCC tissues (represented by rows) arranged in order of high RAP2A expression (top) to low RAP2A expression (bottom).

Discussion

HCC is one of the leading causes of cancer-related deaths worldwide. Significant advancement has been made in the treatment of this malignancy over the past decade, however, clinical response is highly heterogeneous. Further, treatment strategies have been highly adapted to be based on the progression of the disease at the time of diagnosis. Nevertheless, several molecular biomarkers have been determined with high prognostic value and future studies are required to determine novel molecular features as therapeutic targets and prognostic biomarkers. In the current study, we uncovered distinct genomic and epigenomic features of RAP family genes in HCC. Our study revealed that among five RAP genes, RAP2A expression is highly altered in HCC and is associated with multiple oncogenic features in HCC.

Little is known about the specific roles of RAP2A; in its active form RAP2A interacts with several effectors including MINK1, TNIK, and MAP4K4 and activates various signaling pathways involved in cytoskeletal rearrangements, cell migration, cell adhesion, and cell proliferation (Mittal and Linder, 2006). RAP2A interacts directly with upstream MAPK signaling element MAP4K4, and thus, increased RAP2A activity can enable downstream signaling (Machida et al., 2004). So far, the role of RAP2A in human malignancies remains controversial, with some suggesting it as a tumor suppressor gene while other studies refer to it as an oncogene. Upregulation of RAP2A has been observed in several human malignancies such as follicular thyroid cancer (Prabakaran et al., 2011), prostate cancer (Bigler et al., 2007), renal cancer (Wu et al., 2017), gastric cancer (Zhang J. et al., 2020) and bladder cancer (Wang et al., 2020).

In prostate cancer cells, RAP2A promotes androgen hypersensitivity and cell growth (Bigler et al., 2007). In lung cancer cells, ectopic expression of RAP2A enhances the migration and invasion of the cells (Wu et al., 2014). In bladder cancer cells, the expression of RAP2a was found significantly higher as compared to normal cells. The proliferation and invasion of cells were repressed by miR-3127 through directly targeting the 3′-UTR of RAP2A and associated with poor overall survival in bladder cancer patients (Wang et al., 2020). In gastric cancer, the role of RAP2A was also observed in drug resistance where expression of RAP2A increased the viability, migration, and metastasis of cells by suppressing apoptosis and DNA damage (Zhang J. et al., 2020). In renal cancer, overexpression of RAP2A enhances the protein levels of p-Akt and promotes migration and invasion of cells by increasing p-Akt expression (Wu et al., 2017). Contrary to these, RAP2A seems to play tumor suppressor functions in glioma as its downregulation is associated with glioma progression and its inhibition in the glioma cell line reduces migration and invasion (Wang et al., 2014). Results of the current study indicate that in hepatocellular carcinoma, RAP2A may act as an important oncogene and its mRNA expression is strongly associated with patient prognosis in HCC. Furthermore, other RAP genes also exhibit a strong positive correlation with RAP2A expression. This might be due to the conservation of regulatory sequences during evolution. We were further interested in whether RAP genes share common features for association with molecular characteristics in HCC.

It was recently demonstrated that RAP2A expression is regulated by p53 and RAP2A mediated cell migration and invasive properties are driven by downstream activation of the matrix metalloproteinases (MMP) MMP2 and MMP9 via phosphorylation of AKT (Wu et al., 2015). Consistent with this, we observed higher expression of multiple RAP genes, including RAP2A in p53 mutant HCC. Further, we also observed that expression of RAP1A, RAP1B, RAP2A, and RAP2B were reduced in HCC tissues which harbor a mutation in CTNNB1, the gene encoding for beta-catenin protein. This is contrary with the previous report where RAP1B has shown to activate Wnt/beta-catenin signaling in esophageal squamous cell carcinoma (Jia et al., 2017). Further, RAPGEF2, a guanine nucleotide exchange factor for RAP1, was shown to regulate adherence junction (AJ) formation in radial glial cells through ERK-mediated upregulation of β-catenin (Farag et al., 2017). While CTNNB1 mutations in HCC are associated with higher activity of Wnt-beta catenin signaling (Tornesello et al., 2013), its association with RAP signaling appears to be negatively related in this case. Therefore, our results suggested potential crosstalk of Wnt-beta catenin signaling in RAP signaling in HCC tissues.

In light of its aberrant overexpression in HCC, we explored whether the expression of RAP2A is driven by copy number alteration and DNA methylation in HCC. Our results collectively demonstrated that the RAP2A harbors alterations in both of the abovementioned features. Our results highlighted a specific intragenic region in the RAP2A where DNA methylation was highly reduced in tumor tissues compared to normal liver tissues. Further, DNA methylation at this region is negatively correlated to RAP2A gene expression in both tumor and normal tissues. DNA methylation of RAP2A has not been previously studied in cancer, therefore, epigenetic regulation of RAP signaling requires detailed exploration.

While our study is based on mRNA expression, a recent study by, Zheng et al. has also demonstrated that RAP2A protein expression is associated with oncogenic features in HCC (Zheng et al., 2017). Therefore, our findings further provide a detailed understanding of the role of all five members of this gene family involvement in HCC. Among all five RAPs, RAP2A expression exhibited a strong ability to differentiate tumor tissues from normal tissues. Further, its higher expression also exhibited association with higher tumor grade, metastasis, increased AFP levels, and poor patient prognosis. Furthermore, our multivariate survival analysis including major clinical and pathological features revealed that the RAP2A expression is independently associated with poor overall survival, disease-specific survival, and progression-free interval in HCC.

Pathway analysis revealed strong associations of RAP2A expression in HCC with several HCC relevant pathways, including cell cycle-related pathways and metabolic pathways. Interestingly, RAP1A expression has previously been shown to be regulated during the cell cycle (Cruise et al., 1997). The causal relationship between RAP2A expression and these pathways requires further validation. We also analyzed the immunological association of RAP2A expression in HCC, which revealed that its expression is highly associated with the immune composition of HCC tumors. While, the role of RAP2A has been previously demonstrated in the regulation of lipopolysaccharide induced innate cell functions (Carvalho et al., 2019a; Carvalho et al., 2019b), detailed role of RAP2A in the modulation of tumor immunity remains to be studied in detail. Conclusively, the current study provides detailed molecular and clinical features associated with the expression of RAP genes in HCC, however, some of these associations require further exploration for the causal relationships. Further, these results support the potential of RAP2A as a therapeutic target and prognostic biomarker in this malignancy.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

Conceptualization, SC and MA; supervision and resources, SC, RY, and AC; experimentation, MA, SK, and JS; Manuscript writing, MA and LK; manuscript editing, SC, RY, and AC. All authors have approved the final version of the submitted manuscript.

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.

Acknowledgments

MA acknowledges financial support as a fellowship from the Council of Scientific and Industrial Research, Government of India. SK acknowledges financial support as a fellowship from the Department of Health Research, Government of India. JS acknowledges financial support as a fellowship from the Department of Biotechnology, Government of India.

Supplementary Material

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

References

Bartha, Á., and Győrffy, B. (2021). TNMplot.com: A Web Tool for the Comparison of Gene Expression in Normal, Tumor and Metastatic Tissues. Int. J. Mol. Sci., 22. doi:10.3390/ijms22052622

CrossRef Full Text | Google Scholar

Bigler, D., Gioeli, D., Conaway, M. R., Weber, M. J., and Theodorescu, D. (2007). Rap2 Regulates Androgen Sensitivity in Human Prostate Cancer Cells. Prostate 67, 1590–1599. doi:10.1002/pros.20644

PubMed Abstract | CrossRef Full Text | Google Scholar

Bokoch, G. M. (1993). Biology of the RAP Proteins, Members of the Ras Superfamily of GTP-Binding Proteins. Biochem. J. 289 (Pt 1), 17–24. doi:10.1042/bj2890017

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: A Cancer J. Clinicians 68, 394–424. doi:10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruix, J., and Sherman, M.American Association for the Study of Liver Diseases (2011). Management of Hepatocellular Carcinoma: an Update. Hepatology 53, 1020–1022. doi:10.1002/hep.24199

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, B. C., Oliveira, L. C., Rocha, C. D., Fernandes, H. B., Oliveira, I. M., Leãõ, F. B., et al. (2019a). Data in Support of RAP2A GTPase Expression, Activation and Effects in LPS-Mediated Innate Immune Response and NF-κB Activation. Data in Brief 24, 103965. doi:10.1016/j.dib.2019.103965

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, B. C., Oliveira, L. C., Rocha, C. D., Fernandes, H. B., Oliveira, I. M., Leão, F. B., et al. (2019b). Both knock-down and Overexpression of RAP2A Small GTPase in Macrophages Result in Impairment of NF-κB Activity and Inflammatory Gene Expression. Mol. Immunol. 109, 27–37. doi:10.1016/j.molimm.2019.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data: Figure 1. Cancer Discov. 2, 401–404. doi:10.1158/2159-8290.CD-12-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

Che, Y.-L., Luo, S.-J., Li, G., Cheng, M., Gao, Y.-M., Li, X.-M., et al. (2015). The C3G/RAP1 Pathway Promotes Secretion of MMP-2 and MMP-9 and Is Involved in Serous Ovarian Cancer Metastasis. Cancer Lett. 359, 241–249. doi:10.1016/j.canlet.2015.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, B., Khodadoust, M. S., Liu, C. L., Newman, A. M., and Alizadeh, A. A. (2018). Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol. Biol. Clifton NJ 1711, 243–259. doi:10.1007/978-1-4939-7493-1_12

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruise, J. L., Rafferty, M. P., and Riehle, M. M. (1997). Cell-cycle Regulated Expression of RAP1 in Regenerating Liver. Biochem. Biophysical Res. Commun. 230, 578–581. doi:10.1006/bbrc.1996.6003

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, W., Wang, Y., Liu, Z., Cheng, H., and Xue, Y. (2014). HemI: A Toolkit for Illustrating Heatmaps. PLoS One 9, e111988. doi:10.1371/journal.pone.0111988

PubMed Abstract | CrossRef Full Text | Google Scholar

Di, J., Huang, H., Qu, D., Tang, J., Cao, W., Lu, Z., et al. (2015a). RAP2B Promotes Proliferation, Migration and Invasion of Human Breast Cancer through Calcium-Related ERK1/2 Signaling Pathway. Sci. Rep. 5, 12363. doi:10.1038/srep12363

PubMed Abstract | CrossRef Full Text | Google Scholar

Di, J., Huang, H., Wang, Y., Qu, D., Tang, J., Cheng, Q., et al. (2015b). p53 Target Gene RAP2B Regulates the Cytoskeleton and Inhibits Cell Spreading. J. Cancer Res. Clin. Oncol. 141, 1791–1798. doi:10.1007/s00432-015-1948-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Díez-Villanueva, A., Mallona, I., and Peinado, M. A. (2015). Wanderer, an Interactive Viewer to Explore DNA Methylation and Gene Expression Data in Human Cancer. Epigenet. Chromatin 8, 22. doi:10.1186/s13072-015-0014-8

CrossRef Full Text | Google Scholar

Ehrhardt, A., Ehrhardt, G., Guo, X., and Schrader, J. (2002). Ras and Relatives-Job Sharing and Networking Keep an Old Family Together. Exp. Hematol. 30, 1089–1106. doi:10.1016/s0301-472x(02)00904-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Farag, M. I., Yoshikawa, Y., Maeta, K., and Kataoka, T. (2017). RAPGEF2, a Guanine Nucleotide Exchange Factor for RAP1 Small GTPases, Plays a Crucial Role in Adherence junction (AJ) Formation in Radial Glial Cells through ERK-Mediated Upregulation of the AJ-Constituent Protein Expression. Biochem. Biophys. Res. Commun. 493, 139–145. doi:10.1016/j.bbrc.2017.09.062

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, G., Liu, Y., Yuan, J., Zheng, H., Shi, T., Lei, W., et al. (2009). [Identification and Functional Analysis of A Novel Candidate Oncogene RAP2B in Lung Cancer.]. Zhongguo Fei Ai Za Zhi 12, 273–276. doi:10.3779/j.issn.1009-3419.2009.04.03

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signaling 6, pl1. doi:10.1126/scisignal.2004088

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldman, M. J., Craft, B., Hastie, M., Repečka, K., McDade, F., Kamath, A., et al. (2020). Visualizing and Interpreting Cancer Genomics Data via the Xena Platform. Nat. Biotechnol. 38, 675–678. doi:10.1038/s41587-020-0546-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, Z., Yang, Y., Dengyan, Z., Chunyang, Z., Donglei, L., Kai, W., et al. (2017). RAP1B, a DVL2 Binding Protein, Activates Wnt/beta-Catenin Signaling in Esophageal Squamous Cell Carcinoma. Gene 611, 15–20. doi:10.1016/j.gene.2017.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, A., De Meyer, T., Jeschke, J., and Van Criekinge, W. (2015). MEXPRESS: Visualizing Expression, DNA Methylation and Clinical TCGA Data. BMC Genomics 16, 636. doi:10.1186/s12864-015-1847-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, A., Jeschke, J., Van Criekinge, W., van Engeland, M., and De Meyer, T. (2019). MEXPRESS Update 2019. Nucleic Acids Res. 47, W561–W565. doi:10.1093/nar/gkz445

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77, e108–e110. doi:10.1158/0008-5472.CAN-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Li, S., and Huang, L. (2018). Knockdown of RAP2B, a Ras Superfamily Protein, Inhibits Proliferation, Migration, and Invasion in Cervical Cancer Cells via Regulating the ERK1/2 Signaling Pathway. Oncol. Res. 26, 123–130. doi:10.3727/096504017X14912172235777

PubMed Abstract | CrossRef Full Text | Google Scholar

Liberzon, A., Birger, C., Thorvaldsdóttir, H., Ghandi, M., Mesirov, J. P., and Tamayo, P. (2015). The Molecular Signatures Database Hallmark Gene Set Collection. Cel. Syst. 1, 417–425. doi:10.1016/j.cels.2015.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Lichtenberg, T., Hoadley, K. A., Poisson, L. M., Lazar, A. J., Cherniack, A. D., et al. (2018). An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell 173, 400–e11. doi:10.1016/j.cell.2018.02.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, L., Wang, J., Wu, Y., Wan, P., and Yang, G. (2016). RAP1A Promotes Ovarian Cancer Metastasis via Activation of ERK/p38 and Notch Signaling. Cancer Med. 5, 3544–3554. doi:10.1002/cam4.946

PubMed Abstract | CrossRef Full Text | Google Scholar

Machida, N., Umikawa, M., Takei, K., Sakima, N., Myagmar, B.-E., Taira, K., et al. (2004). Mitogen-activated Protein Kinase Kinase Kinase Kinase 4 as a Putative Effector of RAP2 to Activate the C-Jun N-Terminal Kinase. J. Biol. Chem. 279, 15711–15714. doi:10.1074/jbc.C300542200

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, Z., Qiu, Y., Lin, K. C., Kumar, A., Placone, J. K., Fang, C., et al. (2018). RAP2 Mediates Mechanoresponses of the Hippo Pathway. Nature 560, 655–660. doi:10.1038/s41586-018-0444-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal, V., and Linder, M. E. (2006). Biochemical Characterization of RGS14: RGS14 Activity towards G-Protein α Subunits Is Independent of its Binding to RAP2A. Biochem. J. 394, 309–315. doi:10.1042/BJ20051086

PubMed Abstract | CrossRef Full Text | Google Scholar

Mo, S.-J., Hou, X., Hao, X.-Y., Cai, J.-P., Liu, X., Chen, W., et al. (2018). EYA4 Inhibits Hepatocellular Carcinoma Growth and Invasion by Suppressing NF-κB-dependent RAP1 Transactivation. Cancer Commun. 38, 9. doi:10.1186/s40880-018-0276-1

CrossRef Full Text | Google Scholar

Nagy, Á., Lánczky, A., Menyhárt, O., and Győrffy, B. (2018). Validation of miRNA Prognostic Power in Hepatocellular Carcinoma Using Expression Data of Independent Datasets. Sci. Rep. 8, 9227. doi:10.1038/s41598-018-27521-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Y.-G., Zhang, Z.-Q., Chen, Y.-b., and Huang, J.-A. (2016). RAP2b Promotes Proliferation, Migration, and Invasion of Lung Cancer Cells. J. Receptors Signal Transduction 36, 459–464. doi:10.3109/10799893.2015.1122044

CrossRef Full Text | Google Scholar

Prabakaran, I., Grau, J. R., Lewis, R., Fraker, D. L., and Guvakova, M. A. (2011). RAP2AIs Upregulated in Invasive Cells Dissected from Follicular Thyroid Cancer. J. Thyroid Res. 2011, 1–6. doi:10.4061/2011/979840

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, D., Huang, H., Di, J., Gao, K., Lu, Z., and Zheng, J. (2016). Structure, Functional Regulation and Signaling Properties of RAP2B. Oncol. Lett. 11, 2339–2346. doi:10.3892/ol.2016.4261

PubMed Abstract | CrossRef Full Text | Google Scholar

Sequera, C., Manzano, S., Guerrero, C., and Porras, A. (2018). How RAP and its GEFs Control Liver Physiology and Cancer Development. C3G Alterations in Human Hepatocarcinoma. Hepatic Oncol. 5, HEP05. doi:10.2217/hep-2017-0026

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, Y., Ding, S., Chen, K., Chen, J., Wang, S., Zou, C., et al. (2014). Functional Analysis of miR-101-3p and RAP1B Involved in Hepatitis B Virus-Related Hepatocellular Carcinoma Pathogenesis. Biochem. Cel. Biol. 92, 152–162. doi:10.1139/bcb-2013-0128

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Z., Peng, H., Chen, J., Liu, Y., Yan, S., Yu, G., et al. (2018). RAP1B Enhances the Invasion and Migration of Hepatocellular Carcinoma Cells by Up-Regulating Twist 1. Exp. Cel. Res. 367, 56–64. doi:10.1016/j.yexcr.2018.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorsson, V., Gibbs, D. L., Brown, S. D., Wolf, D., Bortone, D. S., Ou Yang, T. H., et al. (2018). The Immune Landscape of Cancer. Immunity 48, 812–e14. doi:10.1016/j.immuni.2018.03.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Tornesello, M. L., Buonaguro, L., Tatangelo, F., Botti, G., Izzo, F., and Buonaguro, F. M. (2013). Mutations in TP53, CTNNB1 and PIK3CA Genes in Hepatocellular Carcinoma Associated with Hepatitis B and Hepatitis C Virus Infections. Genomics 102, 74–83. doi:10.1016/j.ygeno.2013.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Villanueva, A. (2019). Hepatocellular Carcinoma. N. Engl. J. Med. 380, 1450–1462. doi:10.1056/NEJMra1713263

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhan, W., Xie, S., Hu, J., Shi, Q., Zhou, X., et al. (2014). Over-expression of RAP2A Inhibits Glioma Migration and Invasion by Down-Regulating P-AKT. Cell Biol. Int. 38, 326–334. doi:10.1002/cbin.10213

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhu, B., Wang, S., Wu, Y., Zhan, W., Xie, S., et al. (2017). Regulation of Glioma Migration and Invasion via Modification of RAP2A Activity by the Ubiquitin Ligase Nedd4-1. Oncol. Rep. 37, 2565–2574. doi:10.3892/or.2017.5572

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Meng, R., and Hu, Q.-M. (2020). LINC00319-Mediated miR-3127 Repression Enhances Bladder Cancer Progression through Upregulation of RAP2A. Front. Genet. 11, 180. doi:10.3389/fgene.2020.00180

PubMed Abstract | CrossRef Full Text | Google Scholar

Wheeler, D. A., and Roberts, L. R. (2017). Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma. Cell 169, 1327–e23. doi:10.1016/j.cell.2017.05.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Sang, M., Cao, W., Zheng, J., and Pei, D. (2014). [Identification Analysis of Eukaryotic Expression Plasmid RAP2A and its Effect on the Migration of Lung Cancer Cells]. Zhongguo Fei Ai Za Zhi 17, 643–648. doi:10.3779/j.issn.1009-3419.2014.09.01

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J.-X., Zhang, D.-G., Zheng, J.-N., and Pei, D.-S. (2015). RAP2A Is a Novel Target Gene of P53 and Regulates Cancer Cell Migration and Invasion. Cell Signal. 27, 1198–1207. doi:10.1016/j.cellsig.2015.02.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J.-X., Du, W.-Q., Wang, X.-C., Wei, L.-L., Huo, F.-C., Pan, Y.-J., et al. (2017). RAP2A Serves as a Potential Prognostic Indicator of Renal Cell Carcinoma and Promotes its Migration and Invasion through Up-Regulating P-Akt. Sci. Rep. 7, 6623. doi:10.1038/s41598-017-06162-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, X., Liu, H., Wang, M., Ding, F., Xiao, H., Hu, F., et al. (2015). miR-342-3p Targets RAP2B to Suppress Proliferation and Invasion of Non-small Cell Lung Cancer Cells. Tumor Biol. 36, 5031–5038. doi:10.1007/s13277-015-3154-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Duan, H.-b., and Yang, Y.-s. (2017). Knockdown of RAP2B Inhibits the Proliferation and Invasion in Hepatocellular Carcinoma Cells. Oncol. Res. 25, 19–27. doi:10.3727/096504016X14685034103914

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Wei, Y., Min, J., Wang, Y., Yin, L., Cao, G., et al. (2020). Knockdown of RAP2A Gene Expression Suppresses Cisplatin Resistance in Gastric Cancer Cells. Oncol. Lett. 19, 350–358. doi:10.3892/ol.2019.11086

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Wu, J., Yang, Y., Xia, D., Li, J., Quan, H., et al. (2020). Donor Polymorphisms of RAP1A Rs494453 Contribute to a Higher Risk of Hepatocellular Carcinoma Recurrence Following Liver Transplantation. J. Cancer 11, 3082–3088. doi:10.7150/jca.39712

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, X., Zhao, W., Ji, P., Zhang, K., Jin, J., Feng, M., et al. (2017). High Expression of RAP2A Is Associated with Poor Prognosis of Patients with Hepatocellular Carcinoma. Int. J. Clin. Exp. Pathol. 10, 9607–9613.

PubMed Abstract | Google Scholar