MDA-MB-231 breast cancer cells overexpressing single VEGF isoforms display distinct colonisation characteristics.

BACKGROUND: Vascular endothelial growth factor (VEGF) is a multifunctional cytokine that has important roles in angiogenesis. Our knowledge of the significance of VEGF isoforms in human cancer remains incomplete.
METHODS: Bioluminescence imaging and transcriptomic analysis were used to study the colonisation capacity of the human breast cancer cells MDA-MB-231 controlling or overexpressing the VEGF165 or VEGF189 isoform (named cV-B, V165-B and V189-B, respectively) in nude mice.
RESULTS: When injected into the bloodstream, V189-B cells induced less metastasis in the lungs and bone than V165-B and cV-B control cells, consistent with longer survival of these mice and delay in tumour uptake in the mice injected with a V189-B clone. Histological analysis confirmed that there were less αSMA-positive cells in the lungs of the mice injected with V189-B. In vitro V189-B cells decreased both cell invasion and survival. Using transcriptomic analysis, we identified a subset of 18 genes expressed differentially between V189 and V165 cell lines and in 120 human breast tumours. V165 was associated with poor prognosis, whereas V189 was not, suggesting a complex regulation by VEGF isoforms. Our results showed a negative correlation between the expression pattern of VEGF189 and the levels of expression of seven genes that influence metastasis.
CONCLUSION: Our findings provide the first evidence that VEGF isoforms have different effects on breast cancer cell line colonisation in vivo.

The vascular endothelial growth factor-A (VEGF-A) is a key regulator of developmental, physiological and pathological neovascularisation, including tumour growth and metastasis (Folkman, 1995; Ferrara and Davis-Smyth, 1997). This growth factor consists of a family of proteins generated from a single gene by alternative splicing of the primary transcript (Tischer ). Processing results in human VEGF isoforms consist of proteins 121, 145, 165, 189 and 206 amino acids. The presence or absence of genomic exons 6 and 7 determines whether these isoforms are secreted in soluble forms (VEGF121 and VEGF165) or remains associated with a cell or extracellular matrix (VEGF189, VEGF206 and to some extent VEGF165) (Gitay-Goren ; Park ). The VEGF189 has the largest capacity for binding to extracellular matrix components. The most abundant isoforms, VEGF121 and VEGF165, have been the focus of intense studies. Interaction between VEGF isoforms and their endothelial receptors, VEGF-R1 (Flt-1), VEGF-R2 (KDR/Flk-1) and neuropilin-1 (NRP-1), stimulates receptor-associated kinase activity and initiates signalling pathways leading to angiogenesis (Ferrara and Davis-Smyth, 1997; Neufeld ; Gille ). By contrast, little is known about the role of the cell-associated VEGF189 isoform. In vitro, VEGF189 binds preferentially to VEGF-R1 and NRP-1 (Neufeld ; Hervé ). The VEGF189 can be cleaved to shorter bioactive forms by proteases, such as plasmin, urokinase plasminogen activator and MMPs (Plouët ; Lee ). The VEGF189 was recently shown to stimulate the proliferation and migration of endothelial cells (Hervé ). Using mouse embryos with low levels of heparin-binding VEGF-A isoform expression (VEGFwt/120), or expressing only VEGF120 that cannot bind to heparin (VEGF120/120), Ruhrberg demonstrated the involvement of heparin-binding VEGFs in vascular branching complexity at the earliest stages of angiogenic invasion in several organs . Besides, VEGF189 expression increases in the human endometrium during the secretory phase (Ancelin ), suggesting that this isoform also has a role in physiological vascular remodelling during the reproductive process.

Our knowledge of the significance of VEGFA isoform expression in human cancer remains incomplete. Only a few studies have considered the role of the VEGF189 isoform in tumour growth and metastasis. In mouse models, expression of the VEGF165 or VEGF189 isoforms has been associated with differences in tumour growth (Oshika ; Tomii ; Hervé ). Xenograft studies have shown that VEGF189-overexpressing colon cancer cells (SW-480-transfected cells) grow more slowly than VEGF165-overexpressing cells (Hervé ); this effect was interpreted as reflecting particular features of VEGF189 molecules. In breast cancers, VEGF165 and VEGF121 were shown to accelerate tumour development (Tomii ). We recently showed that VEGF189 contributes to mammary tumour growth through both angiogenic and non-angiogenic functions associated with an autocrine effect on breast cancer cells (Hervé ). However, VEGF189-expressing breast cancer cells that delay tumour uptake (Hervé ) are less angiogenic and form less rapidly growing xenograft tumours compared with VEGF121-expressing cells (Zhang ). Different kinds of results were obtained with VEGF189 overexpression, as dilated (Hervé ) or thin (Yuan ) microvessels were observed, depending on the studies. Several other groups have published on VEGF188-overexpressing fibrosarcomas. Grunstein showed that VEGF188 was associated with hypervascular density but not tumour growth. Tozer showed that fibrosarcomas derived from transgenic mice expressing only VEGF188 under constitutive promoter control developed highly vascularised well-defined tumours. Melanoma cells transfected with VEGF189 remain non-tumourigenic and dormant (Yu ). One study, in which VEGF isoform profile was classified into three groups, VEGF121 only (type 1), VEGF121+165 (type 2) and VEGF121+VEGF165+VEGF189 (type 3) (Tokunaga ), showed the higher incidence of liver metastasis in colon cancer patients for the type 3. Higher levels of cell-associated VEGF189 expression were also observed in lung cancers (Yuan ). In addition, high levels of the VEGF189 isoform are associated with poor prognosis in lung cancer and osteosarcoma (Lee ; Matsuyama ), whereas VEGF189 levels are inversely correlated with tumour stage and size in papillary renal cell carcinomas (Jacobsen ). Contradictory, results suggest that VEGF189 has a complex role in tumour development and progression.

