Expression of DOCK1 Is Associated with an Adverse Clinical Outcome for HER2⁺ and Basal Breast Cancer Patients.

Immunohistochemistry (IHC) was performed on a panel of breast tumors to assess if DOCK1 is expressed in this disease. DOCK1 protein was detectable in the tumor epithelial cells of 144/145 samples (Fig. S1A). We examined the relationship between DOCK1 mRNA expression and the probability of survival of breast cancer patients according to clinical subtypes using a large set of microarray data linked to clinical outcome (n = 6,327, including 3,466 patients for whom survival data were available) (14). DOCK1 was broadly expressed but higher in estrogen receptor (ER)⁺/HER2⁻ high proliferation (luminal B) subtype (Fig. S1B). As shown in Kaplan–Meier survival curves of patients classified based on DOCK1 expression, high levels of its expression was significantly associated with poor disease-free survival in HER2⁺ (*P = 0.048) and ER⁻/HER2⁻ (*P = 0.03) subtypes (Fig. 1A); association with survival was not observed for the luminal subtypes (ER⁺/HER2⁻ low and high proliferation tumors; Fig. S1C). These analyses indicated that DOCK1 expression is linked to survival of patients afflicted with aggressive breast cancers and provided the impetus for investigating the mechanism(s) whereby this GEF contributes to breast cancer.

• Download figure
• Open in new tab
• Download powerpoint

Fig. 1.

DOCK1, a negative prognostic factor for human breast cancer survival, activates Rac and promotes cell migration downstream of HER2. (A) High levels of DOCK1 expression are associated with a poor prognostic for HER2⁺ and basal-like breast cancer patients. (B) DOCK1 is found in a complex with HER2 upon treatment of T47D cells with HRG (n = 3). (C) DOCK1 is phosphorylated on Y¹⁸¹¹ in T47D breast cancer cells upon treatment with 20 ng/mL HRG for 15 min (n = 5). (D) DOCK1 is phosphorylated on Y¹⁸¹¹ by SRC kinase. In vitro kinase assay was performed using recombinant SRC and purified GST or GST-DOCK1^1228–1865 proteins (n = 3). (E and F) Pharmacological inhibition of DOCK1 impairs HRG-mediated Rac activation and migration in T47D cells. (E) Cells were treated with 20 ng/mL HRG for 15 min in the presence of 100 μM CPYPP and active Rac was measured by precipitation with the purified p21-Binding Domain of PAK protein kinase expressed as a GST fusion (GST-PAK-PBD) (n = 5). (F) Migration of T47D cells toward 20 ng/mL HRG was measured in a Boyden chamber assay in the presence or absence of 100 μM CPYPP (n = 7). (G) Migration of T47D cells transfected with 100 µM NON-targeting or ON-target smart pool DOCK1 siRNA toward 20 ng/mL HRG was measured in Boyden chamber assay (n = 7). (H) Expression profiling uncovers DOCK1 as a highly expressed RacGEF in HER2-induced mouse mammary tumors. Data were expressed as mean ± SD (n = 4). One-way ANOVA followed by a Bonferroni test calculated the P values; *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, and **** P ≤ 0.0001.

DOCK1, a Target of HER2, Promotes Heregulin-Mediated Rac Activation.

Because DOCK1 binds RTKs and activate Rac during cell motility (10, 12, 13), we examined whether it can enter in a protein complex with HER2. In human ductal breast epithelial tumor cell line T47D, endogenous HER2 coimmunoprecipitated with DOCK1 only when cells were treated with Heregulin β1 (HRG) (Fig. 1B). We also expressed Flag-DOCK1 and an oncogenic and activated form of rat HER2 (NeuNT) and performed coimmunoprecipitation assays. A robust coprecipitation of NeuNT was detected in Flag-DOCK1 immunoprecipitates (Fig. S2A). Treatment of T47D cells with HRG, or overexpression of NeuNT, resulted in the phosphorylation of DOCK1 on a regulatory site, Y¹⁸¹¹, known to promote its RacGEF activity (Fig. 1C and Fig. S2B) (12). Phosphorylation on Y¹⁸¹¹ was not directly regulated by HER2 but instead mediated by SRC protein kinases (Fig. 1D and Fig. S2 C and D). HRG-mediated activation of HER2 is reported to activate Rac and increase cell motility (8). To examine whether DOCK1 can facilitate HRG-induced Rac activation and cell migration, its GEF activity was inhibited via a small molecule, 4-[3′-(2″-chlorophenyl)-2′-propen-1′-ylidene]-1-phenyl-3,5-pyrazolidinedione (CPYPP) (15). Treatment of T47D cells with HRG induced a ninefold increase in Rac activation that was blocked in the presence of the inhibitor (Fig. 1E). Whereas treatment of T47D cells with HRG induced cell migration, this effect was abrogated either in the presence of CPYPP or when DOCK1 expression levels were reduced using RNAi (Fig. 1 F and G). These results demonstrate that DOCK1 is activated and mediates Rac activation and migration downstream of HER2.

