Impact Factor 4.848 | CiteScore 3.5
More on impact ›

Original Research ARTICLE

Front. Oncol., 09 September 2020 | https://doi.org/10.3389/fonc.2020.543562

VEGFR2 Promotes Metastasis and PD-L2 Expression of Human Osteosarcoma Cells by Activating the STAT3 and RhoA-ROCK-LIMK2 Pathways

Bingxin Zheng1, Chuanli Zhou2, Guojian Qu3, Chongmin Ren1, Peng Yan1, Wei Guo4* and Bin Yue1*
  • 1Department of Orthopedic Oncology, The Affiliated Hospital of Qingdao University, Qingdao, China
  • 2Department of Spinal Surgery, The Affiliated Hospital of Qingdao University, Qingdao, China
  • 3Department of General Surgery(adult), Qingdao Women and Children’s Hospital, Qingdao, China
  • 4Musculoskeletal Tumor Center, Peking University People’s Hospital, Beijing, China

The survival rate of osteosarcoma, the most prevalent primary bone tumor, has not been effectively improved in the last 30 years. Hence, new treatments and drugs are urgently needed. Antiangiogenic therapy and immunotherapy have good antitumor effects in many kinds of tumors. It is hypothesized that there may be a synergistic effect between immune checkpoint inhibitors and antiangiogenic therapy. Nevertheless, its potential mechanism is still unclear. Vascular endothelial growth factor receptor-2 (VEGFR2) expression was detected by immunohistochemistry in 18 paired osteosarcoma tissues. Moreover, we investigated the effects of apatinib treatment and VEGFR2 knockdown on osteosarcoma as well as the relevant underlying mechanism. Immunohistochemistry assays showed that, compared with that in primary osteosarcoma, VEGFR2 expression was higher in lung metastases. VEGFR2 was positively correlated with PD-L2 expression in osteosarcoma lung metastasis. Transwell assays indicated that VEGFR2 inhibition reduced osteosarcoma cell metastatic abilities in vitro. We also demonstrated that VEGFR2 inhibition downregulated the STAT3 and RhoA-ROCK-LIMK2 pathways, thereby attenuating migration and invasion. Additionally, VEGFR2 inhibition targeted STAT3, through which it reduced PD-L2 expression in osteosarcoma cells. VEGFR2 inhibition markedly attenuated osteosarcoma lung metastatic ability in vivo. In this study, we presented the pro-metastatic functional mechanism of VEGFR2 in osteosarcoma. VEGFR2 inhibition exhibits antitumor effects through antiangiogenic effects and inhibition of immune escape, which possibly provides potential clinical treatment for metastatic osteosarcoma.

Introduction

As a primary malignant bone tumor that is commonly seen in children and teenagers, osteosarcoma exhibits an obvious tendency toward locally aggressive and pulmonary metastasis (13). Doctors mainly rely on integrating neoadjuvant chemotherapy with surgery to treat osteosarcoma. Nevertheless, chemotherapy for osteosarcoma reaches a plateau period, and the 5-year survival of pulmonary metastases patients is only 28% (4). Thus, it is urgent to develop novel and effective treatments for osteosarcoma, particularly for patients with metastatic osteosarcoma.

Both programmed death ligand-1 (referred to as PD-L1) and programmed death ligand-2 (referred to as PD-L2) are the ligands of programmed death-1 (PD-1), and their combination leads to T-cell exhaustion and immune escape. Recently, immune checkpoint inhibitors (ICIs) have achieved good curative effects in varying solid tumors. However, most studies mainly center on the role of PD-1 ligands in the reciprocity between T-cells and tumor cells, and few studies have been performed on the intrinsic effects of PD-1 ligands in tumor cells, which may be associated with ICI treatment effects (5, 6). According to some studies, PD-L1/PD-L2 are associated with cellular biological behaviors such as autophagy, metastasis, propagation, and epithelial-mesenchymal transition (EMT) (711).

More recently, VEGFR2 inhibition with apatinib has shown remarkable effects in the treatment of multiple tumors (1215). Moreover, some researchers have hypothesized that the integration of antiangiogenic therapy with ICI immunotherapies is likely to produce a collaborative result (16). However, the deeper and specific underlying mechanism of VEGFR2 affecting the metastasis and tumor immunity of osteosarcoma remain poorly understood.

We mainly addressed three aspects in this study: (1) VEGFR2 expression in primary and metastatic osteosarcoma; (2) significance of VEGFR2 inhibition to migration and invasion of osteosarcoma cells both in vivo and in vitro; and (3) possible mechanism of osteosarcoma metastasis mediated by VEGFR2.

Materials and Methods

Osteosarcoma Cell Samples and Patient Data

Primary osteosarcoma specimens fixed with formalin and embedded in paraffin and paired pulmonary metastasis specimens were acquired from the Musculoskeletal Tumor Center at Peking University People’s Hospital (Beijing, China) with relevant informed patient consent. To carry out this study, we also obtained approval from the Ethics Committees of The Affiliated Hospital of Qingdao University and Peking University People’s Hospital.