Therefore, we studied the role of VEGF189 in tumour cell colonisation capacity, which is an important step in the formation of metastasis. To this end, we generated stable human breast carcinoma cell lines (MDA-MB-231) overexpressing the VEGF165 or VEGF189 isoform, named V165 and V189, respectively. The effect of VEGF189 overexpression on the extravasation potential was determined in vivo. We find that VEGF189 expression delays and decreases the colonisation behaviour of V189 cell lines in the lungs and bone. Moreover, analysis of the transcriptomic profiles of our clones identifies genes regulated differently in V189 and V165 cell lines. This study thus provides evidence of a role of VEGF189 in delaying breast cancer cell colonisation.

Materials and methods

Generation of stable bioluminescent MDA-MB-231 clones overexpressing VEGF isoforms

Human MDA-MB-231 breast cancer cells (ATCC, Molsheim, France) were used to generate stable transfected clones overexpressing VEGF165, VEGF189 or control vector (referred to as V165, V189 and cV, respectively) as previously described (Hervé ). Cancer cells were cultured in Dulbecco's modified Eagle's medium (DMEM, Invitrogen, France) supplemented with 10% fetal bovine serum (FBS; DMEM-10% FBS) and 100 μg μl−l sodium pyruvate in the presence of antibiotics (100 units per ml penicillin, 100 μg ml−1 streptomycin, 500 μg ml−1 geneticin G418, Life Technologies). Two clones of each V189, V165 and cV were used to generate stable bioluminescent cells. V165, V189 and cV clones were transfected with a plasmid expressing the luciferase gene and the hygromycin resistance gene (a gift from P. Balaguer, France), in the presence of Fugene 6 reagent, according to the manufacturer's instructions (Roche Diagnostics, Meylan, France). A parallel experiment was performed with a plasmid expressing the GFP gene. Successful gene transfer was confirmed 72 h after infection, by the visualisation of GFP on fluorescence microscopy. Cells transfected with the luciferase gene plasmid were selected on hygromycin (700 μg ml−1; Sigma-Aldrich) in DMEM-10% FBS. Resistant clones were isolated and maintained in DMEM-10% FBS supplemented with G418 (500 μg ml−1) and hygromycin (300 μg ml−1). All cell lines were cultured at 37 °C in a humidified atmosphere containing 5% CO2. Luciferase activity was confirmed in vitro by dispensing 1 × 105 cells into complete DMEM in each well of a 96-well plate. After 24 h, D-luciferin (150 μg ml−1, Caliper) was added to the medium in each well and the plate was incubated for 5 min. Luciferase activity was measured with the IVIS Imaging System (Xenogen, Caliper Life Science, Hopkinton, MA, USA) and analysed with Living Image software, as previously described (Abdelkarim ). The clones generated were named V165-B, V189-B and cV-B, according to the construct used for transfection. V189-B31, V165-B45 and cV-B9 cells were selected on the basis of their having equivalent levels of bioluminescence (104, 97 and 81 photons second per cell, respectively) (Supplementary Data 1A). The sh-NRP-1 stable clones were performed as described before (Vintonenko ).

Intracardiac experimental colonisation model

All animal experiments were carried out with the local and IUH Saint Louis ethics committees' approval. The ethical guidelines that were followed meet the standards required by the UKCCCR guidelines (Workman ). Female nude mice (8–10 weeks old, Janvier, Legenest St Ile, France) were anaesthetised by intraperitoneally (i.p.) injection of 120 mg kg−1 ketamine and 6 mg kg−1 xylazine (2% Rompun, Bayer SAS, Lyon, France), and by exposure to 1–3% isoflurane on subsequent imaging days. Intracardiac injection was performed as previously described (Jenkins ; Abdelkarim ). On day 0, bioluminescent tumour cells (105 in 100 μl PBS) were injected into the left ventricle of the heart. D-luciferin (150 mg kg−1) was injected i.p. and the animals were transferred to a light–tight chamber for whole-body optical imaging, with the acquisition of both dorsal and ventral views with an IVIS Imaging System (Xenogen) and an intensified charge-coupled device camera, over a period of 5 min (Abdelkarim ). Successful intracardiac injection was indicated on day 0 by systemic bioluminescence throughout the animal. During the whole experiment, colonisation sites were identified by persistent bioluminescence signal in the experiment time span using bioluminescence IVIS Imaging System (Xenogen). Both dorsal and ventral sites were observed. Colonised sites were determined by the presence of bioluminescence signal on the acquisition images without limitation of size. Colonisation was monitored in vivo, at least twice a week, after a period of 10 weeks. Killing was performed at the end of the study, or before in accordance with ethical guidelines, in particular weight loss (Workman ). As bioluminescent regions corresponded to the presence of tumour cells, we excised organs from the whole mice in each positive bioluminescent area. These organs were immediately visualised by ex vivo imaging (IVISTM Imaging System) to verify the persistence of bioluminescent signal and to avoid confusion in localisation (between thorax bone and lung for example); further, histological examination was performed to determine the presence of tumour cells. In each case, when a bioluminescent signal was detected in an organ, tumour cells were detected by histological observation. For all mice, all lung lobes were taken out.

Histopathology

Tissues were fixed in 4% paraformaldehyde overnight and embedded in paraffin (Hervé ). Supplementary decalcification was applied for the bone. Each sample was cut into 5-μm sections and stained with haematoxylin and eosin (H&E). All lung lobes were analysed for the presence of tumour cells.

Immunocytochemistry using a specific smooth muscle (SM)-alpha actin antibody (A2547, Sigma, St Quentin Flavier, France; 1/400) was performed on the mice's lung lobes in which the presence of tumour cells was detected by bioluminescence signal and confirmed by H&E. (All these correlations were not shown and ex vivo images were in Supplementary data 2C.) For each animal, two sections of lungs positive for tumour cell were immunostained with SM-alpha actin antibody. The quantification of positively stained cells was based on a grid-supported manual count. The results were expressed as the following percentage: (number of SM-alpha positive cells/total tumour cells in the lung) × 100.