Quantitative Expression Profile of Rho GTPases and Their Regulators in Tumors.

Based on the clinical data and mechanistic insights above suggesting that DOCK1 is operating downstream of HER2, we aimed to investigate its function in an in vivo model of HER2 breast cancer. We first examined the major Rho GTPases, GEFs, GTPase Activating Proteins (GAPs), and guanine nucleotide dissociation inhibitors (GDIs) expressed in HER2 mouse breast tumors. mRNA deep sequencing was used to generate a quantitative transcriptome of mouse mammary tumor virus (MMTV)-Neu (variant NDL2-5) transgenic tumors (n = 4). The most expressed GEFs included LARG, DOCK9, DOCK8, TIM, α-PIX, p190RHOGEF, FARP1, and DOCK1 (Fig. S3A). α−PIX, DOCK1, β−PIX, TIAM1, and DOCK7 are the highest expressed RacGEFs (Fig. 1H). The partner of DOCK1, ELMO1, was highly expressed in comparison with ELMO2/3 in these tumors (Fig. S3B). Moreover, Rac1 was found to be the most expressed Rho GTPase in tumors (Fig. S3C). In term of negative regulators, RhoGDIs are minimally expressed (RhoGDIβ/γ; RhoGDIα was not detected) and the RhoGAP ARHGAP31/CDGAP stood out as the majorly expressed member (Fig. S3 D and E). These data provide a snapshot of the Rho regulation system and confirmed that DOCK1 is highly express in HER2 tumors.

DOCK1 Contributes to Tumor Growth.

To examine the role of DOCK1 in HER2-induced tumorigenesis in vivo, a conditional floxed (flx) mutant mouse line of DOCK1 was generated (Fig. S4 A–E). Before carrying out tumorigenesis studies, we investigated whether deletion of DOCK1 in mammary glands would have adverse effects on the development of the tissue. Whole-mount outgrowth analyses at 9, 12, and 15 wk established that DOCK1 is dispensable for mammary development (Fig. S4F). Additionally, MMTV-Cre⁺DOCK1^flx/flx mice were capable of nursing their offsprings to maturity, suggesting that loss of DOCK1 does not impair mammary gland function.

To assess the role of DOCK1 in a HER2 model of breast cancer, Mouse Mammary Tumor Virus (MMTV)-NeuNDL2-5-Internal Ribosome Entry Site (IRES)-Cre (NIC) transgenics were intercrossed with DOCK1^flx mice and tumor progression was examined in females. In this model, expression of NeuNDL2-5 and Cre recombinase are coupled and females develop ductal carcinoma in situ that progress to invasive carcinoma and lung metastasis (16). Three cohorts of mice were generated (NIC⁺DOCK1^wt/wt, NIC⁺DOCK1^wt/flx, and NIC⁺DOCK1^flx/flx), and animals were monitored weekly for apparition of tumors that were next allowed to grow for 5 wk. Mice developed mammary tumors irrespective of the genotype and a small delay of tumor onset in NIC⁺DOCK1^wt/flx and NIC⁺DOCK1^flx/flx animals was observed (Fig. 2A). Histological analyses demonstrated that tumors from all genotypes displayed a solid adenocarcinoma appearance (Fig. 2B). We further confirmed that Cre-mediated recombination of the DOCK1 floxed allele was efficient in tumors (Fig. S4E). Furthermore, genetic ablation of DOCK1 in mammary tumors reduced the levels of DOCK1 protein as verified by Western blotting and IHC, whereas it had no impact on the expression of Neu (Fig. 2 C and D). Whereas mammary tumor initiation is normal in NIC⁺DOCK1^flx/flx animals, they presented fewer tumor nodules as well as a significant decrease in their cumulative tumor burden after 5 wk (Fig. S5A and Fig. 2E). We observed fewer large tumors in DOCK1-null mammary glands, whereas the small and medium formed equally in all genotypes (Fig. S5B). A similar coverage of CD31⁺ blood vessels in DOCK1-deficient and wild-type NIC tumors was observed, ruling out a defect in angiogenesis to explain the growth defect (Fig. S5C). We next investigated an early time point of tumor onset by quantifying mammary intraepithelial neoplastic lesions (MINs) on inguinal mammary glands that did not develop tumors and noted a significant reduction in the DOCK1-deficient mammary glands (Fig. 2 F and G).