Cell Growth and Reagents

KHOS and U2OS cells were obtained from American Type Culture Collection (ATCC), and we used short tandem repeat (STR) analysis for the authentication of the two cells at Beijing Microread Genetics Co., Ltd. Cells were cultured in RPMI 1640 medium (HyClone) supplemented with 1% penicillin and streptomycin (Invitrogen) and 10% fetal bovine serum (Gibco) in an environment with a temperature of 37°C and 5% CO2.

The antibodies we used in relevant experiments were as follows: antibodies against p-STAT3, STAT3, p-cofilin, cofilin and VEGFR2, which we purchased from Cell Signaling Technology in United States; anti-PDL2, p-LIMK2, LIMK2 and the RhoA activation assay, which we obtained from Abcam in United States; and anti-GAPDH, which we purchased from Santa Cruz Biotechnology. Y-27632 2HCl was purchased from Selleck.

Gene Knockdown With siRNA/shRNA and Ectopic Expression

The lentiviruses targeting VEGFR2 and control shNC were obtained from GenePharma (Suzhou, China), and two pairs of shVEGFR2 primers were designed: #1, sense strand 5′-CTCG GTCATTTATGTCTAT-3′ and antisense strand 5′-ATAGACA TAAATGACCGAG-3′ and #2, sense strand 5′-CTTCGAAGCA TCAGCATAA-3′ and antisense strand 5′-TTATGCTGATG CTTCGAAG-3′. After infection with lentivirus, puromycin (2 μg/ml) was used to select cells stably expressing shVEGFR2 and shNC.

The source of siSTAT3 and control siRNA was GenePharma (Suzhou, China), and the detailed sequences targeting STAT3 were as follows: #1, 5′-AAAGAAUCACAUGCCACU UTT-3′ (sense) and 5′-AAGUGGCAUGUGAUUCUUUGC-3′ (antisense); #2, 5′-ACAAUCUACGAAGAAUCAATT-3′ (sense) and 5′-UUGAUUCUUCGUAGAUUGUGC-3′ (antisense); and #3, 5′-CGUCCAGUUCACUACUAAATT-3′ (sense) and 5′-UUUAGUAGUGAACUGGACGCC-3′ (antisense). SiRNAs were transfected into tumor cells by using Lipofectamine 3000 obtained from Invitrogen, United States, as per the protocols of the manufacturer. Further analyses were carried out on the cells collected after 48 h of transfection.

The plasmid containing VEGFR2 cDNA or negative control was transfected into KHOS cells with Lipofectamine 3000 (Invitrogen). The medium was replaced after 24 h incubation, and then the cells were treated with Y-27632 2HCl.

Cell Counting Kit 8 Assay

After 2000 cells were seeded in each well of 96-well plates, the cells were treated with apatinib at various concentrations. The viability analysis was carried with CCK-8 (Dojindo Laboratories, Kumamoto, Japan) reagent as per the protocols of the manufacturer.

Cell Colony Formation Assay

Cells (1 × 103) in the logarithmic growth phase were coated into six-well plates. After forming colonies for 10 days, cells were fixed with methanol and stained with gentian violet staining solution.

Transwell Assay

A total of 6 × 104 cells were plated into the upper part of the Transwell insert (8 μm pore size; Corning) in the absence of coating for the migration assay or in the presence of a Matrigel coat (BD Bioscience, 354234) for the invasion assay. After 24 h of culture, the cells were fixed with methanol and stained with 0.1% gentian violet staining solution. The quantities of migrated cells were calculated with a microscope in 5 fields per well.

Quantitative RT-PCR

Total RNA was obtained from KHOS and U2OS cells with TRIzol (Invitrogen). Then, complementary DNA (cDNA) was generated with isolated RNA and oligo dT primers by using SuperScript III First-Strand Synthesis SuperMix (Invitrogen). Real-time quantitative PCR was conducted with the SYBR Green PCR Master Mix from Applied Biosystems, United States and was run on Bio-Rad CFX96 following the manufacturer’s instructions.

Bioinformatics Assay

The osteosarcoma dataset from GEO was used for the data mining and bioinformatics studies (Nos. GSE21257 and GSE33382). Cluster analysis and heat map generation were performed by MeV software (version 4.9). The association of the genes of interest and confirmed gene sets was demonstrated by using gene set enrichment analysis (GSEA). REVIGO was used for the gene annotation network analysis.

Western Blotting and GTPase Assay

Western blotting was conducted as described previously (12). Briefly, the same volume of protein obtained from different cell lysates was loaded onto 7.5–15% SDS-PAGE gels and then transferred onto PVDF membranes. The membranes were cultured with suitable antibodies overnight at 4°C after blocking in non-fat milk for 1 h. The western blot bands were observed with a Bio-Rad detection system (Hercules, CA, United States). The GTPase assay was performed as per the manufacturer’s instructions.