Viability assay

Cells were plated on 96-well plates (5 × 103 per well) and grown in DMEM-10% FBS for 3 days. The growth of bioluminescent cells was measured using the XTT assay (Roche Diagnostics), as previously described (Hervé ). Absorbance was recorded at 450 nm, with a reference at 620 nm serving as a blank, with a microplate reader. The experiments were performed three times, each in triplicate.

Invasion assay

Cell migration experiments were performed in Boyden invasion chambers equipped with 8-μm pores filters and coated with Matrigel (Falcon, MA, USA). Equal volumes (0.5 ml per chamber) of cell suspension (3 × 104 cells per ml) were used to seed the Matrigel invasion chambers, which were incubated at 37 °C for 16 h. We used 10% FCS in DMEM as a chemoattractant. After 16 h of incubation, cells in the upper chambers that had not migrated were removed by wiping with a cotton swab. The cells on the lower surface of the filter were fixed, stained with crystal violet and counted (Abdelkarim ). The number of cells invading the Matrigel was determined by counting cells from 10 microscopic fields. Experiments were performed in triplicate.

Apoptosis/cell death detection

MDA-MB-231 cells transfected with a VEGF isoform or control construct were starved for 72 h in a serum-depleted medium (DMEM 1 g l−1 D-glucose, 0.5% FBS). Adherent cells were rinsed with PBS and subjected to double staining with fluorescein isothiocyanate-labelled annexin V and propidium iodide (PI), according to the manufacturer's instructions for the Annexin V apoptosis assay (Beckman Coulter, Fullerton CA, USA). Annexin V detects a phosphatidylserine translocation specific to early (API) and late (APII) stages of apoptosis, whereas APII and the first stage of necrosis are revealed by the incorporation of PI, which enters the cells when the cell membrane is damaged. Immediate analysis was carried out in a FACScan (Becton Dickinson, Sunnyvale, CA, USA) on 10 000 cells gated on the basis of forward and side scatter. Cell debris was excluded from the analysis on the basis of scatter signals. The position of the quadrant gate was set, so that the principal populations (>70%) of relevant untreated cells were viable (Annexin V–/PI–).

Isolation of cancer cells from mouse lungs

Human metastases were isolated from the lungs with collagenase I and cultured in DMEM (Gibco, St Aubin, France) supplemented with 10% FBS, 1% penicillin–streptomycin, 1% sodium pyruvate supplemented with geneticin G418 (0.5 mg ml−l; Gibco) in dishes coated with 0.5% gelatin. The human metastases were selected according to their resistance to geneticin used at 1 mg ml−1. All cell lines were cultured in a humidified atmosphere containing 5% carbon dioxide, at 37 °C, with a total medium replacement every other day.

RT-Q-PCR analysis

Reverse transcription of 1 μg of total RNA from each sample was performed with 100 units of Superscript II RNase H reverse transcriptase (Invitrogen, CergyPontoise, France) and 3 μM random hexamers (Pharmacia, Uppsala, Sweden). Differences in gene expression were assessed by calculating the fold change in gene expression with real-time quantitative RT–PCR (ABI Prism 7900 sequence Detection System, Perkin-Elmer Applied Biosystems, Foster City, SA, USA), with the TBP gene (TATA box-binding protein) as an endogenous RNA control.

Transcriptomic analyses

Whole-genome expression profiling of MDA-MB 231 overexpressing VEGF189 (V189), MDA-MB 231 overexpressing VEGF165 (V165) or control cells (cV) was carried out with AffymetrixGeneChip Human Gene 1.0 ST arrays. Similarly, MDA-MB-231-derived lung metastases obtained after the intracardiac injection of V165, V189 and cV clones, respectively, V165M, V189M or control clones (cVM) were also analysed with AffymetrixGeneChip Human Gene 1.0 ST arrays. These microarrays have 764 885 features corresponding to 28 869 known unique human genes. In all those clones, we isolated total RNA from 18 cell cultures (four biological replicates corresponding to four different passages of V189, V165 and cV, and two biological replicates corresponding to two different passages of 189M, 165M and cM clones) with the RNeasy mini kit (Qiagen, Valencia, CA, USA), by direct cell lysis in 10-cm culture dishes containing 600 μl RLT Buffer according to the manufacturer's instructions. Total RNA was quantified with a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Hudson, NH, USA) and its quality was assessed by running the eukaryotic total RNA nano assay on a Bioanalyzer 2100 (Agilent Technologies, Chandler, AZ, USA) with the Agilent RNA 6000 Nano Kit. Probes (300 ng of total RNA) were prepared for hybridisation to AffymetrixGeneChip Human Gene 1.0 ST arrays (Affimetrix, Cleveland, OH, USA), according to the manufacturer's instructions.

Data analyses

The scanned raw data were processed with the apt-probe-set summarise program from Affymetrix Power Tools (v1.10.2) for background correction and normalisation. The RMA (Robust Multichip Analysis) algorithm (Irizarry ), with full quintile normalisation, was applied, with array design library files V1 release 4, and normalised values were expressed on a log base 2 scale. Transcripts and genes were annotated with the Affymetrix release 29 annotation file based on the human genome version 18 assembly. Data were analysed with the R Bioconductor LIMMA package (Smyth, 2005). The microarray data have been deposited in the NCBI Gene Expression Omnibus (Edgar ) and are accessible through the GEO Series accession number GSE37911 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE37911). Probe sets corresponding to control genes or with a low-intensity signal (median log2 of intensity <4) or low variability were not considered for further analyses. Only genes with a false discovery rate <0.05 and a log fold change value ⩾±1 in paired comparison between treatments and their corresponding controls were considered to be differentially regulated. Data for V189 and V165 were compared by calculating the difference in log fold change for each gene (log fold change in the V165 data set−log fold change in the V189 data set). Only genes with a differential log fold change ⩾±1 and a false discovery rate <0.05 were considered to be differentially regulated.

