Skip to main content

ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 07 July 2021
Sec. Molecular and Cellular Oncology
This article is part of the Research Topic Tumor Micro-environment and Drug Resistance View all 33 articles

Transcriptomic Analysis Identified ARHGAP Family as a Novel Biomarker Associated With Tumor-Promoting Immune Infiltration and Nanomechanical Characteristics in Bladder Cancer

  • 1Department of Urology, Huashan Hospital, Fudan University, Shanghai, China
  • 2Fudan Institute of Urology, Huashan Hospital, Fudan University, Shanghai, China
  • 3National Clinical Research Center for Aging and Medicine, Fudan University, Shanghai, China
  • 4School of Mechatronics Engineering and Automation, Shanghai University, Shanghai, China

Bladder cancer (BCa) is a common lethal urinary malignancy worldwide. The role of ARHGAP family genes in BCa and its association with immuno-microenvironment remain largely unknown. ARHGAP family expression and immune infiltration in BCa were analyzed by bioinformatics analysis. Then, we investigated cell proliferation, invasion, and migration in vivo and in vitro of the ARHGAP family. Furthermore, atomic force microscopy (AFM) was employed in measuring cellular mechanical properties of BCa cells. The results demonstrated that ARHGAP family genes correlate with a tumor-promoting microenvironment with a lower Th1/Th2 cell ratio, higher DC cell infiltration, higher Treg cell infiltration, and T-cell exhaustion phenotype. Silencing ARHGAP5, ARHGAP17, and ARHGAP24 suppressed BCa cell proliferation, migration, and metastasis. Knocking down of ARHGAPs in T24 cells caused a relatively higher Young’s modulus and lower adhesive force and cell height. Taken together, ARHGAP family genes promote BCa progressing through establishing a tumor-promoting microenvironment and promoting cancer progression.

Introduction

Bladder cancer (BCa) is the ninth leading diagnosed malignancy, which causes 17,980 mortality worldwide in 2020 (Siegel et al., 2020). The majority of BCa cases are classified as non-muscle invasive BCa representing the curable pathological type, while the other 25% are muscle-invasive BCa characterized by rapid progression and high recurrence rate, and eventually progress into metastatic disease (Smith et al., 2012; Nekolla et al., 2019). In recent years, various attempts to detect BCa at an early stage or explore the potential mechanism of lethal BCa have been made (Chen et al., 2020; Xie et al., 2021), so treatment strategies based on underlying molecular mechanisms in the metastasis and progression of BCa is in urgent need (Duggan et al., 2000).

Cancer cells are endowed with unique biological capabilities allowing them for constitutive survival, proliferation, and invasion ability during carcinogenesis (Chaffer and Weinberg, 2011). Adaptive immune response, a key regulator in regulating oncogenesis, has been well studied in the past decades. Several cancer-related hallmarks have been found responsible for the interaction of tumor and immune cells. Some of them can impair the differentiation and maturation of various immune cell subpopulations, which eventually leads to weakened antitumor functions in the tumor microenvironment, helping cancer cells to escape from immune surveillance (Chen and Mellman, 2017). In several tumor types, a better prognosis is mainly associated with infiltration of CD8+ T cells, type 1 T helper (Th1) cells, natural killer (NK) cells, and M1 macrophages. Conversely, poor prognosis usually correlates with high infiltration of Treg cells, Th2 cells, MDSCs, and M2 macrophages, and neutrophils (Fridman et al., 2012; Becht et al., 2016). These findings implied a distinctive mechanism of BCa-related immunological regulation, which could bring a radical revolution in BCa treatment other than the classical cisplatin-based chemotherapy.

The Rho GTPase family, a part of the Ras superfamily, consists of some highly conserved genes in regulating biological processes like cytoskeleton organization, vesicle trafficking, cell cycle, cell polarity, cell invasion, and cell migration (Haraguchi et al., 2019). Rho GTPases catalyze the conversion between active form and inactive form of Ras superfamily, henceforth, suppressing Rho GTPase downstream cellular biological processes. GTPase-activating proteins (RhoGAPs), which negatively regulates Rho GTPases are known as ARHGAP family (Post et al., 2013). While the roles of several ARHGAP family members have been identified in several types of cancer (Gen et al., 2009; Wang et al., 2014a, b; Luo et al., 2016; Hashimoto et al., 2018), the function of ARHGAP family in BCa, especially in immunological microenvironment of BCa, has not been elucidated.

Through comprehensive assessment of ARHGAP genes, we found a lower ARHGAP family expression in BCa, of which ARHGAP5, ARHGAP8, ARHGAP11A, ARHGAP17, ARHGAP24, ARHGAP37 (STARD13), and ARHGAP38 (STARD8) were highly related with prognosis. Axon guidance, focal adhesion, and leukocyte transendothelial migration were associated with the ARHGAP family genes, while focal adhesion, chemokine signaling pathway, and T-cell receptor signaling pathway were enriched in ARHGAP coexpression genes. We demonstrated that the ARHGAP family genes correlate with a tumor-promoting microenvironment. Silencing ARHGAP5, ARHGAP17, and ARHGAP24 suppressed BC cell proliferation, migration, and metastasis. Meanwhile, we validated that knocking down of ARHGAPs in T24 cells caused a relatively higher Young’s modulus and lower adhesive force and cell height, suggesting that ARHGAP could maintain a malignant cellular mechanical property of BCa. Our work offers novel insights into a potential mechanism of ARHGAP family genes’ correlation with immune microenvironment in BCa.

Materials and Methods

cBioPortal Database Analysis

The cBioPortal database1 offers visualized analysis of cancer genomics based on available high-throughput sequencing data (Cerami et al., 2012). To analyze the mutation of the ARHGAP family genes, a network of the mutation profile of ARHGAP family was generated based on data of 413 samples in The Cancer Genome Atlas database (TCGA, Bladder Cancer, Cell, 2017) (Robertson et al., 2017) in cBioPortal.

Oncomine Database Analysis

The Oncomine database (Rhodes et al., 2004)2 compiled 86,733 cancer samples within 715 gene expression into construction of a comprehensive data-mining database. We assessed the ARHGAP family expression in various kinds of tumor types.

TIMER Database Analysis

The TIMER3 database contributes to analyzing immune cell infiltration level in multiple cancers through pathological validation and improved statistical methodology to evaluate tumor immune infiltration (Li et al., 2017). We employed this database to assess differences in expression level of several candidates in the ARHGAP family in BCa to explore the association between ARHGAP family expression and infiltration level of particular immune cell subsets including B cells, CD4+ T and CD8+ T cells, neutrophils, macrophages, and dendritic cells. The Kaplan–Meier curve is also applied to analyze patient survival of the differential expression of the ARHGAP family gene-associated immune cell infiltration. The correlation between candidate ARHGAP family genes expression and gene markers of immune cells through relevant modules were analyzed.

GEPIA Database Analysis

The GEPIA database (Tang et al., 2017)4 facilitates the standardized analysis of RNA-seq data from 9,736 tumors and 8,587 normal samples originating from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) data sets. We therefore assessed the link between the ARHGAP family gene expression and patients’ prognoses in BCa.

Gene Set Enrichment Analysis

To further investigate the biological concepts associated with the ARHGAP family genes, we performed Gene Set Enrichment Analysis (GSEA) on the RNA-seq data downloaded from TCGA database. The mRNA expression of BCa was divided into high and low groups according to expression level. The carcinoma-related pathways gene sets (h.all.v6.0.symbols.gmt) were downloaded from the Molecular Signatures Database-MsigDB5. Enrichment analysis was performed by randomly repeating 1,000 times. Gene sets with a p-value < 0.05, normalized enrichment score (NES) > 1 or < -1, and false discovery rate (FDR) < 0.25 were considered as significantly enriched.

Correlation and Protein Interaction Analysis

We both performed correlation and protein–protein interaction (PPI) analysis on the ARHGAP family genes. Herein, we generated correlation analysis by using BCa’s mRNA expression of the ARHGAP family genes downloaded from TCGA database. The results of correlation and relative p-value were constructed by using “corrplot” package in R (version 4.0.2). Then we applied GeneMANIA6 to analyze the ARHGAP family gene interaction and showed the PPI network including coexpression, pathway, predicted, colocalization, and genetic interactions. We also explored 50 most potentially relevant proteins by STRING7 database to serve the PPI network and showed the MCODE components identified by Metascape8.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Functional Enrichment Analysis