• Download figure
• Open in new tab
• Download powerpoint

Fig. 2.

DOCK1 contributes to HER2 tumorigenesis in a mouse model of breast cancer. (A) Kaplan–Meier analysis of tumor onset. NIC⁺DOCK1^wt/flx and NIC⁺DOCK1^flx/flx mice display a minor delay in tumor onset. (B) H&E staining showing that all animals develop tumors with an adenocarcinoma histology. (Scale bar: 100 µm, 20×.) (C) Western blot analyses demonstrating the absence of DOCK1 and equal Neu expression in tumor lysates. (D) IHC analyses showing Cre (Upper) and DOCK1 (Lower) expression in NIC⁺DOCK1^wt/wt and NIC⁺DOCK1^flx/flx tumors. (Scale bar: Upper, 250 µm, 10×; Lower, 50 µm, 40×.) (E) DOCK1 contributes to HER2-induced tumor growth. Average cumulative tumor burden per animal (in cubic centimeters). (F) NIC⁺DOCK1^flx/flx mice develop fewer mammary intraepithelial neoplasic lesions (MINs). Quantification of mammary MINs per tumor-free mammary gland 5 wk after tumor onset is presented. (G) Representative carmine red staining of mammary glands showing MINs in all genotypes. (Scale bar: 1 mm, 2.5×.) (H) DOCK1 promotes cell proliferation in HER2 tumors. Average percentage of Ki67-positive cells in tumors. (I) DOCK1 exerts an antiapoptotic role in HER2 tumors. Average percentage of activated caspase-3–positive cells in tumors. (J) Impaired c-JUN activation in DOCK1-null HER2 tumors. IHC analyses showing pc-JUN levels in MINs. (Scale bar: 100 µm, 20×.) One-way ANOVA followed by a Bonferroni test calculated the P values; *P ≤ 0.05 and **P ≤ 0.01.

The proliferative and apoptotic status of DOCK1-null mammary tumor cells was assessed to gain mechanistic insights into the tumor growth phenotype. A statistically significant decrease of Ki67-positive cells in the DOCK1^flx/flx tumors was noted compared with heterozygote or wild-type tumors (Fig. 2H and Fig. S5D). Conversely, a significant increase in caspase-3 positive cells was detected in the DOCK1-mutant tumors, suggesting that DOCK1 provides a survival signal during Neu oncogenic stress (Fig. 2I and Fig. S5D). One proliferative signal downstream of HER2/Integrin β4 (Itgβ4) crosstalk is c-JUN (17). Interestingly, its activator JNK is a target of DOCK1–Rac signaling (18). We examined the levels of phosphorylated c-JUN (pc-JUN) in MINs and found a decrease in DOCK1-mutant lesions in comparison with wild-type counterparts (Fig. 2J and Fig. S6). Notably, this difference was only noticeable in MINs as pc-JUN levels were similar in fully developed wild-type and mutant tumors (Figs. S6 and S7). These results indicate that DOCK1 contributes to HER2-mediated tumor growth.

DOCK1 Is Essential for HER2-Induced Metastasis.