Patients and samples

Human breast tumours were excised from women at the Institut Curie-Centre René Huguenin in St Cloud, France. All patients who entered this institution before 2007 were informed that their tumour samples might be used for scientific purposes and had the opportunity to decline. Since 2007, patients entering the institution have given their approval also by signed informed consent. This study was approved by the local ethics committee (Breast Group of the René Huguenin Hospital). We investigated the relationships between mRNA levels for genes of interest and metastasis-free survival by analysing 120 (60 which metastasised to distant sites and 60 which did not) human breast tumours excised from women at Institut Curie-Centre René Huguenin (St Cloud, France). This study was approved by the local ethics committee (Breast Group of the René Huguenin Hospital). Samples containing >70% tumour cells were considered suitable for this study. Immediately after surgery, the tumour samples were frozen in liquid nitrogen, in which they were stored until RNA extraction. The characteristics of the patients are shown in Table 1. We determined the oestrogen receptor, progesterone receptor and ERBB2 status of the tumours at the protein level, by biochemical methods (enzymatic immunoassay), with confirmation by real-time PCR. The malignancy of infiltrating carcinomas was scored according to Scarff Bloom and Richardson's histoprognostic system (Bloom and Richardson, 1957).

Table 1

Characteristics of the 120 human breast tumours

	Number of patients (%)	Number of relapse (%)	P-valuea
Total	120 (100.0)	60 (50)
Age
⩽50	20 (16.7)	11 (55)	0.60 (NS)
>50	100 (83.3)	49 (49)
SBR histological gradeb,c
I	14 (11.9)	5 (35.7)	0.14 (NS)
II	58 (49.2)	29 (50)
III	46 (39.0)	25 (54.3)
Lymph node statusd
0	30 (25.2)	14 (46.7)	0.48 (NS)
1–3	62 (52.1)	30 (48.4)
>3	27 (22.7)	16 (59.3)
Macroscopic tumour sizec
⩽25 mm	56 (47.5)	22 (39.3)	0.013
>25 mm	62 (52.5)	37 (59.7)
ER status
Negative	36 (30)	21 (58.3)	0.0019
Positive	84 (70)	39 (46.4)
PR status
Negative	56 (46.7)	31 (55.4)	0.041
Positive	64 (53.3)	29 (45.3)
ERBB2 status
Negative	98 (81.7)	49 (50)	0.72 (NS)
Positive	22 (18.3)	11 (50)
Molecular subtypes
RH− ERBB2−	25 (20.8)	15 (60)	0.016
RH− ERBB2+	11 (9.2)	6 (54.5)
RH+ ERBB2−	73 (60.8)	34 (46.6)
RH+ ERBB2+	11 (9.2)	5 (45.5)
PIK3CA mutation status
Wild type	79 (65.8)	42 (53.2)	0.070 (NS)
Mutated	41 (34.2)	18 (43.9)

Abbreaviation: NS=non significant.

Log-rank: survival distributions were estimated by the Kaplan–Meier method, and the significance of differences in survival rate (P) was assessed in log rank tests.

Scarff Bloom Richardson classification.

Information available for 118 patients.

Information available for 119 patients.

Bold values: P<0.05.

In silico analyses

We identified genes involved in breast carcinogenesis by analysing genes expressed differentially in the V165 and V189 clones (V165 vs V189) or in 165M and 189M lung metastases (165M vs 189M), on Oncomine (www.oncomine.org), with the Richardson breast 2 cancer data set in particular. This data set corresponds to 40 ductal breast carcinomas and seven normal breast samples that were analysed on Affymetrix U133 Plus 2.0 microarrays (Richardson ) (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE3744s). For each gene, we checked the expression level (fold change) and the corresponding P-value (P<0.05 considered significant) in the Richardson breast 2 data set. We took into account only the genes considered to be significantly misregulated. We then analysed those misregulated genes with BRB-ArrayTools, version 3.3β3 (Rockville, MD, USA), developed by Dr R. Simon and Dr A. Peng Lam on a series of human breast tumours (n=344) described in detail elsewhere (Wang ). We used a filter to refine the metastasis-related gene set. The Cox proportional hazards model and Wald statistics were used for univariate analysis. Metastasis-free survival and P-values<0.05 were considered significant in each univariate test. Only significantly misregulated genes were taken into account in subsequent real-time PCR analysis.

Statistical analyses

The graphs show the means and the standard errors of the means. We used a suitable two-way ANOVA multiple comparison test including the Kruskal–Wallis test followed by the Dunn post-test. The statistical analyses were performed using GraphPad Prism 5.01 software (La Jolla, CA, USA). *P<0.05. To analyse the prognostic value of the genes tested by RT–PCR in the series of 120 human breast tumours, the optimal cutoff for each gene was determined by calculating the area under the curve for a comparison of tumours without metastases (n=60) and tumours giving rise to metastases (n=60). We used these optimal cutoffs to determine survival for each gene. Survival distributions were estimated by the Kaplan–Meier method, and the significance of differences in survival rate was assessed in log rank tests. The relationships between mRNA levels for the various target genes were investigated in Spearman's rank correlation tests (assessing the correlation between two quantitative parameters). Differences were considered significant if P<0.05.

Results

VEGF189 overexpression reduced the number of mice with MDA-MB-231 colonisation sites delaying the lung and bone tumour uptakes