To explore the biological functions and novel pathways of the ARHGAP family genes, we applied Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes database (KEGG) enrichment. All GO and KEGG functional enrichment analysis were carried out by “clusterProfiler” and “enrichplot” package in R (version 4.0.2). Results of GO enrichment analysis including “Biological Process,” “Cellular Component,” and “Molecular Function” were visualized as a dot plot. Several most representative GO and KEGG analysis were also shown by R (version 4.0.2). A p–value < 0.05 was set as the cutoff criterion and considered as statistically significant.

Tumor-Infiltrating Immune Cell Portraying

Two methods, xCell (Aran et al., 2017) and TIMER (Tumor Immune Estimation Resource,9 were used to analyze the tissue-infiltrating cell-type abundance from bulk RNA-seq data. Pre-calculated TCGA data by xCell was downloaded from the xCell web tool10. Cell type enrichment scores of TCGA-BLCA samples were extracted as reflection of certain cell type abundance. Metrics including immune score, stroma score, and microenvironment score were defined according to the source code of xCell on GitHub. Among 64 cell types xCell could output, cell types including CD8+ T cells, NK cells, dendritic cells (DCs), CD4+ T cells, regulatory T cells (Tregs), macrophage M1, macrophage M2, type 2 T-helper cells (Th2 cells), and type 1 T-helper cells (Th1 cells) were chosen to perform further analysis. R package “ggpubr” was used to visualize enrichment scores between groups. Samples with certain gene expression level (FPKM normalized expression level) higher than the median of that of all selected samples were defined as highly expressed, while the remaining samples as lowly expressed. Wilcoxon signed rank test was used to compare means between two groups. “Gene module” of TIMER was used to explore the correlation between gene expression and abundance of immune infiltrates, with BLCA (bladder urothelial carcinoma) selected as cancer types.

Patient Samples

The study protocol was approved by the Ethics Committee of Huashan Hospital (Shanghai, China; approval no. KY2011-009) and conducted in accordance with the tenets of the Declaration of Helsinki. All patients consented to the use of resected tissues for research purposes.

A total of 90 pairs of BC tissues and adjacent tissues were collected for tissue microarray (TMA) construction from BC patients after surgical treatment in Huashan Hospital, Fudan University, between January 2007 and January 2013 including a 5-year follow-up. Two experienced pathologists confirm the pathological diagnoses of BC according to the 7th edition of the TNM classification of the Union for International Cancer Control (Wittekind et al., 2019).

Another seven pairs of high-grade and seven pairs of low-grade bladder cancer with normal samples were obtained from Huashan Urology Tissue Bank under an approval from the Ethics Committee of Huashan Hospital. Samples were harvested and immediately snap frozen in liquid nitrogen. Meanwhile, 20 fresh bladder cancer tumor tissues were obtained for flow cytometry.

Immunohistochemical Staining

Single immunohistochemical staining was performed on TMA slides and 20 FFPE tissues from fresh resected tumors. The immunohistochemical staining procedure is illustrated as below. The slides were deparaffinized and rehydrated in advance with dimethylbenzene and ethanol. Heated sodium citrate buffer (0.01 M, pH = 6) was then applied for the slides for antigen repair. The slides were incubated in normal equine serum for 1 h at 37°C for blocking non-specific binding. The slides were washed and incubated with CD8 (Abcam, ab17147), CD4 (Abcam, ab67001), FOXP3 (Abcam, ab22510), ARHGAP5 (Abcam, ab32328), ARHGAP17 (Abcam, ab229221), and ARHGAP24 (Abcam, ab203874) overnight at 4°C in a wet chamber. The slides were stained with horseradish peroxidase-conjugated secondary antibody for 1 h at room temperature and then developed with DAB and hematoxylin. For the quantification of immunohistochemical staining, three randomized high-power fields of each sample were quantitated for the positive-staining cells, and the mean value was adopted. We quantitated FOXP3+ cells in each whole section since its relatively low infiltration in tumor. The 90 samples of the TMA slides were evenly divided into the ARHGAP family high/low groups according to the median value of ARHGAP5+, ARHGAP17+, and ARHGAP24+ cells. Then the three cutoff values were applied for the 20 fresh tumor tissues to define the ARHGAP family high/low groups.

Cell Culture and Infection

Bladder transitional cell carcinoma cell lines (RT4, UM-UC-3, T24, 5,637, and J82), immortalized uroepithelium cell line (SV-HUC-1), and human embryonic kidney cells (HEK-293) were received from Shanghai Yuanye Bio-Technology Co. (Shanghai, China). All cells were maintained at 37°C with 5% CO2 in Dulbecco’s modified Eagle’s medium (DMEM; Gibco, New York, NY, United States) supplemented with 10% fetal bovine serum (FBS; Gibco, New York, NY, United States). To knock down ARHGAP5, 17, 24, the recombinant plasmid vector or control was purchased from Vigene (Jinan, China). HEK-293T cells were seeded into a 10-cm culture dish 24 h before transfection with PSPAX2 and PMD2G. The medium containing virus particles was collected at 72 h post-transfection. After filtration (0.45-μm filter), the medium was added to T24 and UM-UC-3 for viral infection with fresh medium replenished 24 h later. After 48 h postinfection, selection was done with 1.0 μg/ml of puromycin (Sangon, China).

Cell Proliferation Assay

For colony formation assay, cells were seeded into six-well plates at a density of 600 cells per well. After 7–9 days of incubation, cells were fixed with 4% (w/v) paraformaldehyde (PFA) and stained with crystal violet solution.

5-Ethynyl-2′-Deoxyuridine Assay

5-Ethynyl-2′-deoxyuridine (EdU) assay kit (RiboBio, Guangzhou, China) was used to quantitatively investigate cell proliferation through DNA synthesis. A Nikon microscope (Nikon, Tokyo, Japan) was used to measure DNA synthesis reflected by red signals.

Atomic Force Microscopy

The mechanical properties of cells were quantified using a commercial atomic force microscopy (AFM) instrument (BioScope Resolve; Bruker Corporation, Billerica, MA, United States). For characterization, dishes containing cultured primary cells were placed on a stage equipped with a vacuum pump with low noise and drum effects. A PeakForce Quantitative Nanomechanical Mapping—Live Cell probe (Bruker Corporation), with a tip length of 17 μm, tip radius of 65 nm, tip half angle of 18°, and spring constant of 0.076 N/m, was applied to probe the cell surface of the contact model. For characterization, a constant loading force of 1 nN was applied. Deflection images (32 × 32 pixels) were acquired at a ramp rate of 10 Hz and a ramp area of 5 μm × 5 μm. Automatic gain control was used to improve the feedback for surface tracking. YM and AF signal channels were used to map the cell force. All experiments were completed within 2 h to ensure cellular health, and the nanomechanical properties were determined from the AFM force map using Nanoscope Analysis software v1.80 (Bruker Nano Surfaces, Goleta, CA, United States).

Transwell Assay

Cell migration was analyzed using Transwell chambers (Corning, New York, NY, United States) in accordance with the manufacturer’s protocol. After incubation for 24 h, the cells on the upper surfaces of the Transwell chambers were removed, and the cells located on the lower surfaces were fixed with 4% PFA, followed with crystal violet staining. The stained cells were photographed and counted in five randomly selected fields.

Western Blot Analysis

Lysates from cells and tumor tissues were prepared to determine protein levels using the Bradford assay (Bio-Rad). Proteins were separated by 10% SDS-PAGE and transferred to polyvinylidene difluoride transfer membranes. The blots were blocked with freshly prepared 5% non-fat milk in PBST for 2 h at room temperature. Then the blots were incubated at 4°C overnight with primary antibodies. After washing with PBST, the blots were incubated with horseradish peroxidase-conjugated (HRP-conjugated) donkey anti-rabbit IgG or sheep anti-mouse IgG (Invitrogen, Shanghai, China) at room temperature for 2 h. ECL substrate (CLiNX, Shanghai, China) was used for detecting HRP-conjugated antibody.

Xenograft Mouse Model and Metastasis Model

Male BALB/c nude mice (5 weeks old) were obtained from SLACOM (Shanghai, China) and used as xenograft hosts. Nude mice were maintained under a specific pathogen-free condition with approval of the Animal Care Committee of Fudan University.

For subcutaneous tumor formation assay, 1 × 107 transfected T24 cells were suspended in 0.2 ml of PBS and subcutaneously injected (n = 6 for each group). From day 4 on, tumor size was measured every 5 days by a caliper and was recorded as volume (1/6 × Π × L × W × H). The mice were sacrificed for dissection 21 days after injection.

For metastasis analyses, 1 × 105 transfected T24-luc cells were intravenously injected into mice tails (n = 6 mice per group). After 28 days, mice were intraperitoneally injected with 150 mg D-luciferin/kg body weight (Beyotime, Shanghai, China), then tumor development was monitored by bioluminescence imaging with the IVIS Spectrum in vivo Imaging System and Living Image software (PerkinElmer, Waltham, MA, United States).

Flow Cytometry

Twenty fresh resected bladder cancer tissues were obtained from Huashan Hospital. The flow cytometry procedure is illustrated below. Tumor tissues were digested with 1 mg/ml of collagenase V (Sigma, St. Louis, MO, United States, C9263) and 0.2 mg/ml of DNase I (Sangon Biotech, Shanghai, China, A610099) for 12 h, incubated in RBC lysis buffer for 5 min and then filtered through a 70-μm strainer to collect single-cell populations. The single cell suspensions were stained with membrane antibodies including CD45 (BD Biosciences, Franklin Lakes, NJ, United States, 557833), CD3 (BD Biosciences, 552852), CD8 (BD Biosciences, 564526), PD-1 (BD Biosciences, 561272), Tim-3 (BD Bioscience, 565558), and LAG-3 (BD Biosciences, 565616) for 30 min at 4°C after Fc receptor blocking. The single-cell suspensions were then treated with BD Fixation/Permeabilization Solution Kit with BD Golgistop (BD Biosciences, 554715) to stain intracellular proteins including IFN-γ (BD Biosciences, 557643), GZMB (BD Biosciences, 560212), and perforin (BD Biosciences, 556437). Cells were collected and then analyzed on a BD FACSVerse Flow Cytometer (BD Biosciences, Franklin Lakes, NJ, United States).

Statistical Analysis

The databases in our research were used for generating survival plots to analyze HR and p-value through log-rank test. Spearman’s correlation was used to gauge correlation between particular variables.

Results

Dysregulated Mutation and Expression of Rho-GTPase-Activating Proteins Genes Correlate With RFS and OS in Bladder Cancer Patients

With the discovery and validation of Rho GTPase-activating proteins, we merely included ARHGAP1 to ARHGAP49 in our work with the summary of subfamily in the ARHGAP family genes (Table 1 and Figure 1A).

TABLE 1
www.frontiersin.org

Table 1. Summary of sub-family in ARHGAP family genes.

FIGURE 1
www.frontiersin.org

Figure 1. Identification of expression profile and prognosis related to the Rho-GTPase-activating proteins (ARHGAP) family gene in bladder cancer (BCa). (A) ARHGAP mutation level in BCa according to the analysis of cBioPortal. (B) Expression level of ARHGAP5, ARHGAP17, and ARHGAP24 in several tumor and normal tissues and Expression, OS and DFS of ARHGAP5 in BCa.

Overall, ARHGAP genes harbored a 10% mutation level in BCa according to the analysis of cBioPortal (Figure 1A). The expression of ARHGAP family genes was analyzed in several cancer types in Oncomine, and a relatively lower ARHGAP gene expression was observed in BCa (Supplementary Figure 1A). Moreover, we found that higher expression of ARHGAP mainly predicted unpreferable prognosis in BCa. We further attempted to ascertain the BCa prognosis-related members of ARHGAP family genes and identified that ARHGAP5, ARHGAP8, ARHGAP11A, ARHGAP17, ARHGAP24, ARHGAP37 (STARD13), and ARHGAP38 (STARD8) correlated with either or both RFS and OS in BCa (Figure 1B and Supplementary Figure 1B).

Correlation and Functional Enrichment Analysis Reveals Rho-GTPase-Activating Proteins Family Genes Correlate With Immune-Related Biological Processes in BC Patients

We first identified the correlation of the ARHGAP family genes in BC patients’ expression levels. As Supplementary Figure 2A showed, colors in blue or red revealed positive or negative correlation, respectively. For example, ARHGAP37 (STARD13) and coexpression with ARHGAP38 (STARD8), positively correlated with ARHGAP31 and ARHGAP24, meaning these ARHGAP genes may exert a synergism biological effect. Then we explored the coexpression genes of the ARHGAP family genes by GeneMANIA (Supplementary Figure 2B). We also identified 50 relevant altered genes by STRING database and constructed PPI network (Supplementary Figure 2C). GO and KEGG functional enrichment analysis were applied to explore the functional process and novel pathways of the ARHGAP family genes. Figure 2A showed 10 most relevant functional processes of ARHGAP family genes in biological process (BP), cellular component (CC), and molecular function (MF), respectively, most of which enriched as regulating GTPase activity and downstream actin organization. As shown in Figure 2B, axon guidance, focal adhesion, and leukocyte transendothelial migration were associated with the ARHGAP family genes. Regulation of small GTPase-mediated signal transduction, positive regulation of GTPase activity, regelation of GTPase activity, etc., were the three most according functional processes (Figure 2C). To further distinguish the correlation of reciprocal genes of the ARHGAP family, we applied GO and KEGG analysis of 50 altered neighboring genes identified by STRING (Figures 3A,B). As we expected, focal adhesion, chemokine signaling pathway, and T-cell receptor signaling pathway were enriched. Then the protein–protein interaction enrichment analysis network was constructed by Metascape following the previous manuscript. As shown in Figure 3C, the subsets that contain the close genes of enriched terms have been selected and presented as a network plot. Besides, we applied Molecular Complex Detection (MCODE) algorithm to identify hub genes among the correlating genes. Pathways and functional enrichment analysis were performed by each MCODE component independently. Each three most significant pathways of MCODE are shown in Supplementary Table 1. All results above implied an association between the ARHGAP family genes and tumorigenesis, migration, and tumor-immunologic microenvironment. ARHGAP regulates GTPase-mediated cell cytoskeleton organization including actin or filament-related focal adhesion and cell motility, which is crucial for immune cells, especially tumor-infiltrating cells for transendothelial migration from the bloodstream to participate in the tumor-immunologic microenvironment. As focal adhesion, leukocyte transendothelial migration from the ARHGAP family genes and focal adhesion, the chemokine signaling pathway and T-cell receptor signaling pathway were enriched from the ARHGAP family genes, we made the hypothesis that the ARHGAP genes regulate immune-related tumor microenvironment in BCa (Figure 3D).

FIGURE 2
www.frontiersin.org

Figure 2. Functional enrichment analysis of the ARHGAP family genes in BCa. (A) Top 10 functional enrichment of biological process (BP), cellular component (CC), and molecular function (MF) of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway of the ARHGAP family genes in BCa. (B) The KEGG pathway analysis of the ARHGAP family genes in The Cancer Genome Atlas (TCGA) cohorts, p < 0.05. (C) Top 10 enrichment of the Gene Ontology (GO) pathway enrichment.

FIGURE 3
www.frontiersin.org

Figure 3. Network and potent construction analysis of 50 ARHGAP family neighboring genes in BCa. (A) KEGG pathway analysis of 50 most neighboring correlated genes, p < 0.05. (B) GO analysis of 50 most neighboring correlated genes, p < 0.05. (C) Nodes contain differential genes, which are typically close to each other are displayed with the same color. Terms with a similarity > 0.3 are connected by edges. (D) Protein–protein interaction network and Molecular Complex Detection (MCODE) components identified of the hub correlating genes.

Expression of Prognosis-Related Rho-GTPase-Activating Proteins Family Genes Correlated With Immune Cell Infiltration in Bladder Cancer

Ras homolog family member A (RhoA) and CDC42 are key downstream components of ARHGAP correlating with leukocyte transendothelial, chemokine signaling pathway, and T-cell receptor signaling pathway in enriched ARHGAP family genes PPI network, indicating that the ARHGAP family genes are involved in several immune pathways in immune-related tumor microenvironment in BCa (Figure 4A). As such, we next explored the relationship between prognosis-related ARHGAP family gene expression and the degree of immune cell infiltration in BCa using xCell and TIMER to analyze the tissue-infiltrating cell-type abundance from bulk RNA-seq data in TCGA database. First, all prognosis-related ARHGAP family gene correlates with tumor purity (Figure 4B and Supplementary Figure 3A). Interestingly, the ARHGAP family genes correlated with immune cell infiltration in different levels. ARHGAP5 and ARHGAP11A showed no significant correlation with immune infiltration score, while ARHGAP8, ARHGAP17, ARHGAP24, ARHGAP37 (STARD13), and ARHGAP38 (STARD8) significantly correlated with immune infiltration score (Figure 4C and Supplementary Figure 3B). For example, a higher ARHGAP5 correlated with lower Th1/Th2 cell ratio, higher Treg cell, and lower M1 macrophage infiltration, indicating a relatively tumor-promoting microenvironment. A higher ARHGAP8 correlated both with lower Th1/Th2 cell ratio and lower M1 macrophage infiltration, which indicated a tumor-promoting microenvironment and lower DC cell infiltration, which contributed to a tumor suppression microenvironment. Lower Th1/Th2 cell ratio, higher DC cell infiltration, and higher Treg cell infiltration were observed in ARHGAP17, ARHGAP24, and ARHGAP37 (STARD13), indicating a tumor-promoting microenvironment (Figures 4B,D, and Supplementary Figure 3C). We further explored the correspondence between the ARHGAP family gene expression and the type of immune infiltration based on immuno-related markers in BCa, especially Tfh cells, Th17 cells, and exhausted T cells in addition to our previous analysis. Besides further validation of the infiltration of monocyte, TAM, M1 macrophage, M2 macrophage, neutrophils, NK cell, DC, Th1, and Th2 cell in prognosis-related ARHGAP family genes in BCa, we further found that higher ARHGAP17, ARHGAP24, ARHGAP37 (STARD13), andARHGAP38 (STARD8) expression positively correlated with T-cell exhaustion markers: PD-1, CTLA4, TIM3, GZMB, and LAG3 (Damo and Joshi, 2019), while ARHGAP5 and ARHGAP8 showed negative or absence of correlation with T-cell exhaustion markers (Table 2). CD8 + T cells exert antitumor activity through regulation of Th1/Th2 cells, Treg cells, DC cells, and macrophages (Xiao et al., 2020) and higher T-cell exhaustion markers indicating a dysfunction state of CD8 + T-cell-related tumor-suppressing microenvironment. In conclusion, through comprehensive analysis of immune cell infiltration in BCa, we identified that the ARHGAP family genes correlate with a tumor-promoting microenvironment.

FIGURE 4
www.frontiersin.org

Figure 4. Correlations between prognosis-related ARHGAP family gene-associated immune cell infiltration. (A) Immune-related ARHGAP family gene protein–protein interaction (PPI) network was constructed in STRING. (B) Immune infiltration correlation with ARHGAP5, ARHGAP17, and ARHGAP24 in TCGA-BCa samples were carried out by TIMER. (C) Boxplots were used to visualize ARHGAP family-associated certain cell type enrichment scores and logarithm-transformed immune scores, stroma scores, and microenvironment scores of different groups through xCell-analyzed TCGA samples in BCa. (D) ARHGAP family-associated infiltration of CD8+ T cells, CD4+ T cells, type 1 T-helper (Th1) cells, type 2 T-helper (Th2) cells, natural killer (NK) cells, Treg cells, M1 macrophages, M2 macrophages, and dendritic cell were analyzed through xCell-analyzed TCGA samples in BCa. Grouping was done according to the expression level of prognosis-related ARHGAP family gene (*p < 0.05; **p < 0.01; ***p < 0.001; ns, not statistically significant, Wilcoxon signed rank test).

TABLE 2
www.frontiersin.org

Table 2. Correspondence between ARHGAP family gene expression and type of immune infiltration based on immuno-related markers in BCa.

Knockdown of ARHGAP5, ARHGAP17, and ARHGAP24 Suppressed BC Cell Proliferation, Migration, and Metastasis in vitro and in vivo

We further began to study the biological role of prognosis-related ARHGAP genes in BCa. Through system analysis of ARHGAP family gene expression in BCa cell lines through CCLE (Barretina et al., 2012; Figure 5A), we selected ARHGAP5, ARHGAP17, and ARHGAP24 for further exploration because of their relatively high expression and immuno-related role in BCa. ARHGAP5 and ARHGAP24 were significantly downregulated in BCa cell lines compared with normal SV-HUC-1 cells and BCa tumor compared with adjacent normal tissue, while ARHGAP17 showed no significant difference (Figures 5B,C and Table 3). We designed three shRNAs targeting ARHGAP5, ARHGAP17, and ARHGAP24, respectively, and the most significant reduction in parental gene expression of shRNA was selected. Cell cloning formation assays showed that knockdown of ARHGAP5, ARHGAP17, and ARHGAP24 significantly inhibited proliferation of T24 and UM-UC-3 (Figures 5D,E). The EdU assay presented similar results based on detecting cellular DNA synthesis (Figures 5F,G). T24 cells transfected with scramble shRNA (sh-NC) and shRNA targeting ARHGAP5, ARHGAP17, and ARHGAP24 were injected into nude mice subcutaneously, and knockdown of these genes significantly reduced BC tumor growth (Figures 5H–J). To verify the metastasis function of the ARHGAP family in vivo, we chose the tail vein injection lung metastasis mouse model. The metastatic foci were significantly visible from mice injected with T24 cells in comparison with sh-T24 cells (Figure 6A). Through correspondence with the micropathological results, H&E-stained tissues showed that sh-T24 cells formed up fewer nodule areas in the lungs (Figures 6A,B). Bioluminescence imaging revealed that silencing ARHGAP5, ARHGAP17, and ARHGAP24 reduced the metastasis potential of T24 cells (Figures 6C,D). Transwell assay suggested that ARHGAP5, ARHGAP17, and ARHGAP24 knockdown inhibited migration of T24 and UM-UC-3 (Figures 6E,F). Taken together, our data showed that, in addition to the immuno-related role in BCa, silencing ARHGAP5, ARHGAP17, and ARHGAP24 suppressed BC cell proliferation, migration, and metastasis.

FIGURE 5
www.frontiersin.org

Figure 5. Knockdown of ARHGAP5-, ARHGAP17-, and ARHGAP24-suppressed BC cell proliferation in vivo and in vitro. (A) System analysis of expression of ARHGAP family gene in BCa cancer cell lines in the CCLE database. (B) ARHGAP5, ARHGAP17, and ARHGAP24 expression in seven pairs of low-grade BCa and seven pairs of high-grade BCa; GAPDH serves as a loading control. (C) ARHGAP5, ARHGAP17, and ARHGAP24 expression in four bladder transitional cell carcinoma cell lines (RT4, UM-UC-3, T24, and 5637) compared with immortalized epithelium cell line (SV-HUC-1). (D,E) Cloning formation of T24, UMUC-3, and ARHGAP5, ARHGAP17, and ARHGAP24 knocking down T24, UMUC-3 (***p < 0.001 versus shNC). (F,G) 5-Ethynyl-2′-deoxyuridine (EdU) assay of T24 and shARHGAP5-, ARHGAP17-, and ARHGAP24-transfected T24 (***p < 0.001 versus shNC). (H) Tumor growth size of subcutaneous T24 tumors transfected with shNC or shARHGAP5, ARHGAP17, and ARHGAP24 21 days after tumor implantation. (I) Tumor growth curve of subcutaneous T24 tumors transfected with shNC or shARHGAP5, ARHGAP17, and ARHGAP24 in 21 days of tumor implantation. (J) Average weights of tumors from shNC and shARHGAP5, ARHGAP17, and ARHGAP24 are shown in the histogram. Error bars indicate SD.

TABLE 3
www.frontiersin.org

Table 3. Clinicopathological features in 14 pairs of bladder cancer frozen in liquid nitrogen for analysis of ARHGAP isoforms.

FIGURE 6
www.frontiersin.org

Figure 6. Knockdown of ARHGAP5, ARHGAP17, and ARHGAP24 suppressed the metastasis and migration of BC cell in vivo and in vitro. (A,B) Representative hematoxylin and eosin-stained taken from shNC-T24 and shARHGAP5, ARHGAP17, and ARHGAP24-T24-bearing mice (***p < 0.001 versus shNC). (C,D) Bioluminescence imaging in mice T24 metastasis model and knockdown ARHGAP5, ARHGAP17, and ARHGAP24 T24 cell metastasis model (***p < 0.001 versus shNC). (E,F) Transwell assay showed that ARHGAP5, ARHGAP17, and ARHGAP24 knockdown hindered migration of T24 and UM-UC-3 (***p < 0.001 versus shNC).

Discovering Underlying Molecular Mechanism of ARHGAP5, ARHGAP17, and ARHGAP24 in Bladder Cancer

We further performed gene set enrichment analysis (GSEA) of ARHGAP5, ARHGAP17, and ARHGAP24 in BCa based on TCGA database to analyze downstream pathway. As for immune-related pathway, ARHGAP5 (Figure 7A and Supplementary Table 2) only correlates with TGF-β signaling, while ARHGAP17 (Figure 7B, Supplementary Table 3) and ARHGAP24 (Figure 7C, Supplementary Table 4) correlate with TGF-β, TNF-α, IL-2/STAT5, IL-6/JAK/STAT3, and the inflammatory response pathway, which, in line with our previous result, showed an immune-related role of ARHGAP17 and ARHGAP24 rather than ARHGAP5 (Figure 4C). Also, we observed that ARHGAP5 correlates with HALLMARK_PI3K_AKT_MTOR_SIGNALING, HALLMARK _KRAS_SIGNALING, and HALLMARK_HEDGEHOG_SIGNALING, while ARHGAP17 correlates with HALLMARK_APOPTOSIS, HALLMARK_EPITHELIAL_MESENCHYMAL_TRANSITION, HALLMARK_PI3K_AKT_MTOR_SIGNALING and HALLMARK_ KRAS_SIGNALING, and ARHGAP24 correlates with HALLMARK_KRAS_SIGNALING, HALLMARK_ HEDGEHOG_SIGNALING, HALLMARK_ EPITHELIAL_ MESENCHYMAL_ TRANSITION, and HALLMARK_ APOPTOSIS. We further validated the key pathway indicator protein via particular knocking down of ARHGAP5, ARHGAP17, and ARHGAP24 in T24 and UMUC-3. Significant downregulation of p-s6k was observed in sh-ARHGAP17-T24 or UMUC-3, not in ARHGAP5 or ARHGAP24 knocking down cells. Inhibition of vimentin was specifically seen after knocking down ARHGAP17 and ARHGAP24, indicating that ARHGAP17 and ARHGAP24 could activate EMT in BCa cells. Activation of YAP was abolished in ARHGAP5 knocking down BCa cells T24 and UMUC-3, while inhibition of YAP was only observed in ARHGAP17, and ARHGAP24 knocked down UMUC-3. No significant cleaved form of PARP1 or caspase 9 was observed in ARHGAP knocked down BCa, indicating little contribution of ARHGAP in cellular apoptosis. Furthermore, significant deficiency of p-STAT3 and p-STAT5 in shARHGAP17 and shARHGAP24 BCa cells were validated according to GSEA enrichment (Figure 7D). GSEA analysis based on TCGA database and validation experiment indicated that ARHGAP5, ARHGAP17, and ARHGAP24, representing the ARHGAP family, showed, to some extent, difference, rather than similarity, in downstream pathway, performing diverse biological activities while serving as the activator of Rho GTPase.

FIGURE 7
www.frontiersin.org

Figure 7. Molecular mechanism of representative ARHGAP family genes in BC biological process. Gene set enrichment analysis (GSEA) plot of (A) ARHGAP5, (B) ARHGAP17, and (C) ARHGAP24, p < 0.05. False discovery rate (FDR) < 0.25. (D) Western bolt analysis of p-STAT3, STAT3, p-STAT5, STAT5, p-S6K1, S6K1, p-YAP1, YAP1, vimentin, caspase 9, and Parp1 in T24, UMUC-3, and ARHGAP5, ARHGAP17, and ARHGAP24 knocking down T24, UMUC-3.

ARHGAP17 and ARHGAP24 Correlated With CD8+ T Cells and Treg Infiltration and T-Cell Exhaustion Markers in Bladder Cancer Microenvironment

We have demonstrated that the ARHGAP family genes correlated with an immuno-related tumor-promoting microenvironment through bioinformatic analysis, and our in vivo and in vitro experiments validated that ARHGAP5, ARHGAP17, and ARHGAP24 promoted the proliferation, migration, and metastasis phenotype of BCa cells. However, these in silico results, which indicated ARHGAP expression, correlating with a tumor-supporting microenvironment, need further validation to elucidate the relative contribution of tumor and leukocyte ARHGAP in BCa. We performed immunohistochemistry staining on 90 tissues of TMA (Table 4) and conducted flow cytometry for the tumor-infiltrating CD8 + T cells to explore the association between ARHGAP family and tumor immune microenvironment. ARHGAP5, ARHGAP17, and ARHGAP 24 were detected in bladder cancer tumor sections as well as CD8 + T cells, CD4 + T cells, and Treg cells (Figures 8A,B). We found that decreased CD8 + T cells were infiltrated in ARHGAP17 and ARHGAP 24 high groups than the low ones. As in in silico analysis, expression of ARHGAP5, ARHGAP17, and ARHGAP24 shows no correlation with CD4 + T-cell infiltration (Figure 8C). Moreover, the Treg cells showed an increased infiltration in the ARHGAP24 high group. The function of CD8 + T cells were detected through flow cytometry (Table 5 and Supplementary Figure 4). Interestingly, less IFN-γ and GZMB were expressed on CD8 + T cells in high vs. low ARHGAP17 and ARHGAP24 groups, and the CD8 + T cells in the high ARHGAP 17 group expressed declined perforin, which indicate a dampened antitumor function of CD8 + T cells with high ARHGAP17/24 expression (Figure 8D). Moreover, PD-1 was highly expressed on CD8 + T cells in the high ARHGAP5/17/24 groups, and more CD8 + T cells in the high ARHGAP17/24 groups expressed Tim-3 (Figure 8E). These results demonstrated that ARHGAP17/24 may lead to an immunosuppressive tumor microenvironment and an impaired antitumor state of CD8 + T cells in bladder cancer. Considering that higher T-cell exhaustion markers indicating a CD8 + T-cell-related tumor-promoting microenvironment and CD8 + T cells exert antitumor activity through Treg cell regulation (Xiao et al., 2020), our data demonstrated that ARHGAP17 and ARHGAP24 correlated with CD8 + T cells and Treg cell infiltration and function in BCa microenvironment. We noticed that ARHGAP5 showed less correlation with leukocyte infiltration and function. Interestingly, ARHGAP17 and ARHGAP24 influenced immuno-related IL-2/STAT5 and IL-6/JAK/STAT3 pathways, while ARHGAP5 mainly correlated pathways, which have less interaction with immune related environment. Also, in Figure 4C, ARHGAP5 showed no correlation with immuno-microenvironment, while ARHGAP17 and ARHGAP24 correlated with immuno-microenvironment in BCa through xcell analysis. In conclusion, our results indicated that through crosstalk of tumor and leukocyte in BCa, ARHGAP17 and ARHGAP24 correlate with a tumor-promoting microenvironment through regulating CD8 + T cells and Treg infiltration and T cell function.

TABLE 4
www.frontiersin.org

Table 4. Clinicopathological features in Tissue Microarray of 90 pairs of bladder cancer.

FIGURE 8
www.frontiersin.org

Figure 8. Immunosuppressive tumor microenvironment was associated with ARHGAP family expression. (A,B) Representative immunohistochemistry staining of the ARHGAP family [(A); ARHGAP5, ARHGAP17, and ARHGAP24] and immune cells [(B); CD8 + T cells, CD4 + T cells, and Treg cells]. (C) Comparison of CD8 + T cell, CD4 + T cell, and Treg cell infiltration in bladder cancer according to ARHGAP5, ARHGAP17, and ARHGAP24 high/low groups. (D,E) Comparison of effector molecules [(D); IFN-γ, GZMB, and perforin] and immune checkpoint [(E); PD-1, Tim-3, and LAG-3] in bladder cancer according to ARHGAP5, ARHGAP17, and ARHGAP24 high/low groups.

TABLE 5
www.frontiersin.org

Table 5. Clinicopathological features of 20 fresh bladder cancer for Flow cytometry analysis.

Rho-GTPase-Activating Proteins Maintain Malignancy-Related Cellular Mechanical Properties of Bladder Cancer

Cell motility, mostly determined by cytoskeleton mediated by ARHGAP-activated Rho GTPase, played a vital role in tumor invasion and metastasis as well as immune cell chemotaxis and infiltration. To further validate ARHGAP-related cell motility, we measured cellular mechanical properties of BCa (Hecht et al., 2012). Generally, cancer cells are typically less “stiff” than normal cells with a trend in a reduced Young’s modulus (YM) value and an increase in adhesive force (AF) via AFM quantifying cellular nanomechanical properties (Plodinec et al., 2012; Chen et al., 2013). Cells with a low YM are more likely to undergo transendothelial migration or translocation into the tissues under bloodstream pressure (Dufrêne and Pelling, 2013; Li et al., 2014). A higher AF contributes to cell adhesion in the microenvironment of targeted tissues (Park and Lee, 2014). Both topography and mechanical properties of T24 and T24-shARHGAP5, ARHGAP17, and ARHGAP24 were measured. T24 cells showed an isotropic well-spread shape, while sh-ARHGAP cells showed a more spindle-like or irregular shape with relative ARHGAP knocking down. As shown in the elastic maps, the cell edge of each cell showed a larger Young’s modulus value than the cell center, which corresponded to the area of cell nucleus. This high Young’s modulus of the cell edge is due to the substrate and is consistent with a previous report (Park and Lee, 2014; Weder et al., 2014; Figure 9A). We also measured irregular phalloidin, which indicates dysfunction of cytoskeleton system and reduction in cell motility, in ARHGAP5 and ARHGAP24 knocked down T24 cells. The T24 cells possess the lowest Young’s modulus, the highest adhesive force and cell height, while the T24-shARHGAP cells possess relatively higher Young’s modulus and lower adhesive force and cell height, suggesting that ARHGAP maintained malignancy-related cellular mechanical properties of BCa (Figure 9B).

FIGURE 9
www.frontiersin.org

Figure 9. Identification the nanomechanical parameters via knocking down of ARHGAP5, ARHGAP17, and ARHGAP24 in BCa. (A) High-magnification images of the distinctive contrast patterns observed for T24 and knocking down ARHGAP5, ARHGAP17, and ARHGAP24 T24 cells. First-row images correspond to Young’s modulus (YM) maps. Second-row images correspond to cell height. Third-row images correspond to AF maps. (B) All cellular nanomechanical properties were analyzed in GraphPad Prism5 (*p < 0.05; **p < 0.01; ***p < 0.001; ns, not statistically significant, Wilcoxon signed rank test).

Discussion

Rho-GTPases are GTP-binding cytoskeletal proteins that were thought to construct mechanical support to the cell membrane for preserving biological processes (Kumar and Epstein, 2011). Regulation of Rho-GTPases can be highly related to cancer progression (Li et al., 2020). Rho-GTPase-activating proteins (ARHGAP) serve as negative regulators of Rho-GTPases, and alteration in the genes and transcriptomes of the ARHGAP family can exert a carcinogenesis effect through Rho-like GTPase dysregulation in several cancer types. ARHGAP7 downregulation correlates with unfavorable prognosis in breast cancer patients especially in metastatic lesions (Chen et al., 2019). In contrast, ARHGAP18 upregulation correlates with favorable prognosis in breast cancer (Humphries et al., 2017). ARHGAP15, induced via androgen, plays a tumor suppressor role (Takagi et al., 2018). ARHGAP regulates several pathways of cytoskeleton dynamic remodeling and assembling, which may exert an important role in immune cell migration, thus, eventually influencing tumor immune activity and infiltration. However, the role of the ARHGAP family has not been clearly identified in BCa, especially in immune-related BCa microenvironment, which has not been elucidated.

Through integrative analysis of ARHGAP genes in BCa via cBioPortal, Oncomine, and Gepia, we found that the ARHGAP family mutates at 10% level in BCa, and a relatively lower ARHGAP gene expression was observed in BCa. ARHGAP5, ARHGAP8, ARHGAP11A, ARHGAP17, ARHGAP24, ARHGAP37 (STARD13), and ARHGAP38 (STARD8) were identified as prognosis-related ARHGAP genes in BCa.

GeneMANIA, STRING, and MCODE analysis of neighboring and coexpression genes of the ARHGAP family followed by GO and KEGG functional enrichment identified axon guidance, focal adhesion, and leukocyte transendothelial migration to be associated with the ARHGAP family genes, and focal adhesion, chemokine signaling pathway, and T-cell receptor signaling pathways were enriched in ARHGAP coexpression genes. Surprisingly, the ARHGAP family genes and coexpression genes were found to influence several non-cytoskeleton-related pathways including the chemokine signaling pathway and T-cell receptor signaling pathway in the tumor immune microenvironment.

Immune cell infiltration in a cancer-associated immune microenvironment is closely related with prognosis and immunotherapy efficiency in several cancer types. High TILs often relate to a better prognosis than low TILs. Tumor-infiltrating CD4 + T cells, CD8 + T cells, tumor-associated macrophages (TAMs), and neutrophils were associated with patient prognosis and tumor chemosensitivity. We analyzed the correlation between prognosis-related ARHGAP family gene expression and the degree of immune cell infiltration in BCa using xCell and TIMER. Higher ARHGAP5 correlates with lower Th1/Th2 cell ratio, higher Treg cell, and lower M1 macrophage infiltration, which indicates a relative tumor-promoting microenvironment. Lower Th1/Th2 cell ratio, higher DC cell infiltration, and higher Treg cell infiltration were observed in ARHGAP17, ARHGAP24, and ARHGAP37 (STARD13), indicating a tumor-promoting microenvironment. Higher ARHGAP17, ARHGAP24, ARHGAP37 (STARD13), and ARHGAP38 (STARD8) expression positively correlated with T-cell exhaustion markers: PD-1, CTLA4, TIM3, GZMB, and LAG3, while ARHGAP5 and ARHGAP8 shows negative or absence correlation with T-cell exhaustion markers. The CD8 + T cell serves as the most crucial member of the immune-related tumor microenvironment. They mediate tumor cell-specific immune responses through themselves or interaction with other immune cells. High infiltration and low T-cell exhaustion marker of CD8 + T cells predict a higher response to chemotherapy and an overall better prognosis. We also observed a higher DC cell infiltration and a higher Treg cell infiltration in tumor-suppressing ARHGAP microenvironment, which indicates that the ARHGAP family gene is closely linked with tumor DC penetration. DCs increase tumor cell metastasis via enhancing Treg cells and suppressing CD8 + T-cell antitumor cytotoxicity. Collectively, we demonstrated that immune cell infiltration in prognosis-related ARHGAP family genes correlates with a tumor-promoting microenvironment.

Furthermore, we validated that silencing ARHGAP5, ARHGAP17, and ARHGAP24 suppressed BC cell proliferation, migration, and metastasis. Immune-related TGF-β, TNF-α, IL-2/STAT5, IL-6/JAK/STAT3, and the inflammatory response pathway and other signal pathways were related in ARHGAP5-, ARHGAP17-, and ARHGAP24-related BC carcinogenesis. Finally, since the different nanomechanical properties between normal and cancerous living cells enable cellular mechanical properties to be efficient markers for grading malignant-level cancer cells, we validated that T24-shARHGAP cells possess a relatively higher Young’s modulus and lower adhesive force and cell height, suggesting that ARHGAP maintains malignancy-related cellular mechanical properties of BCa. We still need to test the nanomechanical properties of immune cell mediated via knocking down ARHGAP. Taken together, we found that the ARHGAP family genes promote BC progressing through establishing a tumor-promoting microenvironment and promote BC cell proliferation, migration, and metastasis through cellular mechanical property-mediated cell motility (Figure 10).

FIGURE 10
www.frontiersin.org

Figure 10. Scheme of ARHGAP family gene promotes BC progressing through establishing a tumor-promoting microenvironment and promotes BC cell proliferation, migration, and metastasis though cellular mechanical property-mediated cell motility.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The animal study was reviewed and approved by Animal Care Committee of Fudan University.

Author Contributions

HJ, YC, NL, and ZM were responsible for the conception and design, and study supervision. SW, QZ, and ZZ were responsible for the development of the methodology, analysis, and experiments. XC, YO, CY, CX, and XD performed the statistical and bioinformatic analysis. All authors read and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 81802569) and Clinical Research and Training Project (SHDC12017X10) provided by Shanghai Shen Kang Hospital Development Center.

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.

