Abstract
Purpose: The Wnt/β-catenin (β-cat) signaling cascade is a key regulator of development, and dysregulation of Wnt/β-cat contributes to selected cancers, such as colorectal, breast, and hepatocellular carcinoma, through abnormal activation of Wnt target genes. To identify novel modulators of the Wnt/β-cat pathway that may emerge as therapeutic targets, we did an unbiased high-throughput RNA interference screen.
Experimental Design: A synthetic oligonucleotide small interfering RNA library targeting 691 known and predicted human kinases was screened in Wnt3a-stimulated human cells in a live cell luciferase assay for modulation of Wnt/β-cat–dependent transcription. Follow-up studies of a selected high-confidence “hit” were conducted.
Results: A robust quartile-based statistical analysis and secondary screen yielded several kinases worthy of further investigation, including Cdc2L1, Lmtk3, Pank2, ErbB3, and, of note, vascular endothelial growth factor receptor (VEGFR)1/Flt1, a receptor tyrosine kinase (TK) with putative weak kinase activity conventionally believed to be a negative regulator of angiogenesis. A series of loss-of-function, genetic null, and VEGFR TK inhibitor assays further revealed that VEGFR1 is a positive regulator of Wnt signaling that functions in a glycogen synthase kinase-3β (GSK3β)–independent manner as a potential synthetic lethal target in Wnt/β-cat–addicted colon carcinoma cells.
Conclusions: This unanticipated non-endothelial link between VEGFR1 TK activity and Wnt/β-cat signaling may refine our understanding of aberrant Wnt signaling in colon carcinoma and points to new combinatorial therapeutics targeted to the tumor cell compartment, rather than angiogenesis, in the context of colon cancer. (Clin Cancer Res 2009;15(24):7529–37)
Activation of Wnt signaling, which is observed constitutively in most colon cancers, targets >50 genes by a transactivator complex comprising β-catenin (β-cat) and T-cell factor. To uncover novel kinases that modulate Wnt signal transduction and are essential to the survival of cells addicted to Wnt/β-cat signaling, we systematically inactivated, one at a time, all human genes that direct the synthesis of kinases in cells stimulated with Wnt. Our screen identified an unanticipated mode of cross-talk between vascular endothelial growth factor receptor (VEGFR)1/Flt1 and the Wnt/β-cat pathway, two pathways of tremendous importance in cancer. Initial mechanistic studies provided evidence for VEGFR1/Flt1-linked tyrosine phosphorylation of β-cat, thereby modifying transcriptional activity. Thus, a mechanistically rational strategy for synthetic lethal therapy in the treatment of colon carcinoma is targeted suppression of VEGFR1-linked tyrosine kinase activity within the Wnt/β-cat–addicted tumor cell compartment per se, not endothelia.
The Wnt signal transduction cascade relays progrowth signals through β-catenin (β-cat), a key transcriptional coactivator whose cellular levels are tightly regulated. In nonactive states, serine/threonine phosphorylation of β-cat by glycogen synthase kinase-3β (GSK3β), the component kinase of the adenomatous polyposis coli–mediated destruction complex, results in β-cat polyubiquitination mediated by the E3 ligase β-TrCP and consequent proteasomal degradation. In the presence of activating Wnt ligands, the destruction complex is inhibited, and an active β-cat population accumulates in the cytoplasm and translocates into the nucleus, where it forms a transactivator complex with lymphoid-enhancing factor or T-cell factor, resulting in the transcription of >50 target genes (1, 2). Aberrant activation of the Wnt pathway resulting in ligand-independent signaling and constitutively high levels of Wnt target gene expression is common in colon cancers (1, 2). Although such Wnt/β-cat activation contributes to certain forms of cancer, it may also create therapeutically exploitable vulnerabilities within the context of synthetic lethality (3, 4). Two genes are classically synthetic lethal if cells with a mutation in either one of the genes are viable, but simultaneous mutation of both genes causes death or impairs cellular fitness (5, 6). The idea can be extended to common mechanisms of cancer transformation, such as gain-of-function oncogenes exemplified by aberrant Wnt/β-cat activation. In this regard, inhibition or mutation of a second gene may be lethal or antiproliferative only in the context of specific oncogenic upregulation (3). Thus, it may be possible to use the contextual difference of aberrant Wnt/β-cat activation to find selective targets for therapeutic intervention. Kinases, an important family of functional regulators, have emerged as therapeutically alluring targets because single kinase inhibition can often eliminate an entire stimulatory or inhibitory arm of a signal transduction cascade. We have generated a screening strategy that combines the study of Wnt pathway activation with loss of individual kinases to uncover candidate kinase modulators that may emerge as therapeutic targets in the context of Wnt-addicted colon cancer.
Materials and Methods
Cell lines and reagents
pCMV-FLuc has been previously described (7). pTOPFLASH and pcDNA3.1 were purchased from Upstate Products (Billerica, MA) and Invitrogen Corporation (Carlsbad, CA), respectively. STF293 reporter cells were kindly provided by Dr. Jeremy Nathans (Johns Hopkins University School of Medicine, Baltimore, MD) and cultured in STF293 medium. HEK293, HEK293T, and HeLa cells (ATCC, Manassas, VA) were cultured in DMEM supplemented with 10% heat-inactivated fetal bovine serum and 1% glutamine. SW480 and KM12L4a colon carcinoma cells (kindly provided by Dr. Loren Michel, Washington University School of Medicine, St. Louis, MO) were cultured in DMEM supplemented with 10% fetal bovine serum and 1% glutamine. L cells and L Wnt-3a cells were purchased from ATCC, and conditioned medium was collected according to the ATCC protocol. A HEK293 cell line stably expressing FLuc (293Luc) was established by transfecting 293 cells with pCMV-FLuc with the use of FuGENE 6 transfection reagent (Roche, Indianapolis, IN) for 48 h followed by selection of single colonies after continuous exposure to G418 (0.4 mg/mL). SB216763 (Sigma-Aldrich, St. Louis, MO), and VEGFR tyrosine kinase (TK) inhibitors II and III (EMD Biosciences, Gibbstown, NJ) were reconstituted in DMSO. Plasmids expressing short hairpin RNA (shRNA) sequences were provided by the Genome Sequencing Center, Washington University School of Medicine. Each sequence was provided in a pLKO.1 expression vector.
High-throughput screening
Small interfering RNA (siRNA) screening was done in black, clear-bottomed, 96-well culture plates (Corning 3904) with the use of a Beckman Coulter (Fullerton, CA) Core Robotics system, including an FX liquid handler, controlled by the Sagian graphical method development tool (SAMI scheduling software). STF293 cells (10,000 per well) in STF293 medium (DMEM/F12 supplemented with 10% heat-inactivated fetal bovine serum and 1% glutamine) at 100 μL/well were seeded 1 d before transfection. Plates were maintained in an environmentally controlled Cytomat incubator until use, thereby optimizing health and ensuring uniform treatment of all plates. Forward transfection was done with a 96-channel head on the FX liquid handler, adding 0.2 μL/well of media-complexed DharmaFECT 1 reagent (Dharmacon Research Inc., Valencia, CA) to the aliquoted siRNA library (Kinase siRNA set v2; Qiagen Inc., Valencia, CA) in a 96-well reaction plate. Experimental siRNA oligos were arrayed in columns 2 to 11 of each plate, and individual controls comprising mock-transfected wells, a nontargeting siControl sequence (Dharmacon Research Inc., Lafayette, CO), as well as a firefly luciferase–targeting siRNA sequence (Dharmacon Research Inc.) were placed manually in columns 1 and 12. After incubation of siRNA-DharmaFECT 1 complexes for 20 min at room temperature, 100 μL of the complexed siRNA was added to each well of a plate with cells (×3 plates) with the use of the FX liquid handler, yielding a final concentration of ∼25 nmol/L siRNA per well. Plates were maintained in the Cytomat for 24 h, after which each plate was aspirated with the use of the BioTek (Winooski, VT) ELx405 Select microplate washer. Wnt3a-conditioned media (200 μL/well) containing 150 μg/mL d-luciferin (Biosynth, Naperville, IL) was added with the use of the FX liquid handler, and cells were incubated for 10 min. The luminescent signal was measured under ultrasensitive detection mode on an EnVision plate reader (PerkinElmer, Waltham, MA) 10 min, and 6, 12, and 18 h post Wnt3a addition. Cell viability was then determined with resazurin dye (Sigma R7017; final concentration, 44 μmol/L) after 90-min incubation at 37°C, as monitored on a FLUOstar OPTIMA fluorescence reader (BMG Labtech, Cary, NC; excitation, 544 nm; emission, 590 nm).
Secondary shRNA screen
STF293 cells in DMEM/F12 media supplemented with 10% heat-inactivated fetal bovine serum and 1% glutamine were plated in 24-well or 96-well plates 1 d before transfection such that cells would be ∼70% confluent on the day of transfection. In quadruplicate, cells were transfected with 800 ng (for 24-well plates) or 300 ng (for 96-well plates) of shRNA expression vector. Four individual shRNA sequences (except three for Pank2) that targeted nonoverlapping segments of the coding region of each original high stringency “hit” gene were used. Also included were a nontargeting scrambled sequence (negative control) and a sequence targeting firefly luciferase (positive control). Plasmids were transfected with the use of Lipofectamine 2000 at a 2.5:1 ratio (volume of reagent to micrograms of DNA) in 100 μL (for 24-well plates) or 50 μL (for 96-well plates) of serum-free DMEM/F12 media. Transfection reagent was added to the serum-free media and incubated for 5 min, followed by addition of DNA. After incubation for 20 min, the transfection mixture was added dropwise to wells of assay plates. Cells were incubated in the transfection mixture for 4 h, and then the media was replaced with fresh DMEM/F12 media containing serum. After transfection, cells were incubated for 24 h followed by doubling of the volume in each well by addition of Wnt3a-conditioned media that contained 150 μg/mL of d-luciferin. Cells were incubated for another 24 h and then imaged on a bioluminescence imaging system (IVIS 100, Caliper Life Sciences, Hopkinton, MA; filter, open; f-stop, 1; exposure, 5 min; binning, 8). After bioluminescence imaging, an MTS cell viability assay was done (CellTiter 96 AQueous One Solution Cell Proliferation Assay, Promega, Madison, WI) by measuring the formation of the formazan product (absorbance at 490 nm, 1 h, 37°C) with the use of a FluoStar Optima plate reader for 24-well plates or an EL 311 plate reader (BioTek Instruments) for 96-well plates.
RT2 Profiler PCR Wnt Array
Embryonic stem cells from wild-type R1 or VEGFR1/Flt1(−/−) mice were provided by Dr. Guo-Hua Fong (University of Connecticut Health Center, Farmington, CT) and Dr. Kyunghee Choi (Washington University School of Medicine), and maintained at the Murine Embryonic Stem Cell Core at Washington University. Cells were cultured in DMEM containing FCS, l-glutamine, nonessential amino acids, HEPES, leukemia inhibitory factor, and 2-mercaptoethanol, and plated on gelatinized 10-cm tissue culture dishes to eliminate mouse embryonic fibroblasts before experimental assays. After incubation of embryonic stem cells in the presence or absence of rWnt3a, RNA was obtained with the use of the RNeasy RNA extraction kit (Qiagen), including the optional DNase digestion protocol. Reverse transcriptase-PCR array profiling was done with the use of the pathway-specific RT2 Profiler PCR Array system (SuperArray, Frederick, MD) according to the manufacturer's instruction, in an ABI 7300 real-time PCR system (Applied Biosystems, Foster City, CA). Data analysis with the use of the ΔCt method was done according to a data analysis template available online.1
Briefly, ΔCt values were calculated by subtracting the average cycle threshold (Ct) values for housekeeping genes (glyceraldehyde-3-phosphate dehydrogenase, β-actin, and Hprt1) from the Ct values for T (Brachyury). Fold differences were calculated as 2−ΔCt (rWnt3a treated)/2−ΔCt (untreated) for both wild-type and VEGFR1(−/−) embryonic stem cells.Semiquantitative reverse transcriptase-PCR
siRNA transfections were done as described above and allowed to proceed for 48 h. RNA was extracted with the use of a RNeasy RNA extraction kit (Qiagen). One-step reverse transcriptase-PCR was done according to the manufacturer's protocol with the use of the Access RT-PCR system (Promega) on a GeneAmp PCR system 2700 (Applied Biosystems). The primer sequences were as follows: VEGFR1 extracellular domain, forward primer TGGGACAGTAGAAAGGGCTT, reverse primer GGTCCACTCCTTACACGACAA; glyceraldehyde-3-phosphate dehydrogenase, forward primer GTGAAGGTCGGAGTCAACGG, reverse primer TGATGACAAGCTTCCCGTTCTC. Amplified products were subjected to agarose gel electrophoresis, and visualized as a 384-bp product for VEGFR1 and a 200-bp product for glyceraldehyde-3-phosphate dehydrogenase with the use of a ChemiImager (Alpha Innotech Corp., San Leandro, CA).
In cellulo transfection assays
In STF293 cells, follow-up siRNA validation studies were done with the use of a forward transfection protocol with DharmaFECT 1 reagent (Dharmacon Research Inc.) at a final concentration of 50 nmol/L siRNA (ON-TARGETplus SMARTpool, L-003136-00; Dharmacon Research Inc.), and 20 nmol/L siRNA (individual duplexes; Dharmacon Research Inc.). In other experiments, cells were transfected with pTOPFLASH with the use of Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. Bioluminescence was captured with the use of an IVIS 100 or IVIS 50 imaging system, and Living Image (Caliper Life Sciences) and Igor (WaveMetrics, Lake Oswego, OR) software. The acquisition variables were as follows: no filter; exposure time, 1 to 2 min; binning, 8; f-stop, 1; field of view, 12 (IVIS 50) or 10 cm (IVIS 100).
MTS assays
STF293, HeLa, KM12l4a, and SW480 (1,000 cells per well) were plated in appropriate media in clear 96-well plates. Numerous wells containing medium only were also plated as background controls. Each cell line was treated with increasing drug concentrations in triplicate, whereas STF293 cells were also treated with Wnt3a together with increasing drug concentrations. Drug treatment proceeded for 72 h, after which MTS reagent (Promega) was added and incubated for 1 to 3 h depending on the cell line, and the absorbance was assayed at 490 nm in an EL 311 plate reader (BioTek Instruments). Background subtraction was done from medium-only wells that had each been treated with the corresponding drug concentrations in duplicate to calculate the corrected absorbance. Data were then normalized as fold drug-untreated values for each cell type to facilitate plotting on the same graph.
Immunofluorescence and confocal microscopy
Cells were transfected (for siRNA assays) or treated (for Inh II assays) in 35-mm glass coverslip dishes (MatTek Cultureware), and treated with Wnt3a or control media for 2 h before fixation in 4% paraformaldehyde solution. Samples were then washed, treated with ice-cold methanol for 10 min at −20°C, washed, and incubated with anti–β-cat antibody (BD Biosciences; 1:500 dilution) overnight at 4°C on a rocking platform. After three washes, a 2-h incubation in Alexa 488–conjugated secondary antibody (Invitrogen Corporation; 1:500 dilution) was done at room temperature in the dark. Samples were incubated with 10 μg/mL of 4′, 6-diamidino-2-phenylindole for 15 min in the dark and stored in PBS at 4°C protected from light. Images (10 μm slice thickness; ×40 magnification) were acquired with the use of an inverted Zeiss Axiovert 200 laser scanning confocal microscope coupled to a Zeiss LSM 5 PASCAL fitted with excitation lasers at 488 nm (for β-cat) and 405 nm (for 4′, 6-diamidino-2-phenylindole), and long pass emission filters at 505 nm (for β-cat) and 420 nm (for 4′, 6-diamidino-2-phenylindole).
Immunoprecipitation and Western blotting
Cells were transfected (for siRNA assays) or treated (for Inh II assays) in 10-cm tissue culture–treated dishes, and harvested in immunoprecipitation buffer [50 mmol/L HEPES (pH 7.5), 150 mmol/L NaCl, 5% glycerol, 1% Triton X-100, 1 mmol/L EDTA, 10 mmol/L NaF, and 30 mmol/L glycerol phosphate] freshly combined with a 1:100 dilution of protease inhibitor cocktail (Sigma-Aldrich) and a 200-μmol/L final concentration of sodium orthovanadate. Protein amounts were quantified with the use of a bicinchoninic acid assay (Pierce, Rockford, IL) according to the manufacturer's protocol. Protein (500 μg) from each sample was incubated for 1 h with anti–β-cat antibody (BD Biosciences; 0.5 μg) on a rocking platform and then complexed with Protein A agarose beads (Calbiochem, La Jolla, CA) overnight at 4°C on a rotator. Samples were washed thrice after centrifugation at 6,000 rpm and boiled for 10 min before separation by SDS-PAGE. After transfer to a polyvinylidene difluoride membrane, blocking was allowed to proceed for 2 h at room temperature, followed by incubation with anti-phosphotyrosine antibody (clone 4G10; Upstate Products; 1:1,000 dilution) overnight at 4°C. Immune complexes were detected by horseradish peroxidase–labeled secondary antibodies and enhanced chemiluminescence reagent (Amersham Biosciences). Exposure and densitometric analysis were done with the use of the IVIS imaging system, and Living Image (Caliper Life Sciences) and Igor (WaveMetrics) software. A similar SDS-PAGE and immunoblotting protocol was followed for detection of β-cat (1:200 dilution), VEGFR1 (1:200 dilution), and actin (1:1,000 dilution) with the use of antibodies from Santa Cruz (Santa Cruz, CA), Chemicon Corporation (Billerica, MA), and Sigma-Aldrich, respectively.
Statistical analysis and hit selection
In the primary siRNA screen, the luminescence data acquired 18 h post Wnt3a treatment were normalized to cell viability data (expressed as value = ×). Median (Q2), first (Q1), and third (Q3) quartile values were calculated for all normalized values and subjected to plate-by-plate analysis. To facilitate experiment-wide analysis, further normalization was done with a control nontargeting siRNA placed on each individual assay plate. Q1, Q2, and Q3 values were calculated as described above. Averages and SE calculations were done for triplicates of log(x/median) values. Upper and lower boundaries were calculated as Q3 + 2c(Q3 − Q2) and Q1 − 2c(Q2 − Q1), respectively, for c = 1.7239, corresponding to a high stringency targeted error rate (α = 0.0027), and c = 0.9826, corresponding to a lower stringency targeted error rate (α = 0.046; ref. 8). High stringency hits were chosen as those targets that scored α ≤ 0.0027 in both plate-by-plate and experiment-wide analyses. Selected targets clustered within the first two plates (numbers 1-150; Fig. 1, red circles) represented hits scored as significant based on the high stringency criteria in an experiment-wide analysis only and not in a plate-by-plate analysis (listed in Supplementary Table S3). This ensued from the alphabetical clustering and nonrandom placement of siRNA sequences on the plates by the vendor, resulting in over-representation of hits in the first two plates, thereby infringing on the sparse-hit hypothesis, which assumes that hits are randomly spaced throughout the dataset.
Experimental statistical analyses
Student's t tests, two-way ANOVA, curve fitting, box plots, and IC50 calculations were done with the use of GraphPad software (GraphPad, La Jolla, CA).
Results
Kinase-targeted high-throughput siRNA screen of Wnt/β-cat–dependent transcription
We used a synthetic oligonucleotide siRNA library targeting known and predicted human kinases in a high-throughput screening format for systematic yet unbiased target discovery in contexts that were both physiologically relevant and therapeutically alluring (3, 4). HEK293 cells stably expressing a luciferase transcriptional reporter of Wnt/β-cat target gene activation (STF293 cells; ref. 9) were transfected with 1,382 siRNA duplexes targeting 691 distinct genes as a pool of 2 duplexes per target in a 96-well high-throughput screen (Supplementary Fig. S1), and subsequently treated with Wnt3a-conditioned media (Wnt3a; ref. 10) to emulate pathway activation. Human cells, together with a functionally targeted synthetic oligonucleotide library, were used to maximize the discovery of potential disease-relevant candidates and minimize off-target effects, such as those that can be observed in genome-wide screens conducted in model organisms (11, 12). A highly stringent targeted error rate in a quartile-based analysis (ref. 8; see Materials and Methods) was used to uncover candidates of statistical relevance (Fig. 1). Eleven kinases (Cdc2L1, Pank2, Lmtk3, Flt1, ErbB3, Fgr, Mapk10, Phkg1, Phkg2, Pbk, and Stk24), representing both serine/threonine and TKs, were recovered as hits, all of which were kinases not previously associated with the Wnt pathway (Supplementary Table S1). Importantly, analysis at lower stringency revealed additional kinases, several of which have known Wnt pathway regulatory functions (transforming growth factor-β receptor, NLK, CK1/2, CDKN2a, AKT2/3, and GSK3α/β), which served to authenticate the screen (Supplementary Fig. S2; Supplementary Table S2). Also, a subset of kinases enriched in cell cycle regulators emerged from an experiment-wide analysis, including the noteworthy genes Cdk8 (13), ATM, and Plk3 (Supplementary Table S3). The loss of one subset of kinases (Cdc2L1, VEGFR1/Flt1, Pank2, and Lmtk3) resulted in reduced β-cat–dependent transcription (Fig. 1), representing potential targetable therapeutic genes.
Validation of hit genes with the use of shRNA knockdown
To validate the changes in Wnt/β-cat signaling that were observed in the primary screen, we rescreened in STF293 cells the 11 hit genes with the use of four separate shRNA sequences corresponding to each gene of interest. Plasmids bearing shRNA were chosen for the secondary screen because of their potential for higher potency, fewer off-target effects, and stability. The shRNA sequences were selected for those that target nonoverlapping segments of the coding regions of hit genes and were assessed individually as opposed to being pooled, as was done in the primary siRNA screen. Genes were qualified as reproducible if, in at least two of three separate experiments, ≥2 shRNA sequences for a target gene increased or decreased luciferase signals above or below the median of a nontargeting negative control shRNA in concordance with the trend previously observed with siRNA (see Fig. 2, which illustrates a representative experiment). Overall, shRNA against 5 of the original 11 hit genes showed the same effect on luciferase signal as was observed in the primary screen, and were deemed reproducible: Cdc2L1, Lmtk3, Pank2, ErbB3, and VEGFR1/Flt1. Only VEGFR1/Flt1 was concordant in all three independent experiments.
Characterization of VEGFR1/Flt1 as a “high-confidence hit”
Direct therapeutic applications may be found for candidates whose loss-of-function by RNA interference results in attenuated Wnt/β-cat signaling, especially in the context of “oncogene-addicted” colon cancer cells, in which the loss of aberrantly active Wnt/β-cat signaling results in decreased cell survival (14). In this regard, the specific identification of VEGFR1/Flt1, but not VEGFR2/KDR or VEGFR3/Flt4, as a “reproducible hit” was striking. VEGFR1, a receptor TK with putative weak kinase activity, is conventionally thought to function as a negative regulator of angiogenesis (15, 16), an endothelial cell process that is actively targeted in the treatment of colorectal cancers (17). To extend the cellular link between VEGFR1 function and Wnt/β-cat signal transduction, we further validated the effects of VEGFR1 knockdown on the β-cat–dependent transcription observed in the primary screen. With the use of STF293 cells, a pool of four independently designed sequences against VEGFR1, as well as three of these four when tested as individual siRNA sequences, resulted in knockdown of VEGFR1 (Fig. 3A) and again attenuated β-cat–dependent transcriptional activation (Fig. 3B and C), confirming the primary screen. No effect was observed in control HEK293-Luc cells (293Luc), a line expressing luciferase from a constitutive non–β-cat–responsive cytomegalovirus promoter (Fig. 3D).
To independently confirm that the absence of VEGFR1 resulted in attenuated Wnt signaling at the level of transcriptional activation of Wnt/β-cat targets, we obtained embryonic stem cells derived from a VEGFR1(−/−) mouse (18). Importantly, quantitative reverse transcriptase-PCR–based profiling for genes implicated in Wnt signaling was done in both wild-type and VEGFR1(−/−) embryonic stem cells. Embryonic stem cells lacking VEGFR1 and treated with rWnt3a showed a mere 5-fold increase in expression of T (Brachyury), a specific transcriptional target of Wnt3a in embryonic stem cells [ref. 19; 2−ΔCt values: 4.5e-02 (treated) versus 8.6e-03 (untreated)]. By comparison, in rWnt3a-treated wild-type embryonic stem cells, a robust 44-fold increase in T was observed [2−ΔCt values: 6.2e-02 (treated) versus 1.4e-03 (untreated)]. These data provided strong evidence that VEGFR1 enhanced Wnt signaling in embryonic stem cells.
Inhibition of VEGFR1 TK activity inhibits Wnt/β-cat–dependent transcription and is synthetic lethal in colon cancer cells
Recent data indicate that the expression of VEGFR1 in primary colon carcinoma tumor samples correlates with disease progression (20). However, the mechanistic basis of this observation has not been elucidated. Thus, to directly test the effect of abrogating VEGFR TK activity on cell survival, SW480 and KM12L4a cells, colon cancer lines that show constitutive β-cat activation (21), were treated with VEGFR TK inhibitor II (ref. 22; Inh II) for 72 hours (Fig. 4A). An Inh II–induced decrease in the MTS cell viability assay was seen in these Wnt/β-cat–addicted colon cancer cells in a concentration-dependent manner. The observed IC50 of 336 nmol/L (for SW480) and 120 nmol/L (for KM12L4a) matched the known IC50 for inhibition of VEGFR1 TK activity (180 nmol/L), a value 9-fold greater than the IC50 for inhibition of VEGFR2 TK activity (20 nmol/L; ref. 22). In contrast, a 72-hour MTS assay with STF293 cells, a line that exhibits normal Wnt signaling, was unaffected by Inh II even in the presence of Wnt stimulation. Furthermore, cervical carcinoma–derived HeLa cells, a line that does not exhibit aberrant Wnt/β-cat pathway activation (23), also showed no response to Inh II (Fig. 4A). These data suggested that VEGFR TK inhibition was synthetic lethal to cells with aberrant Wnt/β-cat signaling, a process essential for colon cancer cell survival (24), but not to cells with normal Wnt/β-cat signaling. However, further analysis of the luciferase transcriptional reporter showed that treatment with Inh II still resulted in a concentration-dependent inhibition of canonical β-cat–dependent transcription in STF293 cells (Fig. 4B) but only in the presence of Wnt3a, as confirmed by transient transfection of HEK293T cells with pTOPFLASH, a similar β-cat–dependent bioluminescence reporter to that used in STF293 cells (Fig. 4B, inset). Finally, in SW480 and KM12L4a colon cancer cells transiently transfected with the pTOPFLASH reporter, constitutively active β-cat–dependent transcription was also attenuated upon treatment with Inh II (Fig. 4C). These effects were recapitulated with VEGFR TK inhibitor III, a second inhibitor targeting VEGFR TK activity (ref. 25; Supplementary Fig. S3A and B), but no effects were observed in off-target transcriptional reporter cells (293Luc) as expected (data not shown). Thus, whereas inhibition of VEGFR1 TK activity blocked Wnt-induced β-cat–dependent transcription in normal cells as well as constitutively active β-cat–dependent transcription in colon cancer cells, it was only in colon cancer cells that VEGFR1 TK blockade was lethal. Overall, with the caveat that the pharmacologic specificity of the tested inhibitors is not fully resolved, these data showed that by targeting VEGFR1 TK activity, the therapeutically exploitable vulnerability of Wnt/β-cat–addicted colon carcinoma was exposed.
Wnt pathway modulation by VEGFR1 is independent of GSK3β-dependent β-cat stabilization
Mutations in APC or β-cat are prevalent in colon carcinoma, resulting in constitutively stabilized β-cat that is no longer a substrate for GSK3β (1, 2). Concordantly, in STF293 cells in which β-cat was stabilized by pharmacologic treatment with the GSK3β inhibitor SB216763 (10 μmol/L; ref. 26), knockdown of VEGFR1 with duplexes of siRNA yielded a similar attenuation of β-cat–dependent transcription (Fig. 5A) as VEGFR1 knockdown in the context of Wnt3a stimulation (Fig. 3C). Thus, the effect of VEGFR1 was independent of GSK3β.
Downstream processing
VEGFR1 did not seem to impact the processing of total β-cat protein levels. For example, Western blot analysis of STF293 cells showed competent stabilization of β-cat protein upon stimulation with Wnt3a, even when effective knockdown of VEGFR1 protein was achieved with siRNA (Fig. 5B) or VEGFR TK activity was blocked (Fig. 5C). Overall, these data would suggest that therapeutic interventions directed at VEGFR TK activity may be mechanistically impacting downstream processing nodes of the Wnt/β-cat pathway within cells, such as nuclear translocation or transcriptional activation.
We thus explored the effect of VEGFR1 on nuclear translocation of β-cat with the use of immunofluorescence microscopy. In STF293 cells, both VEGFR1 knockdown with siRNA (Fig. 6A) or pharmacologic inhibition of VEGFR TK activity with Inh II (Fig. 6B) resulted in normal Wnt-induced β-cat accumulation and nuclear translocation. These data concurred with observations in ApcMin/+ mice that showed no change in nuclear β-cat upon treatment with the VEGFR TK inhibitor AZD2171 (27). Together, the above observations suggested that regulation of Wnt/β-cat signaling by VEGFR1 occurred at the level of post-translational modification of β-cat as a transcriptional coregulator and/or subsequent transcriptional activation of a selected subset of β-cat targets.
Preliminary support for VEGFR1-linked effects on post-translational modification of β-cat was provided by analysis of KM12L4a and SW480 colon cancer cells. Inhibition of VEGFR TK activity by treatment with Inh II, followed by immunoprecipitation and western blot analysis of β-cat in both lines, showed a dramatic decrease in tyrosine phosphorylated β-cat. However, β-cat total protein levels, which are canonically governed by phospho-serine/threonine–dependent proteasomal degradation, were found to be unchanged (Fig. 6C). Thus, inhibition of VEGFR1 TK activity resulted in decreased tyrosine phosphorylation of β-cat, thereby abrogating Wnt/β-cat–dependent pro-proliferative transcriptional activity in colon carcinoma cells per se.
Discussion
A systematic yet unbiased screening approach in the context of activated Wnt signaling in human cells has enabled the discovery of novel candidate kinase regulators of the Wnt/β-cat pathway. Whereas the Wnt pathway has been extensively pursued with the use of model organism-based genetic screens, recent observations have found newer positive functions of two classic negative regulators of the Wnt pathway, APC (28) and GSK3β (29), highlighting the inherent epistatic bias of conventional assays (30). By interrogating Wnt signaling at the level of transcriptional activation in human STF293 cells with the use of targeted RNA interference, we hoped to identify novel candidate repressors as well as activators of the Wnt/β-cat pathway, which could be further investigated in cell-specific contexts. Additionally, when establishing a screening strategy, there may be concern that the multitude of altered functional pathways and additional dysregulated kinases in colon cancer cell lines could confound the discovery of Wnt pathway–specific regulators. For a conventional synthetic lethal assay, we reasoned that the likelihood was high that a siRNA library screen would identify candidates that simply killed colon cancer cells through other pathways and hence would give rise to off-target effects instead of direct effects on Wnt-dependent transcription. Thus, starting with STF293 cells, we showed that one identified candidate, VEGFR1/Flt1, was a regulatory target with clinical relevance by examining Wnt/β-cat–dependent transcription in a wide variety of contexts, including genetic-null embryonic cells from a VEGFR1−/− mouse, loss-of-function by siRNA and shRNA, VEGFR TK inhibition in Wnt-responsive cells, as well as survival of Wnt/β-cat–addicted colon carcinoma cells.
Our data provide potential mechanistic linkages in support of prior observations. For example, VEGFR1 is expressed in various colon cancer lines that show activated Wnt signaling, including SW480 and KM12L4a cells (21). In addition, preclinical studies of a hexapeptide specifically targeting VEGFR1 (31), as well as the VEGFR TK inhibitors SU11248 (32), AZD2171 (27), and CHIR-258 (33), have shown efficacy in various models of colon cancer; however, whether the mechanism involved therapeutic “anti-angiogenesis” remains unclear (21). Furthermore, we found that VEGFR1-linked regulation of β-cat was found to be independent of GSKβ activity and β-cat nuclear translocation while impacting β-cat tyrosine phosphorylation, thus providing preliminary insight into a mechanism for the antiproliferation action of VEGFR1 inhibition. It is acknowledged that the exact site of tyrosine phosphorylation remains to be mapped, and the potential differential roles of soluble VEGFR1, a splice variant (34), versus full-length VEGFR1 need to be resolved. Nonetheless, the implication that VEGFR1 links TK activity, directly or indirectly, to a node of the Wnt pathway downstream of commonly occurring oncogenic mutations tremendously enhances its potential as a therapeutic target in the context of aberrant Wnt signaling in cancer.
Our siRNA-based high-throughput screen, which uncovered an unanticipated tumor cell linkage between VEGFR1 and Wnt/β-cat, extends the context-dependent cross-talk identified between the oncogenic Flt3 and Wnt pathways in acute myeloid leukemia (35), to VEGFR1/Flt1 and Wnt/β-cat in colon carcinoma. Indeed, comparative genomic hybridization analysis of metastatic human colon cancer shows gain of chromosomes 1p, 8q, 13q, and 20q, including a distinct amplicon ranging from 13q11 to 13q21 (36), a region that overlaps with the exact cytogenetic location of VEGFR1. Thus, our results provide an alternative tumor cell–centric basis for the partial successes of anti-angiogenic therapy being delimited to only a few cancer types, such as colorectal cancer, a Wnt/β-cat–addicted tumor type. Future studies will be pursued to confirm this in clinical specimens in the context of targeted therapy. Potentially, a mechanistically rational strategy for synthetic lethal therapy in the treatment of colon carcinoma is selective targeted suppression of VEGFR1 TK activity in the Wnt/β-cat–addicted tumor cells themselves.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Acknowledgments
We thank Dr. Dustin Maxwell for assistance with immunofluorescence microscopy, Dr. James Hsieh for assistance with reverse transcriptase-PCR profiling, and our colleagues at the Washington University Molecular Imaging Center for discussion.
References
Competing Interests
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.