Bioluminescent breast cancer cell lines can be used for a rapid and sensitive in vivo detection of cell dissemination throughout the body of immune-deficient mice (Supplementary Data 1A and B; Jenkins ; Minn ; Abdelkarim ). Bioluminescent V189-B and V165-B clones and cV-B control cells were injected into the left heart ventricle of female nude mice (n=10 per group). Fourteen days later, the bioluminescent tumour cells were clearly detectable and their numbers increased over a period of 6–10 weeks. In contrast to cV and V165 clones, V189 resulted in a significantly longer time for tumour colonies to appear in the mice (Figure 1A, left panel). The number of colonisation sites appeared in a smaller proportion after the injection of VEGF189-overexpressing cells than after injection of the controls or VEGF165-overexpressing cells (V189-B vs cV orV165-B) (Figure 1A; right panel; Supplementary Data 2A). After 16 days, only 10% of the mice receiving injections of the V189-B clone had developed tumour colonised sites, against 50% of the mice receiving injections of cV-B and V165-B cells (Figure 1A, left panel). After 38 days, almost all the mice (90–100%) receiving cV-B or V165-B clones injections had developed tumour sites, vs 70% of the mice injected with V189-B clones (Figure 1A, left). V189-B-injected mice had a significantly smaller total number of colonisation sites than those injected with the V165-B and cV-B clones after 16 days (P=0.04) and this difference is lower and only significant for V189 vs cV-B cells (P=0.037) after 38 days of the time course, underlying the delaying effect of VEGF189 overexpression (Figure 1A, right panel). Moreover, the lower colonisation observed with the V189 clone was associated with greater mouse survival (Supplementary Data 2B).

Figure 1

VEGF189 overexpression in MDA-MB-231 breast cancer cells decreases the number of mice with colonisation sites (A) Tumour formation for the various bioluminescent clones injected into mice (n=10). Left panel: number of mice with colonisation sites. Right panel: number of colonisation sites per mouse. (B) Number of mice with colonisation present in the bone and lungs after injection with the various clones (V165, V189 and cV). At the indicated times, the total numbers of mice with metastasis in the bone (legs) and lungs were counted in each group. *PV189 vs PC-V and PV165<0.05 and **PV189 vs PC-V<0.05. (C) Haematoxylin- and eosin-stained sections of breast cancer tumour in lung tissue (T). Different patterns of colonisation are shown: more or less invasive (middle and right, respectively). Original objective magnification, × 16 (top), × 40 (bottom), scale bars are 500 and 200 μm, respectively.

The intracardiac injection of V189-B cells into nude mice resulted in fewer mice with lung or bone metastasis, than observed with injection of V165-B and cV-B control cells (Figure 1B; left and right panels, respectively). In both cases (bone and lung), we observed a delay in tumour uptake with V189-B cells (25 days), compared with V165-B cells (9 days) and cV-B controls (18 days). Histological examination of the lungs excised from the mice confirmed cancer cell invasion and showed heterogeneity in the lung parenchyma between VEGF189-B, V165-B and cV-B clones (Figure 1C). The presence of tumour burden was correlated to the presence of bioluminescence signal but no significant differences were observed between the mean of ROI quantification in each group (data not shown).

We observed that bone colonisation in the mice also occurred later with V189-B clones (25 days) than with V165-B (9 days) or cV-B clones (9 days) (Figure 1B). Salivary and adrenal glands (kidney) were also colonisation sites in addition to the lungs, bone and brain but to a lesser extent (Supplementary Figure 2C). The number and pattern of tumour formation in the brain differed from that in the lungs and bone (data not shown). V189-B clones induced fewer mice with colonisation in the brain before day 32 than V165-B and cV-B cells (data not shown), but this difference subsequently disappeared and was not observed at day 40.

VEGF as an autocrine factor and regulation via NRP-1

We previously demonstrated that VEGF can affect non-angiogenic aspects of tumour cell progression, through autocrine loop acting on tumour cells (Bachelder , 2003; Lee ; Hervé ; Vintonenko ; Perrot-Applanat and Di Benedetto, 2012). VEGF165 has been shown to increase the invasiveness and survival of breast cancer cells (Orimo , Oskarsson ). We further analysed the invasiveness of VEGF189-overexpressing clones, in the Matrigel invasion chamber assays. As NRP-1 was demonstrated to interact with VEGF189 and to modulate VEGF effect on cell survival (Vintonenko ), we performed the silencing of NRP-1 to investigate the modulation of invasion induced by VEGF 189 and 165 isoforms (Figure 2A). The number of V189-overexpressing MDA-MB-231 cells invading the Matrigel was significantly decreased by half compared with the more invasive control cV-B and V165-B clones (Figure 2A). We observed that silencing the expression of NRP-1 had no effect on VEGF165-B clones or the cV-B control invasion, whereas it reversed the decrease of VEGF189-expressing cell invasion. In addition, the clones were analysed in cell culture in limiting conditions (0.5% FBS and 1 g l−i glucose) by FACS analysis (Figure 2B and C). As expected, the number of apoptotic V189-B cells (AnnV+/PI and AnnV+/PI+) was increased about two-fold as compared with cV-B and V165-B cells (Figure 2B and C). Thus, VEGF189 may act as an autocrine factor via NRP-1, resulting in lower levels of cell invasion than with the VEGF165 isoform.

Figure 2

Lower invasive capacity and survival of V189-B clones in limiting serum conditions. (A) Histogram showing invasive capacities of cV-B, V189-B and V165-B clones and corresponding sh-NRP-1. Results are presented as % invasion with respect to control cells, for which the value was set at 100% (means of four independent experiments ±s.d.). ***P<0.05. Significant values were shown for V189 vs cV or V165 vs cV. **P<0.05. Significant values were shown for V189 vs sh-NRP-1 of VEGF189 clones (B and C). Survival/apoptotic assay. (B) Histogram showing the mean±s.d. of apoptotic cells in the cV-B, V189-B and V165-B clones. *P<0.01. Significant values were shown for V189 and V165 vs cV. (C) Representative flow cytometry profiles for cV-B, V189-B and V165-B clones incubated in medium supplemented with 0.5% FBS.

VEGF189 and VEGF165 lung metastasis exhibit different stromal compartment