Supplementary Material

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

Supplementary Figure 1 | Identification of expression profile and prognosis related ARHGAP family gene in BCa. (A) Expression level of ARHGAP family gene in several tumor and normal tissues in Oncomine. (B) Expression level of ARHGAP8, ARHGAP11A, ARHGAP37, and ARHGAP38 in several tumor and normal tissues and Expression, OS and DFS of ARHGAP8, ARHGAP11A, ARHGAP37, and ARHGAP38 in BCa.

Supplementary Figure 2 | Functional enrichment analysis of ARHGAP family genes in BCa. (A) Correlation analysis of ARHGAP family genes. (B,C) Protein-protein interaction network of ARHGAP family genes and relevant co-expression genes through GeneMANIA and STRING analysis.

Supplementary Figure 3 | Correlations between prognosis related ARHGAP family gene associated immune cell infiltration. (A) Immune infiltration correlation with ARHGAP8, ARHGAP11A, ARHGAP37, and ARHGAP38 in TCGA-BCa samples were carried out by TIMER. (B) Boxplots were used to visualize ARHGAP family associated certain cell type enrichment scores and logarithm transformed immune scores, stroma scores, microenvironment scores of different groups through xCell analyzed TCGA samples in BCa. (C) ARHGAP family associated infiltration of CD8+ T cells, CD4+ T cells, type 1 T helper (Th1) cells, type 2 T helper (Th2) cells, natural killer (NK) cells, Treg cells, M1 macrophages, M2 macrophages and Dendric cell were analyzed through xCell analyzed TCGA samples in BCa. Grouping was done according to the expression level of prognosis related ARHGAP family gene. (p < 0.05; ∗∗p < 0.01; ∗∗∗p < 0.001; ns, not statistically significant, Wilcoxon signed rank test).