NIC mice develop lung metastasis and we investigated whether inactivation of DOCK1 can protect from this cancer progression step. The incidence of metastasis was decreased in lungs from NIC⁺DOCK1^flx/flx (39%) in comparison with lungs from NIC⁺DOCK1^wt/wt (52%) and NIC⁺DOCK1^wt/flx (55%) mice (Fig. 3A). For animals afflicted with lung lesions, a fivefold reduction of the number of metastases was noted in NIC⁺DOCK1^flx/flx compared with NIC⁺DOCK1^wt/wt mice (Fig. 3B). Notably, 100% of the lung metastases in the NIC⁺DOCK1^flx/flx mice were confined to the lung vasculature, whereas close to half the metastases were extravascular in control animals, suggesting a role for DOCK1 in lung invasion (Fig. 3C). In NIC⁺DOCK1^wt/wt animals, lung metastases expressed high levels of DOCK1 (Fig. S8) with respect to the expression detected in the primary tumors (Fig. 2D). Notably, we also noted residual expression of DOCK1 in metastases from NIC⁺DOCK1^flx/flx animals, suggesting that lung lesions are enriched in a small population of cells that resisted Cre-mediated gene inactivation (Fig. S8). We also examined whether immune cells are differentially recruited to breast tumors in the mutant mouse. Equal amounts of CD45⁺ hematopoietic cells were present in tumors from all genotypes (Fig. S9), ruling out a major contribution of the immune stroma to the observed metastatic phenotype. We conclude that the loss of DOCK1 is likely to have a cell intrinsic effect in NIC transformed cells.

• Download figure
• Open in new tab
• Download powerpoint

Fig. 3.

DOCK1 regulates HER2-mediated lung metastasis. (A) Incidence of lung metastases is reduced in NIC⁺DOCK1^flx/flx animals. Quantification of the percentage of mice that develop lung metastases for the indicated genotypes. (B) NIC⁺DOCK1^flx/flx mice have a reduction in the amount of lung lesions. (C) H&E staining showing that metastases in NIC⁺DOCK1^wt/wt mice can be found residing inside blood vessels (Left) and invading the lung parenchyma (Center). Metastases in NIC⁺DOCK1^flx/flx were found inside blood vessels (Right). (Scale bar: 250 µm, 10×.) (D) Western blot analysis demonstrating two independent and stable shRNA-mediated knockdown of DOCK1 in NMuMG-NeuNT cells. (E) Representative pictures of the lungs in F. (F) Average lung lesions from animals injected with NMuMG-NeuNT, NMuMG-NeuNT^sh1 DOCK1, and NMuMG-NeuNT^sh2 DOCK1 cells. (G) DOCK1 regulates STAT3 phosphorylation in HER2-driven MINs. IHC analyses showing pSTAT3 staining in mice lesions. (Scale bar: Upper, 500 µm, 5×; Lower, 100 µm, 20×.) One-way ANOVA followed by a Bonferroni or a Student t test was used to calculate the indicated P values; *P ≤ 0.05 and **P ≤ 0.01.

Because a decrease in the primary tumor volume in NIC⁺DOCK1^flx/flx mice could alter metastasis efficiency, an experimental metastasis assay was performed to validate our in vivo findings. NeuNT-transformed Normal Murine Mammary Gland epithelial (NMuMG) cells form lung metastases upon tail vein injection (19) and they were modified to express two independent shRNAs targeting DOCK1 (Fig. 3D). An equal number of NMuMG-NeuNT, NMuMG-NeuNT^sh1 DOCK1, and NMuMG-NeuNT^sh2 DOCK1 cells were injected in nude mice and metastasis to lungs was quantified. We observed a 9.6- and 17.3-fold decrease in metastases in mice injected with cells where DOCK1 was down-regulated by shRNAs in comparison with control cells (Fig. 3 E and F and Fig. S10). Collectively, these results demonstrate that DOCK1 is a central mediator of HER2-driven breast cancer metastasis.

As DOCK1 and Rac act in the integrin-signaling pathway, we verified whether upstream components were activated normally in DOCK1-null tumors. In agreement with the position of DOCK1 in the signaling cascade, we tested whether intermediates were activated correctly (Fig. S11). An important mediator of metastasis downstream of HER2, and a partner of activated Rac, is STAT3 (17). We investigated the phosphorylation status of STAT3 (pSTAT3) in MINs and found that it is highly phosphorylated and distributed in a polarized manner in control MINs, whereas only a faint and evenly distributed pSTAT3 signal is detectable in DOCK1-deficient MINs (Fig. 3G and Fig. S6). This difference in STAT3 activation was only noted in MINs because a similar pSTAT3 signal was observed in all fully developed tumors (Figs. S6 and S7). These data identify DOCK1 as a critical component to promote STAT3 phosphorylation during the early stages of HER2 tumorigenesis.

Identification of a Gene Signature Associated with DOCK1 Expression in HER2 Mammary Tumors.