Immunohistochemistry Assay

Immunohistochemistry (IHC) staining was performed as previously described (12). Briefly, paraffin slices were incubated with corresponding primary antibodies or with non-immune serum in PBS overnight at 4°C. Then, the sections were stained with secondary antibody (ZSGB-BIO, China) for 30 min at 37°C. The percentage of positive cells was defined as follows: 0: 0% positive; 1: less than 5% positive; 2: between 5 and 50% positive; and 3: more than 50% positive. In addition, the staining intensity was scored as 4 follows: 0: no staining; 1: poor; 2: medium; and 3: intense. The total score was calculated by multiplying the intensity result by the percentage result. More than 10 representative fields (x400 magnification) were scored between the primary osteosarcoma and lung metastasis specimens for IHC analyses. The immunostaining was scored by two pathologists who were blinded to the clinical features of the clinical specimens.

Immunofluorescence Assay

Cells (2 × 105) were plated onto coverslips in six-well plates with or without the corresponding inhibitor treatment. The coverslips were cultured with 100 nM rhodamine phalloidin (manufactured in the United States) at room temperature for 45 min. The sections were observed by a confocal microscope manufactured in Japan.

Tumor Xenografts

To confirm the effect of apatinib treatment on the metastatic ability of KHOS cells, 3 × 106 KHOS cells were injected intravenously into the tail vein of 42-day-old maternal BALB/c nude mice from Vitalriver, China. On the 5th day after injection, the mice were randomly separated into two group (n = 5 per group)and fed DMSO or apatinib 2 mg/kg/day for 1 month total.

To evaluate the effect of VEGFR2 knockdown on the metastatic ability of KHOS cells, 42-day-old maternal BALB/c nude mice were intravenously injected with 3 × 106 KHOS-shVEGFR2 or -shNC cells (n = 5 per group) via the tail vein. After 30 days, all mice were sacrificed.

At the end of the experiment, the lungs were collected, imaged, stained with haematoxylin-eosin dye, and the quantity of pulmonary metastases was confirmed. In addition, animal care and processing procedures were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.

Statistical Analysis

The statistical analysis was conducted by SPSS software (version 21.0). Statistical evaluation was conducted by Student’s t-test. Data are shown as the mean ± S.D. A p-value below 0.05 showed that there were significant differences.

Results

VEGFR2 Expression Is Elevated in Osteosarcoma Lung Metastases

An IHC assay for VEGFR2 was conducted on paired primary osteosarcoma specimens and pulmonary transfer specimens. VEGFR2 exhibited cytoplasmic and nuclear expression (Figure 1A), and VEGFR2 expression was increased in pulmonary metastases compared with local osteosarcoma samples (Figure 1B), implying that VEGFR2 may play a significant role in osteosarcoma metastasis.

FIGURE 1
www.frontiersin.org

Figure 1. Elevated VEGFR2 expression in osteosarcoma lung metastasis. (A) VEGFR2 and PD-L2 expression in 18 pairs of primary osteosarcoma specimens and pulmonary metastasis specimens was investigated via immunohistochemistry. Representative images are presented (magnification at 200×). (B) IHC scores of VEGFR2 and PD-L2 were analyzed between the primary osteosarcoma group and the pulmonary metastasis group. (C) The correlation between PD-L2 and VEGFR2 expression levels in 18 osteosarcoma pulmonary metastases was investigated. (D) VEGFR2 was overexpressed in osteosarcoma in comparison to normal control. Data are shown as the mean ± S.D. *p < 0.05, **P < 0.01.

PD-L2 expression was also investigated in these paired osteosarcoma specimens, and the results indicated that PD-L2 expression was also elevated in lung metastasis compared with primary osteosarcoma (Figure 1A). Thus, we investigated whether VEGFR2 and PD-L2 were related to each other in osteosarcoma lung metastasis, and an important positive relationship between VEGFR2 and PD-L2 expression was confirmed (Figure 1C). VEGFR2 was overexpressed in osteosarcoma in comparison to normal control (Figure 1D).

Hence, the results reveal that VEGFR2 may play a crucial role in metastasis and immune escape of osteosarcoma.

VEGFR2 Inhibition Suppresses the Migration and Invasion of Osteosarcoma Cells in vitro

The proliferation activity of KHOS and U2OS cells was not appreciably influenced by apatinib at concentrations of 0.5, 1.0, and 1.5 μM using the CCK-8 assay (Figure 2A). The 1.5 μM concentration was chosen as the upper limit for subsequent studies. There were few differences between the DMSO and apatinib groups in the cell colony formation analysis (Figure 2A). The two cell lines were transferred with lentivirus targeting VEGFR2, and VEGFR2 expression was later investigated using WB and real-time PCR assays (Figure 2B). Furthermore, we investigated the effect of VEGFR2 reduction on the migration and invasion of osteosarcoma cells using a Transwell assay.