Supplementary Figure 4 | Gating strategy for CD8+ T cells.

Abbreviations

AFM, atomic force microscopy; ARHGAP, Rho GTPase-activating protein; BCa, bladder cancer; BP, biological process; CC, cellular component; CCK8, Cell Counting Kit-8; DC, dendritic cell; DMEM, Dulbecco’s modified Eagle’s medium; EdU, 5-ethynyl-2′-deoxyuridine; EMT, epithelial–mesenchymal transition; FDR, false discovery rate; FBS, fetal bovine serum; FPKM, fragments per kilobase million; KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology; GTEx, Genotype-Tissue Expression; GSEA, gene set enrichment analysis; IL-2, interleukin 2; IL-6, interleukin 6; MCODE, Molecular Complex Detection; MF, molecular function; NES, normalized enrichment score; NK, natural killer; PFA, paraformaldehyde; PPI, protein–protein interaction; RhoA, Ras homolog family member A; TAMs, tumor-associated macrophages; TCGA, The Cancer Genome Atlas; TGF-β, transforming growth factor beta; TNF-α, tumor necrosis factor alpha; Th1, type 1 T helper; Th1, type 2 T helper; WB, Western blotting; YM, Young’s modulus.

Footnotes

  1. ^ http://www.cbioportal.org
  2. ^ https://www.oncomine.org/resource/login.html
  3. ^ https://cistrome.shinyapps.io/timer/
  4. ^ http://gepia.cancer-pku.cn/
  5. ^ https://www.gsea-msigdb.org/gsea/msigdb
  6. ^ https://genemania.org/
  7. ^ https://string-db.org/
  8. ^ https://metascape.org/
  9. ^ https://cistrome.shinyapps.io/timer/
  10. ^ http://xCell.ucsf.edu/

