Purpose:

FGFR1 overexpression has been associated with endocrine resistance in ER+ breast cancer. We found FGFR1 localized in the nucleus of breast cancer cells in primary tumors resistant to estrogen suppression. We investigated a role of nuclear FGFR1 on gene transcription and antiestrogen resistance.

Experimental Design:

Tumors from patients treated with letrozole were subjected to Ki67 and FGFR1 IHC. MCF7 cells were transduced with FGFR1(SP-)(NLS) to promote nuclear FGFR1 overexpression. FGFR1 genomic activity in ER+/FGFR1-amplified breast cancer cells ± FOXA1 siRNA or ± the FGFR tyrosine kinase inhibitor (TKI) erdafitinib was examined by chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq). The nuclear and chromatin-bound FGFR1 interactome was investigated by mass spectrometry (MS).

Results:

High nuclear FGFR1 expression in ER+ primary tumors positively correlated with post-letrozole Ki67 values. Nuclear FGFR1 overexpression influenced gene transcription and promoted resistance to estrogen suppression and to fulvestrant in vivo. A gene expression signature induced by nuclear FGFR1 correlated with shorter survival in the METABRIC cohort of patients treated with antiestrogens. ChIP-Seq revealed FGFR1 occupancy at transcription start sites, overlapping with active transcription histone marks. MS analysis of the nuclear FGFR1 interactome identified phosphorylated RNA-Polymerase II and FOXA1, with FOXA1 RNAi impairing FGFR1 recruitment to chromatin. Treatment with erdafitinib did not impair nuclear FGFR1 translocation and genomic activity.

Conclusions:

These data suggest nuclear FGFR1 contributes to endocrine resistance by modulating gene transcription in ER+ breast cancer. Nuclear FGFR1 activity was unaffected by FGFR TKIs, thus supporting the development of treatment strategies to inhibit nuclear FGFR1 in ER+/FGFR1 overexpressing breast cancer.

Translational Relevance

We show herein evidence to support an unconventional role for FGFR1 in gene transcription regulation in ER+ breast cancer. Nuclear FGFR1 induces a gene expression profile promoting resistance to antiestrogens. Furthermore, nuclear FGFR1 translocation and activity are not impaired by FGFR1 TKIs that inhibit canonical FGFR1 signaling pathway. These findings support development of pharmacological inhibitors to repress this novel receptor function in ER+/FGFR1–overexpressing breast cancer with a therapeutic intent.

FGFR1 belongs to the family of fibroblast growth factor receptors, comprising four highly conserved transmembrane receptor tyrosine kinases (RTK) FGFR1–4, and another membrane-associated receptor lacking the intracellular domain (FGFR5 or FGFRL1; ref. 1). FGFRs are activated upon binding of ligands (FGFs) to their extracellular domain (2) followed by receptor dimerization and phosphorylation of C-terminal tyrosines. These phosphorylated tyrosines dock several adaptor proteins that induce activation of downstream signaling pathways, including RAS/RAF/MEK/ERK, PI3K/AKT, PLCγ, and STATs (3–6).

The genomic locus of FGFR1, 8p11–12, is amplified in various cancer types, including breast, lung, ovarian, and bladder tumors (7–10). FGFR1 gene amplification has been identified in about 15% of patients with ER+ breast cancers (11). We have previously shown that FGFR1 gene amplification is associated with resistance to estrogen suppression in a cohort of patients with ER+ breast cancer treated with the aromatase inhibitor letrozole (12).

The oncogenic role of FGFR1 is generally the result of genomic aberrations, such as gene amplification, activating mutations, gene fusions, or by dysregulated autocrine/paracrine signaling, involving FGF ligands (13). In addition to the canonical intracellular signaling function of membrane-bound FGFR1, there is experimental evidence that FGFR1 can localize in the nucleus of cancer cells (14, 15). The role of nuclear FGFR1 has been investigated in the context of neuronal development. In embryonic stem cells (ESC), nuclear FGFR1 has been shown to interact with chromatin, alone or in association with the nuclear receptors RXR (Retinoid X Receptor) or NR4A1 (also known as Nur77 or nerve growth factor IB), thereby promoting gene transcription associated with developmental pathways (16). Furthermore, in human medulloblastoma cells, FGFR1 interacts with CREB-Binding Protein (CBP) and ribosomal S6 kinase 1 (RSK1) in the nucleus. In this ternary complex, FGFR1 induces CBP-mediated transcription by releasing CBP from RSK1 inhibition and directly activating CBP (17).

There is a growing body of evidence showing that RTKs, in addition to their signal transduction function as membrane-bound receptors, can also be found in the nucleus (18). For example, it was recently reported that the insulin receptor (IR) associates with several gene promoters, thus regulating gene expression profiles involved in classic insulin-related functions (19). Genomic binding studies recently identified a direct involvement of HER2 (ERBB2) in gene transcription regulation in HER2-amplified breast cancers (20). Other studies showed that EGFR binds to promoters of genes involved in cancer cell proliferation and regulates their expression (21).

FGFR1 overexpression has been associated with resistance to antiestrogen therapy in ER+ breast cancers (22, 23). We previously showed that estradiol deprivation in patients resulted in a significant increase in both cytosolic and nuclear FGFR1 levels in ER+/FGFR1-amplified breast cancer cells in primary tumors (24). In this study, we now report abundant levels of FGFR1 in nuclear and chromatin-bound fractions from ER+/FGFR1-amplified PDXs as well as a correlation between nuclear localization of FGFR1 and resistance to estrogen suppression in ER+ primary tumors in patients. Hence, we have examined whether nuclear FGFR1 has a potential causal role in endocrine resistance. To investigate a transcriptional function of nuclear FGFR1, we analyzed genome-wide chromatin occupancy of FGFR1 by chromatin immunoprecipitation followed by next-generation sequencing (ChIP-Seq) and FGFR1-associated gene expression profiles by RNA sequencing (RNA-seq) in ER+/FGFR1-amplified breast cancer cells. Taken together, these analyses identified a functional role for nuclear FGFR1 in transcriptional regulation, driving gene expression associated with antiestrogen resistance and poor patient outcome.

Cell lines and inhibitors

MCF-7 (ATCC HTB-22), CAMA1 (ATCC HTB-21), and MDA-MD-134-VI (ATCC HTB-23) human breast cancer cells were obtained from the American Type Culture Collection (ATCC) and maintained in ATCC-recommended media supplemented with 10% FBS (Gibco) and 1× antibiotic/antimycotic (Gibco) at 37°C in a humidified atmosphere of 5% CO2 in air. Cell lines were authenticated by ATCC prior to purchase by the short tandem repeat (STR) method. All experiments were performed <2 months after thawing early passage cells. Mycoplasma testing was conducted for each cell line before use. Fulvestrant, tamoxifen and erdafitinib were purchased from Selleck Chemicals.

Xenograft studies