FIGURE 2
www.frontiersin.org

Figure 2. VEGFR2 inhibition decreases osteosarcoma cell invasion and migration. (A) Apatinib did not obviously influence KHOS and U2OS cell viability at concentrations of 0.5, 1.0, and 1.5 μM. Proliferation of these two cell lines in response to apatinib treatment was confirmed via cell colony formation analysis. (B) Both shNC and shVEGFR2 lentiviruses were transferred into KHOS and U2OS cells, following estimation of VEGFR2 expression via WB and RT-PCR. (C,D) Both apatinib treatment and VEGFR2 knockdown significantly attenuated the KHOS and U2OS cell migratory and invasive abilities, as demonstrated by a Transwell assay. (Magnification at 100×), (mean ± S.D., n = 3), ***p < 0.001.

The results showed that apatinib reduced the abilities of KHOS and U2OS cells to migrate and invade (Figure 2C). Moreover, both the KHOS-shVEGFR2 and U2OS-shVEGFR2 groups showed lower migratory and invasive abilities than the control groups (Figure 2D).

In short, these results indicate that VEGFR2 is crucial for the migratory and invasive abilities of osteosarcoma cells, and its specific molecular mechanism needs further study.

Potential Mechanisms That Underlie VEGFR2 Associations With Metastasis

A review of the expression of VEGFR2-related genes and their relationships with overall survival and HUVOS grade in osteosarcoma specimens from the GEO dataset is shown in Figure 3A. To investigate the potential mechanisms of the relationship of VEGFR2 with metastasis and prognosis, functional enrichment and gene expression of VEGFR2-related genes were calculated and visualized. The heat map of the top differentially expressed genes between the VEGFR2 high-expression group and the VEGFR2 low-expression group is shown in Figure 3B.

FIGURE 3
www.frontiersin.org

Figure 3. Potential mechanisms that underlie the VEGFR2 associations with metastasis and survival time. (A) VEGFR2-related genes, HUVOS grade and overall survival were analyzed by heat map observation and cluster analysis. (B) Differentially expressed genes were analyzed by visual analysis of heat maps in the VEGFR2 high- and low-expression groups. (C,D) Differentially expressed VEGFR2-related genes were investigated by analysis of the gene annotation network and GSEA.

The GSEA results indicated that VEGFR2 was positively associated with STAT3 signaling pathway, activation of immune response pathways, cytoskeletal arrangement and relevant metastasis characteristics (Figure 3C). Gene annotation network analysis was conducted on the differentially expressed genes (Figure 3D: VEGFR2 high against VEGFR2 low).

These results show that VEGFR2 is related to immunosuppression, cytoskeletal arrangement and relevant metastasis pathways. Next, we investigated the specific mechanism by which VEGFR2 affects metastasis and immune escape in vitro.

VEGFR2 Inhibition Inactivates RhoA-ROCK-LIMK2 Signaling and Suppresses PD-L2 Expression in Osteosarcoma Cells by Targeting STAT3

As Figure 4A shows, apatinib treatment reduced PD-L2 expression in KHOS and U2OS cells. Combined with the above IHC results, these findings showed that VEGFR2 significantly affected PD-L2 expression in osteosarcoma. A previous study suggested that PD-L2 knockdown inhibited osteosarcoma pulmonary metastasis by the RhoA-ROCK-LIMK2 pathway (11); thus, we further investigated the effect of apatinib on the RhoA-ROCK-LIMK2 pathway. As Figure 4A shows, apatinib treatment also attenuated the activation of p-cofilin and p-LIMK2 in KHOS and U2OS cells.

FIGURE 4
www.frontiersin.org