References

Aran, D., Hu, Z., and Butte, A. J. (2017). xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 18:220. doi: 10.1186/s13059-017-1349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Barretina, J., Caponigro, G., Stransky, N., Venkatesan, K., Margolin, A. A., Kim, S., et al. (2012). The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 483, 603–607. doi: 10.1038/nature11003

PubMed Abstract | CrossRef Full Text | Google Scholar

Becht, E., Giraldo, N. A., Germain, C., de Reyniès, A., Laurent-Puig, P., Zucman-Rossi, J., et al. (2016). Immune Contexture, Immunoscore, and Malignant Cell Molecular Subgroups for Prognostic and Theranostic Classifications of Cancers. Adv. Immunol. 130, 95–190. doi: 10.1016/bs.ai.2015.12.002

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. Cancer Discov. 2, 401–404. doi: 10.1158/2159-8290.cd-12-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaffer, C. L., and Weinberg, R. A. (2011). A perspective on cancer cell metastasis. Science 331, 1559–1564. doi: 10.1126/science.1203543

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C.-L., Mahalingam, D., Osmulski, P., Jadhav, R. R., Wang, C.-M., Leach, R. J., et al. (2013). Single-cell analysis of circulating tumor cells identifies cumulative expression patterns of EMT-related genes in metastatic prostate cancer. Prostate 73, 813–826. doi: 10.1002/pros.22625

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, D. S., and Mellman, I. (2017). Elements of cancer immunity and the cancer-immune set point. Nature 541, 321–330. doi: 10.1038/nature21349

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W.-X., Lou, M., Cheng, L., Qian, Q., Xu, L.-Y., Sun, L., et al. (2019). Bioinformatics analysis of potential therapeutic targets among ARHGAP genes in breast cancer. Oncol. Lett. 18, 6017–6025. doi: 10.3892/ol.2019.10949

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Zhang, J., Ruan, W., Huang, M., Wang, C., Wang, H., et al. (2020). Urine DNA methylation assay enables early detection and recurrence monitoring for bladder cancer. J. Clin. Invest. 130, 6278–6289. doi: 10.1172/JCI139597