As the presence of SM-alpha actin-positive cells is strongly associated with tumour aggressiveness (Orimo ), we performed immunostaining with the SM-alpha actin antibody to detect myofibroblasts in the compartment of lungs positive for V165-B and V189-B tumour cells (Figure 3A). The percentage of positive cells was higher in the lungs colonised by V165-B cells than in ones colonised by V189-B or cV-B cells (Figure 3A). The overexpression of VEGF165 and VEGF189 in isolated lungs positive for tumour cells (called V165, V189 and cV, respectively) was confirmed by RT–PCR (Figure 3C). The correlation between the expression evaluation by RT–PCR of the VEGF isoforms and their protein levels was previously demonstrated (Hervé ). Besides, V189 cells also expressed significantly lower levels of MMP1mRNA than V165 cells (Figure 3C). We further investigated the expression in those clones of tenascin-C (TNC), a well-known marker of stroma associated with epithelial malignancy (Mackie ; Dandachi ). V189 cells expressed significantly lower levels of TNC mRNA than cV and V165 cells (Figure 3D).

Figure 3

Stromal compartment and migration of gene expression. (A and B) Immunostaining (A) and quantification (B) for SM-alpha actin. For each animal, two sections per lung positive for tumour cells were immunostained with SM-alpha actin. Original objective magnification, × 32; T (tumor), L (Lung), **P<0.01, significant for V165, respectively, vs cV. (C and D) RT-Q-PCR analyses of VEGF isoforms, MMP-1 (C) and TNC (D) in isolated cVM, V189M and V165M lung tumours cells. TBP was used as an internal control. The relative mRNA levels of each gene in the control are set to 1, and values are expressed as means±s.e.m. *P<0.05, **P< 0.01. Scale bar, 250 μm.

Transcriptomic profiles of the V189 vs the V165 cell line and of the 189M vs the 165M lung metastatic cell line

We compared the transcriptomic profiles of the V165 and V189 cell lines (Supplementary Data S4). We identified a subset of 149 genes with significantly different expression levels in the two cell lines (V165 vs V189) (P<0.05 and |Log FC|>1) (Supplementary Data S4). We reduced the number of genes for further analysis, by comparing our 149 deregulated genes with a list of human genes known to be specifically involved in mammary carcinogenesis. We used the Richardson breast 2 database (oncomine.org), which compared the expression of 19 514 human genes between 40 breast tumours and seven normal breast tissue samples. We found that 36 genes expressed differentially between V165 and V189 cells (V165 vs V189 list) were also included in the Richardson breast 2 database (Figure 4A). We then analysed the expression profiles of the 36 genes identified with BRB-ArrayTools and a metastasis potential event selection filter, and identified 11 genes expressed significantly (Figure 4A).

Figure 4

(A) In silico analysis comparing V165 with V189 as described in Materials and methods, and genes significantly expressed in the BRB data set and associated with invasive potential in the comparison of V165 with V189. (B) In silico analysis comparing V165M with V189M as described in Materials and methods, and genes significantly expressed in the BRB data set and associated with invasive potential in the comparison of V165M with V189M.

As described above, animals developed distant lung colonisation after intracardiac injection of the V165, V189 and cV clones. We applied the same protocol to these lung metastatic cell lines. We first identified a subset of 150 genes displaying significant differential expression between 165M and 189M cells (165M vs 189M) (P<0.05 and |Log FC| >1) (Figure 4B; Supplementary Data S5). After the specific selection of genes involved in metastatic processes, from the Richardson breast 2 database and BRB ArrayTools, respectively (Figure 4B), we identified seven genes expressed significantly (Figure 4B).

We thus identified a subset of 18 genes both regulated differentially between V165 (or 165M) and V189 (or 189M) and involved in breast cancer and metastasis potential (Figures 4A and B).

Overexpression of 18 selected genes in 120 human breast cancers

We further analysed the expression of these 18 genes and four other genes as controls (TNC, VEGFA, VEGF165 and VEGF189) in a series of 120 human breast tumours (Table 1), which we compared with 10 normal breast tissue samples by quantitative RT–PCR (Table 2). Values of 3 (or more) were considered to indicate marked overexpression, and values of 0.33 (or less) were considered to indicate marked underexpression. Overexpression in breast tumours was the change in expression profile most frequently found in the 22 selected genes compared with normal tissues (Table 2). Changes in gene expression were detected in 13 of the 22 genes. We found that 13 of the 22 genes were frequently overexpressed in breast tumours compared with normal tissues (% overexpression >10%), none of these genes was exclusively underexpressed in breast tumours (% underexpression >10%) and 2 of these 18 genes (i.e., HRASLS and TNC) were overexpressed in some breast tumours and underexpressed in others (Table 2). A higher expression of VEGFA was also observed in breast tumours than in normal breast tissues (27.5%, median mRNA level: 1.82, range: 0.42–13.37), and this was more marked for VEGF165 isoforms, which was associated with poor prognosis (26.7%, median mRNA level: 1.89, range: 0.35–16.43) than for VEGF189 isoforms (19.2%, median mRNA level: 1.72, range: 0–15.54) (Table 2).

Table 2

Quantitative PCR analyses of the 22 genes of interest in 120 human breast tumours