Figure 4. VEGFR2 inhibition inactivates the RhoA-ROCK-LIMK2 pathway and suppresses PD-L2 expression in osteosarcoma cells by targeting STAT3. (A) Apatinib reduced the expression of PD-L2, STAT3, p-cofilin and p-LIMK, as demonstrated by western blotting. (B) VEGFR2 knockdown reduced the expression of PD-L2, STAT3, p-cofilin and p-LIMK, as indicated by western blotting. (C) STAT3 expression in KHOS cells was downregulated by three STAT3 siRNA sequences (top). STAT3 knockdown decreased the expression of PD-L2, p-LIMK and p-cofilin (bottom). (D) Overexpression of VEGFR2 increased Y-27632 2HCl-induced decline of p-LIMK and p-Cofilin by western blot. (E) A coexistence phenomenon was observed between RhoA and STAT3 by bioinformatics analysis (http://genemania.org/). (F) The change in RhoA activation in U2OS cells was confirmed by the GTPase test. ***P < 0.001.

Furthermore, we transfected shRNA targeting VEGFR2 into these two cell lines to confirm the effects of apatinib on the expression of PD-L2 and activation of RhoA-ROCK-LIMK2 signaling. As shown in Figure 4B, VEGFR2 knockdown attenuated the expression of PD-L2, p-LIMK2 and p-cofilin in KHOS and U2OS cells.

Previous studies showed that STAT3 could regulate PD-L1/PD-L2 expression (1720). VEGFR2 inhibition could induce inactivation of STAT3 in this study (Figures 4A,B). Therefore, we transfected siRNA targeting STAT3 into the two cell lines and tested the effects of STAT3 knockdown on the expression of PD-L2 and activation of RhoA-ROCK-LIMK2 signaling. As shown in Figure 4C, STAT3 knockdown decreased the expression of PD-L2, p-LIMK2 and p-cofilin in KHOS and U2OS cells. KHOS cells that were transfected by VEGFR2 overexpression plasmid increased Y-27632 2HCl-induced decline of p-LIMK and p-Cofilin (Figure 4D), indicating that VEGFR2 promoted migration and invasion by activating RhoA-ROCK-LIMK2 pathway. Additionally, bioinformatics analysis indicated that there may be a coexistence phenomenon between STAT3 and RhoA (Figure 4E). A reduction in RhoA activation in the U2OS-siSTAT3 group compared with the controls was observed via the GTPase assay (Figure 4F).

In conclusion, these data strongly suggest that VEGFR2 inhibition suppresses PD-L2 expression and the RhoA-ROCK-LIMK2 pathway by targeting STAT3.

VEGFR2 and STAT3 Inhibition Affects Cytoskeletal Rearrangement in Osteosarcoma Cells

The actin cytoskeleton is very important in cell motility, and LIMK and cofilin are key regulators in this process. Phosphorylation of cofilin suppresses its ability to adhere to actin filaments, which inhibits filament breakdown. In addition, cofilin phosphorylation is mediated by LIMK (2123). Since the phosphorylation of cofilin and LIMK2 was influenced by VEGFR2 or STAT3 expression, as mentioned above, we hypothesized that some morphological alterations occurred when VEGFR2 or STAT3 expression was inhibited.

As Figure 5A shows, obvious lamellipodial prominences were observed near the submembranous part in the control group, while the opposite result was observed in the apatinib treatment group, which indicated rearranged cytoskeleton and well-arranged F-actin in cells. In the same manner, compared with those in the KHOS-shVEGFR2 group, the lamellipodial prominences and colored F-actin filaments gathered in the margin of the KHOS-shNC group (Figure 5B). Additionally, F-actin was mostly arranged in KHOS-siSTAT3 cell cytoplasm; however, F-actin and lamellipodial prominences were concentrated near the edge in the KHOS-siNC group (Figure 5C).

FIGURE 5
www.frontiersin.org

Figure 5. VEGFR2 and STAT3 expression influences osteosarcoma cell cytoskeletal arrangement. The effects of (A) apatinib treatment, (B) VEGFR2 knockdown and (C) STAT3 knockdown on cytoskeletal rearrangement were observed via confocal microscopy in KHOS cells. Representative images are shown. Cell nuclei were stained with DAPI. Scale bar represents 10 μm.

These results suggest that the inhibition of VEGFR2 or STAT3 expression has a distinct effect on the cytoskeletal rearrangement of osteosarcoma cells.

VEGFR2 Inhibition Attenuates Osteosarcoma Cell Metastasis in vivo

To confirm the effects of apatinib treatment or VEGFR2 knockdown on the metastatic ability of KHOS cells, corresponding KHOS cell lines were injected intravenously via the tail vein of BALB/c nude mice as described in the tumor xenograft section.

The apatinib treatment group showed distinctly fewer pulmonary metastases than the control group (Figure 6A). Similarly, lung metastases were more often observed in the control group than in the shVEGFR2 group (Figure 6B). Representative H&E images are shown (Figures 6A,B).

FIGURE 6
www.frontiersin.org

Figure 6. VEGFR2 inhibition attenuates osteosarcoma cell metastasis in vivo. (A,B) The quantities of pulmonary metastases in the apatinib and DMSO groups or shNC and shVEGFR2 groups are presented. Representative images are shown. (C,D) Local primary osteosarcoma was further estimated by immunoblotting and IHC (magnification at 400×) for PD-L2, p-STAT3 and p-LIMK expression. (E) A schematic drawing of VEGFR2 in osteosarcoma metastasis. (mean ± S.D., n = 4) **p < 0.01.

We previously found that there was no significant difference in the tumor size of primary osteosarcoma between the apatinib and DMSO groups (12), and here the primary tumors were further assessed via IHC and WB assays. Consistent with the above results, decreased p-STAT3, PD-L2 and p-LIMK2 expression was observed in the apatinib group compared with the DMSO group (Figures 6C,D).

In conclusion, VEGFR2 inhibition attenuates human osteosarcoma cell metastasis.

Discussion

Osteosarcoma is a local invasive malignancy, and the incidence of pulmonary metastasis is high. Metastasis is the dominant element influencing the treatment effect in osteosarcoma patients, which is a very complex course involving various molecular biological mechanisms.

In this study, we demonstrate that VEGFR2 expression is higher in pulmonary metastases than in primary specimens using IHC. We also investigated PD-L2 expression in these paired osteosarcoma specimens, and the results indicated that PD-L2 expression is also enhanced in pulmonary metastasis compared to the primary specimen. Moreover, a significant positive relationship between PD-L2 and VEGFR2 expression was confirmed. In addition to its external immune escape effects of tumor cells, PD-L2 also has a few internal effects in tumor cell processes, including EMT, metastasis and autophagy (10, 11). These results indicate that both VEGFR2 and PD-L2 may play important roles in the pulmonary metastasis of osteosarcoma cells. In this study, we show that both apatinib treatment and VEGFR2 knockdown inhibit KHOS and U2OS cell migration and invasion. In the future, we will continue investigating their underlying mechanisms.

To determine the potential elements that support VEGFR2 associations with metastasis, the mRNA expression of 10 VEGFR2-relevant genes in osteosarcoma specimens was measured and observed via bioinformatics analysis. Gene expression methods, GSEA and gene annotation network analysis of VEGFR2-relevant genes were performed. For angiopoiesis, the results show that VEGFR2 may be related to the immune response, cytoskeletal rearrangement and relevant metastasis pathways.

To verify the bioinformatics analysis results, we investigated the effect of VEGFR2 inhibition on the signaling pathway in osteosarcoma. Apatinib (1.5 μM) and shRNA targeting VEGFR2 were used. Similar to the previous results, both apatinib treatment and VEGFR2 knockdown attenuated PD-L2 expression in these two cell lines, consistent with the IHC results. These data indicate that VEGFR2 actively regulates PD-L2 expression in osteosarcoma.

Previous studies have demonstrated that the cell actin cytoskeleton plays a considerable role in adhesion and motility, and cofilin and LIMK are critical regulators (2123). A previous study indicated that PD-L2 knockdown inhibits the ability of osteosarcoma cells to metastasize by the RhoA-ROCK-LIMK2 signaling pathway (11); thus, we further investigated the effect of VEGFR2 inhibition on RhoA-ROCK-LIMK2 signaling. As shown in Figure 4B, both apatinib treatment and VEGFR2 knockdown decreased p-cofilin and p-LIMK2 expression in KHOS and U2OS cells.

Previous studies suggest that STAT3 could regulate PD-1 ligand expression (1720). In previous research, we demonstrated that apatinib inhibits PD-L1 presentation by targeting STAT3 in osteosarcoma. Recently, with the promising curative effect of ICIs, it has been proposed that ICI therapy together with antiangiogenic therapy may show a collaborative impact (16). Nevertheless, the potential mechanism remains unclear. In this study, we reveal that STAT3 knockdown attenuates the expression of PD-L2, p-LIMK2 and p-cofilin in osteosarcoma, suggesting that VEGFR2 inhibition suppresses migration, invasion and PD-L2 expression by targeting the STAT3 and RhoA-ROCK-LIMK2 pathways.

It is well known that angiogenesis is related to tumor metastasis (2426), but its relationship with immune escape remains poorly understood. At present, there are few studies investigating the connection between PD-L2 and VEGFR2. In this sense, we take the lead in revealing, through this study, that VEGFR2 inhibition decreases PD-L2 expression in osteosarcoma cells. Thus, through the tumor cell extrinsic and intrinsic functions of PD-L2, VEGFR2 can affect both osteosarcoma cell metastasis and immune escape through the deactivation of STAT3 and the RhoA-ROCK-LIMK2 pathway (Figure 6E).

Overall, the pro-metastatic functional mechanism of VEGFR2 in osteosarcoma is presented herein. VEGFR2 inhibition targeted the RhoA-ROCK-LIMK2 pathway and STAT3 both in vitro and in vivo, through which migration, invasion and PD-L2 expression were reduced. It is important to verify the efficacy and mechanism of antiangiogenesis therapy combined with immunotherapy in a humanized mouse model, which may be one of the shortcomings of our study, and we will verify the above results in our further study. Based on the above, our understanding of the regulatory characteristics of angiogenesis in the immune escape of osteosarcoma cells was improved, which offers a feasible strategy for treating metastatic osteosarcoma.

Data Availability Statement

Publicly available datasets were analyzed in this study, these can be found in the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) (GSE21257 and GSE33382).

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committees of The Affiliated Hospital of Qingdao University and Peking University People’s Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. The animal study was reviewed and approved by the Ethics Committees of The Affiliated Hospital of Qingdao University and Peking University People’s Hospital.