PubMed Abstract | CrossRef Full Text | Google Scholar

Damo, M., and Joshi, N. S. (2019). Treg cell IL-10 and IL-35 exhaust CD8+ T cells in tumors. Nat. Immunol. 20, 674–675. doi: 10.1038/s41590-019-0389-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Dufrêne, Y. F., and Pelling, A. E. (2013). Force nanoscopy of cell mechanics and cell adhesion. Nanoscale 5, 4094–4104. doi: 10.1039/c3nr00340j

PubMed Abstract | CrossRef Full Text | Google Scholar

Duggan, B., Kelly, J., Keane, P. F., Williamson, K., and Johnston, S. R. (2000). Bcl-2 expression identifies patients with advanced bladder cancer treated by radiotherapy who benefit from neoadjuvant chemotherapy. BJU Int. 86:757. doi: 10.1046/j.1464-410x.2000.00895-6.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fridman, W. H., Pagès, F., Sautès-Fridman, C., and Galon, J. (2012). The immune contexture in human tumours: impact on clinical outcome. Nat. Rev. Cancer 12, 298–306. doi: 10.1038/nrc3245

PubMed Abstract | CrossRef Full Text | Google Scholar

Gen, Y., Yasui, K., Zen, K., Nakajima, T., Tsuji, K., Endo, M., et al. (2009). A novel amplification target, ARHGAP5, promotes cell spreading and migration by negatively regulating RhoA in Huh-7 hepatocellular carcinoma cells. Cancer Lett. 275, 27–34. doi: 10.1016/j.canlet.2008.09.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Haraguchi, M., Fukushige, T., Kanekura, T., and Ozawa, M. (2019). E-cadherin loss in RMG-1 cells inhibits cell migration and its regulation by Rho GTPases. Biochem. Biophys. Rep. 18:100650. doi: 10.1016/j.bbrep.2019.100650