Genes	Ct median of normal breast tissues (n=10)	Normal breast tissues (n=10)	Breast tumours (n=120)	% of underexpressed tumours (mRNA level <0.33)	% of normal expressed tumours	% of overexpressed tumours (mRNA level >3)
CAMSAP1L1	27.84 (25.16–29.6)a	1.00 (0.53–1.55)b	1.02 (0.26–27.10)b	1.67%	97.50%	0.83%
CARS	28.53 (26.57–30.21)	1.00 (0.63–1.72)	1.23 (0.46–14.24)	0.00%	93.33%	6.67%
CPNE3	26.38 (24.37–27.75)	1.00 (0.77–1.36)	1.43 (0.34–17.15)	0.00%	89.17%	10.83%
CTSL1	27.18 (25.58–27.95)	1.00 (0.42–2.46)	1.11 (0.35–14.87)	0.00%	94.17%	5.83%
FAM82B	28.00 (25.63–29.59)	1.00 (0.65–1.35)	1.17 (0,29–3.96)	6.80%	85.80%	7.40%
GGH	28.99 (27.43–30.56)	1.00 (0.52–1.45)	1.79 (0.19–13.94)	2.50%	68.33%	29.17%
GINS2	30.15 (28.30–32.89)	1.00 (0.33–1.79)	2.55 (0,17–50.25)	0.83%	56.67%	42.50%
GPRC5A	28.08 (25.99–30.82)	1.00 (0.27–1.92)	3.52 (0,06–84.28)	3.33%	40.00%	56.67%
HIST1H2BK	23.97 (21.94–26.18)	1.00 (0.59–4.55)	1.46 (0.44–10.02)	0.00%	83.33%	16.67%
HRASLS	31.89 (29.19–34.00)	1.00 (0.49–1.43)	1.27 (0–12.94)	19.17%	58.33%	22.50%
IMPA1	29.23 (27.24–30.25)	1.00 (0.73–1.31)	1.22 (0.35–12.06)	0.00%	97.50%	2.50%
MMP1	34.89 (34.17–49.02)	1.00 (0–19.2)	75.31 (0–71.89)	2.50%	0.00%	97.50%
PTPN7	31.47 (28.81–34.17)	1.00 (0.33–5.20)	2.07 (0.36–16.93)	0.00%	69.17%	30.83%
PXDN	27.49 (25.86–29.16)	1.00 (0.64–5.73)	2.13 (0.35–33.53)	8.10%	90.40%	1.50%
SFT2D2	26.22 (23.96–27.93)	1.00 (0.68–1.22)	0.74 (0.27–7.67)	2.50%	95.83%	1.67%
TBC1D8	28.44 (26.23–29.85)	1.00 (0.90–1.30)	1.19 (0.33–13.40)	8.10%	90.40%	1.50%
TPD52	27.03 (24.13–28.93)	1.00 (0.56–1.77)	3.73 (0.83–47.30)	0.00%	35.83%	64.17%
ZDHHC13	29.27 (26.83–30.98)	1.00 (0.56–1.34)	1.31 (0.18–6.04)	2.50%	91,67%	5.83%
TNC	27.38 (25.73–31.20)	1.00 (0.16–1.80)	0.74 (0.07–30.81)	25.00%	60.83%	14.17%
VEGF165	29.17 (26.76–30.20)	1.00 (0.61–1.68)	1.89 (0.35–16.43)	0.00%	73.33%	26.67%
VEGF189	30.55 (29.46–31.93)	1.00 (0.60–2.52)	1.72 (0–15,54)	1.67%	79.16%	19.17%
VEGFA	28.21 (26.62–29.71)	1.00 (0.69–2.15)	1.82 (0.42–13.37)	0.00%	72.50%	27.50%

Abbreviation: TNC=tenascin-C.

TBP was used as an internal control. Relative mRNA levels for each gene in the control are set at 1. The median Ct for normal breast tissue is indicated as median Ct (min–max) and median gene expression level is indicated as median (min–max). Relative mRNA levels of 3 or more were considered to indicate marked overexpression and values of 0.33 or less were considered to indicate marked underexpression.

Median (range) of genes Ct values.

Median (range) of gene mRNA levels; the mRNA values of the samples were normalised such so that the median of the 10 normal breast tissues mRNA value was 1. Bold values: under or overexpression >10%.

Relationship between overexpression of VEGF189 and the 18 selected genes

We then carried out Spearman's rank correlation tests to assess the extent to which the expression of the 19 genes identified (18+TNC) was correlated with that of VEGF189/VEGFA ratio (Figure 5A). We found a statistically negative correlation between the expression pattern of VEGF189 and the levels of expression of seven genes: CARS (r=−0.381, P=0.000035), CPNE3 (r=−0.226, P=0.012), GINS2 (r=−0.273, P=0.0027), GPRC5A (r=−0.219, P=0.016), HIST1H2BK (r=−0.338, P=0.00023), TBC1D8 (r=−0.181, P=0.045) and ZDHHC13 (r=−0.270, P=0.003). Levels of VEGF165 expression were found to be correlated with the levels of expression of CTSL-1, GGH, HIST1H2BK, HRASLS, SFT2D2, TNC and MMP1 (data not shown).

Figure 5

Statistical analysis of mRNA levels. (A) Spearman's rank correlation test for the genes of interest and the proportion of total VEGF-A in the form of VEGF189. Differences were considered significant if P<0.05. (B) Survival curves with an optimal cutoff point. Survival was estimated by the Kaplan–Meier method, and groups were compared in log-rank tests. The optimal cutoff point was calculated for each gene and used to distinguish two groups as described in Materials and methods.

We also tested the prognostic value of those seven genes correlated with expression pattern of VEGF189. As in our series of 120 human breast tumours (60 metastasising to distant sites and 60 displaying no metastasis; Table 1), we used a log rank test to determine the relationship between metastasis-free survival (MFS) and the mRNA levels for the seven selected genes. The differences in MFS between patients with low and high levels of expression of the selected genes were significant for two genes (Figure 5B). A significant difference in MFS between patients with low and high levels of GPRC5A, TBC1D8 and VEGF expression was observed (P=0.000012, P=0.009 and P=0.044, respectively). Thus, high expression of GPRC5A, TBC1D8 and VEGFA expressions were strongly correlated with longer MFS (Figure 5B), suggesting that the overexpression of those genes in tumours is associated with decreased metastasis potential with GPRC5A and TBC1D8 as opposed to high metastasis potential with VEGF165.

Discussion

