The impact of KRAS mutations on VEGF-A production and tumour vascular network
© Figueras et al.; licensee BioMed Central Ltd. 2013
Received: 18 October 2012
Accepted: 8 March 2013
Published: 18 March 2013
Skip to main content
© Figueras et al.; licensee BioMed Central Ltd. 2013
Received: 18 October 2012
Accepted: 8 March 2013
Published: 18 March 2013
The malignant potential of tumour cells may be influenced by the molecular nature of KRAS mutations being codon 13 mutations less aggressive than codon 12 ones. Their metabolic profile is also different, with an increased anaerobic glycolytic metabolism in cells harbouring codon 12 KRAS mutations compared with cells containing codon 13 mutations. We hypothesized that this distinct metabolic behaviour could be associated with different HIF-1α expression and a distinct angiogenic profile.
Codon13 KRAS mutation (ASP13) or codon12 KRAS mutation (CYS12) NIH3T3 transfectants were analyzed in vitro and in vivo. Expression of HIF-1α, and VEGF-A was studied at RNA and protein levels. Regulation of VEGF-A promoter activity was assessed by means of luciferase assays using different plasmid constructs. Vascular network was assessed in tumors growing after subcutaneous inoculation. Non parametric statistics were used for analysis of results.
Our results show that in normoxic conditions ASP13 transfectants exhibited less HIF-1α protein levels and activity than CYS12. In contrast, codon 13 transfectants exhibited higher VEGF-A mRNA and protein levels and enhanced VEGF-A promoter activity. These differences were due to a differential activation of Sp1/AP2 transcription elements of the VEGF-A promoter associated with increased ERKs signalling in ASP13 transfectants. Subcutaneous CYS12 tumours expressed less VEGF-A and showed a higher microvessel density (MVD) than ASP13 tumours. In contrast, prominent vessels were only observed in the latter.
Subtle changes in the molecular nature of KRAS oncogene activating mutations occurring in tumour cells have a major impact on the vascular strategy devised providing with new insights on the role of KRAS mutations on angiogenesis.
Ras proteins have been the subject of intense research as signalling molecules in normal and neoplastic cells . Yet, a complete understanding of their exact mode of action is still to come. Among the three RAS genes (H-RAS, KRAS and N-RAS) KRAS is the most commonly activated in human tumours. Several lines of evidence suggest that not only the presence or absence of a KRAS mutation but its molecular nature influences tumour cell behaviour [2, 3]. A reduced transforming capacity of codon 13 mutation as compared with codon 12 is observed in vitro and in vivo, with short latency times to tumour-appearance for codon 12 KRAS overexpressing cells [4–6]. Moreover, our previous results indicate that distinct mutations associate with specific metabolic phenotypes, an increased anaerobic glycolytic metabolism in cells containing codon 12 KRAS compared with cells containing codon 13 mutations. Switching to a glycolytic metabolism is a rapid adaptation to hypoxia that can be related to HIF1α expression .
Perpetual blood vessel formation and remodelling (angiogenesis) is a hallmark of cancer and a prerequisite for three-dimensional tumour growth, invasion, and metastasis . Hypoxia, by inducing HIF-1α, promotes the expression of VEGF-A, the main pro-angiogenic hypoxia-induced gene . However, oncogenes are also per se potent inductors of angiogenesis . Ras proteins are a paradigm for oncogene-dependent induction of tumour angiogenesis due to their involvement in the regulation of key pro and anti angiogenic factors [11–14]. However, its cross-talk with hypoxia-dependent signals is not so clear.
To gain further insight into the metabolic potential and distinct aggressiveness of different activating KRAS mutations, we examined the expression levels of HIF-1α and VEGF-A in stable mutated 12 and 13 NIH3T3 transfectants. Our results in vivo and in vitro indicate that the distinct KRAS mutations generated different normoxic HIF-1α responses. Moreover, different VEGF-A expression patterns were observed that are independent of the HIF-1α status but dependent upon ERKs stimulation. These alterations associated with distinct tumoral angiogenic profiles.
NIH3T3 cells were produced as previously described [4, 15], with plasmid DNA containing a KRAS minigene with a G:C A:T mutation (CYS12) at the first position of codon 12 (pMLK12), a G:C A:T mutation (ASP13) at the second position of codon 13 (pMLK13), and a control plasmid containing the expression vector alone (pMLneo). pMLK12, pMLK13, and pMLKwt plasmids were a gift of Dr. Manuel Perucho of the Burnham Institute at La Jolla, CA. Levels of expression of the KRAS protein in the selected clones used were similar .
Clones were cultured in DMEM supplemented with 20% Fetal Calf Serum and 500 μg/ml of neomycin G418. Mutations were verified by direct sequencing prior to the initiation of every experiment.
Transfected cells cultured 12 hours in FCS deprivation were incubated 15 minutes with the corresponding kinase inhibitor maintaining FCS deprivation. PI3K inhibitor LY294002 (15 μM), p44/42 ERKs inhibitors PD98859 (0.06 mM) or U0126 (20 μM) were obtain by Calbiochem, Ca. Afterwards, next fifteen minutes cells were in contact with FBS and without inhibitors. At the end of incubations, transfected cells were removed from the dishes and we obtained proteins or mRNA as convinced.
Athymic male nu/nu Swiss mice (Charles River Laboratory, Sta Perpetua, Spain) were injected subcutaneously (s.c.) as previously described , according to the protocols approved by the Institutional Animal Care and Use Committee. Tumours were measured periodically with a calliper, and the volume was calculated as length × width2 × 1/2 . Tumours were surgically removed and analysed when they reached a diameter of ~ 1 cm.
Cells and tissue samples were lysed with RIPA buffer (1% SDS 10%, 1% NP40, 0.5% de Sodium deoxycholate) plus protease inhibitors. Forty μg of each protein sample were subjected to 10% SDS/PAGE under reducing conditions, and transferred to polyvinylidene fluoride membranes (PVDF-Bio-Rad, Hercules, CA). Membranes were blocked in TBST buffer (0.9% NaCl, 0.02 M Tris (pH 7.5), 0.05% Tween 20; 5% skimmed milk; 1 h, RT), and probed with primary antibodies: anti-pan-Ras (Oncogene Research Products San Diego CA mouse monoclonal clone RAS10 cat: OP40), anti-HIF-1α (a generous gift from Dr. Edurne Berra, CICBiogune, Bilbao, Spain), anti-GLUT-1 (Abcam, Cambridge, UK, ref. Ab652), anti-VEGF-A (Neomarkers, CA, Ref. MS-35), anti-Sp1 (Santa Cruz, CA, ref. SC-59), anti-p-ERKs (Sigma-Aldrich Inc. Monoclonal Anti-Map Kinase activated clone MAPK-YT), anti-p-Akt (Cell Signaling Rabbit Phospho-Akt (Ser473) antibody #9271) and anti-Tubulin (Sigma-Aldrich Inc). Detection was performed using peroxidase-conjugated secondary antibodies. The resulting complexes were visualized by enhanced chemiluminiscence autoradiography (Amersham Life Science,Chicago,Il). Autoradiographs were quantified by scanning densitometry Quantity One Quantitation Software™ (Bio-Rad, Hercules, CA).
Expression levels of culture medium cells and tissue associated VEGF were also examined by enzyme linked immunosorbent assay (ELISA: Quantikine immunoassay kits; R&D Systems) according to the manufacturer’s instructions.
It was performed on paraffin-embedded tissues with VEGF (C-1) mouse monoclonal antibody (Santa Cruz Biotechnology Cat#sc-7269). We used anti-mouse DakoCytomation EnVision System HRP to visualize the reaction.
Trizol Reagent (Life Technologies Gibco) according to manufacturer’s instructions was used to total mRNA extraction. One μg of RNA was reverse-transcribed into cDNA using pdN6 primers using High Capacity Reverse Transcriptase (Applied Biosystems, Foster city CA). Subsequent Real-Time PCR reaction for Vegf-A mRNA levels was performed in duplo in the LightCycler® System SYBGreen480 (Roche) using primers: Rev: CACCTTCATGGGACTTCTGCT; Fwd: GCACTGGACCCTGGCTTTA. Angpt2 mRNA levels were assessed in duplo using inventoried TAQMAN gene expression assays Mm00545822 (Applied Biosystem, Foster city CA). Every gene expression quantification was corrected using three housekeeping genes: mitochondrial ribo-somal protein L19 Mn00452754; Hypoxanthine guanine phosphoribonyl transferase 1 Mn01545399 and Pepti-dylpropyl isomerase A Mn02342430. Threshold cycle data were analyzed using the following formula: ratio = ((E target ) ΔCP target (control–sample)) .
Cellular clones were cultured in 12 well plates and incubated 15, 30 and 180 minutes with Actinomycin D (5 μg/ml) just before RNA total extraction was performed using Trizol Reagent (Life Technologies Gibco) and following manufacturer’s protocol.
NIH3T3 stable KRAS clones were transfected with three different Vegf promoter constructions and a plasmid containing multiple HRE inserts that were a kind gift of Dr. Pages and Dra Berra [18–21]. Construct 1 contained the complete Vegf Promoter, construct 2 lacked AP2 and Sp1 response elements and construct 3 lacked HRE elements. Lipofectamine 2000 was used as the transfection agent according to manufacturer’s instructions. (Invitrogen Life Technologies S.A.). At 48 h post transfection, luciferase expression was determined using the Dual Luciferase Reporter System (Promega), and the relative luciferase value was determined after normalizing to the β-Gal control. Results are expressed as mean of three independent experiments.
HUVE cells (obtained from Advancell company) were deprived of growth factors for 24 h in 199 medium containing 0.1% Fetal Calf Serum. Cells were then stimulated 48 h with conditioned growth media from the different NIH3T3 clones, containing 0.25 μCi/ml (6-3H) thymidine (Amersham Pharmacia Biotech) (3 mM final concentration). After 20 hours of incubation, the incorporated radioactivity was counted by liquid scintillation. Results are expressed as a percentage of the (6-3H) thymidine incorporation in the presence of the medium of control NIH3T3.
One mg of the protein lysate from ASP13 or CYS12 NIH3T3 cell clones was added to 30 μg of glutathione S-transferase (GST)- RBD fusion protein. RBD, encompassing amino acids 51 to 131 of Raf-1 protein, is the minimal domain required for the binding of Ras-GTP. Presence of active Ras was detected by Western blotting with monoclonal anti-pan-Ras (Oncogene Research Products mouse monoclonal clone RAS10 cat: OP40). Autoradiography bands were quantified by scanning densitometry using Quantity One Quantitation Software™ (Bio-Rad Laboratories, S.A., Alcobendas, Madrid).
It was assessed using rat anti-PECAM (BD Pharmingene,NJ USA 550274/MEC13.3) staining on cryopreserved sections. MVD was evaluated as mean of the number of vessels in 5 hot spot fields at 400X. Only endothelial PECAM staining was considered .
Double immunofluorescence was used to evaluate Desmin(+)/PECAM(+). Rat anti-mouse monoclonal PEC-AM (BD Pharmingene,NJ USA 550274/MEC13.3) and rabbit anti-Desmin antibody (Lab Vision/NeoMarkers, CA, USA; RB-9014) were simultaneously incubated overnight at 4°C. Secondary goat anti-rat Alexa 546 and anti-rabbit Alexa 488 were incubated for 1 h. After washing, staining was evaluated in five hot spot fields at 400X. LEICA software counted the areas with Desmin and PECAM overlapping staining and also measured major diameters of the lumen of vessels. Correlation was performed with the presence of necrosis.
Rabbit anti α-Smooth muscle actin (Lab Vision/NeoMarkers, CA, USA; RB-9010, ready to use) staining was performed on paraffin-embedded tissue sections. Envision anti-rabbit (Dako Denmark A/S, Glostrup, DK) plus diaminobenzidine were used to visualize the reaction. Immunoreactivity was semiquantitatively evaluated by comparing the level of expression between distinct groups. Carbonic anhydrase IX (Ab15086,Abcam plc,. Cambridge, UK) and GLUT-1 (Ab652; Abcam plc. Cambridge UK) staining was also performed. Cytoplasmic positive cells were expressed as a percentage of total cells counted. For all antibodies, no staining was observed with negative control samples (absence of primary antibody, or incubation with an irrelevant antibody or IgG).
Next, we decided to evaluate the impact of this differential expression on two HIF-1α-dependent genes, GLUT-1 the ubiquitous glucose transporter protein, and VEGF-A . As observed in Figure 1C, and as expected from its more glycolytic phenotype, CYS12 mutant cells presented higher total levels of GLUT-1 as well as an increase in the glycosylated forms , when compared with ASP13 cells. Surprisingly, VEGF-A protein levels were higher in ASP13 cells than in CYS12 (Figure 1C).
Deletion of HRE within the VEGF-A promoter in all clones did not affect its activity. (Figure 3B Construct 2 vs Construct 1). These results further confirm the HIF-1α independent regulation of VEGF-A expression. In contrast, the selective deletion of SP1/AP2 response elements resulted in a significant decrease of VEGF promoter activity in both transfectants that was only significant to ASP 13 mutants (Figure 3B Construct 3 vs Construct 1).
Differences in latency period in vivo between ASP13 and CYS12 transfectants
Clone 4 (n = 10)
Clone 6 (n = 10)
Clone 2 (n = 10)
Clone 5 (n = 10)
Clone 4.3 (n = 10)
Clone 4.6 (n = 10)
Days ± SD
30,1 ± 6,5
24,5 ± 6,1
42 ± 0
9,7 ± 7,9
8.4 ± 0.9
13,2 ± 1,7
Mean ± SD
32.2 ± 8.9
10.4 ± 2.4
P < 0.04
In the context of KRAS-driven tumourigenesis, mutations located at codon 12 and 13 display distinct malignant potential and differentially regulate apoptosis, cell cycle [4, 15], or metabolic profiles . Here we show that minor differences in the molecular nature of KRAS mutations stimulate distinct intracellular signalling pathways in normoxic conditions with different impact in basal levels of HIF-1α VEGF-A production and generation of a distinct vascular network in tumours.
Upregulation of VEGF by the KRAS pathway has been previously shown . Here we show that cells expressing ASP13 KRAS mutant present higher levels of VEGF-A, the main pro-angiogenic gene induced by hypoxia, in the absence of high HIF-1α levels . In contrast, CYS12 mutants present a high glycolytic phenotype  through HIF-1α-dependent induction of glycolytic enzymes including GLUT-1 glucose transporter supporting the role of HIF-1α in switching to a glycolytic metabolism .
We have attempted to gain insight into the molecular mechanisms underlying the differential VEGF-A overexpression, apparently independent of HIF-1α in ASP13 clones, Our data support a direct transcriptional effect of ASP13 acting on VEGF-A promoter. This effect is mediated by a distinct activation of Raf-ERKs pathway and AP2/Sp1 elements within the proximal VEGF-A promoter. Of note it is independent of hypoxia-dependent elements and of PI3K activity. Extracellular signals that induce VEGF-A through this proximal region include, among others, growth factors such as EGF, insulin and PDGF in fibroblasts , prostaglandin E2 in human muscle cells , M-CSF in monocytes  and lysophosphatidic acid (LPA) in ovarian cancer cells . All of them affect promoter activity through modulation of at least Sp1 transcriptional activity. Noteworthy, Sp1 is also regulated by different signalling pathways including ERKs, PKA and PI3K-Akt [18, 30]. We have not detected changes in total Sp1 protein levels between ASP13 and CYS12 mutants, but other mechanisms with impact in the activity of this transcription factor could be implicated, such as acetylation, sumoylation, glycosylation or phosphorylation .
In our xenograft model, ASP13 xenografts consistently develop angiogenic sprouts of large diameter, invested by mural cells. These structures seem to be sufficient to support the increased utilization of the oxidative pentose phosphate pathway observed in the more benign ASP13 tumours . While development of these complex vascular structures may account for the initial delay observed in tumour growth, we speculate that they are able to support the very rapid growth occurring later . Nonetheless, the presence of significant tumour necrosis and less Carbonic anhydrase IX to hypoxic adaptation, observed in established ASP13 tumours may depict the relative insufficiency of this vascular tree . In contrast, histological analysis reveals that the more aggressive CYS12 tumours educe a dense endothelial-lined microvascular network that allows an early, steady and sustained tumour growth. This vascular strategy appears to be effective for these tumour cells that are more resistant to hypoxia, do not proliferate fast  and have relatively low energetic requirements associated with an increased anaerobic glycolysis.
The vascular pattern observed in ASP13 xenografts is in line with previous observations linking high VEGF-A levels with an increased diameter of newly forming vessels [32, 33]. The prominent stimulation of DNA synthesis in primary HUVECs by whole ASP13 conditioned medium, and in a less conspicuous manner by CYS12 supernatants, propose significant paracrine effects of tumour cell-derived VEGF-A in neovascularization . Also, ASP13 tumours vessels are covered with α-Sma(+)/Desmin(+) cells  further highlighting the contribution of VEGF-A to vessel maturation and tumour growth.
The retarded growth of ASP13 tumours harbouring elevated VEGF-A levels is consistent with reports challenging the concept that VEGF is just a positive angiogenic regulator. While angiopoietin2 levels did not show differences between transfectants, we cannot exclude a role of other angiogenic factors in differences observed between ASP13 and CYS12 tumoral vessels .
The impact of the genetic background of tumour cells on the angiogenic phenotype is relevant since they may have consequences regarding efficacy of specific antiangiogenic strategies. An evolving tumour with an ever-changing genetic background likely educes a dynamic vascular strategy that may escape to specific antiangiogenic treatment such as those targeting VEGFRs or its ligand . This is of importance now that more antiangiogenic drugs are being introduced to the clinical setting and there is a need for biomarkers that help in the selection of patients to be treated. KRAS mutations are used as negative predictors of antiEGFR therapies in colorectal cancer . The role of KRAS mutation as a predictive marker of bevacizumab-based treatment has been also explored. Indeed, better response rates to bevacizumab can be observed in KRAS wt colorectal tumors when compared to KRAS mutant . Of note, some authors have explored a potential differential behaviour of codon 13 mutant tumors with no conclusive results . It will of interest to explore in the adequate clinical setting whether our experimental observations correlate with clinical outcome in other tumor types such as colorectal cancer.
Mutations in the KRAS gene are among of the most prevalent in human tumours and they are known to have pleiotropic effects on tumour biology. The less aggressive ASP13 mutation, through Raf-Ras-ERKs activation of the VEGF-A promoter, creates a prominent VEGF-A-associated vascular network in the absence of high HIF-1α levels. This vascularisation is less effective than the dense microvascular network observed in CYS12 tumours. In our model system, we have shown that the molecular nature of KRAS mutations clearly influences the vascular strategy devised by the tumour cell. These observations provide us with a deeper insight of the complex role of major angiogenic regulators such as VEGF on tumour vasculature development and their relationship with oncogene activation.
This study was supported by grants from SAF 2009–07319, SAF2012-3363, FIS 01/1264 to Gabriel Capellà, FIS 03/0290 to Maria Antonia Arbós and from the Ministerio de Ciencia y Tecnología (SAF2010-20859), the Spanish Ministry of Health (RTICC RD2006-0092) and Generalitat de Catalunya (2009SGR283) to Francesc Viñals. The research team belongs to the Network of Cooperative Research on Cancer RD06/0020/0150 RD12/0036/0031 funded by the Instituto Carlos III, Ministerio de Sanidad y Consumo of Spain, F05-01 from the Fundació Gastroenterologia Dr. Francisco Vilardell, and Acción Transversal contra el Cáncer and Fundación Científica de la AECC.
Role of the funding source
Funding sponsors had no involvement in any scientific area to this work.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.