PubMed Abstract | CrossRef Full Text | Google Scholar

Hashimoto, K., Ochi, H., Sunamura, S., Kosaka, N., Mabuchi, Y., Fukuda, T., et al. (2018). Cancer-secreted hsa-miR-940 induces an osteoblastic phenotype in the bone metastatic microenvironment via targeting ARHGAP1 and FAM134A. Proc. Natl. Acad. Sci. U. S. A. 115, 2204–2209. doi: 10.1073/pnas.1717363115

PubMed Abstract | CrossRef Full Text | Google Scholar

Hecht, E., Thompson, K., Frick, M., Wittekindt, O. H., Dietl, P., Mizaikoff, B., et al. (2012). Combined atomic force microscopy-fluorescence microscopy: analyzing exocytosis in alveolar type II cells. Anal. Chem. 84, 5716–5722. doi: 10.1021/ac300775j

PubMed Abstract | CrossRef Full Text | Google Scholar

Humphries, B., Wang, Z., Li, Y., Jhan, J.-R., Jiang, Y., and Yang, C. (2017). ARHGAP18 Downregulation by miR-200b Suppresses Metastasis of Triple-Negative Breast Cancer by Enhancing Activation of RhoA. Cancer Res. 77, 4051–4064. doi: 10.1158/0008-5472.can-16-3141

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, J., and Epstein, D. L. (2011). Rho GTPase-mediated cytoskeletal organization in Schlemm’s canal cells play a critical role in the regulation of aqueous humor outflow facility. J. Cell Biochem. 112, 600–606. doi: 10.1002/jcb.22950

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Liu, L., Xi, N., and Wang, Y. (2014). Research progress in quantifying the mechanical properties of single living cells using atomic force microscopy. Chin. Sci. Bull. 59, 4020–4029. doi: 10.1007/s11434-014-0581-2

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, Z., Wang, Q., Peng, S., Yao, K., Chen, J., Tao, Y., et al. (2020). The metastatic promoter DEPDC1B induces epithelial-mesenchymal transition and promotes prostate cancer cell proliferation via Rac1-PAK1 signaling. Clin. Transl. Med. 10:e191. doi: 10.1002/ctm2.191

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, N., Guo, J., Chen, L., Yang, W., Qu, X., and Cheng, Z. (2016). ARHGAP10, downregulated in ovarian cancer, suppresses tumorigenicity of ovarian cancer cells. Cell Death Dis. 7:e2157. doi: 10.1038/cddis.2015.401