Author Contributions

BZ, BY, and WG planned and designed the study. BZ and CZ were carried out the relevant experiments. BZ and GQ analyzed the data. CR and PY were material contributors. BZ prepared the manuscript. All authors unanimously approved the final draft.

Funding

The source of funding for the study was grants from the A Project of Shandong Province Higher Educational Science and Technology Program (No. J16LL06). The funder did not partake in any steps of the study from study design to manuscript preparation.

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.

Abbreviations

ICI, immune checkpoint inhibitor; LIMK2, LIM domain kinase 2; PD-1, programmed death-1; PD-L1, programmed death ligand-1; PD-L2, programmed death ligand-2; STAT3, signal transducer and activator of transcription 3; VEGFR2, vascular endothelial growth factor receptor-2.

References

1. Nagao-Kitamoto H, Setoguchi T, Kitamoto S, Nakamura S, Tsuru A, Nagata M, et al. Ribosomal protein S3 regulates GLI2-mediated osteosarcoma invasion. Cancer Lett. (2015) 356:855–61. doi: 10.1016/j.canlet.2014.10.042

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sevelda F, Mayr L, Kubista B, Lotsch D, van Schoonhoven S, Windhager R, et al. EGFR is not a major driver for osteosarcoma cell growth in vitro but contributes to starvation and chemotherapy resistance. J Exp Clin Cancer Res. (2015) 34:134. doi: 10.1186/s13046-015-0251-5

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Judson I, Verweij J, Gelderblom H, Hartmann JT, Schoffski P, Blay JY, et al. Doxorubicin alone versus intensified doxorubicin plus ifosfamide for first-line treatment of advanced or metastatic soft-tissue sarcoma: a randomised controlled phase 3 trial. Lancet Oncol. (2014) 15:415–23. doi: 10.1016/S1470-2045(14)70063-4