DOCK1 is generally assumed to exert its function at the membrane via Rac activation and cytoskeleton remodeling (9). To further characterize the defects in tumor growth and metastasis occurring in DOCK1-deficient HER2 tumors, we established the transcriptome of NIC⁺DOCK1^wt/wt (n = 4) and NIC⁺DOCK1^flx/flx (n = 4) mammary tumors by next generation RNA sequencing and uncovered significant variations in the expression of 45 genes between our experimental conditions (more than twofold change; P value <0.05; of ∼16,000 transcripts sequenced) (Fig. 4A and Table S1). Globally, this represents a constrained change of 0.28% in gene expression. These differential expression data were confirmed by quantitative (Q)-PCR for a subset of the genes (Fig. S12 and Table S2). Among up-regulated genes in the DOCK1-knockout tumors were Casein α/β (CSN1S1 and CSN2), Lactalbumin (LALB1), and Estrogen receptor 1 (ESR1) (Fig. 4A), characteristic markers of differentiated mammary epithelial cells, suggesting that DOCK1-mutant tumors are less dedifferentiated.

• Download figure
• Open in new tab
• Download powerpoint

Fig. 4.

Identification of a DOCK1-null gene signature enriched in IFN response genes. (A) Heat map of the 45 genes differentially expressed between NIC⁺DOCK1^wt/wt and NIC⁺DOCK1^flx/flx tumors. Red, elevated; green, decreased; intensity of color represents relative change. Genes in red are bona fide players of IFN signaling according to the literature. (B) Gene ontology analyses. (C) Expression levels of a subset of genes from the NIC⁺DOCK1^flx/flx signature correlate with disease-free survival in human HER2⁺ breast cancer patients. Correlation coefficient >1 means high expression of the gene correlate with a poor prognostic and <1 with a good prognostic. (D) Expression levels of RTP4 and STAT1 independently predict disease-free survival in HER2⁺ cancer patients.

We performed gene ontology (GO) analyses to gain insights into the biological processes regulated by the differentially expressed genes. Unexpectedly, 34/45 differentially expressed genes, and specifically 34/37 of the down-regulated genes in DOCK1-null tumors, are enriched for IFN response genes (Fig. 4 A and B; genes shown in red). These genes can be broken down further into four functional categories (Fig. 4B). A number of the identified genes in our DOCK1-null tumors have previously been linked to a STAT1-governed cancer gene signature that is predictive of risk of resistance to radiotherapy (20) (Fig. 4A). Likewise, we identified STAT1 as a down-regulated gene in our DOCK1-null gene signature and we could confirm lower levels of STAT1 protein in our DOCK1-null tumors both by IHC and Western blot (Figs. S6 and S7). As all of the other down-regulated genes are targets of STAT1, we envision that our DOCK1-null gene signature could reflect variation in STAT1 activity.

STAT1 is known to regulate transcription of a large set of response genes when IFN signaling is engaged. Unexpectedly, oncogenic transformation by Ras was shown to induce the expression of IFN-responsive genes including Interferon Stimulated Gene 15 (ISG15) that is also linked to poor prognosis and increased cell migration (21⇓–23). We investigated whether such a gene regulation profile could be reproduced ex vivo upon NeuNT transformation. For the genes tested (RTP4, ISG15, STAT1, and UBA7), higher expression could be detected in NMuMG-NeuNT in comparison with parental NMuMG-EV explants (Fig. S13). These genes are also regulated by DOCK1 as they were less expressed in NMuMG-NeuNT^sh1 DOCK1 cells (Fig. S13).

We examined the importance of our DOCK1-null gene signature in breast cancer and found that expression levels of five genes (UBA7, RTP4, CSN1S1, DDX60, and IFI44) are individually predictive of disease-free survival in HER2⁺ patients (Fig. 4C). Furthermore, we examined the survival curves and identified six genes (STAT1, RTP4, OASL, PARP12, LGALS9, and LGALS3BP) for which high levels of expression correlated with a worse survival prognostic in HER2⁺ patients (Fig. 4D and Fig. S14A). We tested whether the DOCK1-associated gene signature could predict breast cancer outcome in a cohort of patients irrespective of the molecular subtype and found that our DOCK1^flx/flx-like signature correlates with good patient survival (Fig. S14B). These data demonstrate that DOCK1 controls the expression of a gene signature that predicts the outcome of breast cancer patients.