PubMed Abstract | CrossRef Full Text | Google Scholar

Nekolla, K., Brieu, N., Gavriel, C. G., Widmaier, M., Budco, A., Medrikova, D., et al. (2019). Prognostic immunoprofiling of muscle invasive bladder cancer (MIBC) patients in a multicentre setting. Ann. Oncol. 30, v49–v50. doi: 10.1093/annonc/mdz239.067

CrossRef Full Text | Google Scholar

Park, S., and Lee, Y. J. (2014). AFM-based dual nano-mechanical phenotypes for cancer metastasis. J. Biol. Phys. 40, 413–419. doi: 10.1007/s10867-014-9353-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Plodinec, M., Loparic, M., Monnier, C. A., Obermann, E. C., Zanetti-Dallenbach, R., Oertle, P., et al. (2012). The nanomechanical signature of breast cancer. Nat. Nanotechnol. 7, 757–765.

Google Scholar

Post, A., Pannekoek, W.-J., Ross, S. H., Verlaan, I., Brouwer, P. M., and Bos, J. L. (2013). Rasip1 mediates Rap1 regulation of Rho in endothelial barrier function through ArhGAP29. Proc. Natl. Acad. Sci. U. S. A. 110, 11427–11432. doi: 10.1073/pnas.1306595110

PubMed Abstract | CrossRef Full Text | Google Scholar

Rhodes, D. R., Yu, J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., et al. (2004). ONCOMINE: a Cancer Microarray Database and Integrated Data-Mining Platform. Neoplasia 6, 1–6. doi: 10.1016/s1476-5586(04)80047-2

CrossRef Full Text | Google Scholar

Robertson, A. G., Kim, J., Al-Ahmadie, H., Bellmunt, J., Guo, G., Cherniack, A. D., et al. (2017). Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer. Cell 171, 540–556.e25.

Google Scholar

Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer statistics, 2020. CA Cancer J. Clin. 70, 7–30. doi: 10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, A., Nix, J. W., Nielsen, M. E., Wallen, E., and Pruthi, R. (2012). Prospective randomized controlled trial of robotic versus open radical cystectomy for bladder cancer: median 3-year follow-up. J. Clin. Oncol. 30:284–284. doi: 10.1200/jco.2012.30.5_suppl.284

CrossRef Full Text | Google Scholar

Takagi, K., Miki, Y., Onodera, Y., Ishida, T., Watanabe, M., Sasano, H., et al. (2018). ARHGAP15 in Human Breast Carcinoma: a Potent Tumor Suppressor Regulated by Androgens. Int. J. Mol. Sci. 19:804. doi: 10.3390/ijms19030804

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 45, W98–W102.

Google Scholar

Wang, J., Qian, J., Hu, Y., Kong, X., Chen, H., Shi, Q., et al. (2014a). ArhGAP30 promotes p53 acetylation and function in colorectal cancer. Nat. Commun. 5:4735.

Google Scholar

Wang, J., Tian, X., Han, R., Zhang, X., Wang, X., Shen, H., et al. (2014b). Downregulation of miR-486-5p contributes to tumor progression and metastasis by targeting protumorigenic ARHGAP5 in lung cancer. Oncogene 33, 1181–1189. doi: 10.1038/onc.2013.42

PubMed Abstract | CrossRef Full Text | Google Scholar

Weder, G., Hendriks-Balk, M. C., Smajda, R., Rimoldi, D., Liley, M., Heinzelmann, H., et al. (2014). Increased plasticity of the stiffness of melanoma cells correlates with their acquisition of metastatic properties. Nanomedicine 10, 141–148. doi: 10.1016/j.nano.2013.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Wittekind, C., Asamura, H., and Sobin Leslie, H. (2019). TNM Atlas. Hoboken, NJ: Wiley-Blackwell.

Google Scholar

Xiao, Z., Hu, L., Yang, L., Wang, S., Gao, Y., Zhu, Q., et al. (2020). TGFβ2 is a prognostic-related biomarker and correlated with immune infiltrates in gastric cancer. J. Cell. Mol. Med. 24, 7151–7162.

Google Scholar

Xie, R., Chen, X., Cheng, L., Huang, M., Zhou, Q., Zhang, J., et al. (2021). NONO Inhibits Lymphatic Metastasis of Bladder Cancer via Alternative Splicing of SETMAR. Mol. Ther. 29, 291–307. doi: 10.1016/j.ymthe.2020.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: bladder cancer, ARHGAP, immune infiltration, tumor microenvironment, cellular mechanical properties

Citation: Yang C, Wu S, Mou Z, Zhou Q, Zhang Z, Chen Y, Ou Y, Chen X, Dai X, Xu C, Liu N and Jiang H (2021) Transcriptomic Analysis Identified ARHGAP Family as a Novel Biomarker Associated With Tumor-Promoting Immune Infiltration and Nanomechanical Characteristics in Bladder Cancer. Front. Cell Dev. Biol. 9:657219. doi: 10.3389/fcell.2021.657219

Received: 22 January 2021; Accepted: 24 May 2021;
Published: 07 July 2021.

Edited by:

Liwu Fu, Sun Yat-sen University, China

Reviewed by:

Alessandra Ghigo, University of Turin, Italy
Tianxin Lin, Sun Yat-sen Memorial Hospital, China

Copyright © 2021 Yang, Wu, Mou, Zhou, Zhang, Chen, Ou, Chen, Dai, Xu, Liu and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chen Yang, YangC_Huashan@163.com; Haowen Jiang, haowj_sh@fudan.edu.cn; Na Liu, liuna_sia@shu.edu.cn

ORCID: Chen Yang, orcid.org/0000-0001-6011-4067; Siqi Wu, orcid.org/0000-0001-6464-2772

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.