CrossRef Full Text | Google Scholar

4. Kempf-Bielack B, Bielack SS, Jürgens H, Branscheid D, Berdel WE, Exner GU, et al. Osteosarcoma relapse after combined modality therapy: an analysis of unselected patients in the cooperative osteosarcoma study group (COSS). J Clin Oncol. (2005) 23:559–68. doi: 10.1200/JCO.2005.04.063

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Brower, V. Hyperprogressive disease with anti-PD-1 and anti-PD-L1. Lancet Oncol. (2016) 17:e527. doi: 10.1016/S1470-2045(16)30590-3

CrossRef Full Text | Google Scholar

6. Champiat S, Dercle L, Ammari S, Massard C, Hollebecque A, Postel-Vinay S, et al. Hyperprogressive disease (HPD) is a new pattern of progression in cancer patients treated by anti-PD-1/PD-L1. Clin Cancer Res. (2017) 23:1920–8. doi: 10.1158/1078-0432.CCR-16-1741

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Cao Y, Zhang L, Kamimura Y, Ritprajak P, Hashiguchi M, Hirose S, et al. B7-H1 overexpression regulates epithelial-mesenchymal transition and accelerates carcinogenesis in skin. Cancer Res. (2011) 71:1235–43. doi: 10.1158/0008-5472.CAN-10-2217

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kim S, Koh J, Kim MY, Kwon D, Go H, Kim YA, et al. PD-L1 expression is associated with epithelial-to-mesenchymal transition in adenocarcinoma of the lung. Hum Pathol. (2016) 58:7–14. doi: 10.1016/j.humpath.2016.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Clark CA, Gupta HB, Sareddy G, Pandeswara S, Lao S, Yuan B, et al. Tumor-Intrinsic PD-L1 signals regulate cell growth pathogenesis, and autophagy in ovarian cancer and melanoma. Cancer Res. (2016) 76:6964–74. doi: 10.1158/0008-5472.CAN-16-0258

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lou Y, Diao L, Cuentas ER, Denning WL, Chen L, Fan YH, et al. Epithelial-mesenchymal transition is associated with a distinct tumor microenvironment including elevation of inflammatory signals and multiple immune checkpoints in lung adenocarcinoma. Clin Cancer Res. (2016) 22:3630–42. doi: 10.1158/1078-0432.CCR-15-1434

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ren T, Zheng B, Huang Y, Wang S, Bao X, Liu K, et al. Osteosarcoma cell intrinsic PD-L2 signals promote invasion and metastasis via the RhoA-ROCK-LIMK2 and autophagy pathways. Cell Death Dis. (2019) 10:261. doi: 10.1038/s41419-019-1497-1

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zheng B, Ren T, Huang Y, Guo W. Apatinib inhibits migration and invasion as well as PD-L1 expression in osteosarcoma by targeting STAT3. Biochem Biophys Res Commun. (2018) 495:1695–701.

Google Scholar