Transformed cells must acquire additional abilities to surmount the natural barriers against metastasis such as survival and migration/invasion. Thus, after entering the bloodstream, cancer cells must infiltrate, survive and then overtake new tissue. This study is the first to investigate the effects of VEGF isoforms in an important step involved in metastasis that is the colonisation of tissues. Previous results from our group for breast cancer (Hervé ) and other cancers (Tokunaga ; Yuan ) have suggested that VEGF189 has an important role locally in aggressive tumours, especially by delaying tumour growth potential, invasion and migration of xenografted cells, and tumour angiogenesis (Perrot-Applanat and Di Benedetto, 2012). Here our results show that VEGF189 and VEGF165 differently influence the circulating tumour cell ability to colonise new tissues. The main colonisations occurred in the lungs and bone and to a lesser extent in other organs. V189-injected mice had less colonisation sites in vivo than mice receiving injections of the V165 and cV clones (Figure 1), particularly in the bone and lungs. Differences in colonisation were less important when considering the brain, suggesting that the lungs or bone are the first sites colonised by breast tumour cells (Yoneda ). Also, we observed that the intracardiac injection of MDA-MB-231 cells overexpressing VEGF induced colonisation of the lymph node, skeletal tissues, lungs and brain, with all clones, as described for breast cancer cells, including MDA-MB-231 (Minn ; Steeg, 2006; El-Mabhouh ; Bos ; Price, 2010). Although injecting cells into the bloodstream is limited compared to the orthotopic model and that measurement of total tumour lung burden would have improved the study, our data as a whole suggest distinct effects of VEGF189 compared with VEGF165 isoforms on colonisation uptake, which could be associated with different effects on tumour cell functions.

Lungs with VEGF189 tumours were associated with significantly lower levels of αSMA-positive cells than VEGF165 tumours. Although we could not conclude on the stromal or tumour origin of this staining, lungs with VEGF189 cell colonisation expressed lower levels of TNC and MMP1 than those containing VEGF165 cells. Joan Massagué's laboratory recently showed that TNC is associated with aggressive lung metastasis as described before (Mackie ; Oskarsson ). Indeed, TNC enhances the expression of stem cell signalling components, thereby promoting the survival and outgrowth of micro-metastases in the lung, a finding highlighting the relevance of TNC as an ECM component of the metastatic niche (Oskarsson ). The metalloproteinase MMP1 was shown to facilitate angiogenesis and to increase ECM degradation, two crucial processes for the invasive and migratory phenotype of tumour cells (Gialeli ). MMP-1 level in tumour and stromal cells is significantly associated with tumour progression (Foley ). When expressed in breast cancer cells, MMP-1 increases the assembly of new tumour blood vessels, the release of tumour cells into the bloodstream and the breaching of lung capillaries by circulating tumour cells to initiate lung metastasis (Gupta and Massagué, 2006; Foley ).

Little is known about the differential expression and function of VEGF-A isoforms in human breast tumours. Transcriptomic analysis identified new genes expressed differentially in V165, V189 and cV clones. On the basis of the breast cancer metastasis gene signature (www.oncomine.org), we showed that those genes were specifically involved in metastasis formation. In addition, we carried out quantitative PCR to quantify the transcripts of those genes of interest in 120 breast tumours from patients with and without distant metastases. VEGF165 overexpression was detected in 26.7% of human breast cancer tumours, whereas VEGF189 overexpression was detected in 19.2% in our series of 120 tumours. High VEGF165 expression was correlated with poor prognosis in our series of 120 breast tumours (data not shown). VEGFA, GPRC5A and TBC1D8 low expression is also associated with the decrease of patient survival (Figure 5). However, although high VEGF189/VEGFA expression was significantly correlated with low expression of GPRC5A and TBC1D18, we could not significantly associate VEGF189 overexpression with poor prognosis, suggesting a complex regulation of these genes by VEGF189.

We found that VEGF189 increased cell adhesion to ECM (as previously described by Hervé ) and reduced both invasion in the Matrigel assay and survival in stress conditions, as shown by comparison with VEGF165. Furthermore, the higher levels of apoptosis observed with VEGF189 in limiting culture conditions probably prevented many cells from surviving in the bloodstream or from remaining in the metastatic niche. The molecular basis of those differences is unknown. Neuropilin 1, a coreceptor for VEGF (Soker ), is strongly expressed in several tumour cells, in which it has been implicated in cell migration and survival (Bachelder , 2003; Vintonenko ; Perrot-Applanat and Di Benedetto, 2012). We recently showed that VEGF189 binds to NRP1 (Hervé ), in addition to VEGF165. The affinity of VEGF189 for NRP1 appears to be stronger than that of VEGF165 (Vintonenko ). NRP1 was also shown to interact with α5β1 integrin, enhancing pancreatic cancer cell survival and invasion (Fukasawa ). In contrast in VEGF189 overexpressing cells, we found that VEGF189 bound to NRP-1 decreased breast cancer cell survival. (Vintonenko ). Herein, we also demonstrated that VEGF189 autocrine effect on invasion was modulated through the NRP1 pathway as it was the case for survival (Vintonenko ). Furthermore, these results were in accordance with those of our previous work demonstrating that sh-NRP-1 of VEGF189 clones resulted in the inhibition of UPA expression and activity.

Survival is more important for colonisation than proliferation, particularly in the bloodstream. These properties may account for the lower levels of colonisation observed with V189 than with V165 clones in the first steps. Cell migration and invasion result from a balance between cell adhesion and detachment, both of which are required for motility. If the adhesion force is too strong, the cells attach tightly to the substrate and are unable to disperse. The greater adhesion of cells to the ECM with the VEGF189 clone, previously reported by Hervé , is thus also consistent with a lower capacity of the cells to invade distant organs and to colonise. Furthermore, the higher levels of apoptosis observed with VEGF189 in limiting culture conditions probably prevented many cells from surviving in the bloodstream or from remaining in the metastatic niche.

In conclusion, we found that overexpression of VEGF189 results in fewer mice with colonisation in the lungs and bone than VEGF165 overexpression. Our results are important as they suggest that the nature of overexpressed VEGF isoform influences tissue colonisation by tumour cell, in a different manner.