Mouse experiments were approved by the UT Southwestern Institutional Animal Care and Use Committee (IACUC). Female ovariectomized athymic mice (Envigo) were implanted with a 14-day-release 17β-estradiol pellet (0.17 mg; Innovative Research of America). The following day, 1 × 107 MCF-7EV or MCF-7FGFR1(SP-)(NLS) cells suspended in PBS and Matrigel at 1:1 ratio were injected subcutaneously (s.c.) into the right flank of each mouse. Approximately 4 weeks later, mice bearing tumors measuring ≥150 mm3 were randomized to treatment with vehicle (control) or fulvestrant (5 mg per week; s.c.). Tumor diameters were measured with caliper every four days and tumor volume in mm3 was calculated with the formula: volume  =  width2 × length/2. After 4 weeks, tumors were harvested and snap-frozen in liquid nitrogen or fixed in 10% neutral buffered formalin followed by embedding in paraffin. Four-micron paraffinized sections were used for IHC using ERα (Santa Cruz Biotechnology #8002) antibody. Sections were scored by an expert pathologist (P.G.E) blinded to treatment arm.

Immunohistochemistry

For IHC analysis, samples were prepared as previously described (24). Briefly, formalin-fixed paraffin-embedded (FFPE), 4-μmol/L tissue sections were deparaffinized. Antigen retrieval was performed with citrate buffer pH 6. Endogen peroxidase was blocked and protein block was applied. Sections were then incubated with FGFR1 ab10646 (Abcam) antibody at 1:2000 overnight at 4°C. The visualization system Envision (Agilent Technologies, Santa Clara, CA). DAB was used as chromogen (Agilent Technologies) and hematoxylin was applied as counterstain. Slides were digitally acquired using an AxioScan Z1 slide scanner (Carl Zeiss) at 20x. Automated semiquantitative scoring was performed on whole slide images using QuPath software (25). Color deconvolution stains were set for hematoxylin and DAB. Cell segmentation was determined on hematoxylin. An object classifier was trained to define tumor and stroma compartments. Percentage of FGFR1 positive cells over total tumor cell and H-score were calculated with the positive cell detection algorithm according to the nuclear DAB optical densities mean. Each selected region was visually assessed for correct performance of the quantification algorithm. Ki67 expression in biopsies from patients treated with letrozole prior to surgery was measured by automated quantitative immunofluorescence, as previously described (12, 26). Studies involving patients' samples were conducted in accordance with the Declaration of Helsinki and approved by the Vanderbilt-Ingram Cancer Center Institutional Review Board. All patients provided a written informed consent.

ChIP and reChIP assays

Cells cultured in 15-cm dishes were washed with PBS once and incubated with serum free DMEM media containing 1% formaldehyde, for 10 min at room temperature. Crosslinking was quenched with glycine. Plates were next washed with ice-cold PBS three times and cells scraped off the plate with ice-cold PBS with protease inhibitor tablet (Mini, EDTA-free Protease Inhibitor Cocktail – ROCHE – Cat no. 04693159001). Cells were centrifuged at 720 × g for 10 min at 4°C and resuspended in 1.5 ml sonication buffer (1% SDS, 10 mmol/L EDTA, 50 mmol/L Tris-HCl, pH 8.0–8.1) with protease inhibitor. After rotation at 4°C for 10 min, chromatin was sonicated using a Bioruptor Plus sonication device (Diagenode) at 4°C. Sonicated chromatin was adjusted to a final concentration of 200 mmol/L NaCl and incubated at 65°C overnight. The following day, the chromatin was incubated with RNAse A for 30 min at 37°C, followed by Proteinase K for 1 h at 42°C. DNA was purified using QIAquick PCR Purification Kit (Cat 28106). DNA shearing was checked, with an average fragment size of 200–750 bp. To conduct chromatin immunoprecipitation (ChIP), the sonicated chromatin was eluted with ChIP dilution buffer (0.01% SDS buffer, 1.1% Triton X-100, 1.2 mmol/L EDTA, 167 mmol/L NaCl, 16.7 mmol/L Tris-HCl, pH 8.0–8.1) with protease inhibitors and precleared with Gammabind G Sepharose beads (GE Healthcare, Cat no. 17–0885–01) previously washed three times with ChIP-dilution buffer and preblocked in 0.5% BSA at 4°C for 1 h. The precleared chromatin was incubated at 4°C overnight with primary antibody. The following day, Gammabind G Sepharose beads were added to the antibody pulldowns for 2 h at 4°C. Next, beads were washed once in Buffer I (0.1% SDS, 1% Triton X-100, 2 mmol/L EDTA, 20 mmol/L Tris-HCl, 150 mmol/L NaCl and protease inhibitors, pH 8.0–8.1), once in Buffer II (0.1% SDS, 1% Triton X-100, 2 mmol/L EDTA, 20 mmol/L Tris-HCl, 500 mmol/L NaCl, and protease inhibitors), once in Buffer III (0.25 M LiCl, 1% NP40, 1% Na-deoxycholate, 1 mmol/L EDTA, 10 mmol/L Tris-HCl, pH 8.0–8.1) and twice in TE buffer, pH 8.0–8.1. Complexes were eluted off the beads in Elution Buffer (0.1 M NaHCO3, 1%SDS) at 65°C for 10 min, twice, in a thermomixer. Eluates from the antibody pulldowns and sonicated chromatin (as input controls) were adjusted to 200 mmol/L NaCl and incubated at 65°C overnight. The following day, sequential incubations with RNAse A for 30 min at 37°C, and Proteinase K for 1 h at 42°C were performed. DNA was purified using QIAquick PCR Purification Kit. ChIPed and input DNA was analyzed by qRT-PCR or used to construct sequencing libraries. For ChIP-qPCR, the enrichment was shown as percent of input. ChIP-qPCR results were reproduced in two or more independent experiments. We used the following primary antibodies for ChIP experiments: rabbit anti-FGFR1 (Abcam, ab10646), normal rabbit IgG (Cell Signaling Technology, 2729), rabbit anti-histone H3 acetyl K27 (Abcam, ab4729), rabbit anti-histone H3 trimethyl K4 (Abcam, ab8580), mouse anti-Rpb1 CTD 4H8 (Cell Signaling Technology, 2629), mouse anti-HNF-3α (Santa Cruz Biotechnology, sc-101058). For ChIP-qPCR, primer pair sequences were designed based on the nucleotide sequences identified by ChIP-Seq peaks and derived from UCSC Genome Browser. The identified sequence was input in BLAST (Basic Local Alignment Search Tool) to design valuable pair sequences. Primer pair sequences are indicated in Supplementary Table S1.

For reChIP, ChIP with rabbit anti-FGFR1 or normal rabbit IgG was performed as described above but after the last TE wash, beads were incubated with reChIP elution buffer (10 mmol/L Tris pH 7.6, 1 mmol/L EDTA, 2% SDS, 20 mmol/L DTT) and heated at 37°C for 30 min (19). Samples were centrifuged at 1000 × g × 1 min and the supernatant was diluted 20-fold in ChIP dilution buffer containing protease inhibitors. Next, samples were split and incubated with rabbit IgG or mouse anti-HNF-3α or mouse anti-Rpb1 CTD 4H8 at 4°C overnight. The reChIP samples were collected, washed, and eluted as described above for ChIP. The enrichment value was measured by qRT-PCR and shown as percent (%) of input of the first round of ChIP.

ChIP-sequencing analysis

Cells were cultured in 15-cm dishes and kept for 48 h in IMEM/10% Charcoal Stripped Serum (CSS) prior to apply the ChIP protocol described above. To test the effect of FGF2 and erdafitinib on FGFR1 genomic distribution, CAMA1 cells were cultured for 24 h in serum free DMEM media, prior to treatment with FGF2 (R&D Systems, 233-FB) (5 ng/mL, 3 h) ± erdafitinib (250 nmol/L, 3 h). All ChIP-Seq in CAMA1 cells experiments were conducted on two biological replicates. ChIP-Seq experiments in MCF7FGFR1(SP-)(NLS) cells were conducted on three biological replicates. Libraries from 5–10 ng ChIP DNA were prepared using KAPA Hyper Library Preparation Kit. ChIPed DNA was quantified on the Qubit 4 Fluorometer (Invitrogen). Samples were end repaired, 3′ ends adenylated and barcoded with multiplex adapters, followed by size selection with Ampure XP beads and PCR amplification. Samples were next quantified by Qubit and validated on the Agilent Tapestation 4200, normalized and pooled, and then run on the Illumina NextSeq 500 using 75 cycle SBS v2.5 reagents. Reads were aligned to human reference genome (hg38) using Bowtie2(v2.2.3) (27) with default parameters. Low quality reads and duplicate reads were removed from aligned files using “samtools view -bh -F 0 × 04 -q 10” (28) and “Picard MarkDuplicates.jar” (https://github.com/broadinstitute/picard) commands. Model-based Analysis of ChIP-Seq (MACS) (29) software tool (v.2.1.2) was used to identify TF bound regions from the ChIP-seq data. MACS uses 200-bp as the fragment length and each read was shifted by 100-bp to identify candidate peaks with significant tag enrichment. Default parameters of MACS were applied using an FDR cutoff of 5% to select the peaks for further analysis. For ChIP-Seq experiments in CAMA1 cells, only peaks reproduced in both replicates were considered relevant and used for further analysis. For experiments in MCF7FGFR1(SP-)(NLS) cells only peaks reproduced in at least two replicates were considered for further analysis. For motif analysis, 200-bp around the peak region was used. De novo motif discovery analysis was performed using findMotifsGenome module in homer (v4.9). We used bdgcmp executable available in MACS software tools to calculate the background signal from the input sample and subtract it from the ChIP signal. We used the “ppois” method to calculate the normalized signal. To plot sequencing read densities relative to specific positions in the genome, we used annotatePeaks function available in Homer tools (30). This function takes the peak coordinates, tags directories as input and extends each tag by their estimated ChIP-fragment length, and calculates ChIP fragment coverage represented in per base pairs per peak. Heatmaps were plotted using signal density in the 5000-bp region around each peak. Density plots around TSS were plotted using ngsplot tool (31). GREAT (http://great.stanford.edu/public/html/) was used to identify the potential biological functions of the FGFR1 DNA binding sites with default parameters. All the processed ChIP-Seq data generated for this manuscript are included in Supplementary Table S2.

RNA-sequencing analysis

CAMA1 cells were cultured in 10-cm dishes in either IMEM/10% CSS for 48 h or 24 h in serum free DMEM media, prior to treatment with DMSO or FGF2 (5 ng/mL, 6 h). Cells were harvested and RNA was purified using a RNA purification kit (Maxwell, Promega). Each experiment was conducted in three biological replicates. Samples were run on the Agilent Tapestation 4200 to ensure use of only high quality RNA (RIN Score ≥8). The Qubit fluorometer was used to determine sample concentration prior to preparation of libraries. One μg of total DNAse treated RNA was then prepared with the TruSeq Stranded mRNA Library Prep Kit from Illumina. Poly-A RNA was purified and fragmented before strand specific cDNA synthesis. cDNAs were then a-tailed and indexed adapters were ligated. After adapter ligation, samples were PCR amplified and purified with AmpureXP beads, then validated again on the Agilent Tapestation 4200. Before being normalized and pooled, samples were quantified by Qubit then run on the Illumina NextSeq 500 using V2.5 reagents. RNA-seq reads were mapped to human genome (hg38) and junctions were identified using tophat(v2.1.2). RNA-seq expression counts for each gene were quantified using featureCounts module from Subread package (v.1.6.3). Differential expression analysis was performed using DESeq2 (v.1.24.0) R package. Binding and Expression Target Analysis (BETA) v1.0.7, with the BETA-plus protocol, was used to evaluate significant association between FGFR1-bound genomic loci in MCF7FGFR1(SP-)(NLS) and differentially expressed genes in MCF7EV vs. MCF7FGFR1(SP-)(NLS) cells. BETA was also used to evaluate significant association between FGFR1-bound genomic loci and FGF2-induced DEGs in CAMA1 cells. GSEA analysis were performed using a JAVA GSEA 3.0 program. All the processed RNA-Seq data generated for this manuscript are included in Supplementary Table S3.

Gene expression analysis

RNA was extracted and purified from cells or xenografts using Maxwell RSC simplyRNA Cells Kit or simplyRNA Tissue Kit (Promega Corporation). cDNA was generated using iScript cDNA Synthesis Kit (Bio-Rad). qPCR was performed with a cDNA equivalent of 50 ng RNA, 1 mmol/L each of the forward and reverse primers, and SYBR Green PCR Master Mix (Applied Biosystems), using a QuantStudio 3 Real-Time PCR System Machine (Applied Biosystems). All primers were purchased from QIAGEN: GAPDH (PPH00150F) PGR (PPH01007F), CCND1 (PPH00128F), VEGFA (PPH00251C), TFF1 (PPH00998C), PDZK1 (PPH08038E), CDK4 (PPH00118F), CDK12 (PPH05712A), DUSP1 (PPH00406A), FTO (PPH16000B), JUNB (PPH00179A), GREB1 (PPH20761F), MYC (PPH00100B), FOS (PPH00094A), BRD2 (PPH09948A). Ct (threshold cycle) values were determined in triplicate samples by subtracting the target gene Ct from the GAPDH Ct; 2–ΔΔCt was used to determine the expression of selected mRNAs relative to GAPDH.

Immunoprecipitation and Western blot analysis

Cells were seeded in 15-cm dishes and after three washes with ice-cold PBS, harvested in lysis buffer (50 mmol/L Tris pH 7.5, 150 mmol/L NaCl, 0.5 mmol/L EDTA, and 1% Triton X-100) supplemented with protease inhibitors (Mini, EDTA-free Protease Inhibitor Cocktail – ROCHE) and phosphatase inhibitor (PhosSTOP; ROCHE). Lysates were sonicated using a Bioruptor Plus sonication device (Diagenode) at 4°C, at maximum power for 10 min with 30 sec on and 30 sec off twice. Lysates were incubated with DNase I (Thermo Fisher Scientific) and Benzonase (Sigma-Aldrich) at 37°C for 10 min, followed by centrifugation at 20,000 xg at 4°C for 10 min. Protein concentration in the supernatants was determined using Pierce BCA Protein Assay Reagents (Thermo Fisher). Lysates were incubated overnight with primary antibodies at 4°C and antibody pulldowns were captured by incubation with Protein G dynabeads for 2 h at 4°C. Beads were washed and the immune complexes eluted from the beads following manufacturer's instructions (Dynabeads Protein G Immunoprecipitation Kit – Thermo Fisher). Eluted proteins were separated by SDS-PAGE and transferred to nitrocellulose membranes. Membranes were blocked with 5% milk at room temperature for 1 h and then incubated overnight with primary antibodies at 4°C, followed by incubation with HRP-conjugated rabbit or mouse secondary antibodies for 1–2 h at room temperature. Protein bands were detected with an enhanced chemiluminescence substrate (Pierce ECL Western Blotting substrate, Thermo Fisher). The following primary antibodies were used for immunoprecipitation: rabbit anti-FGFR1 (Abcam, ab10646), mouse anti-FLAG M2 (Sigma-Aldrich), mouse anti-Rpb1 CTD 4H8 (Cell Signaling Technology, 2629), rabbit Phospho-Rpb1 CTD (Ser2) (E1Z3G) (Cell signaling technology, 13499), mouse anti-CDK7 (CST, 2916), mouse anti-CDK9 (Santa Cruz Biotechnology, sc-13130), mouse anti-HNF-3α (Santa Cruz Biotechnology, sc-101058), normal rabbit IgG (CST, 2729), mouse IgG2b isotype control (CST, 53484). The following primary antibodies were used for immunoblotting: rabbit anti-FGFR1 (Abcam ad76464), mouse anti-FLAG M2 (Sigma-Aldrich), mouse anti-Rpb1 CTD 4H8 (Cell Signaling Technology, 2629), rabbit Phospho-Rpb1 CTD (Ser2) (E1Z3G) (Cell signaling technology, 13499), mouse anti-CDK7 (CST, 2916), mouse anti-CDK9 (Santa Cruz Biotechnology, sc-13130), mouse anti-HNF-3α (Santa Cruz Biotechnology, sc-101058), rabbit anti-α/β Tubulin (CST, 2148), rabbit anti-AIF (CST, 5318), mouse anti-Lamin A/C (CST, 4777), rabbit anti-Histone H3 (CST, 4499), rabbit anti-phospho-FGF Receptor (Tyr653/654) (CST, 3471), rabbit anti-phospho-FRS2α (Tyr436) (CST, 3861), rabbit anti-phospho-p44/42 MAPK (Erk1/2) (Thr202/204) (CST, 4370), rabbit anti-p44/42 MAPK (Erk1/2)(137F5) (CST, 4695), mouse anti-β-actin (CST, 3700). For immunoprecipitation followed MS, CAMA1FGFR13XFLAG cells underwent subcellular fractionation according to the manufacturer's protocol (Subcellular Protein Fractionation Kit for Cultured Cells, CST #78840). Soluble nuclear and chromatin-bound fractions were combined, sonicated and nuclease-treated as described above and used for immunoprecipitation with mouse anti-FLAG M2 (Sigma-Aldrich) and mouse IgG2b isotype control (CST, 53484) antibodies. For subcellular fractionation assays, CAMA1 and MDA-MB-134 cells were cultured for 24 h in serum free DMEM media, prior to treatment with FGF2, 5 ng/mL ± erdafitinib 250 nmol/L for 3 h.

Mass spectrometry

Immunoprecipitated proteins were separated by SDS-PAGE. Excised Gel band samples were digested overnight with trypsin (Promega) following reduction and alkylation with DTT and iodoacetamide (Sigma–Aldrich). Digests underwent solid-phase extraction cleanup with an Oasis MCX plate (Waters) followed by run on an Orbitrap Fusion Lumos mass spectrometer coupled to an Ultimate 3000 RSLC-Nano liquid chromatography system. Samples were injected onto a 75 um i.d., 75-cm long EasySpray column (Thermo). As previously described (32), the mass spectrometer operated in positive ion mode. We used source voltage of 2.2 kV and ion transfer tube temperature of 275°C. Orbitrap acquired MS scans at 120,000 resolution; MS/MS spectra, up to 10, were obtained in the ion trap using higher-energy collisional dissociation (HCD) for ions with charges 2–7. Dynamic exclusion was set for 20 sec after an ion was selected for fragmentation. Proteome Discoverer v2.2 (Thermo) was used to analyze raw MS data. Sequest HT was used to identify peptides from UniProt human protein database. Carbamidomethylation of Cys and oxidation of Met were set as a fixed and variable modifications, respectively. We applied a false-discovery rate (FDR) cutoff of 1% for all peptides. We used three replicates of FLAG-immunoprecipitates and two replicates for IgG-immunoprecipitates. Only the proteins present in all three FLAG pulldown replicates and absent in IgG pulldown replicates or with a fold enrichment in FLAG over IgG >10 were considered as positive hits. Proteins with less than two peptide hits were not included.

Immunofluorescence

Cells were seeded in MatTek 35 mm dishes (MatTek Life Sciences) and after 48 h, fixed with 4% PFA for 30 min. Cells were next incubated in 0.25% Triton X-100 in PBS for 5 min, then in blocking buffer (10% BSA in PBS) for 60 min, followed by an overnight incubation with an FGFR1 antibody (Abcam cat #76464) at 4 °C. The following day, cells were incubated with a goat anti-rabbit IgG (H+L) cross-adsorbed secondary antibody, Alexa Fluor 488 (Thermo Fisher Scientific, cat #A-11008) for 60 min, followed by DAPI staining. Vectashield antifade mounting medium (Vector Laboratories, cat #H-1000–10) was used as mounting solution. Images were collected with a fully automated inverted microscope (DMI8, Leica microsystems). Quantification of nuclear fluorescence was performed with the software ImageJ.

FGFR1 fluorescent in situ hybridization (FISH) analysis

Four-μmol/L tissue sections were mounted on charged slides and hybridized overnight with the SPEC FGFR1:CEN8 Dual Color Probe (ZytoVision, catalog# Z-2072–200). Briefly, deparaffinization, protease treatment and washes were performed as per standard protocols. After this pretreatment, the slides were denatured in the presence of 10 μL of the probe for 6 min at 72°C, and hybridized at 37°C overnight in StatSpin (Thermobrite, Abbott Molecular, Inc.). Post-hybridization saline-sodium citrate washes were performed at 72°C and the slides were then stained with DAPI before analysis. Tumor tissue was scanned for amplification hot spots under 20× magnification (Olympus BX60 Fluorescent microscope). If the FGFR1 signals were homogeneously distributed, then random areas were used for counting the signals. Images for cell counting were captured with a 100x oil immersion objective using Cytovision software. Twenty to sixty tumor cell nuclei were individually evaluated by counting green FGFR1 and orange centromere 8 signals. The FGFR1:CEN8 ratio and the average FGFR1 copy number for each specimen were calculated next. Based on HER2 guidelines (33), samples were considered to be FGFR1 amplified under one of the following conditions: a) FGFR1:CEN8 ratio ≥2.0; b) average number of FGFR1signals ≥6.

Clonogenic assay

Cells (5 × 104/well) were seeded in triplicate in 10% DMEM-FBS in 6-well plates and treated with 0.5–10 nmol/L fulvestrant or with increasing doses of erdafitinib or tamoxifen. Seven days later, when control wells reached 50–70% confluence, monolayers were fixed and stained with 20% methanol/80% water/0.5% crystal violet for 20 min, washed with water, and dried. Photographic images of the plates were obtained. The crystal violet was solubilized with 20% acid acetic and the image intensity of the monolayers was quantified by spectrophotometric detection at 490 nm using a plate reader (GloMax-Multi Detection System, Promega).

PDX-derived cell sorting

HCI-011 PDXs established in ovariectomized SCID/beige mice were harvested and dissociated using digestion buffer [125 μg/mL DnaseI (#LS002138, Worthington), 10 μg/mL Dispase (#LS02109, Worthington), 500 μg/mL Collagenase 3 (#LS004182, Worthington), and 5× triple antibiotics (#15240–062, Invitrogen) in CnT-PR media (Cellntec)]. Dissociated single cells were washed with ice-cold PBS, stained with APC anti-human CD298 Antibody (BioLegend, #341706) for 30 minutes at 4°C, and fixed with 1% formaldehyde for 15 minutes at 4°C. FACS analysis was performed on the FACSAria Fusion flow cytometer (BD Biosciences). APC+ cells were collected, washed with ice-cold PBS twice and then resuspended in sonication buffer (1% SDS, 10 mmol/L EDTA, 50 mmol/L Tris-HCl, pH 8.0–8.1) with protease inhibitors.

Plasmid and viral transfection

3X FLAG-FGFR1–expressing lentiviral constructs were generated in the pLX302 Gateway vector (Open Biosystems). To generate stably transduced CAMA1 cells, 4 μg of the 3XFLAG-FGFR1-pLX302 construct were cotransfected with 3 μg psPAX2 and 1 μg pMD2G envelope plasmid into 293FT cells using FuGENE 6 (Promega). 293FT growth media was changed 24 hours posttransfection; virus-containing supernatants were harvested 48 and 72 hours posttransfection, passed through a 0.45-mm filter, diluted 1:4, and applied to target cells with 8 μg/mL polybrene (Sigma Aldrich). Virus-producing cells were selected in 1 μg/mL puromycin. pCDNA3.1(-)FGFR1(SP-)(NLS) has been described previously (34). FGFR1(SP-)(NLS) genomic sequence was cloned from pCDNA3.1(-)FGFR1(SP-)(NLS) into pLVX-M-puro plasmid by restriction enzyme cloning, to generate the pLVX-M-puro-FGFR1(SP-)(NLS) plasmid. pLVX-M-puro empty vector and pLVX-M-puro-FGFR1(SP-)(NLS) plasmid were cotransfected with 3 μg psPAX2 and 1 μg pMD2G envelope plasmid into 293FT cells as described above. Virus-containing supernatants from 293FT cells were collected as described above and applied to MCF-7 cells with 8 μg/mL polybrene (Sigma Aldrich). Virus-producing cells were selected in 1 μg/mL puromycin. The mutant FGFR1 D432N was generated from wild-type human FGFR1 in pLX302 using site-directed mutagenesis (Q5 Site-Directed mutagenesis kit, New England Biolabs) using the primer 5′-GGTGTCTGCTAACTCCAGTGC-3′, and its reverse complement. Plasmid was cotransfected with psPAX2 and pMD2G envelope plasmids into 293FT to generate lentiviral particles as described above. FGFR1 was knocked out using pX458 (Addgene plasmid # 48138). Briefly, two sgRNAs that target FGFR1 were separately cloned into pX458 as described previously (35). The plasmids were then cotransfected into MCF-7. Forty-eight hours posttransfection, top 10% of the GFP-positive single MCF-7 cells were sorted into 96-well plates to grow FGFR1 knockout clones. The knockout of FGFR1 was validated by Western blot analysis. The primers used are reported in Supplementary Table S1.

siRNA transfection

Silencer Select predesigned siRNAs targeting FGFR1, FOXA1, and SP1 were purchased from Ambion. Cells were reverse-transfected with siRNAs of interest using Lipofectamine RNAiMAX transfection reagent (Thermo Fisher Scientific) as per manufacturer's instructions. Twenty-four hours posttransfection, the culture media was replaced with Iscove's modified Dulbecco medium/10% CSS. Seventy-two hours posttransfection, the ChIP or RNA extraction protocols described above were applied. Simultaneously, transfected cells were harvested in lysis buffer (50 mmol/L Tris pH 7.5, 150 mmol/L NaCl, 0.5 mmol/L EDTA, and 1% Triton X-100) for protein extraction.

Statistical analysis

Log-rank Mantel–Cox test, two-tailed Mann–Whitney test, unpaired t test, multiple t test analysis, two-way ANOVA with Sidak multiple comparisons test or Wilcoxon matched-pairs test were performed as indicated in figure legends. All data quantification and the statistical analysis described above were performed using Prism (v8).

Data availability

The datasets generated and analyzed during this study are available in the GEO repository, code GSE148313 and can be found at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE148313. All the analyzed chromatin immunoprecipitation (ChIP-seq) and RNA sequencing (RNA-seq) data can be found in Supplementary Material as Supplementary Tables S2 and S3.

Nuclear FGFR1 promotes antiestrogen resistance

We performed IHC and FISH of FGFR1 in 155 biopsies from women with operable ER+ breast cancer who received letrozole for 10–21 days before surgery (in trial NCT00651976; Supplementary Table S4; Supplementary Fig. S1A; refs. 12, 24). In this published trial, response to estrogen suppression with letrozole was defined on the basis of posttreatment Ki67 levels [sensitive: natural log (ln) of post-letrozole Ki67 ≤1.0% or ≤2.4% tumor cells; intermediate: ln = 1.1–1.9 or 2.5–7.3% tumor cells; resistant: (ln ≥2.0 or ≥7.4% tumor cells)]. Nineteen (19/155, 12%) tumors exhibited FGFR1 amplification, defined as an FGFR1:CEN8 ratio ≥2.0 (Supplementary Table S4; Supplementary Fig. S1B). IHC revealed that, in addition to faint FGFR1 cytosolic staining, FGFR1 was clearly localized in breast cancer cell nuclei (Fig. 1A). We observed similar FGFR1 localization in three ER+/FGFR1-amplified PDXs (Supplementary Fig. S1C and S1D). Subcellular fractionation of the PDX lysates revealed that, in addition to the expected localization in the plasma membrane, a substantial fraction of FGFR1 was detected in nuclear fractions (Supplementary Fig. S1E). ER+/FGFR1-amplified CAMA1 and MDA-MB-134 human breast cancer cell fractions also showed nuclear and chromatin-bound FGFR1 (Supplementary Fig. S1F).

Figure 1.

Nuclear FGFR1 promotes antiestrogen resistance. A, Sections of FFPE of ER+/FGFR1-amplified breast primary tumors were subjected to IHC with an FGFR1 antibody as described in the Methods. B, Plot showing direct correlation between nuclear FGFR1 H-Score, measured by IHC, and posttreatment Ki67 in 155 ER+/HER2 breast cancer biopsies from women treated with preoperative letrozole for 10–21 days in the NCT00651976 trial (12, 24). In blue, samples exhibiting FGFR1 amplification (n = 19/155), defined as an FGFR1:CEN8 ratio ≥2.0. C, MCF-7EV and MCF-7FGFR1(SP-)(NLS) cells were seeded in 6-well plates in estrogen-deprived medium. Medium was replenished every 72 hours. Monolayers were stained with crystal violet on days 1, 4, 7, and 10. Quantification of the integrated intensity values as fold change normalized to day 1 (Sidak multiple comparisons test, two-way ANOVA, **** P < 0.0001). D, MCF-7EV and MCF-7FGFR1(SP-)(NLS) cells were seeded in 6-well plates in full medium and treated with vehicle (DMSO) or fulvestrant. After seven days, monolayers were stained with crystal violet. E, Quantification of the integrated intensity values as fold change normalized to vehicle-treated controls (multiple Student t test). F, MCF-7EV and MCF-7FGFR1(SP-)(NSL) xenografts were established in ovariectomized athymic mice implanted with a subcutaneous 14-day release, 0.17-mg 17β-estradiol pellet. Once tumors established, mice were randomized to treatment with vehicle or fulvestrant (5 mg/week). Each data point represents the mean tumor volume in mm3 ± SEM (n = 8 per arm, ****P < 0.0001; two-way ANOVA). G, GSEA of RNA-Seq data for MCF7FGFR1(SP-)(NLS) versus MCF7EV, showing enrichment of gene sets associated with antiestrogen resistance in MCF7FGFR1(SP-)(NLS) cells.

Figure 1.

Nuclear FGFR1 promotes antiestrogen resistance. A, Sections of FFPE of ER+/FGFR1-amplified breast primary tumors were subjected to IHC with an FGFR1 antibody as described in the Methods. B, Plot showing direct correlation between nuclear FGFR1 H-Score, measured by IHC, and posttreatment Ki67 in 155 ER+/HER2 breast cancer biopsies from women treated with preoperative letrozole for 10–21 days in the NCT00651976 trial (12, 24). In blue, samples exhibiting FGFR1 amplification (n = 19/155), defined as an FGFR1:CEN8 ratio ≥2.0. C, MCF-7EV and MCF-7FGFR1(SP-)(NLS) cells were seeded in 6-well plates in estrogen-deprived medium. Medium was replenished every 72 hours. Monolayers were stained with crystal violet on days 1, 4, 7, and 10. Quantification of the integrated intensity values as fold change normalized to day 1 (Sidak multiple comparisons test, two-way ANOVA, **** P < 0.0001). D, MCF-7EV and MCF-7FGFR1(SP-)(NLS) cells were seeded in 6-well plates in full medium and treated with vehicle (DMSO) or fulvestrant. After seven days, monolayers were stained with crystal violet. E, Quantification of the integrated intensity values as fold change normalized to vehicle-treated controls (multiple Student t test). F, MCF-7EV and MCF-7FGFR1(SP-)(NSL) xenografts were established in ovariectomized athymic mice implanted with a subcutaneous 14-day release, 0.17-mg 17β-estradiol pellet. Once tumors established, mice were randomized to treatment with vehicle or fulvestrant (5 mg/week). Each data point represents the mean tumor volume in mm3 ± SEM (n = 8 per arm, ****P < 0.0001; two-way ANOVA). G, GSEA of RNA-Seq data for MCF7FGFR1(SP-)(NLS) versus MCF7EV, showing enrichment of gene sets associated with antiestrogen resistance in MCF7FGFR1(SP-)(NLS) cells.

Close modal

We next correlated FGFR1 localization with response to estrogen suppression. Nuclear FGFR1 levels in these tumors positively correlated with post-letrozole Ki67 values (Fig. 1B, P = 0.034 and S2A, P = 0.0055). Although nuclear FGFR1 levels tended to be higher in FGFR1-amplified compared with non-FGFR1–amplified tumors, this difference was not statistically significant (Supplementary Fig. S2B, P = 0.096). Furthermore, analysis of only those tumors without FGFR1 amplification (n = 136) revealed that high nuclear FGFR1 correlated with higher on-treatment Ki67 (Supplementary Fig. S2C, P = 0.027), suggesting that nuclear FGFR1 protein levels are associated with antiestrogens resistance regardless of FGFR1 gene copy number. We evaluated the distribution of nuclear FGFR1 H-score according to the Ki-67 defined groups. Nuclear FGFR1 H-Score values were borderline higher in resistant samples (Supplementary Fig. S2D; P = 0.056). We next divided the cohort in two groups based on a median nuclear FGFR1 H-Score of 156.77. We found that 21 of 78 (27%) tumors with nuclear FGFR1 H-Score above the median were resistant (Supplementary Fig. S2E; χ2 test = 0.04).

Next, to explore a causal association of nuclear FGFR1 with antiestrogen resistance, we transduced MCF-7 cells with an FGFR1(SP-)(NLS) plasmid (ref. 34; Supplementary Fig. S2F). The lack of an N-terminal Signal Peptide (SP-) sequence and the presence of a strong Nuclear Localization Signal (NLS; PKKKRKV) are intended to prevent membrane localization of the resulting protein, while enforcing its nuclear localization. Indeed, the levels of FGFR1 in the nucleus of MCF-7 cells were increased upon transfection with FGFR1(SP-)(NLS), as shown by subcellular fractionation assay (Supplementary Fig. S2G) and immunofluorescence analysis (Supplementary Fig. S2H). MCF7FGFR1(SP-)(NLS) cells exhibited higher estrogen-independent growth (Fig. 1C) and reduced sensitivity to fulvestrant (Fig. 1D and E) compared with control cells. We next implanted MCF7EV and MCF7FGFR1(SP-)(NLS) cells in ovariectomized athymic mice. In agreement with the in vitro findings, MCF7FGFR1(SP-)(NLS) xenografts exhibited higher estrogen-independent growth compared to control tumors (Fig. 1F). Further, while fulvestrant treatment induced tumor regression in control tumors, MCF7FGFR1(SP-)(NLS) xenografts did not regress (Fig. 1F). Gene Set Enrichment Analysis (GSEA) of the 3,014 differentially expressed genes (DEGs) in MCF7FGFR1(SP-)(NLS) vs. MCF7EV cells (shown in Supplementary Fig. S2I) also revealed enrichment of gene sets associated with resistance to endocrine therapy and to tamoxifen (Fig. 1G; Supplementary Table S3), suggesting that nuclear FGFR1 overexpression induces a transcriptional profile causal to antiestrogen resistance in breast cancer. Previous studies reported that Granzyme B cleaves FGFR1 at Asp432, producing a C-terminal fragment that translocates to the nucleus (15). Hence, we transduced MCF7 cells with an FGFR1 D432N mutant, previously reported as resistant to Granzyme B cleavage (15), to investigate a potential effect on antiestrogen resistance. Nuclear fractions of MCF7FGFR1 D432N cells showed markedly higher expression of both full-length (∼140 KDa) and the cleaved form (∼55 KDa) of FGFR1 (Supplementary Fig. S3A), suggesting that Granzyme B does not contribute to FGFR1 cleavage in these cells. Finally, we observed that overexpression of FGFR1 wild-type, FGFR1 D432N and FGFR1(SP-)(NLS) in MCF7 cells resulted in enhanced estradiol-independent growth and reduced sensitivity to fulvestrant and tamoxifen (Supplementary Fig. S3B–S3E).

Nuclear FGFR1 occupies chromatin at transcription start sites of genes associated with antiestrogen resistance

We next examined a role of nuclear FGFR1 in transcriptional regulation in breast cancer. To study chromatin occupancy by FGFR1, we performed ChIP-seq in CAMA1 cells cultured for 48 hours in estradiol-deprived media, which identified 4,412 high-confidence binding sites. (Fig. 2A). De novo consensus motif discovery analysis of the FGFR1 DNA-binding sites revealed GC-rich motifs among the top enriched motifs (Fig. 2B; Supplementary Fig. S4A).

Figure 2.

Nuclear FGFR1 occupies chromatin at TSSs of genes associated with antiestrogen resistance. A, Heatmap showing ChIP-seq read densities around the FGFR1-bound regions in CAMA1 cells, in two replicates. The x-axis represents read densities within 5-kb region around the peak summit; the y-axis represents each predicted binding site. Cells were cultured for 48 hours in IMEM/10% CSS. B, Top consensus motifs, identified by de novo motif analysis, at genomic loci bound by FGFR1 identified in A. Statistical significance expressed as P value for each motif is shown. C, Plot representing the density of the FGFR1 distribution around the TSSs in the two ChIP-seq replicates. D, Genomic annotations for the FGFR1-binding sites identified by ChIP-seq showing enrichment of promoter regions. E, Expression level of genes whose promoter is bound by FGFR1 (red; 2,704 genes) versus gene expression of all other genes (blue; 10,452 genes). Data are derived from RNA-seq on CAMA1 cells cultured for 48 hours in IMEM/10% CSS (**** P < 0.0001, Mann–Whitney test). F, Gene expression of genes whose promoter is bound by FGFR1 in CAMA1 cells upon siRNA-mediated FGFR1 knockdown (n = 3; **** P < 0.0001, two-tailed Wilcoxon test). G, Disease-free survival (DFS) and overall survival (H) of the METABRIC cohort of 950 patients with ER+/HER2 breast cancer treated with antiestrogens as a function of the FGFR1-associated polygenic score (lowest quartile vs. highest quartile). The signature score was calculated by GSVA (77).

Figure 2.

Nuclear FGFR1 occupies chromatin at TSSs of genes associated with antiestrogen resistance. A, Heatmap showing ChIP-seq read densities around the FGFR1-bound regions in CAMA1 cells, in two replicates. The x-axis represents read densities within 5-kb region around the peak summit; the y-axis represents each predicted binding site. Cells were cultured for 48 hours in IMEM/10% CSS. B, Top consensus motifs, identified by de novo motif analysis, at genomic loci bound by FGFR1 identified in A. Statistical significance expressed as P value for each motif is shown. C, Plot representing the density of the FGFR1 distribution around the TSSs in the two ChIP-seq replicates. D, Genomic annotations for the FGFR1-binding sites identified by ChIP-seq showing enrichment of promoter regions. E, Expression level of genes whose promoter is bound by FGFR1 (red; 2,704 genes) versus gene expression of all other genes (blue; 10,452 genes). Data are derived from RNA-seq on CAMA1 cells cultured for 48 hours in IMEM/10% CSS (**** P < 0.0001, Mann–Whitney test). F, Gene expression of genes whose promoter is bound by FGFR1 in CAMA1 cells upon siRNA-mediated FGFR1 knockdown (n = 3; **** P < 0.0001, two-tailed Wilcoxon test). G, Disease-free survival (DFS) and overall survival (H) of the METABRIC cohort of 950 patients with ER+/HER2 breast cancer treated with antiestrogens as a function of the FGFR1-associated polygenic score (lowest quartile vs. highest quartile). The signature score was calculated by GSVA (77).

Close modal

Most mammalian promoters (i.e., regions upstream the transcription start sites; TSS) are enriched with GC-rich sequences, also named CpG islands (CGI; refs. 36, 37). Also, many human active promoters are enriched with SP1, Nrf-1, E2F, and ETS transcription factor–binding sites motifs, all containing CGIs (38, 39). Consistent with the predominance of GC-rich motifs identified (Fig. 2B; Supplementary Fig. S4A), the majority of FGFR1-binding sites were near the TSS (Fig. 2C), with 67% of ChIP-Seq peaks occurring at gene promoters (Fig. 2D). On the basis of the nucleotide sequence of FGFR1-bound genomic loci identified by ChIP-Seq (Supplementary Fig. S5A), we next designed specific primer pair sequences to validate FGFR1 occupancy at nine selected promoter regions in CAMA1 (Supplementary Fig. S5B) and MDA-MB-134 cells (Supplementary Fig. S5C) by ChIP-qPCR. We further examined the binding of FGFR1 to the same promoter regions in tumor epithelial cells isolated by FACS from disaggregated HCI-011 PDXs (Supplementary Fig. S5D). For this, we performed ChIP-qPCR on the isolated tumor epithelial cells and confirmed enrichment of FGFR1 at similar promoter regions examined in CAMA1 and MDA-MB-134 cells (Supplementary Fig. S5E). These results suggest that FGFR1 binds to gene promoters and support a role for FGFR1 in transcriptional regulation.

RNA-Seq analysis of CAMA1 cells cultured for 48 hours in estrogen-deprived media revealed that the average expression of genes whose promoter was bound by FGFR1 was significantly higher than the expression of all other genes (Fig. 2E; Mann–Whitney, P < 0.0001), implying that FGFR1 target genes are actively transcribed. Silencing of FGFR1 with siRNA significantly reduced the expression of genes whose promoters is bound by FGFR1 (Fig. 2F, Wilcoxon test P < 0.0001).

We intersected the list of 2,901 genes whose promoter was bound by FGFR1, as identified by ChIP-seq, with the 2,048 genes downregulated by FGFR1 siRNA in CAMA1 cells (Supplementary Fig. S6A). A polygenic risk score was generated based on the 561 overlapping genes (Supplementary Fig. S6A; Supplementary Table S5). Analysis of 212 ER+/HER2 primary breast tumors in TCGA with known 8p11.23 copy number, where the FGFR1 locus lies, showed that the signature score was enriched in cancers with 8p11.23 amplification, compared with all other samples not harboring 8p11.23 copy gain or amplification (Supplementary Fig. S6B; t test, P = 0.039). In the METABRIC cohort of 1,294 ER+/HER2− breast cancers, the signature score was enriched in tumors exhibiting 8p11.23 amplification, compared with cancers without FGFR1 copy gain or amplification (Supplementary Fig. S6C; t-test, P = 0.038). These results suggest that FGFR1-induced transcriptional signature is statistically enriched in primary breast tumors harboring FGFR1 gene amplification. We next investigated whether the FGFR1-associated risk score correlated with antiestrogen resistance. In the METABRIC cohort of ER+ breast cancers treated with antiestrogens (n = 950), the risk score correlated with a worse DFS (Fig. 2G; log rank P = 0.0103, HR = 1.57; 95% CI, 1.12–2.21, median 278.4 vs. not reached) and OS (Fig. 2H; log rank P < 0.0001; HR = 1.72; 95% CI, 1.35 – 2.18, median 126.5 vs. 188.2 months). Taken together, these findings further suggest a role for nuclear FGFR1 in resistance to endocrine therapy.

FGFR1 associates with RNA polymerase II and promotes gene expression

To further dissect the function of FGFR1 in transcriptional regulation, we sought to uncover which components of the transcriptional machinery physically interact with nuclear FGFR1. For this purpose, we transduced CAMA1 cells with an FGFR1–3XFLAG lentiviral vector. Next we immunoprecipitated FLAG from mixed nuclear and chromatin fractions of CAMA1FGFR1–3XFLAG cells followed by mass spectrometry (MS) of the antibody pulldowns. Among the proteins recovered in the nuclear FGFR1 interactome, we found several RNA polymerase II (Pol II) subunits (Fig. 3A and B; Supplementary Table S6). RNA Pol II plays a pivotal role in gene transcription, which also requires other proteins critical for RNA polymerase recruitment to promoters, transcriptional initiation, and elongation (40). These findings were in line with enrichment of FGFR1 binding near transcription start sites (shown in Fig. 2C), which are also the sites where Pol II is recruited for transcription initiation. We confirmed the FGFR1-Pol II association by coprecipitation with an FGFR1 antibody followed by Pol II Western blot analysis of CAMA1 and MDA-MB-134 cell lysates (Fig. 3C; Supplementary Fig. S7A) and of HCI-011 PDX lysates (Supplementary Fig. S7B). FGFR1-Pol II coprecipitation was observed when we performed sonication and nuclease treatment, allowing the release of more soluble Pol II.

Figure 3.

FGFR1 associates with RNA polymerase II and promotes gene expression. A, FLAG-immunoprecipitation followed by MS was performed on sonicated and nuclease-treated mixed nuclear and chromatin fractions of CAMA1FGFR13XFLAG cells. FGFR1 and RNA Polymerase II subunits are highlighted (red). The Western blot on the right shows FLAG-FGFR1 overexpression in CAMA1FGFR13XFLAG cells. B, List of RNA polymerase II subunits that coprecipitate with FLAG-FGFR1. Coverage (%) indicates the percentage of the protein sequence that was covered by the identified peptides. #PSMs indicates the number of Peptide Spectrum Matches or the number of spectra assigned to peptides that contributed to inference of the protein. C, Coprecipitation of FGFR1 and RNA Polymerase phosphorylated on Serine residue 5 (Pol II S5P), with or without sonication and nuclease treatment. D, Heatmaps of ChIP-seq read densities around the FGFR1 (red) and Pol II S5P bound regions (green) in CAMA1 cells. E, Venn diagram and (F) density plots related to FGFR1- and Pol II S5P-ChIP-seq shown in D. G, Distribution of FGFR1 and Pol II S5P binding peaks at the VEGFA and CCND1 promoters (UCSC genome browser). H, ChIP-reChIP assay performed by sequential ChIP-qPCR with an FGFR1 antibody followed by an antibody against Pol II S5P (Rpb1) or normal rabbit IgG (control), at the VEGFA and CCND1 promoters. Enrichment values expressed as percent (%) of input. I, Venn diagram of FGFR1-bound genes by ChIP-seq in MCF7FGFR1(SP-)(NLS) cells (red circle), and genes upregulated (green), and downregulated (blue) in MCF7FGFR1(SP-)(NLS) versus MCF7EV cells. J, Prediction of activating or repressing transcription function of FGFR1 by Binding and Expression Target Analysis (BETA). ChIP-Seq and RNA-seq data from I were integrated. Genomic regions bound by FGFR1 are predicted to regulate the expression of upregulated genes, but not of downregulated genes. P values indicate the significance of the associations compared with background nonregulated genes.

Figure 3.

FGFR1 associates with RNA polymerase II and promotes gene expression. A, FLAG-immunoprecipitation followed by MS was performed on sonicated and nuclease-treated mixed nuclear and chromatin fractions of CAMA1FGFR13XFLAG cells. FGFR1 and RNA Polymerase II subunits are highlighted (red). The Western blot on the right shows FLAG-FGFR1 overexpression in CAMA1FGFR13XFLAG cells. B, List of RNA polymerase II subunits that coprecipitate with FLAG-FGFR1. Coverage (%) indicates the percentage of the protein sequence that was covered by the identified peptides. #PSMs indicates the number of Peptide Spectrum Matches or the number of spectra assigned to peptides that contributed to inference of the protein. C, Coprecipitation of FGFR1 and RNA Polymerase phosphorylated on Serine residue 5 (Pol II S5P), with or without sonication and nuclease treatment. D, Heatmaps of ChIP-seq read densities around the FGFR1 (red) and Pol II S5P bound regions (green) in CAMA1 cells. E, Venn diagram and (F) density plots related to FGFR1- and Pol II S5P-ChIP-seq shown in D. G, Distribution of FGFR1 and Pol II S5P binding peaks at the VEGFA and CCND1 promoters (UCSC genome browser). H, ChIP-reChIP assay performed by sequential ChIP-qPCR with an FGFR1 antibody followed by an antibody against Pol II S5P (Rpb1) or normal rabbit IgG (control), at the VEGFA and CCND1 promoters. Enrichment values expressed as percent (%) of input. I, Venn diagram of FGFR1-bound genes by ChIP-seq in MCF7FGFR1(SP-)(NLS) cells (red circle), and genes upregulated (green), and downregulated (blue) in MCF7FGFR1(SP-)(NLS) versus MCF7EV cells. J, Prediction of activating or repressing transcription function of FGFR1 by Binding and Expression Target Analysis (BETA). ChIP-Seq and RNA-seq data from I were integrated. Genomic regions bound by FGFR1 are predicted to regulate the expression of upregulated genes, but not of downregulated genes. P values indicate the significance of the associations compared with background nonregulated genes.

Close modal

FGFR1 interacted with Pol II phosphorylated on Ser5 and Ser2 of the C-terminal domain (CTD) of the largest Pol II subunit, RPB1 (Pol II S5P, Pol II S2P; Fig. 3C; Supplementary Fig. S7A–S7C). Phosphorylation of Ser5 and Ser2 on the CTD are posttranslational modifications required for transcription initiation and elongation, respectively (41, 42). MS analysis of CAMA1FGFR1–3XFLAG pulldowns also revealed proteins known to be directly involved in Pol II phosphorylation and subsequent transcriptional activity, such as CDK7 and CDK9 (43) (Supplementary Table S6). We next confirmed the association of FGFR1 with CDK7-CDK9 in CAMA1 and MDA-MB-134 cells by coprecipitation with an FGFR1 antibody (Supplementary Fig. S7D). The binding of FGFR1 with Pol II S5P, Pol II S2P, CDK7 and CDK9 was also confirmed by immunoprecipitation assays with Pol II S5P, Pol II S2P, CDK7, and CDK9 antibodies, followed by Western blot analysis for FGFR1, in CAMA1 cells (Supplementary Fig. S7E). Comparison of the genomic coordinates of FGFR1 occupancy with data from previously published ChIP-seq datasets for transcription factors (44) identified RNA Polymerase II subunit B1, RPB1 (POLR2A), as one of the top hits (Supplementary Fig. S8A). Using ChIP-seq to evaluate the distribution of Pol II S5P across the CAMA1 cells genome, we found that 2867/4412 (65%) FGFR1-binding peaks overlapped with Pol II S5P–binding sites (Fig. 3D and E), with distribution around similar locations near the TSS (Fig. 3F and G; Supplementary Fig. S8B), further suggesting that FGFR1 genomic distribution occurs at actively transcribed loci. ChIP-reChIP assays, with sequential FGFR1 and Pol II S5P antibodies, confirmed FGFR1-Pol II S5P association at selected gene promoters in CAMA1 (Fig. 3H; Supplementary Fig. S8C) and MDA-MB-134 cells (Supplementary Fig. S8D). In addition, by qRT-PCR we confirmed that silencing FGFR1 by siRNA in CAMA1 cells resulted in downregulation of selected FGFR1 target genes (Supplementary Fig. S8E). Taken together, these findings suggest that FGFR1 is part of a chromatin-associated Pol II complex engaged in active transcription.

Using MCF7FGFR1(SP-)(NLS) cells, we next investigated the effect of exogenous nuclear FGFR1 in gene transcription. By ChIP-seq, we identified 24,475 FGFR1 DNA–binding sites in MCF7FGFR1(SP-)(NLS) cells, corresponding to 11,840 genes (Fig. 3I). We intersected the FGFR1-bound genes with the differentially regulated transcripts (FDR<0.05) in MCF7FGFR1(SP-)(NLS) versus MCF7EV cells (shown in Supplementary Fig. S2G), finding that 1,009 of the 1,454 upregulated genes and 1,141 of the 1,560 downregulated genes were also bound by FGFR1 at a genomic level (Fig. 3I). To examine whether nuclear FGFR1 directly affects gene expression in MCF7FGFR1(SP-)(NLS) cells, we used Binding and Expression Target Analysis (BETA) platform, a computational program which integrates ChIP-seq and RNA-seq data (45). Results from this analysis suggested a direct role for genomic-bound FGFR1 in activating gene expression (P = 8.01e-06; Fig. 3J). Conversely, there was no significant association between FGFR1-bound loci and repressed genes (P = 0.139; Fig. 3J). Next, to further demonstrate a direct role for genomic-bound FGFR1 in activating gene transcription, we first silenced the expression of FGFR1 using CRISPR/Cas9 in MCF7 cells (MCF7FGFR1 KO; Supplementary Fig. S8F). Next, we rescued the expression of FGFR1(SP-)(NLS), FGFR1 D432N and wild-type FGFR1 (Supplementary Fig. S8G). By RT-qPCR, we demonstrated that re-expression of FGFR1(SP-)(NLS), FGFR1 D432N and wild-type FGFR1 in MCF7FGFR1 KO cells induced the expression of selected FGFR1 target genes (Supplementary Fig. S8H). GSEA of the 1,009 genes upregulated in MCF7FGFR1(SP-)(NLS) and bound by FGFR1 at a genomic level (shown in Fig. 3I) revealed a strong enrichment of estrogen response early (FDR = 2.29e-44) and late (FDR = 6.43e-33) genes (Supplementary Fig. S9A), suggesting that nuclear FGFR1 induces an ERα-associated transcriptional profile that may contribute to endocrine resistance. By qRT-PCR, we confirmed the upregulation of canonical ERα-regulated genes in MCF7FGFR1(SP-)(NLS) cells compared with control cells (Supplementary Fig. S9B). We hypothesized that FGFR1 may regulate the expression of canonical ERα-regulated genes independent of ERα transcription. To support this, we performed qRT-PCR analysis of xenografts treated with fulvestrant (Fig. 1F). In both MCF7EV and MCF7FGFR1(SP-)(NLS) xenografts, fulvestrant treatment suppressed ERα protein expression by IHC (Supplementary Fig. S9C). Expression of canonical ERα target genes was downregulated in fulvestrant-treated MCF7EV xenografts but not in fulvestrant-treated MCF7FGFR1(SP-)(NLS) xenografts (Supplementary Fig. S9D), suggesting that nuclear FGFR1 induces the expression of an ERα-associated signature independent of ERα transcription.

Nuclear FGFR1 overlaps with active transcription histone marks

To examine if FGFR1 is associated with active gene transcription, we first superimposed the genomic coordinates of FGFR1-binding sites in CAMA1 cells with those from publicly available ChIP-seq datasets for histone marks (44). This analysis revealed a strong enrichment for active transcription histone marks, such as H3K4me3 and H3K27ac (Supplementary Fig. S10A). H3K4me3 is commonly enriched around TSS of actively transcribed genes (46). Similarly, H3K27 acetylation is a marker of active regulatory elements and found at both promoters and enhancers (47). We next analyzed the genomic distribution of H3K27ac in CAMA1 cells by ChIP-seq. Notably, 3961/4175 (94.8%) of FGFR1 peaks overlapped with H3K27ac, suggesting that FGFR1 binding to DNA occurs at genomic loci that are enriched for promoters of actively transcribed genes (Fig. 4A and B). We also investigated the distribution of the H3K4me3 histone mark in CAMA1 cells. In agreement with the finding of FGFR1 enrichment at TSSs/promoter regions (shown in Fig. 2C and D), we identified a strong overlap between FGFR1 and H3K4me3 (4,055/4,265 sites or 95%), implying that FGFR1 is occupying sites engaged in transcription initiation (Fig. 4C and 4D). Taken together, our findings suggest that FGFR1 is recruited to genomic regions where epigenetic modifications, such as H3K27 acetylation and H3K4 trimethylation, cooccur to promote a chromatin state accessible for active transcription complexes at promoter regions (Fig. 4E).

Figure 4.

Nuclear FGFR1 overlaps with active transcription histone marks. A, Heatmaps of ChIP-seq read densities around the FGFR1-bound regions (red) and areas of increased acetylation of H3K27 (H3K27ac, blue) in CAMA1 cells. Two replicates for each antibody are shown. B, Venn diagram of FGFR1 ChIP-Seq peaks and H3K27ac regions. C, Heatmaps of ChIP-Seq read densities around the FGFR1-bound regions (red) and areas of increased trimethylation of H3K4 (H3K4me3; yellow) in CAMA1 cells. Two replicates for each antibody are shown. D, Venn diagram of FGFR1 ChIP-Seq peaks and H3K4me3 regions. E, Distribution of FGFR1-binding peaks, H3K27ac and H3K4me3 histone marks at the CREEBP, KDM4B, and CDK12 promoters (UCSC genome browser).

Figure 4.

Nuclear FGFR1 overlaps with active transcription histone marks. A, Heatmaps of ChIP-seq read densities around the FGFR1-bound regions (red) and areas of increased acetylation of H3K27 (H3K27ac, blue) in CAMA1 cells. Two replicates for each antibody are shown. B, Venn diagram of FGFR1 ChIP-Seq peaks and H3K27ac regions. C, Heatmaps of ChIP-Seq read densities around the FGFR1-bound regions (red) and areas of increased trimethylation of H3K4 (H3K4me3; yellow) in CAMA1 cells. Two replicates for each antibody are shown. D, Venn diagram of FGFR1 ChIP-Seq peaks and H3K4me3 regions. E, Distribution of FGFR1-binding peaks, H3K27ac and H3K4me3 histone marks at the CREEBP, KDM4B, and CDK12 promoters (UCSC genome browser).

Close modal

FOXA1 mediates FGFR1 recruitment to chromatin

We next sought to identify mechanisms regulating FGFR1 recruitment to chromatin. To this end, we further examined the list of FGFR1-interacting proteins identified by MS of FLAG antibody pulldowns from mixed nuclear and chromatin bound CAMA13XFLAG-FGFR1 cell fractions (Supplementary Table S6). Our attention was captured by FOXA1 (Forkhead Box A1; Fig. 5A), a master regulator that modifies chromatin accessibility for transcription factors (48). FOXA1 has been shown to be required for estrogen receptor binding to chromatin and to induce antiestrogen resistance in ER+ breast cancer (49, 50). FOXA1 also influences the genomic binding of other transcription factors, such as the androgen receptor in prostate cancer, contributing to castration resistance in these tumors (51). FOXA1 is required for the survival of ER+ breast cancer cells, as demonstrated by Dependency Map Project (https://depmap.org/portal/; ref. 52; Supplementary Fig. S11A), and has the ability to promote the association of transcription factors with chromatin (53). Thus, we sought to determine whether FOXA1 affects FGFR1 transcriptional function.

Figure 5.

FOXA1 mediates FGFR1 recruitment to chromatin. A, Immunoprecipitation of FLAG-FGFR1 followed by mass spectrometry was conducted on mixed nuclear plus chromatin-bound fractions of CAMA13XFLAG-FGFR1 cells. The plot shows the enrichment of FOXA1 in the FLAG-FGFR1-bound fraction compared to IgG control. B, Coprecipitation of FGFR1 and FOXA1, followed by Western blot analysis, in CAMA1 (left) and MDA-MB-134 (right) cells. C, Venn diagram of FGFR1 peaks (red circle), identified by ChIP-seq (in Fig. 2A) and FOXA1 DNA–binding loci identified by ChIP-Seq (blue circle) in CAMA1 cells. D, ChIP-reChIP assay performed by sequential ChIP-qPCR with a FGFR1 antibody followed by an antibody against FOXA1 or normal rabbit IgG, at the BRD2, CDK12, VEGFA, and JUNB promoters. Enrichment values expressed as percent (%) of input. E, Heatmaps of ChIP-seq read densities around the FGFR1 bound regions in CAMA1 cells transfected with nontargeting (siNT) or FOXA1 siRNAs. Twenty-four hours posttransfection, dishes were replenished with IMEM/10% CSS, and cells were collected 48 hours later for ChIP. F, ChIP-qPCR confirmation of reduced FGFR1 binding at selected genomic loci, VEGFA (t test, P = 0.0071), BRD2 (P = 0.0145), CDK12 (P = 0.0114), JUNB (P = 0.0252) promoters, upon siRNA-mediated FOXA1 knockdown in CAMA1 cells. Enrichment values expressed as percent (%) of input.

Figure 5.

FOXA1 mediates FGFR1 recruitment to chromatin. A, Immunoprecipitation of FLAG-FGFR1 followed by mass spectrometry was conducted on mixed nuclear plus chromatin-bound fractions of CAMA13XFLAG-FGFR1 cells. The plot shows the enrichment of FOXA1 in the FLAG-FGFR1-bound fraction compared to IgG control. B, Coprecipitation of FGFR1 and FOXA1, followed by Western blot analysis, in CAMA1 (left) and MDA-MB-134 (right) cells. C, Venn diagram of FGFR1 peaks (red circle), identified by ChIP-seq (in Fig. 2A) and FOXA1 DNA–binding loci identified by ChIP-Seq (blue circle) in CAMA1 cells. D, ChIP-reChIP assay performed by sequential ChIP-qPCR with a FGFR1 antibody followed by an antibody against FOXA1 or normal rabbit IgG, at the BRD2, CDK12, VEGFA, and JUNB promoters. Enrichment values expressed as percent (%) of input. E, Heatmaps of ChIP-seq read densities around the FGFR1 bound regions in CAMA1 cells transfected with nontargeting (siNT) or FOXA1 siRNAs. Twenty-four hours posttransfection, dishes were replenished with IMEM/10% CSS, and cells were collected 48 hours later for ChIP. F, ChIP-qPCR confirmation of reduced FGFR1 binding at selected genomic loci, VEGFA (t test, P = 0.0071), BRD2 (P = 0.0145), CDK12 (P = 0.0114), JUNB (P = 0.0252) promoters, upon siRNA-mediated FOXA1 knockdown in CAMA1 cells. Enrichment values expressed as percent (%) of input.

Close modal

Coimmunoprecipitation of CAMA1 and MDA-134 cell lysates with an FGFR1 antibody followed by FOXA1 western showed association of FOXA1 with FGFR1 (Fig. 5B; Supplementary Fig. S11B). We next investigated the FOXA1 genomic distribution in CAMA1 cells by ChIP-seq. We found that 1,669 of 4,438 (38%) of FGFR1 DNA–binding sites were shared with FOXA1 (Fig. 5C; Supplementary Fig. S11C). ChIP-reChIP experiments, conducted by sequential chromatin immunoprecipitation of FGFR1 and FOXA1, confirmed that these two proteins are both present in the transcriptional complex binding at selected genomic loci in CAMA1 and MDA-MB-134 cells (Fig. 5D; Supplementary Fig. S11D). Finally, we silenced FOXA1 in CAMA1 cells using siRNA. Notably, 52% of FGFR1-binding peaks were reduced or lost upon FOXA1 suppression (Fig. 5E). Using ChIP-qPCR, we confirmed that FGFR1 recruitment to target genomic loci was significantly reduced upon FOXA1 knockdown in CAMA1 and MDA-MB-134 cells (Fig. 5F; Supplementary Fig. S11E). Expression of FGFR1 target genes was reduced by siFOXA1 as measured by RT-qPCR (Supplementary Fig. S11F). Of note, FOXA1 knockdown did not have any effect on FGFR1 expression or protein levels (Supplementary Fig. S11G). Altogether, these data suggest a role for FOXA1 in regulating FGFR1 recruitment to chromatin and, consequently, its transcriptional activity.

In addition, because the de novo motif analysis of the FGFR1 DNA–binding sites (shown in Fig. 2B) identified SP1 as the top enriched motif, we investigated whether SP1 influenced the genomic distribution of FGFR1. We silenced SP1 expression by siRNA in CAMA1 cells (Supplementary Fig. S12A). ChIP-qPCR assays revealed that siRNA-mediated SP1 knockdown did not inhibit FGFR1 recruitment to gene promoters (Supplementary Fig. S12B) identified by FGFR1 ChIP-seq (shown in Fig. 2A and listed in Supplementary Table S2). Furthermore, SP1 knockdown enhanced the FGFR1 binding at CCND1, CDK12, JUNB, and PPP6R3 promoters (Supplementary Fig. S12B). Concordant with these results, the expression of these genes was increased upon SP1 knockdown (Supplementary Fig. S12C). These results may suggest a scenario in which SP1 competes with FGFR1 to bind to selected genomic loci, with SP1 having a potential repressive action on FGFR1 recruitment to chromatin and consequent transcriptional activity.

Nuclear FGFR1 activity is enhanced by FGF2 but not affected by receptor tyrosine kinase blockade

Previous studies have shown that growth factors such as EGF and insulin can induce nuclear trafficking of EGFR and insulin receptor (IR), respectively, enhancing their transcriptional activity (19, 21). FGF2 and FGFR1 have been shown to colocalize in the nucleus of pancreatic cancer cells, promoting an invasive phenotype (14). In contrast to FGFR1, FGF2 contains an NLS peptide sequence, which may promote importin β–mediated translocation of FGFR1 into the nucleus (54). Thus, we examined whether growth factor stimulation influences FGFR1 transcriptional function in breast cancer cells. Stimulation with FGF2 increased FGFR1 levels in subcellular nuclear fractions of CAMA1 cells (Supplementary Fig. S13A) as well as the FGFR1-Pol II S5P association measured by coimmunoprecipitation followed by immunoblot analysis (Supplementary Fig. S13B). To further characterize the effect of FGF2 on FGFR1-mediated transcription, we first performed RNA-sequencing analysis of FGF2-stimulated CAMA1 cells. This analysis identified 7,923 differentially expressed (DEGs) as a function of FGF2 treatment (FDR <0.05; Fig. 6A). Integration of the RNA-sequencing data and FGFR1 ChIP-seq results from FGF2-stimulated CAMA1 cells (shown in Supplementary Fig. S13C) with the BETA platform (45) suggested that FGFR1 DNA–binding sites have a significant role in promoting FGF2-induced transcription (Fig. 6B, P = 0.029). Conversely, there was not a statistical association between FGFR1 DNA-binding sites and FGF2-repressed genes (P = 0.392). We also employed the Genomic Regions Enrichment of Annotations Tool (GREAT; ref. 55) to interpret the biological functions of FGFR1 DNA–binding sites identified by ChIP-seq in FGF2-treated CAMA1 cells. Consistent with the BETA platform, genes associated with FGFR1 DNA–binding sites were highly enriched for signatures of cellular response to growth factor stimulation and activation of the FGFR pathway (Fig. 6C; Supplementary Fig. S7; Supplementary Table S7).

Figure 6.

Nuclear FGFR1 activity is enhanced by FGF2 but not affected by RTK blockade. A, MA plot of M (log ratio) versus A (log mean average) normalized counts, showing differentially expressed genes with P <1e-5 in response to FGF2 (5 ng/mL, 6 hours). B, Prediction of activating or repressing transcription function of FG2-stimulated FGFR1 by BETA platform. FGFR1 ChIP-seq results from CAMA1 cells stimulated with FGF2 (5 ng/mL, 3 hours) and RNA-seq data from A were integrated. Genomic regions bound by FGFR1 are predicted to modulate expression of genes upregulated upon FGF2 stimulation, but not downregulated genes; P values indicate the significance of the associations compared with background nonregulated genes. C, GREAT analysis the FGFR1 DNA–binding sites identified by FGFR1 ChIP-seq in CAMA1 cells stimulated with FGF2 (5 ng/mL, 3 hours). Top ten enriched Gene Ontology (GO) Biological Processes and binomial FDR values are shown. D, Subcellular fractionation on CAMA1 cells, showing FGFR1 in cytosolic, membrane, soluble nuclear, and chromatin-bound compartments, after treatment with FGF2 (5 ng/mL, 3 hours) ± erdafitinib (250 nmol/L, 3 hours). Red box shows longer exposure. E, Coprecipitation of FGFR1 and Pol II S5P from sonicated and nuclease-treated lysates of CAMA1 cells ± FGF2 (5 ng/mL, 3 hours) ± erdafitinib (250 nmol/L, 3 hours). F, Heatmaps of ChIP-seq read densities around the FGFR1-bound regions in CAMA1 cells treated with FGF2 (5 ng/mL, 3 hours) ± erdafitinib (250 nmol/L, 3 hours).

Figure 6.

Nuclear FGFR1 activity is enhanced by FGF2 but not affected by RTK blockade. A, MA plot of M (log ratio) versus A (log mean average) normalized counts, showing differentially expressed genes with P <1e-5 in response to FGF2 (5 ng/mL, 6 hours). B, Prediction of activating or repressing transcription function of FG2-stimulated FGFR1 by BETA platform. FGFR1 ChIP-seq results from CAMA1 cells stimulated with FGF2 (5 ng/mL, 3 hours) and RNA-seq data from A were integrated. Genomic regions bound by FGFR1 are predicted to modulate expression of genes upregulated upon FGF2 stimulation, but not downregulated genes; P values indicate the significance of the associations compared with background nonregulated genes. C, GREAT analysis the FGFR1 DNA–binding sites identified by FGFR1 ChIP-seq in CAMA1 cells stimulated with FGF2 (5 ng/mL, 3 hours). Top ten enriched Gene Ontology (GO) Biological Processes and binomial FDR values are shown. D, Subcellular fractionation on CAMA1 cells, showing FGFR1 in cytosolic, membrane, soluble nuclear, and chromatin-bound compartments, after treatment with FGF2 (5 ng/mL, 3 hours) ± erdafitinib (250 nmol/L, 3 hours). Red box shows longer exposure. E, Coprecipitation of FGFR1 and Pol II S5P from sonicated and nuclease-treated lysates of CAMA1 cells ± FGF2 (5 ng/mL, 3 hours) ± erdafitinib (250 nmol/L, 3 hours). F, Heatmaps of ChIP-seq read densities around the FGFR1-bound regions in CAMA1 cells treated with FGF2 (5 ng/mL, 3 hours) ± erdafitinib (250 nmol/L, 3 hours).

Close modal

The effects of FGF2 on nuclear translocation and transcriptional regulation suggested a causal link with the FGFR1 tyrosine kinase activity. Thus, we investigated the effect of the pan-FGFR tyrosine kinase inhibitor (TKI) erdafitinib (56). Treatment with erdafitinib abolished FGFR1, FRS2α, and ERK phosphorylation/activation (Supplementary Fig. S13D) but did not abrogate the FGF2-induced increase in FGFR1 levels in CAMA1 (Fig. 6D) and MDA-MB-134 (Supplementary Fig. S13E) cell nuclear fractions. Furthermore, treatment with erdafitinib did not affect the FGF2-induced association of FGFR1 with Pol II S5P association (Fig. 6E; Supplementary Fig. S13F) nor the FGFR1 genomic distribution as shown by ChIP-seq of CAMA1 cells (Fig. 6F; Supplementary Fig. S13G). Finally, treatment with erdafitinib did not impair growth of MCF7FGFR1(SP-)(NLS) cells, compared with MCF7EV cells (Supplementary Fig. S13H), suggesting the FGFR1 activated TK is not causally associated with the receptor's activity in the nucleus.

FGFR signaling has multifaceted roles in normal physiology and development, such as embryogenesis, organ development, tissue repair, remodeling, angiogenesis, and metabolism (57–60). Gene rearrangements, activating tyrosine kinase domain mutations, and fusions in FGFR genes are oncogenic drivers in several cancer types, such as breast, bladder, and biliary duct tumors (61–63). The canonical FGFR pathway is triggered by activation of the receptor's tyrosine kinase at the plasma membrane followed by induction of signaling nodes, primarily RAS/RAF/MEK/ERK and PI3K/AKT. However, recent work has shown the presence of FGFR receptors, such as FGFR2, in the nucleus of cancer cells where it may affect gene transcription (64). In medulloblastoma cells, nuclear FGFR1 has been shown to bind CREB-Binding Protein (CBP) and modulate its transcriptional activity (17). We found FGFR1 abundance in the nucleus of breast cancer cells in ER+ primary tumors and PDXs (Fig. 1A; Supplementary Fig. S1C and S1E). Therefore, we investigated herein the role of nuclear FGFR1, which may shed light on mechanisms by which aberrant FGFR signaling drives cancer progression and its therapeutic implications.

FGFR1 overexpression has been shown to induce antiestrogen resistance (12, 22, 23). In line with these reports, we found that high nuclear levels of FGFR1 are associated with a lower response to estrogen suppression in ER+ tumors in patients (Fig. 1B; Supplementary Fig. S2A). This correlation between nuclear FGFR1 levels and resistance to estrogen suppression was independent of FGFR1 gene amplification (Supplementary Fig. S2C). We acknowledge that other genomic and nongenomic alterations, not investigated in this clinical cohort may have influenced the results of the correlative analysis reported in Fig. 1B. However, in agreement with our clinical findings, MCF7 cells transduced with an FGFR1 expression vector containing a strong nuclear localization sequence exhibited robust growth in the absence of estrogen and, when established as xenografts in ovariectomized mice, were relative resistant to the ER antagonist fulvestrant (Fig. 1C–F). In these cells, genomic-bound FGFR1 induced the expression of estrogen receptor early and late genes sets (Supplementary Fig. S9A), suggesting that nuclear FGFR1 contributes to endocrine resistance promoting the ERα-associated transcriptional profile. Furthermore, a gene signature score based on the gene promoters bound by FGFR1 using an approach previously used to predict genes near ERα-binding sites in breast cancer (65), was associated with antiestrogen resistance in the METABRIC cohort of patients with ER+ breast cancer treated with antiestrogens (Fig. 2G and H). Further suggesting a clinical significance of our findings, the FGFR1-induced polygenic risk score also correlated with FGFR1 copy number alterations in two large datasets of patients, TCGA and METABRIC (Supplementary Fig. S6B and S6C).

Unlike the findings from Chioni and colleagues (15), nuclear fractions of MCF7 cells transduced with an FGFR1 D432N mutant exhibited higher expression of both full-length (∼140 KDa) and the cleaved form (∼55 KDa) of FGFR1 (Supplementary Fig. S3A). In agreement with these results, neither GZMB knockdown by siRNA, nor treatment with recombinant human Granzyme B (rhGrB), altered nuclear levels of full-length and cleaved FGFR1 in MCF7, CAMA1, and MDA-MB-134 cells (data not shown). These results suggest that Granzyme B does not contribute to FGFR1 cleavage in these cells.

Mapping of FGFR1 genomic occupancy in ER+/FGFR1-amplified breast cancer cells by ChIP-seq revealed marked enrichment at TSS/promoter regions (Fig. 2C and D). Supporting a role of FGFR1 recruitment in active gene transcription, FGFR1 coupled with phosphorylated RNA-Polymerase II (Fig. 3C; Supplementary Fig. S7A–S7C). This hypothesis was further strengthened by the overlap of FGFR1 with histone marks H3K27ac and H3K4me3 (Fig. 4A and C), representative of promoters of actively transcribed genes. These findings are in line with previous studies reporting FGFR1 localization to nuclear speckles, which are sites of active transcription, and colocalization with RNA Pol II in TE671 human medulloblastoma cells (66). To the best of our knowledge, this is the first report describing the genome-wide distribution of FGFR1 in cancer.

The interaction of FGFR1 with CDK7, CDK9, Cyclin H, CTR9, and MED family proteins, revealed by MS (Supplementary Fig. S7D; Supplementary Table S6), further supported the notion that FGFR1 is part of a complex that regulates gene transcription. FGFR1 does not harbor a bona fide DNA-binding domain. Thus, it is likely that FGFR1 is tethered to promoters through interaction with other proteins, presumably general transcription factors. We also showed that FOXA1 plays a significant role in the recruitment of FGFR1 to chromatin. FOXA1 orchestrates transcription factor binding to chromatin and drives gene expression that results in endocrine therapy resistance in cancer (67, 68). FOXA1 pioneering function has been shown to be independent on steroids and hormone signaling (53) but can be affected by mitogenic signals (65). Also, because FOXA1 is currently considered undruggable, understanding how extracellular and intracellular signals modulate FOXA1–FGFR1 interaction may provide rationale to test new therapeutic strategies to inhibit FGFR1 transcriptional function.

Our findings provide new insight into the biology and oncogenic functions of RTKs. For example, nuclear EGFR has been associated with resistance to the anti-EGFR antibody cetuximab (69). EGFR nuclear translocation has been associated with resistance to chemotherapy and ionizing radiation in lung cancer cells, primarily by modulating the DNA repair process (70). Nuclear HER2 was reported to be involved in resistance to the HER2 antibody trastuzumab in HER2+ breast cancer cells (71). More recently, Hancock and colleagues reported an unexpected role of the nuclear insulin receptor (INSR) as a transcription factor that induces gene expression associated with insulin-related functions and insulin resistance (19). Previous studies showed that ligand activation of membrane-bound FGFR1 can induce its migration to the nucleus, mainly by importin β family proteins (14, 54). There is also evidence that Sec61 α family proteins can mediate translocation of newly synthesized FGFR1 from the endoplasmic reticulum to the nucleus (72). A complete understanding of the factors regulating these mechanisms could identify new therapeutic targets to suppress FGFR1 nuclear localization and FGFR1-induced gene expression.

FGFR TKIs have been clinically effective against cancers harboring FGFR alterations acting as dominant oncogenic drivers, such as FGFR3 mutations or fusions in urothelial carcinoma, or FGFR2 fusions in cholangiocarcinoma (73, 74). Even though testing has been more limited in patients with FGFR1-amplified breast cancer, FGFR TKIs have not shown any meaningful clinical activity against these tumors (75, 76). The reasons for this failure may be attributed to various reasons. First, the inhibitors that have been tested are not specific and exhibit non-FGFR–associated toxicities that limit their dosing (1). Second, it is not clear to what extent FGFR1-amplified tumors depend on canonical FGFR signaling. Third, there is not agreement regarding the optimal threshold of copy number to be used to stratify patients enrolled in clinical trials. Results from the study showed herein may shed light on an unconventional role of FGFR1 in these tumors, different from the classical RTK-induced signaling and, as such, have translational implications for drug development. For example, treatment with erdafitinib, a pan-FGFR TKI currently in clinical trials in breast cancer, did not affect nuclear FGFR1 levels and its genomic distribution (Fig. 6D and F), nor growth of cells overexpressing nuclear FGFR1 (Supplementary Fig. S13H), ultimately suggesting that FGFR1 transcriptional activity is not influenced by therapeutic inhibition of the FGFR1 tyrosine kinase. We speculate this could explain, at least in part, the low activity of FGFR TKIs against FGFR1-overexpressing tumors. Taken together, these data provide a rationale for the development of targeted therapies that inhibit nuclear FGFR1 function in ER+ breast cancer.

A. Servetto reports grants from University of Texas Southwestern Simmons Cancer Center P30 CA142543, CPRIT RR170061 (to C.L. Arteaga), NCI Breast SPORE P50 CA098131, Vanderbilt-Ingram Cancer Center P30 CA68485, Susan G. Komen Breast Cancer Foundation SAC100013 (to C.L. Arteaga), Breast Cancer Research Foundation (to C.L. Arteaga), NCI R01CA224899 (to C.L. Arteaga and A.B. Hanker), and Susan G. Komen Postdoctoral Fellowship PDF17487926 (to K.-M. Lee) during the conduct of the study. L. Formisano reports grants from AIRC and Eli Lilly during the conduct of the study. D.R. Sudhan reports grants from Takeda, Puma, Pfizer, and Eli Lilly and other from Provista, Y-TRAP, Novartis, Merck, Daiichi Sankyo, Taiho Oncology, OrigiMed, Immunomedics, and AstraZeneca during the conduct of the study. S. Chatterjee reports grants from UTSW Simmons Cancer Center P30 CA142543, CPRIT RR170061 (to C.L. Arteaga), NCI Breast SPORE P50 CA098131, Vanderbilt-Ingram Cancer Center P30 CA68485, Susan G. Komen Breast Cancer Foundation SAC100013 (to C.L. Arteaga), Breast Cancer Research Foundation (to C.L. Arteaga), NCI R01CA224899 (to C.L. Arteaga and A.B. Hanker), Susan G. Komen Postdoctoral Fellowship PDF17487926 (to K. Lee), Associazione Italiana per la Ricerca sul Cancro (AIRC) Investigator grant IG2018–21339 (to R. Bianco), AIRC My First Grant MFGA2018–21505 (to L. Formisano). A.B. Hanker receives research grant support from Takeda and has received a travel grant from Puma Biotechnology, Inc.; receives research grant support from Pfizer; receives travel expenses support from Pfizer and Novartis; has consultory role for Novartis and AstraZeneca; receives or has received research grants from Pfizer, Lilly, and Takeda; holds minor stock options in Provista and Y-TRAP; serves or has served in an advisory role to Novartis, Merck, Lilly, Daiichi Sankyo, Taiho Oncology, OrigiMed, Puma Biotechnology, Immunomedics, AstraZeneca, and Sanofi; and reports scientific advisory board remuneration from the Susan G. Komen Foundation during the conduct of the study. A. Guerrero-Zotano reports grants, personal fees, and non-financial support from Pfizer; personal fees and non-financial support from Novartis; non-financial support from Roche; and personal fees and non-financial support from AstraZeneca outside the submitted work. S. Mendiratta reports grants from UTSW Simmons Cancer Center, CPRIT, NCI Breast SPORE, Vanderbilt-Ingram Cancer Center, Susan G. Komen Breast Cancer Foundation, Breast Cancer Research Foundation, Susan G. Komen Postdoctoral Fellowship, Associazione Italiana per la Ricerca sul Cancro (AIRC) Investigator Grant, AIRC My First Grant, and NCI during the conduct of the study. R. Bianco reports grants from AIRC outside the submitted work. A.B. Hanker reports grants from NIH during the conduct of the study, as well as grants from Takeda and Breast Cancer Research Foundation and non-financial support from Puma Biotechnology, Daiichi Sankyo, Symphogen, and South Texas Accelerated Research Therapeutics outside the submitted work. C.L. Arteaga reports grants from Lilly, Pfizer, and Takeda and personal fees from Novartis, Lilly, Immunomedics, Merck, Daiichi Sakyo, Taiho Oncology, AstraZeneca, OrigiMed, Arvinas, Clovis, and Athenex outside the submitted work, as well as Minor stock options in Provista and Scientific Advisory Board honorarium from the Susan G. Komen for the Cure Breast Cancer Research Foundation. No disclosures were reported by the other authors.

A. Servetto: Conceptualization, data curation, formal analysis, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. R. Kollipara: Resources, data curation, software, formal analysis, visualization, writing–review and editing. L. Formisano: Conceptualization, funding acquisition. C.-C. Lin: Conceptualization, formal analysis, investigation, methodology. K.-M. Lee: Conceptualization, software, formal analysis, investigation, writing–review and editing. D.R. Sudhan: Conceptualization, formal analysis, visualization, methodology. P.I. Gonzalez-Ericsson: Resources, software, formal analysis, investigation, visualization, methodology. S. Chatterjee: Investigation, methodology. A. Guerrero-Zotano: Formal analysis, methodology. S. Mendiratta: Project administration. H. Akamatsu: Investigation, methodology. N. James: Resources, formal analysis, investigation, visualization. R. Bianco: Supervision, funding acquisition, writing–review and editing. A.B. Hanker: Conceptualization, data curation, supervision, funding acquisition, writing–review and editing. R. Kittler: Conceptualization, resources, supervision, methodology, writing–review and editing. C.L. Arteaga: Conceptualization, resources, data curation, supervision, funding acquisition, writing–review and editing.

We thank the assistance of the UTSW Proteomics Core. We acknowledge Professor Yonghao Yu (Department of Biochemistry at UT Southwestern Medical Center, Dallas, TX) for sharing the pCDNA5-zz-FLAG plasmid. We thank Professor Michal Stachowiak (State University of New York, New York, NY) for sharing the pCDNA3.1(-)FGFR1(SP-)(NLS) plasmid. We thank Professor Alana Welm (University of Utah, Salt Lake City) for sharing the HCI-003 and HCI-011 PDXs. We received the following financial support to conduct this work: UTSW Simmons Cancer Center P30 CA142543, CPRIT RR170061 (to C.L. Arteaga), NCI Breast SPORE P50 CA098131, Vanderbilt-Ingram Cancer Center P30 CA68485, Susan G. Komen Breast Cancer Foundation SAB1800010 (to C.L. Arteaga), Breast Cancer Research Foundation (to C.L. Arteaga), NCI R01CA224899 (to C.L. Arteaga and A.B. Hanker), Susan G. Komen Postdoctoral Fellowship PDF17487926 (to K.-M. Lee), Associazione Italiana per la Ricerca sul Cancro (AIRC) Investigator Grant IG2018–21339 (to R. Bianco), and AIRC My First Grant MFGA2018–21505 (to L. Formisano).

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.

Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).

1.
Turner
N
,
Grose
R
. 
Fibroblast growth factor signalling: from development to cancer
.
Nat Rev Cancer
2010
;
10
:
116
29
.
2.
Schlessinger
J
,
Plotnikov
AN
,
Ibrahimi
OA
,
Eliseenkova
AV
,
Yeh
BK
,
Yayon
A
, et al
Crystal structure of a ternary FGF-FGFR-heparin complex reveals a dual role for heparin in FGFR binding and dimerization
.
Mol Cell
2000
;
6
:
743
50
.
3.
Eswarakumar
VP
,
Lax
I
,
Schlessinger
J
. 
Cellular signaling by fibroblast growth factor receptors
.
Cytokine Growth Factor Rev
2005
;
16
:
139
49
.
4.
Peters
KG
,
Marie
J
,
Wilson
E
,
Ives
HE
,
Escobedo
J
,
Del Rosario
M
, et al
Point mutation of an FGF receptor abolishes phosphatidylinositol turnover and Ca2+ flux but not mitogenesis
.
Nature
1992
;
358
:
678
81
.
5.
Klint
P
. 
Signal transduction by fibroblast growth factor receptors
.
Front Biosci
1999
;
4
:
d165
.
6.
Hart
KC
,
Robertson
SC
,
Kanemitsu
MY
,
Meyer
AN
,
Tynan
JA
,
Donoghue
DJ
. 
Transformation and Stat activation by derivatives of FGFR1, FGFR3, and FGFR4
.
Oncogene
2000
;
19
:
3309
20
.
7.
Reis-Filho
JS
,
Simpson
PT
,
Turner
NC
,
Lambros
MB
,
Jones
C
,
Mackay
A
, et al
FGFR1 emerges as a potential therapeutic target for lobular breast carcinomas
.
Clin Cancer Res
2006
;
12
:
6652
62
.
8.
Yang
W
,
Yao
Y-W
,
Zeng
J-L
,
Liang
W-J
,
Wang
L
,
Bai
C-Q
, et al
Prognostic value of FGFR1 gene copy number in patients with non-small cell lung cancer: a meta-analysis
.
J Thorac Dis
2014
;
6
:
803
9
.
9.
Gorringe
KL
,
Jacobs
S
,
Thompson
ER
,
Sridhar
A
,
Qiu
W
,
Choong
DYH
, et al
High-resolution single nucleotide polymorphism array analysis of epithelial ovarian cancer reveals numerous microdeletions and amplifications
.
Clin Cancer Res
2007
;
13
:
4731
9
.
10.
Simon
R
,
Richter
J
,
Wagner
U
,
Fijan
A
,
Bruderer
J
,
Schmid
U
, et al
High-throughput tissue microarray analysis of 3p25 (RAF1) and 8p12 (FGFR1) copy number alterations in urinary bladder cancer
.
Cancer Res
2001
;
61
:
4514
9
.
11.
Comprehensive molecular portraits of human breast tumours
.
Nature
2012
;
490
:
61
70
.
12.
Giltnane
JM
,
Hutchinson
KE
,
Stricker
TP
,
Formisano
L
,
Young
CD
,
Estrada
MV
, et al
Genomic profiling of ER + breast cancers after short-term estrogen suppression reveals alterations associated with endocrine resistance
.
Sci Transl Med
2017
;
9
:
eaai7993
.
13.
Babina
IS
,
Turner
NC
. 
Advances and challenges in targeting FGFR signalling in cancer
.
Nat Rev Cancer
2017
;
17
:
318
32
.
14.
Coleman
SJ
,
Chioni
A
,
Ghallab
M
,
Anderson
RK
,
Lemoine
NR
,
Kocher
HM
, et al
Nuclear translocation of <scp>FGFR</scp>1 and <scp>FGF</scp>2 in pancreatic stellate cells facilitates pancreatic cancer cell invasion
.
EMBO Mol Med
2014
;
6
:
467
81
.
15.
Chioni
A-M
,
Grose
R
. 
FGFR1 cleavage and nuclear translocation regulates breast cancer cell behavior
.
J Cell Biol
2012
;
197
:
801
17
.
16.
Terranova
C
,
Narla
ST
,
Lee
Y-W
,
Bard
J
,
Parikh
A
,
Stachowiak
EK
, et al
Global developmental gene programing involves a nuclear form of fibroblast growth factor receptor-1 (FGFR1)
.
PLoS One
2015
;
10
:
e0123380
.
17.
Fang
X
,
Stachowiak
EK
,
Dunham-Ems
SM
,
Klejbor
I
,
Stachowiak
MK
. 
Control of CREB-binding protein signaling by nuclear fibroblast growth factor receptor-1
.
J Biol Chem
2005
;
280
:
28451
62
.
18.
Carpenter
G
,
Liao
H-J
. 
Receptor tyrosine kinases in the nucleus
.
Cold Spring Harb Perspect Biol
2013
;
5
:
a008979
.
19.
Hancock
ML
,
Meyer
RC
,
Mistry
M
,
Khetani
RS
,
Wagschal
A
,
Shin
T
, et al
Insulin receptor associates with promoters genome-wide and regulates gene expression
.
Cell
2019
;
177
:
722
36
.
20.
Redmond
AM
,
Omarjee
S
,
Chernukhin
I
,
Le Romancer
M
,
Carroll
JS
. 
Analysis of HER2 genomic binding in breast cancer cells identifies a global role in direct gene regulation
.
PLoS One
2019
;
14
:
e0225180
.
21.
Lin
S-Y
,
Makino
K
,
Xia
W
,
Matin
A
,
Wen
Y
,
Kwong
KY
, et al
Nuclear localization of EGF receptor and its potential new role as a transcription factor
.
Nat Cell Biol
2001
;
3
:
802
8
.
22.
Drago
JZ
,
Formisano
L
,
Juric
D
,
Niemierko
A
,
Servetto
A
,
Wander
SA
, et al
FGFR1 amplification mediates endocrine resistance but retains TORC sensitivity in metastatic hormone receptor–positive (HR+) breast cancer
.
Clin Cancer Res
2019
;
25
:
6443
51
.
23.
Turner
N
,
Pearson
A
,
Sharpe
R
,
Lambros
M
,
Geyer
F
,
Lopez-Garcia
MA
, et al
FGFR1 amplification drives endocrine therapy resistance and is a therapeutic target in breast cancer
.
Cancer Res
2010
;
70
:
2085
94
.
24.
Formisano
L
,
Stauffer
KM
,
Young
CD
,
Bhola
NE
,
Guerrero-Zotano
AL
,
Jansen
VM
, et al
Association of FGFR1 with ERα maintains ligand-independent ER transcription and mediates resistance to estrogen deprivation in ER + breast cancer
.
Clin Cancer Res
2017
;
23
:
6138
50
.
25.
Bankhead
P
,
Loughrey
MB
,
Fernández
JA
,
Dombrowski
Y
,
McArt
DG
,
Dunne
PD
, et al
QuPath: open source software for digital pathology image analysis
.
Sci Rep
2017
;
7
:
16878
.
26.
Camp
RL
,
Chung
GG
,
Rimm
DL
. 
Automated subcellular localization and quantification of protein expression in tissue microarrays
.
Nat Med
2002
;
8
:
1323
8
.
27.
Langmead
B
,
Salzberg
SL
. 
Fast gapped-read alignment with Bowtie 2
.
Nat Methods
2012
;
9
:
357
9
.
28.
Li
H
,
Handsaker
B
,
Wysoker
A
,
Fennell
T
,
Ruan
J
,
Homer
N
, et al
The sequence alignment/map format and SAMtools
.
Bioinformatics
2009
;
25
:
2078
9
.
29.
Zhang
Y
,
Liu
T
,
Meyer
CA
,
Eeckhoute
J
,
Johnson
DS
,
Bernstein
BE
, et al
Model-based analysis of ChIP-Seq (MACS)
.
Genome Biol
2008
;
9
:
R137
.
30.
Heinz
S
,
Benner
C
,
Spann
N
,
Bertolino
E
,
Lin
YC
,
Laslo
P
, et al
Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities
.
Mol Cell
2010
;
38
:
576
89
.
31.
Shen
L
,
Shao
N
,
Liu
X
,
Nestler
E
. 
ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases
.
BMC Genomics
2014
;
15
:
284
.
32.
Vilcaes
AA
,
Chanaday
NL
,
Kavalali
ET
. 
Interneuronal exchange and functional integration of synaptobrevin via extracellular vesicles
.
Neuron
2021
;
109
:
971
83
.
33.
Wolff
AC
,
Hammond
MEH
,
Allison
KH
,
Harvey
BE
,
Mangu
PB
,
Bartlett
JMS
, et al
Human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists Clinical Practice Guideline Focused Update
.
Arch Pathol Lab Med
2018
;
142
:
1364
82
.
34.
Peng
H
,
Moffett
J
,
Myers
J
,
Fang
X
,
Stachowiak
EK
,
Maher
P
, et al
Novel nuclear signaling pathway mediates activation of fibroblast growth factor-2 gene by type 1 and type 2 angiotensin II receptors
.
Mol Biol Cell
2001
;
12
:
449
62
.
35.
Ran
FA
,
Hsu
PD
,
Wright
J
,
Agarwala
V
,
Scott
DA
,
Zhang
F
. 
Genome engineering using the CRISPR-Cas9 system
.
Nat Protoc
2013
;
8
:
2281
308
.
36.
Sandelin
A
,
Carninci
P
,
Lenhard
B
,
Ponjavic
J
,
Hayashizaki
Y
,
Hume
DA
. 
Mammalian RNA polymerase II core promoters: insights from genome-wide studies
.
Nat Rev Genet
2007
;
8
:
424
36
.
37.
Fenouil
R
,
Cauchy
P
,
Koch
F
,
Descostes
N
,
Cabeza
JZ
,
Innocenti
C
, et al
CpG islands and GC content dictate nucleosome depletion in a transcription-independent manner at mammalian promoters
.
Genome Res
2012
;
22
:
2399
408
.
38.
Deaton
AM
,
Bird
A
. 
CpG islands and the regulation of transcription
.
Genes Dev
2011
;
25
:
1010
22
.
39.
Landolin
JM
,
Johnson
DS
,
Trinklein
ND
,
Aldred
SF
,
Medina
C
,
Shulha
H
, et al
Sequence features that drive human promoter function and tissue specificity
.
Genome Res
2010
;
20
:
890
8
.
40.
Hahn
S
. 
Structure and mechanism of the RNA polymerase II transcription machinery
.
Nat Struct Mol Biol
2004
;
11
:
394
403
.
41.
Buratowski
S
. 
Progression through the RNA polymerase II CTD cycle
.
Mol Cell
2009
;
36
:
541
6
.
42.
Komarnitsky
P
. 
Different phosphorylated forms of RNA polymerase II and associated mRNA processing factors during transcription
.
Genes Dev
2000
;
14
:
2452
60
.
43.
Jonkers
I
,
Lis
JT
. 
Getting up to speed with transcription elongation by RNA polymerase II
.
Nat Rev Mol Cell Biol
2015
;
16
:
167
77
.
44.
Zheng
R
,
Wan
C
,
Mei
S
,
Qin
Q
,
Wu
Q
,
Sun
H
, et al
Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis
.
Nucleic Acids Res
2019
;
47
:
D729
35
.
45.
Wang
S
,
Sun
H
,
Ma
J
,
Zang
C
,
Wang
C
,
Wang
J
, et al
Target analysis by integration of transcriptome and ChIP-seq data with BETA
.
Nat Protoc
2013
;
8
:
2502
15
.
46.
Heintzman
ND
,
Stuart
RK
,
Hon
G
,
Fu
Y
,
Ching
CW
,
Hawkins
RD
, et al
Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome
.
Nat Genet
2007
;
39
:
311
8
.
47.
Kundaje
A
,
Meuleman
W
,
Ernst
J
,
Bilenky
M
,
Yen
A
,
Heravi-Moussavi
A
, et al
Integrative analysis of 111 reference human epigenomes
.
Nature
2015
;
518
:
317
30
.
48.
Cirillo
LA
,
Lin
FR
,
Cuesta
I
,
Friedman
D
,
Jarnik
M
,
Zaret
KS
. 
Opening of compacted chromatin by early developmental transcription factors HNF3 (FoxA) and GATA-4
.
Mol Cell
2002
;
9
:
279
89
.
49.
Fu
X
,
Pereira
R
,
De Angelis
C
,
Veeraraghavan
J
,
Nanda
S
,
Qin
L
, et al
FOXA1 upregulation promotes enhancer and transcriptional reprogramming in endocrine-resistant breast cancer
.
Proc Natl Acad Sci U S A
2019
;
116
:
26823
34
.
50.
Hurtado
A
,
Holmes
KA
,
Ross-Innes
CS
,
Schmidt
D
,
Carroll
JS
. 
FOXA1 is a key determinant of estrogen receptor function and endocrine response
.
Nat Genet
2011
;
43
:
27
33
.
51.
Robinson
JLL
,
Hickey
TE
,
Warren
AY
,
Vowler
SL
,
Carroll
T
,
Lamb
AD
, et al
Elevated levels of FOXA1 facilitate androgen receptor chromatin binding resulting in a CRPC-like phenotype
.
Oncogene
2014
;
33
:
5666
74
.
52.
Tsherniak
A
,
Vazquez
F
,
Montgomery
PG
,
Weir
BA
,
Kryukov
G
,
Cowley
GS
, et al
Defining a cancer dependency map
.
Cell
2017
;
170
:
564
76
.
53.
Glont
S-E
,
Chernukhin
I
,
Carroll
JS
. 
Comprehensive genomic analysis reveals that the pioneering function of FOXA1 is independent of hormonal signaling
.
Cell Rep
2019
;
26
:
2558
65
.
54.
Reilly
JF
,
Maher
PA
. 
Importin β–mediated nuclear import of fibroblast growth factor receptor
.
J Cell Biol
2001
;
152
:
1307
12
.
55.
McLean
CY
,
Bristor
D
,
Hiller
M
,
Clarke
SL
,
Schaar
BT
,
Lowe
CB
, et al
GREAT improves functional interpretation of cis-regulatory regions
.
Nat Biotechnol
2010
;
28
:
495
501
.
56.
Perera
TPS
,
Jovcheva
E
,
Mevellec
L
,
Vialard
J
,
De Lange
D
,
Verhulst
T
, et al
Discovery and pharmacological characterization of JNJ-42756493 (Erdafitinib), a functionally selective small-molecule FGFR family inhibitor
.
Mol Cancer Ther
2017
;
16
:
1010
20
.
57.
Lan
T
,
Morgan
DA
,
Rahmouni
K
,
Sonoda
J
,
Fu
X
,
Burgess
SC
, et al
FGF19, FGF21, and an FGFR1/β-Klotho-activating antibody act on the nervous system to regulate body weight and glycemia
.
Cell Metab
2017
;
26
:
709
18
.
58.
Oladipupo
SS
,
Smith
C
,
Santeford
A
,
Park
C
,
Sene
A
,
Wiley
LA
, et al
Endothelial cell FGF signaling is required for injury response but not for vascular homeostasis
.
Proc Natl Acad Sci U S A
2014
;
111
:
13379
84
.
59.
Ono
K
,
Kita
T
,
Sato
S
,
O'Neill
P
,
Mak
S-S
,
Paschaki
M
, et al
FGFR1-Frs2/3 signalling maintains sensory progenitors during inner ear hair cell formation
.
PLoS Genet
2014
;
10
:
e1004118
.
60.
Zhao
M
,
Ross
JT
,
Itkin
T
,
Perry
JM
,
Venkatraman
A
,
Haug
JS
, et al
FGF signaling facilitates postinjury recovery of mouse hematopoietic system
.
Blood
2012
;
120
:
1831
42
.
61.
Goyal
L
,
Saha
SK
,
Liu
LY
,
Siravegna
G
,
Leshchiner
I
,
Ahronian
LG
, et al
Polyclonal secondary FGFR2 mutations drive acquired resistance to FGFR inhibition in patients with FGFR2 fusion–positive cholangiocarcinoma
.
Cancer Discov
2017
;
7
:
252
63
.
62.
Pal
SK
,
Rosenberg
JE
,
Hoffman-Censits
JH
,
Berger
R
,
Quinn
DI
,
Galsky
MD
, et al
Efficacy of BGJ398, a fibroblast growth factor receptor 1–3 inhibitor, in patients with previously treated advanced urothelial carcinoma with FGFR3 alterations
.
Cancer Discov
2018
;
8
:
812
21
.
63.
Razavi
P
,
Chang
MT
,
Xu
G
,
Bandlamudi
C
,
Ross
DS
,
Vasan
N
, et al
The genomic landscape of endocrine-resistant advanced breast cancers
.
Cancer Cell
2018
;
34
:
427
38
.
64.
Lee
JE
,
Shin
S-H
,
Shin
H-W
,
Chun
Y-S
,
Park
J-W
. 
Nuclear FGFR2 negatively regulates hypoxia-induced cell invasion in prostate cancer by interacting with HIF-1 and HIF-2
.
Sci Rep
2019
;
9
:
3480
.
65.
Ross-Innes
CS
,
Stark
R
,
Teschendorff
AE
,
Holmes
KA
,
Ali
HR
,
Dunning
MJ
, et al
Differential oestrogen receptor binding is associated with clinical outcome in breast cancer
.
Nature
2012
;
481
:
389
93
.
66.
Somanathan
S
,
Stachowiak
EK
,
Siegel
AJ
,
Stachowiak
MK
,
Berezney
R
. 
Nuclear matrix bound fibroblast growth factor receptor is associated with splicing factor rich and transcriptionally active nuclear speckles
.
J Cell Biochem
2003
;
90
:
856
69
.
67.
Fu
X
,
Jeselsohn
R
,
Pereira
R
,
Hollingsworth
EF
,
Creighton
CJ
,
Li
F
, et al
FOXA1 overexpression mediates endocrine resistance by altering the ER transcriptome and IL-8 expression in ER-positive breast cancer
.
Proc Natl Acad Sci U S A
2016
;
113
:
E6600
9
.
68.
Lupien
M
,
Eeckhoute
J
,
Meyer
CA
,
Wang
Q
,
Zhang
Y
,
Li
W
, et al
FoxA1 translates epigenetic signatures into enhancer-driven lineage-specific transcription
.
Cell
2008
;
132
:
958
70
.
69.
Li
C
,
Iida
M
,
Dunn
EF
,
Ghia
AJ
,
Wheeler
DL
. 
Nuclear EGFR contributes to acquired resistance to cetuximab
.
Oncogene
2009
;
28
:
3801
13
.
70.
Liccardi
G
,
Hartley
JA
,
Hochhauser
D
. 
EGFR nuclear translocation modulates DNA repair following cisplatin and ionizing radiation treatment
.
Cancer Res
2011
;
71
:
1103
14
.
71.
Cordo Russo
RI
,
Béguelin
W
,
Díaz Flaqué
MC
,
Proietti
CJ
,
Venturutti
L
,
Galigniana
N
, et al
Targeting ErbB-2 nuclear localization and function inhibits breast cancer growth and overcomes trastuzumab resistance
.
Oncogene
2015
;
34
:
3413
28
.
72.
Stachowiak
MK
,
Fang
X
,
Myers
JM
,
Dunham
SM
,
Berezney
R
,
Maher
PA
, et al
Integrative nuclear FGFR1 signaling (INFS) as a part of a universal ?feed-forward-and-gate? signaling module that controls cell growth and differentiation
.
J Cell Biochem
2003
;
90
:
662
91
.
73.
Loriot
Y
,
Necchi
A
,
Park
SH
,
Garcia-Donas
J
,
Huddart
R
,
Burgess
E
, et al
Erdafitinib in locally advanced or metastatic urothelial carcinoma
.
N Engl J Med
2019
;
381
:
338
48
.
74.
Goyal
L
,
Shi
L
,
Liu
LY
,
Fece de la Cruz
F
,
Lennerz
JK
,
Raghavan
S
, et al
TAS-120 overcomes resistance to ATP-competitive FGFR inhibitors in patients with FGFR2 fusion–positive intrahepatic cholangiocarcinoma
.
Cancer Discov
2019
;
9
:
1064
79
.
75.
André
F
,
Bachelot
T
,
Campone
M
,
Dalenc
F
,
Perez-Garcia
JM
,
Hurvitz
SA
, et al
Targeting FGFR with dovitinib (TKI258): preclinical and clinical data in breast cancer
.
Clin Cancer Res
2013
;
19
:
3693
702
.
76.
Nogova
L
,
Sequist L
V
,
Perez Garcia
JM
,
Andre
F
,
Delord
J-P
,
Hidalgo
M
, et al
Evaluation of BGJ398, a fibroblast growth factor receptor 1–3 kinase inhibitor, in patients with advanced solid tumors harboring genetic alterations in fibroblast growth factor receptors: results of a global phase I, dose-escalation and dose-expansion Stu
.
J Clin Oncol
2017
;
35
:
157
65
.
77.
Hänzelmann
S
,
Castelo
R
,
Guinney
J
. 
GSVA: gene set variation analysis for microarray and RNA-Seq data
.
BMC Bioinformatics
2013
;
14
:
7
.
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs International 4.0 License.

Supplementary data