13. Hu X, Cao J, Hu W, Wu C, Pan Y, Cai L, et al. Multicenter phase II study of apatinib in nontriple- negative metastatic breast cancer. BMC Cancer. (2014) 14:820. doi: 10.1016/j.bbrc.2017.12.032

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Li J, Qin S, Xu J, Guo W, Xiong J, Bai Y, et al. Apatinib for chemotherapy-refractory advanced metastatic gastric cancer: results from a randomized, placebo-controlled, parallel-arm, phase II trial. J Clin Oncol. (2013) 31:3219–25. doi: 10.1200/JCO.2013.48.8585

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Peng H, Zhang Q, Li J, Zhang N, Hua Y, Xu L, et al. Apatinib inhibits VEGF signaling and promotes apoptosis in intrahepatic cholangiocarcinoma. Oncotarget. (2016) 7:17220–9. doi: 10.18632/oncotarget.7948

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Manegold C, Dingemans AC, Gray JE, Nakagawa K, Nicolson M, Peters S, et al. The potential of combined immunotherapy and antiangiogenesis for the synergistic treatment of advanced NSCLC. J Thorac Oncol. (2017) 12:194–207. doi: 10.1016/j.jtho.2016.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Seki N, Kan-O K, Matsumoto K, Fukuyama S, Hamano S, Tonai K, et al. Interleukin-22 attenuates double-stranded RNA-induced upregulation of PD-L1 in airway epithelial cells via a STAT3-dependent mechanism. Biochem Biophys Res Commun. (2017) 494:242–8. doi: 10.1016/j.bbrc.2017.10.045

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Garcia-Diaz A, Shin DS, Moreno BH, Saco J, Escuin-Ordinas H, Rodriguez GA, et al. Interferon receptor signaling pathways regulating PD-L1 and PD-L2 expression. Cell Rep. (2017) 19:1189–201. doi: 10.1016/j.celrep.2017.04.031

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Xu P, Sun Z, Wang Y, Miao C. Long-term use of indomethacin leads to poor prognoses through promoting the expression of PD-1 and PD-L2 via TRIF/NF-κB pathway and JAK/STAT3 pathway to inhibit TNF-α and IFN-γ in hepatocellular carcinoma. Exp Cell Res. (2015) 337:53–60. doi: 10.1016/j.yexcr.2015.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Horlad H, Ma C, Yano H, Pan C, Ohnishi K, Fujiwara Y, et al. An IL-27/Stat3 axis induces expression of programmed cell death 1 ligands (PD-L1/2) on infiltrating macrophages in lymphoma. Cancer Sci. (2016) 107:1696–704. doi: 10.1111/cas.13065

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Voorneveld PW, Kodach LL, Jacobs RJ, Liv N, Zonnevylle AC, Hoogenboom JP, et al. Loss of SMAD4 alters BMP signaling to promote colorectal cancer cell metastasis via activation of Rho and ROCK. Gastroenterology (2014) 147:196–208.e113. doi: 10.1053/j.gastro.2014.03.052

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ding Y, Milosavljevic T, Lahari SK. Nischarin inhibits LIM kinase to regulate cofilin phosphorylation and cell invasion. Mol Cell Biol. (2008) 28:3742–56. doi: 10.1128/MCB.01832-07

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Yoshioka K, Foletta V, Bernard O, Itoh K. A role for LIM kinase in cancer invasion. Proc Natl Acad Sci U.S.A. (2003) 100:7247–52. doi: 10.1073/pnas.1232344100

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Folkman, J. Angiogenesis in cancer vascular, rheumatoid and other disease. Nat Med. (1995) 1:27–31. doi: 10.1038/nm0195-27

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li Y, Wu Z, Yuan J, Sun L, Lin L, Huang N, et al. Long non-coding RNA MALAT1 promotes gastric cancer tumorigenicity and metastasis by regulating vasculogenic mimicry and angiogenesis. Cancer Lett. (2017) 395:31–44. doi: 10.1016/j.canlet.2017.02.035

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Xu H, Zhang Y, Peña MM, Pirisi L, Creek KE. Six1 promotes colorectal cancer growth and metastasis by stimulating angiogenesis and recruiting tumor-associated macrophages. Carcinogenesis. (2017) 38:281–92. doi: 10.1093/carcin/bgw121

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: VEGFR2, PD-L2, STAT3, metastasis, osteosarcoma

Citation: Zheng B, Zhou C, Qu G, Ren C, Yan P, Guo W and Yue B (2020) VEGFR2 Promotes Metastasis and PD-L2 Expression of Human Osteosarcoma Cells by Activating the STAT3 and RhoA-ROCK-LIMK2 Pathways. Front. Oncol. 10:543562. doi: 10.3389/fonc.2020.543562

Received: 17 March 2020; Accepted: 18 August 2020;
Published: 09 September 2020.

Edited by:

Antonella Zannetti, Institute of Biostructure and Bioimaging (IBB), Italian National Research Council Naples, Italy

Reviewed by:

Jun Yan, Fudan University, China
Andrea Porzionato, University of Padua, Italy

Copyright © 2020 Zheng, Zhou, Qu, Ren, Yan, Guo and Yue. 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: Wei Guo, bonetumor@pku.edu.cn; Bin Yue, bonetumor0532@163.com