- Open Access
Relationship between the expression of ARHGAP25 and RhoA in non-small cell lung cancer and vasculogenic mimicry
BMC Pulmonary Medicine volume 22, Article number: 377 (2022)
Vasculogenic mimicry (VM) is a recently identified pattern of blood supply to tumor tissue. It has long been considered a functional element in the metastasis and prognosis of malignant tumors. Both Rho GTPase-activating protein 25 (ARHGAP25) and Ras homolog family member A (RhoA) are effective predictors of tumor metastasis. In this study, we examined the expression levels of ARHGAP25 and RhoA and the structure of VM in non-small cell lung cancer (NSCLC). At the same time, we used cytology-related experiments to explore the effect of ARHGAP25 on the migration ability of tumor cells. Furthermore, we analyzed the interaction between the three factors and their association with clinicopathological characteristics and the five-year survival time in patients using statistical tools.
A total of 130 well-preserved NSCLC and associated paracancerous tumor-free tissues were obtained. Cell colony formation, wound healing, and cytoskeleton staining assays were used to analyze the effect of ARHGAP25 on the proliferation and migration ability of NSCLC cells. Immunohistochemical staining was used to determine the positivity rates of ARHGAP25, RhoA, and VM. Statistical software was used to examine the relationships between the three factors and clinical case characteristics, overall survival, and disease-free survival.
Cell colony formation, wound healing, and cytoskeleton staining assays confirmed that ARHGAP25 expression affects the proliferation and migratory abilities of NSCLC cells. ARHGAP25 positivity rates in NSCLC and paracancerous tumor-free tissues were 48.5% and 63.1%, respectively, whereas RhoA positivity rates were 62.3% and 18.5%, respectively. ARHGAP25 had a negative relationship with RhoA and VM, whereas RhoA and VM had a positive relationship (P < 0.05). ARHGAP25, RhoA, and VM affected the prognosis of patients with NSCLC (P < 0.05) according to Kaplan–Meier of survival time and Cox regression analyses. Furthermore, lowering ARHGAP25 expression increased NSCLC cell proliferation and migration.
ARHGAP25 and RhoA expression is associated with VM and may be of potential value in predicting tumor metastasis, prognosis, and targeted therapy.
Lung cancer accounts for approximately one in ten (11.4%) confirmed cancer cases and one in five (18.0%) fatalities, according to Global Cancer Statistics 2020, with 2.2 million confirmed cases of cancer and 1.8 million deaths, ranking second and first, respectively, among all cancer types . Non-small cell lung cancer (NSCLC), which includes adenocarcinoma of the lung and squamous cell carcinoma, is the most common histologic subtype of lung cancer according to the World Health Organization. It accounts for approximately 85% of lung malignancies . Despite tremendous advances in the treatment of NSCLC, patients with NSCLC have a dismal five-year survival rate owing to severe side effects and resistance to surgery, radiation, and chemotherapy . Therefore, it is vital to identify new antitumor targets to build on existing treatments. Angiogenesis plays a significant role in the growth, metastasis, and development of various malignancies [4,5,6]. Classical tumor angiogenesis has shown that when the tumor diameter exceeds 1–2 mm, activation of vascular endothelial cells is required for neovascularization to obtain blood supply and allow continued growth . However, Maniotis  first proposed vasculogenic mimicry (VM), a novel microcirculation pattern that can supply blood to tumors, in his study of human melanoma. VM is defined as a tubular structure of neoplastic cells that can contain blood cells through self-deformation and matrix remodeling. The unique structure of VM can cause distant metastasis in the early stage of tumor development and predict a poor prognosis .
Rho GTPases are members of the Ras superfamily and participate in various essential biological processes, including cell cycle progression, cytoskeletal reorganization, and malignant transformation ; they are an essential part of tumor development and progression . Rho GTPase-activating protein 25 (ARHGAP25) belongs to the Rho GTPase-activating protein (RhoGAP) family, which promotes endogenous hydrolysis of GTP and is a negative regulator of Rho GTPases. Thuault et al. revealed for the first time that ARHGAP25 inhibits the invasion of alveolar rhabdomyosarcoma (ARMS) cells . Previous studies have shown that ARHGAP25 overexpression significantly inhibits the growth of many neoplastic cells and suppresses their migration and invasion [13,14,15,16]. However, the association between ARHGAP25 and VM in NSCLC has not been studied. In this study, we analyzed the effect of ARHGAP25 on NSCLC cell motility, and its expression in NSCLC tissue in relation to VM, so as to explore whether ARHGAP25 is beneficial to the prognosis of patients.
Ras homolog family member A (RhoA) is an isoform of the Rho family of small-molecule G proteins, which are involved in regulating various life processes, including cytoskeleton assembly, cell adhesion, motility, cycle progression, cytokinesis, and gene transcription . RhoA is highly expressed in a range of tumor tissues, closely correlated with tumor malignancy, and plays an invaluable role in neoplastic metastasis and development . In this study, we examined RhoA expression in NSCLC, its link to VM and ARHGAP25, and its potential clinical applicability.
Materials and methods
Patients and tissue samples
We collected 130 archived paraffin-embedded NSCLC tissue specimens, without preoperative radiotherapy or chemotherapy, with postoperative immunohistochemically confirmed lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD), and paracancerous tumor-free tissues (> 5 cm from the tumor edge) , from the First Affiliated Hospital of Bengbu Medical College from January 2013 to December 2016. All patients had complete clinical, pathological, and follow-up data, and the enrolled cases were followed until death or up to 2021. To record the patients’ postoperative survival status, patients were followed-up telephonically at six-month intervals. Patient death or the December 2021 endpoint was used to compute overall survival time (OS), while patient death, relapse, or the December 2021 endpoint was used to calculate disease-free survival (DFS). The Ethics Committee of Bengbu Medical College authorized all patients to provide signed informed consent (NO. 2020KY035). The study adhered to the ethical guidelines of the Declaration of Helsinki. The American Joint Committee on Cancer’s 8th edition staging system was used to perform tumor-node-metastasis (TNM) staging of NSCLC. Table 1 shows the clinicopathological characteristics of the patients.
The collected paraffin specimens were serially sectioned at 4-μm thickness, deparaffinized with xylene, and dehydrated with graded alcohols after baking for 2 h. Sections were then washed in distilled water and phosphate-buffered saline (PBS) for 3 min each, respectively, and this was repeated thrice. Tissue slices were submerged in sodium citrate buffer (pH 6.0), pressure-cooked, chilled to 20 °C for antigen retrieval, and rinsed in PBS. To inhibit endogenous peroxidase, we performed tissue antigen retrieval and dripped a 3% hydrogen peroxide solution across all sections, followed by incubation for 15 min at 37 °C. The samples were washed thrice in PBS (pH 7.2) for 3 min each. The sections were incubated overnight at 4 °C with a few drops of rabbit monoclonal anti-ARHGAP25 (1:400, ab192020, Abcam, USA), mouse monoclonal anti-RhoA (1:100, ab54835, Abcam), and anti-CD34 (1:200, ab762, Abcam) antibodies. The secondary antibody was then added in a dropwise manner and the sections were allowed to sit for 30 min at 37 °C. Thereafter, sections were soaked in DAB solution, counterstained with hematoxylin, washed with distilled water, dehydrated with graded alcohols, and sealed with neutral resin.
Evaluation of immunostaining
Two senior pathologists assessed the slides and performed immunohistochemical analysis. The immunohistochemical scores of ARHGAP25 and RhoA protein consisted of two parts: the staining intensity (0, no staining; 1, pale yellow staining; 2, tan staining; 3, brown staining) (Fig. 1)and the number of positive cells in the tumor tissue (0, < 10%; 1, 11–50%; 2, 51–75%; 3, > 75%), with the final score being the sum of the two. A three-point total score was deemed favorable, whereas a score below 3 was considered negative. We used the method of double staining with PAS-CD34 to show the structure of VM, and we considered that VM was present in this tumor when it appeared as CD34-negative but PAS-positive .
Cell culture and SiRNA transfection
The Chinese Academy of Sciences Cell Bank (Shanghai, China) provided the lung cell lines. NCI-H1299 cells were cultivated in RPMI 1640 medium (Invitrogen, Carlsbad, CA, USA) containing 10% fetal bovine serum (FBS), whereas A549 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% FBS. The incubator temperature and CO2 content were set at 37 °C and 5%, respectively. Transfection of small interfering RNAs in NSCLC cells (Reebok Bio, Guangdong, China) according to the Lipofectamine 3000 reagent (Thermo Fisher Scientific) operating manual inhibited ARHGAP25 production when cell fusion rates reached 50–60% . The knockdown efficiency of these siRNAs was assessed by RT-PCR and Western blotting.
Quantitative reverse transcription-polymerase chain reaction analysis of ARHGAP25
Total RNA was extracted using the TRIzol technique. cDNA was prepared using the NovoStart SYBR qPCR SuperMix Plus (Novoprotein, China); quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to measure the relative expression of ARHGAP25 mRNA in cells [21, 22]. β-actin was used as an endogenous reference gene. The following primers were used: ARHGAP25: forward, 5'-GACAAGCGACTCTGATACAA-3', and reverse, 5'-GAAACATTTCCGGTTAGG-3'; β-actin forward, 5'-CATGTACGTTGCTATCCAGGC-3', and reverse, 5'-CTCCTTAATGTCACGCACGAT-3'.
Western blotting of ARHGAP25
Radioimmunoprecipitation assay (RIPA) lysis buffer (Solarbio, China) was used to extract cellular proteins, and protein concentrations were determined using a bicinchoninic acid (BCA) protein quantification kit (Solarbio, China). Absorbance was recorded at 562 nm using a microplate reader. A complete protein (25 µg) was added to each well and proteins were separated using polyacrylamide gel electrophoresis (PAGE) with 8–12% sodium dodecyl sulfate (SDS) gels and transferred to polyvinylidene fluoride (PVDF) membranes [21, 22]. After blocking with 5% skim milk for 1 h, the membranes were incubated overnight at 4 °C with primary antibodies specific for ARHGAP25 (1: 1000, Abcam) and β-actin (1: 1000, Proteintech, China). Membranes were washed with tris-buffered saline with Tween® (TBST) on a shaker before adding the appropriate secondary antibodies (1:10,000, Proteintech, China). An enhanced chemiluminescence kit (Xinsaimei, China) was used to visualize immunoreactive protein bands.
Detection of colony formation
Cells (1000 cells/well) were inoculated into six-well culture dishes and cultured in 2 ml DMEM (Gibco, USA) containing 10% FBS for two weeks. Thereafter, cells were fixed for 10 min in 1 mL of 4% neutral formaldehyde, washed three times with PBS, stained with 1 mL of 1% crystalline violet, photographed, and counted using ImageJ (ACTREC, Navi Mumbai, India) software .
Cell migration assays
Cell migration ability was tested using a wound-healing experiment. A six-well plate was seeded with cells (5.0 × 105 cells/well). When the cells were confluent, three vertical scratches were made in each well using a 100-μL sterile pipette tip. PBS was added slowly along the lateral wall, thereby flushing the cells shed on the wound. At 0 and 24 h, images were taken, samples were obtained, and the scratch lengths in each group were measured using ImageJ software .
Staining of the actin cytoskeleton
Coomassie Brilliant Blue staining solution was used to examine the cytoskeleton of cultivated NSCLC cells. First, 8 × 104 cells per group were seeded into six-well plates. The following day, adherent cells were fixed with 4% paraformaldehyde at room temperature for 20 min, infiltrated with 0.2% Triton-X-100 for 5 min, and stained with 1 mL Coomassie Brilliant Blue R-250 (0.2%). The cells were monitored and imaged using an inverted microscope after five rinses with distilled water .
Statistical scores were obtained from our experimental data using SPSS 26.0 (IBM; Chicago, IL, USA). The correlations between the expression levels of ARHGAP25, RhoA, and VM, as well as clinicopathological data, were analyzed using the chi-square test. Spearman's correlation analysis was used to determine the relationship between ARHGAP25, RhoA, and VM. The association between each component and OS and DFS was compared using Kaplan–Meier analysis. Finally, Cox regression was used to estimate the risk factors that may influence patient prognosis using univariate and multifactor analyses.
Relationship between ARHGAP25, RhoA, VM, and clinicopathological parameters
We used immunohistochemistry and SPSS 26.0 to count the positive rates of ARHGAP25, RhoA, and VM in NSCLC. ARHGAP25 was mainly localized in the nucleus, and ARHGAP25 levels were decreased in lung cancer tissues (48.5%, 63/130; Fig. 2a) compared to paracancerous tissues (63.1%, 82/130; Fig. 2b). The expression levels of ARHGAP25 in NSCLC tissues were significantly associated with tumor size (P = 0.015), TNM stage (P = 0.013), and lymph node metastasis (LNM; P = 0.015), but not with other clinicopathological factors (P > 0.05; Table 2). RhoA was primarily found in the cytoplasm, with a substantially higher expression rate in carcinomas (62.3%, 81/130; Fig. 2c) than in paracancerous tissues (18.5%, 24/130; Fig. 2d). The degree of RhoA expression was linked to tumor TNM stage (P = 0.002) and LNM (P = 0.001), but not to other clinicopathological factors (P > 0.05; Table 2). VM was present among the 130 NSCLC tissues collected in 48 cases (36.9%; Fig. 2e), whereas VM structures were absent in paracancerous tissues (Fig. 2f). The presence of VM structures in NSCLC was positively associated with larger tumor size (P = 0.016), poor grade (P = 0.004), TNM stage (P = 0.011), and LNM (P = 0.001); however, no other clinicopathological factors (P > 0.05) were found to be significantly related (Table 2).
Relationship between ARHGAP25, RhoA, and VM
ARHGAP25 expression was negatively correlated with RhoA expression (r = -0.262, P = 0.003) and VM positivity (r = -0.232, P = 0.008), whereas RhoA expression was positively correlated with VM positivity (r = 0.365, P < 0.001; Table 3).
ARHGAP25 and RhoA expression levels, VM, and clinicopathological factors affect OS
We observed a median survival time of 35.0 months as well found a mean overall survival time of 44.6 ± 2.9 months, and the five-year OS rate was 20.9%, based on the Kaplan–Meier analysis of the five-year overall survival time of the collected NSCLC cases. The ARHGAP25-positive group (56.8 ± 4.1; 32.4%) displayed a considerably higher mean survival time and five-year OS rate than the ARHGAP25-negative samples (29.7 ± 2.6; 8.8%; χ2 = 20.530, P < 0.001; Fig. 3A). Similarly, the RhoA-positive group (31.9 ± 2.8; 12.1%) showed a significantly lower OS and OS rates than the RhoA-negative samples (62.1 ± 4.4; 35.2%; χ2 = 21.604, P < 0.001; Fig. 3B). The VM-positive samples (25.2 ± 2.9; 8.0%) displayed a significantly lower OS and OS rate than the VM-negative samples (53.9 ± 3.5; 28.1%; χ2 = 23.304, P < 0.001; Fig. 3C). OS and OS rates were reduced in patients with tumor size greater than 3 cm (32.9 ± 2.7 vs. 57.3 ± 4.5; 13.1% vs. 32.4%; χ2 = 12.072, P = 0.001; Fig. 3D). The OS and OS rates gradually decreased in patients with TNM stages I, II, and III (73.9 ± 4.1 vs 41.3 ± 3.1 vs. 20.2 ± 2.0; 49.7% vs. 11.4% vs. 2.5%; χ2 = 83.749, P < 0.001; Fig. 3E). Patients with LNM presented with a considerably lower five-year OS rate than those without LNM (22.6 ± 1.8 vs 66.5 ± 3.5; 4.0% vs. 38.4%; χ2 = 74.429, P < 0.001; Fig. 3F). Age, smoking, sex, gross tumor type, tumor grade, and histological type did not significantly affect OS (all P > 0.05).
ARHGAP25 and RhoA expression levels, VM, and clinicopathological factors affect DFS
We observed a median survival time of 31.0 months and a mean disease-free survival time of 41.4 ± 2.9 months, and the five-year DFS rate was 19.9%, based on the Kaplan–Meier analysis of the five-year DFS time of the collected NSCLC cases. The DFS and DFS rates were significantly higher in ARHGAP25-positive samples (53.6 ± 4.1; 30.4%) than in ARHGAP25-negative samples (26.8 ± 2.6; 9.2%; χ2 = 20.012, P < 0.001; Fig. 4A). Similarly, the RhoA-positive group (28.8 ± 2.7; 10.5%) displayed significantly lower DFS and DFS rates than the RhoA-negative samples (59.3 ± 4.6; 35.1%; χ2 = 22.385, P < 0.001; Fig. 4B). The VM-positive samples (22.1 ± 2.8; 7.3%) had significantly lower DFS and DFS rates than the VM-negative samples (51.3 ± 3.6; 27.2%; χ2 = 25.721, P < 0.001; Fig. 4C). DFS and DFS rates were reduced in patients with tumor size greater than 3 cm (30.0 ± 2.6 vs. 54.5 ± 4.6; 12.9% vs. 30.8%; χ2 = 12.206, P < 0.001; Fig. 4D). The DFS and five-year DFS rates decreased in patients with TNM stages I, II, and III (72.2 ± 4.3 vs 38.2 ± 3.1 vs 17.6 ± 2.9; 50.6% vs. 7.9% vs. 2.4%; χ2 = 83.827, P < 0.001; Fig. 4E). Patients with LNM (20.2 ± 1.8; 4.3%) displayed considerably lower DFS and five-year DFS rates than those without LNM (63.6 ± 3.7; 37.1%; χ2 = 72.335, P < 0.001; Fig. 4F). Age, smoking, sex, general tumor type, tumor grade, and histological type did not significantly affect DFS (all P > 0.05).
Cox regression analysis
Using Cox regression analysis, we identified the parameters that influence the prognosis in patients with NSCLC. Univariate analysis suggested that the expression levels of ARHGAP25 and RhoA, VM, tumor size, TNM stage, and LNM were strongly linked with OS and DFS, which are critical factors influencing the future prognosis with NSCLC. Multivariate analysis confirmed that ARHGAP25, RhoA expression, VM, tumor size, TNM stage, and LNM were closely associated with OS and DFS and may be employed as independent prognostic markers of NSCLC (Tables 4 and 5).
ARHGAP25 is downregulated in NSCLC cell lines
The Cancer Genome Atlas (TCGA) database analysis (GEPIA) tool (http://gepia.cancer-pku.cn/) revealed that the level of ARHGAP25 was significantly lower in NSCLC (including 486 squamous lung carcinomas and 483 lung adenocarcinomas) tissues than in normal tissues (P < 0.001; Fig. 5A). To screen for gene expression profiles in NSCLC, we selected a qualified gene expression microarray dataset (TCGA) and determined that the expression level of ARHGAP25 was substantially reduced in NSCLC tissues (n = 1137) compared to normal tissues (n = 108) (Fig. 5B).
Downregulation of ARHGAP25 enhances NSCLC cell proliferation and migration
We first knocked down ARHGAP25 and used qRT-PCR and western blotting to confirm the results of siRNA knockdown in two NSCLC cell lines (A549 cells and H1299 cells). ARHGAP25 mRNA and protein levels were considerably lower in the knockdown cases than in the negative control without siRNA knockdown (Fig. 5C and D). The original images of the Western blotting can be found in Additional files 1, 2, 3, 4. These results indicated that the ARHGAP25 knockdown was successful. Subsequently, cell colony formation (Fig. 5E) and cell scratching assays (Fig. 5F) showed that lung cancer cells produced a higher number of colonies and a remarkable increase in the number of migrated cells after ARHGAP25 expression was reduced. The results of actin cytoskeleton staining supported the enhanced cytoskeleton and increased motility observed in the knockdown group (Fig. 5G). In summary, the ability of A549 and H1299 cell lines to proliferate and migrate was substantially enhanced when ARHGAP25 expression was reduced.
Depending on the TNM stage, histology, genetic modifications, and disease progression, treatments for NSCLC typically involve surgery, radiotherapy, chemotherapy, immunotherapy, and targeted molecular therapy, which can be administered alone or in combination . For patients with early-stage NSCLC (stage I, II, and IIIA [N2 lymph node involvement usually found during surgery]), surgical resection is recommended. For stage II-IIIA disease, platinum-based adjuvant chemotherapy is recommended, with a 5.4% lower risk of death at five years, but with a high recurrence rate and relatively high toxicity . Therefore, this study investigated the expression levels of ARHGAP25 and RHOA in NSCLC, as well as their relationship with VM, hoping to find more accurate and effective immune and molecular targeted therapies.
ARHGAP25 is localized on human chromosome 2p13 and is involved in the regulation of Rho family GTPases . Rho GTPases are involved in cytoskeletal dynamics, cell cycle progression, transcriptional control, cell survival, and vesicle transport, all of which can influence cancer growth . Some Rho GTPases promote cell cycle progression and gene transcription, which may explain their carcinogenic features, such as their ability to facilitate Ras-induced transformation . Angiogenesis must be induced for a tumor to grow beyond a certain size and for malignant cells to release substances that encourage angiogenesis of neighboring pre-existing blood vessels. To promote neovascularization, certain Rho GTPases regulate the release of pro-angiogenic molecules . RhoA is one of the most characteristic and universally highly expressed types of 20 Rho GTPs . Rho GTPases have been demonstrated to influence diseased cell invasion and metastasis by controlling cytoskeletal contraction and cell membrane protrusion . Negative regulators of Rac/Rho-like GTPases, GTPase-activating proteins (GAPs), reduce Rho GTPase activity by boosting the hydrolytic ability of Rho GTPases to convert activated GTP-binding status to inactivated GDP-binding status [29, 30].
ARHGAP25 is a Rac-specific GAP that is primarily expressed in hematopoietic cells. The invasion capability of alveolar rhabdomyosarcoma cells is controlled by the RhoE/ROCK/ARHGAP25 signaling pathway . Furthermore, researchers discovered that ARHGAP25 was downregulated in colorectal cancer (CRC) and that upregulating ARHGAP25 decreased CRC metastasis both in vivo and externally . One study revealed that the abnormal expression of ARHGAP25 reduces lung cancer cell proliferation and migration . ARHGAP25 expression was found to be lower in NSCLC tissues (48.5%) than in nearby normal tissues (63.1%), and was correlated with larger tumor size (P = 0.015), LNM (P = 0.015), and later clinical stage (P = 0.013) in this study. Based on data from the TCGA database, we found that the expression level of ARHGAP25 was decreased in NSCLC compared to normal tissues, which is consistent with immunohistochemical results. Cell colony formation and wound healing assays, as well as actin cytoskeleton staining assays, supported that proliferation and migration were significantly enhanced in ARHGAP25 knockdown NSCLC cells. According to Kaplan–Meier analysis and Cox regression analyses, patients with positive ARHGAP25 expression displayed longer OS than those with negative ARHGAP25 expression. Previous studies have found that ARHGAP25 expression is an independent predictor of NSCLC prognosis. These findings imply that ARHGAP25 may function as a tumor suppressor, thereby slowing tumor growth.
RhoA is a small guanosine triphosphatase (GTPase) that belongs to the Ras superfamily and functions in cytoskeletal reorganization . RhoA activity is controlled by Rho-associated coiled-coil protein kinase (ROCK), which phosphorylates its target protein  and disrupts certain biological functions, such as cell migration, adhesion, proliferation, contraction, and death. The activation of RhoA has been proven to increase the invasion and metastasis of several malignancies [33,34,35]. ROCK has been implicated in the production of VMs in liver cancer cell lines in several investigations [36, 37], and ROCK activation may be important in tumor cell VM regulation. The process by which tumor cells generate highly patterned vascular channels by relocating the F-actin cytoskeleton is characterized as VM, suggesting that RhoA might play a role in progression [21, 38]. In this study, RhoA expression in NSCLC was examined using immunohistochemistry. RhoA positivity was found in 62.3% of NSCLC tissues and 18.5% of normal tissues. RhoA overexpression was associated with LNM (P = 0.001) and later TNM stage (P = 0.013), implying that RhoA is relevant in NSCLC development. RhoA overexpression was associated with shorter OS and DFS, according to Kaplan–Meier survival and Cox regression analyses. Notably, RhoA expression and its impact on prognosis are comparable in NSCLC and other cancers. The expression level of RhoA was positively correlated with VM positivity and negatively correlated with ARHGAP25, according to the correlation analysis, which is consistent with earlier findings.
Differently from the classical model of tumor angiogenesis, VM formation is not dependent on endothelial cells . Tumor cells mimic normal endothelial cells and form tubular VMs containing both erythrocytes and tumor cells. VM formation consists mainly of deformation of tumor cells, remodeling of the extracellular matrix, and vascular-like structures connected to existing blood vessels . Studies have shown that the hypoxic microenvironment is closely associated with the development of VM, and that tumor cells form new blood vessels to obtain the required oxygen and nutrients . Hypoxia promotes the differentiation of cancer stem cells to form endothelial-like structures, and epithelial-mesenchymal transition (EMT) promotes VM formation by reducing the expression of cell adhesion molecules, such as E-cadherin [40, 41]. Previously published studies showed 28.5%, 33.3%, 34%, 22.7%, and 19.2% positive rates of VM in gastric cancer, esophageal mesenchymal tumor, cutaneous melanoma, osteosarcoma, and colorectal cancer, respectively [42,43,44,45,46]. VM occurred in less than 50% of the cancer types studied. In this study, the VM positivity rate was 36.9%, which was linked to larger tumor size, poor differentiation, LNM, and late clinical stage, and the OS and DFS of VM-positive patients were also shorter. This demonstrates that tumor vascularization aided NSCLC invasion and metastasis, resulting in a poor prognosis, which was in line with previous research findings.
In this study, we analyzed the relationship between ARHGAP25, RhoA, and VM in NSCLC and initially explored the effects of ARHGAP25 and RhoA on vascularization, hoping to determine their value in predicting NSCLC metastasis, prognosis, and targeted therapy. Furthermore, through correlation analysis, we observed that VM positivity is negatively correlated with the expression of ARHGAP25, as opposed to RhoA. However, the specific mechanism of ARHGAP25 in VM remains unclear. The RhoA/ROCK signaling pathway has been implicated in the creation of VM in hepatocellular carcinoma, melanoma, carcinoma, and osteosarcoma [36, 38, 47, 48]; therefore, we inferred that RhoA has the same mechanism in NSCLC and that ARHGAP25, as a regulatory protein for RhoA, should also promote VM formation through negative regulation of RhoA. Due to the limitations of the conditions and the relatively simple experimental methods, the specific mechanism of ARHGAP25 and RhoA in promoting VM formation is still unclear and requires further study in future work.
In conclusion, we showed that low ARHGAP25 expression and high RhoA expression are associated with VM and poor prognosis in patients with NSCLC. We believe that ARHGAP25 and RhoA may be used as novel prognostic biomarkers and therapeutic targets.
Availability of data and materials
All data are available from the corresponding author upon reasonable request.
Non-small cell lung cancer
Rho GTPase-activating protein 25
Ras homolog family member A
Lymph node metastasis
Lung squamous cell carcinoma
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Oser MG, Niederst MJ, Sequist LV, Engelman JA. Transformation from non-small-cell lung cancer to small-cell lung cancer: molecular drivers and cells of origin. Lancet Oncol. 2015;16(4):e165–72.
Yuan M, Huang LL, Chen JH, Wu J, Xu Q. The emerging treatment landscape of targeted therapy in non-small-cell lung cancer. Signal Transduct Target Ther. 2019;4:61.
Ayala-Domínguez L, Olmedo-Nieva L, Muñoz-Bello JO, Contreras-Paredes A, Manzo-Merino J, Martínez-Ramírez I, et al. Mechanisms of vasculogenic mimicry in ovarian cancer. Front Oncol. 2019;9:998.
Xie J, Kong X, Wang W, Li Y, Lin M, Li H, et al. Vasculogenic mimicry formation predicts tumor progression in oligodendroglioma. Pathol Oncol Res. 2021;27:1609844.
Andonegui-Elguera MA, Alfaro-Mora Y, Cáceres-Gutiérrez R, Caro-Sánchez CHS, Herrera LA, Díaz-Chávez J. An overview of vasculogenic mimicry in breast cancer. Front Oncol. 2020;10:220.
Zeeshan R, Mutahir Z. Cancer metastasis—tricks of the trade. Bosn J Basic Med Sci. 2017;17(3):172–82.
Maniotis AJ, Folberg R, Hess A, Seftor EA, Gardner LM, Pe’er J, et al. Vascular channel formation by human melanoma cells in vivo and in vitro: vasculogenic mimicry. Am J Pathol. 1999;155(3):739–52.
Viallard C, Larrivée B. Tumor angiogenesis and vascular normalization: alternative therapeutic targets. Angiogenesis. 2017;20(4):409–26.
Jaffe AB, Hall A. Rho GTPases: biochemistry and biology. Annu Rev Cell Dev Biol. 2005;21:247–69.
Heasman SJ, Ridley AJ. Mammalian Rho GTPases: new insights into their functions from in vivo studies. Nat Rev Mol Cell Biol. 2008;9(9):690–701.
Thuault S, Comunale F, Hasna J, Fortier M, Planchon D, Elarouci N, et al. The RhoE/ROCK/ARHGAP25 signaling pathway controls cell invasion by inhibition of Rac activity. Mol Biol Cell. 2016;27(17):2653–61.
Tao L, Zhu Y, Gu Y, Zheng J, Yang J. ARHGAP25: a negative regulator of colorectal cancer (CRC) metastasis via the Wnt/beta-catenin pathway. Eur J Pharmacol. 2019;858: 172476.
Xu K, Liu B, Ma Y. The tumor suppressive roles of ARHGAP25 in lung cancer cells. Onco Targets Ther. 2019;12:6699–710.
Huang WK, Chen Y, Su H, Chen TY, Gao J, Liu Y, et al. ARHGAP25 inhibits pancreatic adenocarcinoma growth by suppressing glycolysis via AKT/mTOR pathway. Int J Biol Sci. 2021;17(7):1808–20.
Zhang Y, Lin Y, Zhu Y, Zhang X, Tao L, Yang M. ARHGAP25 expression in colorectal cancer as a biomarker associated with favorable prognosis. Mol Clin Oncol. 2022;16(4):84.
Aspenström P, Ruusala A, Pacholsky D. Taking Rho GTPases to the next level: the cellular functions of atypical Rho GTPases. Exp Cell Res. 2007;313(17):3673–9.
Zhong Z, Grasso L, Sibilla C, Stevens TJ, Barry N, Bertolotti A. Prion-like protein aggregates exploit the RHO GTPase to cofilin-1 signaling pathway to enter cells. EMBO J. 2018. https://doi.org/10.15252/embj.201797822.
Zhang HY, Li HM, Yu Z, Yu XY, Guo K. Expression and significance of squalene epoxidase in squamous lung cancerous tissues and pericarcinoma tissues. Thorac Cancer. 2014;5(4):275–80.
He X, You J, Ding H, Zhang Z, Cui L, Shen X, et al. Vasculogenic mimicry, a negative indicator for progression free survival of lung adenocarcinoma irrespective of first line treatment and epithelial growth factor receptor mutation status. BMC Cancer. 2021;21(1):132.
Xia Y, Cai XY, Fan JQ, Zhang LL, Ren JH, Li ZY, et al. The role of sema4D in vasculogenic mimicry formation in non-small cell lung cancer and the underlying mechanisms. Int J Cancer. 2019;144(9):2227–38.
Li Q, Lai Q, He C, Fang Y, Yan Q, Zhang Y, et al. RUNX1 promotes tumour metastasis by activating the Wnt/beta-catenin signalling pathway and EMT in colorectal cancer. J Exp Clin Cancer Res. 2019;38(1):334.
Alexander M, Kim SY, Cheng H. Update 2020: management of non-small cell lung cancer. Lung. 2020;198(6):897–907.
Pignon JP, Tribodet H, Scagliotti GV, Douillard JY, Shepherd FA, Stephens RJ, et al. Lung adjuvant cisplatin evaluation: a pooled analysis by the LACE collaborative group. J Clin Oncol. 2008;26(21):3552–9.
Katoh M, Katoh M. Identification and characterization of ARHGAP24 and ARHGAP25 genes in silico. Int J Mol Med. 2004;14(2):333–8.
Benitah SA, Valeron PF, van Aelst L, Marshall CJ, Lacal JC. Rho GTPases in human cancer: an unresolved link to upstream and downstream transcriptional regulation. Biochim Biophys Acta. 2004;1705(2):121–32.
Aspenstrom P. Fast-cycling Rho GTPases. Small GTPases. 2020;11(4):248–55.
Karlsson R, Pedersen ED, Wang Z, Brakebusch C. Rho GTPase function in tumorigenesis. Biochim Biophys Acta. 2009;1796(2):91–8.
Haga RB, Ridley AJ. Rho GTPases: regulation and roles in cancer cell biology. Small GTPases. 2016;7(4):207–21.
Zhang L, Zhou A, Zhu S, Min L, Liu S, Li P, et al. The role of GTPase-activating protein ARHGAP26 in human cancers. Mol Cell Biochem. 2022;477(1):319–26.
Liu J, Wada Y, Katsura M, Tozawa H, Erwin N, Kapron CM, et al. Rho-associated coiled-coil kinase (ROCK) in molecular regulation of angiogenesis. Theranostics. 2018;8(21):6053–69.
Narumiya S, Thumkeo D. Rho signaling research: history, current status and future directions. FEBS Lett. 2018;592(11):1763–76.
Privat M, Cavard A, Zekri Y, Ponelle-Chachuat F, Molnar I, Sonnier N, et al. A high expression ratio of RhoA/RhoB is associated with the migratory and invasive properties of basal-like breast tumors. Int J Med Sci. 2020;17(17):2799–808.
Mei C, Liu C, Gao Y, Dai WT, Zhang W, Li X, et al. eIF3a regulates colorectal cancer metastasis via translational activation of RhoA and Cdc42. Front Cell Dev Biol. 2022;10: 794329.
Liu D, Mei X, Wang L, Yang X. RhoA inhibits apoptosis and increases proliferation of cultured SPCA1 lung cancer cells. Mol Med Rep. 2017;15(6):3963–8.
Zhang JG, Zhang DD, Wu X, Wang YZ, Gu SY, Zhu GH, et al. Incarvine C suppresses proliferation and vasculogenic mimicry of hepatocellular carcinoma cells via targeting ROCK inhibition. BMC Cancer. 2015;15:814.
Zhang JG, Zhou HM, Zhang X, Mu W, Hu JN, Liu GL, et al. Hypoxic induction of vasculogenic mimicry in hepatocellular carcinoma: role of HIF-1 α, RhoA/ROCK and Rac1/PAK signaling. BMC Cancer. 2020;20(1):32.
Xia Y, Cai XY, Fan JQ, Zhang LL, Ren JH, Chen J, et al. Rho kinase inhibitor fasudil suppresses the vasculogenic mimicry of B16 mouse melanoma cells both in vitro and in vivo. Mol Cancer Ther. 2015;14(7):1582–90.
Seftor RE, Hess AR, Seftor EA, Kirschmann DA, Hardy KM, Margaryan NV, et al. Tumor cell vasculogenic mimicry: from controversy to therapeutic promise. Am J Pathol. 2012;181(4):1115–25.
Luo Q, Wang J, Zhao W, Peng Z, Liu X, Li B, et al. Vasculogenic mimicry in carcinogenesis and clinical applications. J Hematol Oncol. 2020;13(1):19.
Wei X, Chen Y, Jiang X, Peng M, Liu Y, Mo Y, et al. Mechanisms of vasculogenic mimicry in hypoxic tumor microenvironments. Mol Cancer. 2021;20(1):7.
Zhao H, Gu XM. Study on vasculogenic mimicry in malignant esophageal stromal tumors. World J Gastroenterol. 2008;14(15):2430–3.
Liu Z, Sun B, Qi L, Li H, Gao J, Leng X. Zinc finger E-box binding homeobox 1 promotes vasculogenic mimicry in colorectal cancer through induction of epithelial-to-mesenchymal transition. Cancer Sci. 2012;103(4):813–20.
van Beurden A, Schmitz RF, van Dijk CM, Baeten CI. Periodic acid Schiff loops and blood lakes associated with metastasis in cutaneous melanoma. Melanoma Res. 2012;22(6):424–9.
Ren K, Yao N, Wang G, Tian L, Ma J, Shi X, et al. Vasculogenic mimicry: a new prognostic sign of human osteosarcoma. Hum Pathol. 2014;45(10):2120–9.
Sun J, Sun B, Sun R, Zhu D, Zhao X, Zhang Y, et al. HMGA2 promotes vasculogenic mimicry and tumor aggressiveness by upregulating Twist1 in gastric carcinoma. Sci Rep. 2017;7(1):2229.
Xia Y, Cai X, Fan J, Zhang L, Li Z, Ren J, et al. RhoA/ROCK pathway inhibition by fasudil suppresses the vasculogenic mimicry of U2OS osteosarcoma cells in vitro. Anticancer Drugs. 2017;28(5):514–21.
Zhang JG, Li XY, Wang YZ, Zhang QD, Gu SY, Wu X, et al. ROCK is involved in vasculogenic mimicry formation in hepatocellular carcinoma cell line. PLoS ONE. 2014;9(9): e107661.
The Key Project of the Anhui Provincial Education Department (no. KJ2021A0725), and the Postgraduate Innovation Program of Bengbu Medical College (No. Byycx21091) provided funding for this study.
The authors declare no competing interests.
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Bengbu Medical College (No. 2020KY035). All patients provided written informed consent. The study adhered to the ethical guidelines of the Declaration of Helsinki.
Consent for publication
The authors declare that there are no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shi, F., Wu, J., Jia, Q. et al. Relationship between the expression of ARHGAP25 and RhoA in non-small cell lung cancer and vasculogenic mimicry. BMC Pulm Med 22, 377 (2022). https://doi.org/10.1186/s12890-022-02179-5
- Rho GTPase
- Vasculogenic mimicry