Abstract
The expression of several members of the FOX gene family is known to be altered in a variety of cancers. We show that in breast cancer, FOXF1 gene is a target of epigenetic inactivation and that its gene product exhibits tumor-suppressive properties. Loss or downregulation of FOXF1 expression is associated with FOXF1 promoter hypermethylation in breast cancer cell lines and in invasive ductal carcinomas. Methylation of FOXF1 in invasive ductal carcinoma (37.6% of 117 cases) correlated with high tumor grade. Pharmacologic unmasking of epigenetic silencing in breast cancer cells restored FOXF1 expression. Re-expression of FOXF1 in breast cancer cells with epigenetically silenced FOXF1 genes led to G1 arrest concurrent with or without apoptosis to suppress both in vitro cell growth and in vivo tumor formation. FOXF1-induced G1 arrest resulted from a blockage at G1-S transition of the cell cycle through inhibition of the CDK2-RB-E2F cascade. Small interfering RNA–mediated depletion of FOXF1 in breast cancer cells led to increased DNA re-replication, suggesting that FOXF1 is required for maintaining the stringency of DNA replication and genomic stability. Furthermore, expression profiling of cell cycle regulatory genes showed that abrogation of FOXF1 function resulted in increased expression of E2F-induced genes involved in promoting the progression of S and G2 phases. Therefore, our studies have identified FOXF1 as a potential tumor suppressor gene that is epigenetically silenced in breast cancer, which plays an essential role in regulating cell cycle progression to maintain genomic stability. Cancer Res; 70(14); 6047–58. ©2010 AACR.
Introduction
The members of FOX gene family encode transcription factors that share a common, structurally related, DNA-binding domain called the forkhead domain (1, 2). Many members of this family have been found to play important roles in embryonic development as well as in the control of cell cycle, cell growth, cellular metabolism, life span, and immune responses (3, 4). The genetic alteration and deregulation of members of several FOX subfamilies has been linked to the progression of certain types of cancers (5–8). A comprehensive understanding of expression profiles, genetic alterations, and epigenetic changes of FOX family genes as well as their binding proteins and target genes is important to develop novel therapeutics and preventives for human diseases, especially for cancer.
The forkhead box-F1 (FOXF1) transcription factor plays a key role in regulating embryonic development. Foxf1−/− embryos exhibited defects in extraembryonic mesoderm development and died by day 8 postcoitum (9). On account of the lethal consequences of a complete loss of Foxf1 alleles, Foxf1+/− mice have been studied extensively for deciphering Foxf1 function. Haploinsufficiency in Foxf1+/− mice resulted in various abnormalities in the mesenchyme of lung, liver, gallbladder, esophagus, and trachea (10–12). In humans, genome-wide array analysis of prostate cancer showed that FOXF1 mRNA expression was downregulated in tumors compared with normal prostate, providing a clue that FOXF1 might behave as a tumor suppressor gene (13).
In this study, we report that FOXF1 expression is aberrantly silenced in breast cancer through epigenetic mechanisms, and show that FOXF1 exerts tumor suppressor activity. Restoration of FOXF1 expression in breast cancer cells suppresses cell growth and tumorigenicity, whereas silencing of FOXF1 expression in FOXF1-expressing breast cancer cells by small interfering RNA (siRNA) results in an increase in DNA re-replication, a sign of a defect in the stringent control of DNA replication. Therefore, we have identified FOXF1 as a putative tumor suppressor gene that is frequently epigenetically inactivated in breast cancer.
Materials and Methods
In silico analysis of FOXF1 gene expression
The Oncomine Cancer Microarray Database (14) was used to perform in silico analysis of the FOXF1 gene expression in cancer versus the normal counterpart tissue.
Identification of the CpG island of FOXF1 gene
We obtained the genomic DNA sequence, including the upstream promoter sequence, of FOX1 gene from the GenBank Database of the National Center of Biotechnology Information (15). We identified the CpG-rich region at the FOXF1 gene locus by using the criteria and algorithm of online CpG Island Searcher (16).
Cell lines and tissue samples
Immortalized, nontumorigenic human mammary epithelial cells (HBL100) and 10 breast cancer cell lines (obtained from the American Type Culture Collection) were cultured according to the online instructions of the manufacturer. Mammary organoids were prepared from reduction mammoplasty specimens of normal women as previously described (17). We obtained fresh reduction mammoplasty specimens and frozen primary breast cancer samples from the Department of Pathology, Johns Hopkins Hospital (Baltimore, MD). All of the human tissue specimens in this study were processed and used with prior approval from the Johns Hopkins Institutional Review Board.
Methylation-specific PCR and bisulfite sequencing analysis
Bisulfite modification of genomic DNA and methylation-specific PCR (MSP) amplification were performed as previously described (18, 19). The sequence information for MSP primers is provided in the Supplemental Experimental Procedures. Bisulfite sequencing analysis of the FOXF1 promoter was performed as previously described (19) by using a pair of external primers (primer sequences are available on request).
Conventional and quantitative reverse-transcription PCR analysis
Detailed methods are provided in the Supplemental Experimental Procedures.
Treatment with 5-aza-2′-deoxycytidine and/or trichostatin A
Breast cancer cell lines were treated with the regimen as previously described (19).
Immunohistochemistry assay
The immunohistochemical analysis was performed using the avidin-biotin-peroxidase complex method with the anti-FOXF1 antibody (Aviva Antibody Corporation). See the Supplemental Experimental Procedures for details.
Statistical analysis
Comparison of clinicopathologic characteristics of patients with invasive ductal carcinoma (IDC) with methylation status of FOXF1 gene, and the association between FOXF1 methylation and FOXF1 protein expression was subjected to χ2 analysis using the SPSS system (version 11.5 for Windows; SPSS, Inc.). Student's t test was used to analyze the significance of the difference between two groups of data. P < 0.05 was regarded as statistically significant.
Construction of the FOXF1 expression plasmid
Full-length cDNA fragments for FOXF1 gene were obtained by reverse-transcription PCR (RT-PCR) of RNA from normal mammary organoids and cloned into the pGEM-T-Easy vector (Promega). Sequence-confirmed FOXF1 cDNA was subcloned into the pcDNA3.1(−) vector (Invitrogen) with the hemagglutinin (HA)-tag coding sequence.
Transfection, Western blotting analysis, and reporter assay
Details are provided in the Supplemental Experimental Procedures.
Colony formation assay
At 24 hours after transfection, transfected cells were plated in six-well plates at a density of 6,000 to 12,000 cells/well in complete medium containing 0.3 to 0.6 mg/mL of geneticin (G418). After 2 to 3 weeks, colonies were fixed and then stained with 0.1% crystal violet. Unbound dye was washed out with water. Bound dye was then eluted with 10% acetic acid and quantitated by measuring the absorbance at 590 nm.
Flow cytometry analysis
Detailed procedures are provided in the Supplemental Experimental Procedures.
Tumorigenicity assay
G418-selected viable HA-FOXF1–transfected Hs578T cells (5 × 106) were injected s.c. into the right rear flanks of six athymic nude mice (The Jackson Laboratory), and empty vector–transfected Hs578T cells were injected into the left rear flanks of the same mice as a control. Weekly measurements of the tumors were used to determine the tumor volume.
Immunofluorescence assay
Transfected cells on coverslips were first fixed and then were sequentially queried with primary (including anti-HA, anti-CDK2, anti-phosphorylated CDK2, anti-RB, anti-phosphorylated RB, and anti-FOXF1 antibodies) and secondary antibodies (Alexa Series, Invitrogen) according to the method detailed in the Supplemental Experimental Procedures.
siRNA transfection
siRNA transfections were performed with 20 nmol/L of each siRNA (obtained from Dharmacon) using OligofectAMINE RNAiMAX (Invitrogen) according to the instructions of the manufacturer. The FOXF1 siRNA sequence is 5′-GAAAGGAGUUUGUCUUCUC-3′.
Bromodeoxyuridine incorporation assay
Cells were bromodeoxyuridine (BrdUrd)-labeled, fixed, denatured, and detected using the anti-BrdUrd antibody as described in the Supplemental Experimental Procedures.
Gene expression profiling of cell cycle regulatory genes
Gene expression profiling of cell cycle regulatory genes was performed using the RT Profiler PCR Array System (SABiosciences). See the Supplemental Experimental Procedures for details.
Results
Epigenetic silencing of FOXF1 gene expression in breast cancer
In silico meta-analysis of the Oncomine Cancer Microarray Database (14) was performed to examine whether FOXF1 gene expression is altered in cancer. The results showed that FOXF1 was underexpressed in various cancer types such as lung, prostate, bladder, ovarian, and breast cancers compared with their respective normal tissues (Supplementary Fig. S1; Supplementary Table S1; refs. 20–32). To validate these findings, we examined FOXF1 expression in a panel of breast cancer cell lines by quantitative RT-PCR (qRT-PCR). An immortalized human mammary epithelial cell line (HBL100) and three epithelial-enriched organoid samples from breast tissues were used as normal controls. Significant loss or downregulation of FOXF1 mRNA expression was observed in 9 out of 10 breast cancer cell lines compared with normal controls (Fig. 1A).
Expression and methylation status of the FOXF1 gene in breast cancer cell lines. A, qRT-PCR analysis of FOXF1 expression in three mammary epithelial-enriched organoids, an immortalized HMEC line (HBL100), and 10 malignant breast cancer cell lines. B, MSP analysis of FOXF1 gene in breast cancer cell lines. A map of CpG dinucleotides in exon 1 of the FOXF1 gene and the genomic region upstream of exon 1 (top). The exon 1, MSP-amplified, and bisulfite sequencing regions are indicated. The MSP results: U, unmethylated; M, methylated (bottom). C, bisulfite sequencing analysis of FOXF1 promoter in normal breast epithelial and breast cancer cell lines. The amplified FOXF1 DNA region for sequencing analysis ranges from −209 to +62 nucleotides relative to the first nucleotide (set as the default +1) of exon 1 based on FOXF1 cDNA sequence (NM_001451). The methylation status of each CpG site is indicated by the open circle (unmethylated), partially filled circle (partially methylated), and completely filled circle (completely methylated). D, RT-PCR analysis of FOXF1 expression in breast cancer cell lines treated with 5-aza-2′-deoxycytidine (5-azaC), TSA alone, or a combination (top). SKBR3 cells were treated with TSA at two different doses as indicated for 24 and 48 h (bottom). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression is shown as an internal control.
As a first step to determine whether epigenetic mechanisms are involved in silencing FOXF1, we performed MSP. MSP detected hypermethylation of the CpG island in the FOXF1 promoter in 8 out of 10 breast cancer cell lines (Fig. 1B). MSP findings were further validated by bisulfite sequencing analysis. The FOXF1 promoter was fully methylated in MDA-MB-231 and MDA-MB-468, partially methylated in MDA-MB-435, MCF7, and T47D, and unmethylated in normal controls and SKBR3 (Fig. 1C). With the exception of SKBR3, hypermethylation of FOXF1 showed an inverse correlation with mRNA expression in the same sample.
To further establish that epigenetic mechanisms play a role in downregulating FOXF1 expression, five hypermethylated cell lines were treated with epigenetic modifiers. Treatment of these lines with the demethylating agent, 5-aza-2′-deoxycytidine, or histone deacetylase inhibitor, trichostatin A (TSA), singly or in combination resulted in re-expression of FOXF1 mRNA, an effect that was enhanced in some cell lines (MDA-MB-231, MDA-MB-468, and T47D) by cotreatment with TSA (Fig. 1D). FOXF1 mRNA transcripts were also re-expressed in T47D cells by a single treatment with TSA (Fig. 1D). Thus, promoter hypermethylation is the predominant mechanism of silencing gene expression in the breast cancer cell lines. TSA treatment revived FOXF1 expression in SKBR3 cells (Fig. 1D), indicating that in this line, histone deacetylation, not DNA methylation, was responsible for silencing of FOXF1 expression.
To examine whether FOXF1 promoter hypermethylation also occurs in primary breast carcinomas, MSP analysis was performed on IDCs (n = 117), normal mammoplasty tissues (Mam; n = 8), organoids from normal reduction mammoplasty specimens (Org; n = 6), and peripheral WBC (n = 9) from normal individuals. The FOXF1 promoter was hypermethylated in 37.6% (44 of 117) of IDC tumors (Table 1), unmethylated in all normal breast tissue samples and organoids, and weakly methylated in two out of nine WBC samples (Fig. 2A). Bisulfite sequencing analysis have further confirmed that IDC cases (IDC-4, IDC-6, IDC-8, IDC-11, and IDC-14) with positive MSP results had extensive CpG methylation throughout the analyzed genomic region compared with normal organoid and WBC DNA samples with an unmethylated CpG island (Fig. 2B). These five IDC samples consistently exhibited downregulation or loss of FOXF1 expression compared with the normal organoid (Fig. 2B). Hypermethylation of FOXF1 gene correlated with tumor grade (P = 0.004; Table 1). Thus, hypermethylation of FOXF1 commonly occurs in primary breast tumors.
Correlation between the promoter hypermethylation of FOXF1 gene and clinicopathologic characteristics of IDC tumors
Epigenetic inactivation of FOXF1 in primary invasive breast tumors. A, MSP analysis was performed on normal mammoplasty tissues (Mam-1 to Mam-8), normal organoids (Org-1 to Org-6), normal peripheral WBC (WBC-1 to WBC-9), and primary invasive ductal carcinomas (IDC-1 to IDC-16). B, bisulfite sequencing analysis of the FOXF1 promoter in invasive ductal carcinomas. Right, FOXF1 expression status in these primary IDC tumors: (↓), underexpression; (+), expression; (−), no expression; nd, not determined. C, qRT-PCR analysis of FOXF1 expression in normal breast organoids versus primary breast IDC tumors with methylated FOXF1 promoter. D, immunohistochemistry analysis of FOXF1 protein expression in normal breast tissue and primary breast IDC tumors with or without methylated FOXF1 gene. Bottom, a summarized result from 26 IDC samples.
To further explore the functional consequence of FOXF1 promoter hypermethylation, FOXF1 expression was examined in primary breast tumors. qRT-PCR analysis of 6 normal breast organoids (Org-1 to Org-6) and 10 IDC tumor tissues with methylated FOXF1 gene showed that FOXF1 mRNA was significantly underexpressed in IDC tumors (Fig. 2C) compared with normal organoids (P < 0.05). To further confirm the RT-PCR data, FOXF1 protein expression was examined by immunohistochemistry assay in paraffin-embedded tissue sections of 26 tumors (Table 1) selected from the 117 IDC cases tested by MSP. Normal breast tissue exhibited positive immunostaining for FOXF1 protein in ductal epithelial cells (Fig. 2D). Of IDC tumors harboring methylated FOXF1 gene, 12 of 18 (66.7%) cases exhibited negative or focal, weak immunoreactivity for FOXF1 protein (<10% of total tumor cells; Fig. 2D). Among the remaining cases, 6 of 18 (33.3%) showed heterogeneous positive immunostaining and were interpreted as immunohistochemistry-positive. In contrast, in the unmethylated group, six of eight (75%) were heterogeneous or homogeneous immunohistochemistry-positive for FOXF1 (Fig. 2D). The inverse correlation between immunohistochemistry reactivity for FOXF1 protein and methylation status for FOXF1 gene was statistically significant (P < 0.05).
Collectively, evidence from promoter methylation, mRNA, and protein expression studies strongly suggest that FOXF1 promoter hypermethylation correlates significantly with reduction or loss of FOXF1 expression.
FOXF1 re-expression in cancer cells inhibits cell proliferation in vitro and tumors in vivo
We elucidated the biological significance of epigenetic silencing of FOXF1 in breast cancer by performing a variety of assays in cells with enforced expression of HA-tagged FOXF1. Protein expression and nuclear localization of the CMV-HA–tagged FOXF1 construct was confirmed by Western blot and immunofluorescence analyses (data not shown). To test whether exogenously introduced FOXF1 functions as a transcription factor, reporter assays were performed using the reporter plasmid containing a minimal promoter with (p4XFREAC-Luc) or without (pApo-Luc) four tandem copies of the FOXF1 binding site (33). HA-FOXF1 protein induced reporter luciferase activity dramatically in p4XFREAC-Luc–transfected T47D and SKBR3 cells, but not in control pApo-Luc–transfected cells (Fig. 3A).
FOXF1 re-expression suppresses the growth and tumorigenicity of breast cancer cells by inducing G1 arrest and/or apoptosis. A, the transcriptional activity of HA-tagged FOXF1 in T47D and SKBR3 cells was measured by using the FOXF-responsive reporter plasmid (p4XFREAC-Luc) versus its empty reporter control pApo-Luc. Columns, mean from three independent experiments; bars, SD (*, P < 0.05, statistical significance). B, FOXF1 inhibits in vitro growth of multiple breast cancer cell lines. G418-resistant colonies were visualized by staining with crystal violet dye (left). Absorbance values at 590 nm of eluted binding dye are shown as the bar graph (right) based on three independent experiments. C, FOXF1 suppresses tumorigenicity of Hs578T cells. Points, mean measurements of tumor volume; bars, SD (n = 6; bottom). D, flow cytometry analysis of cell cycle profile and apoptotic status in FOXF1-transfected MCF7 and SKBR3 cells. Empty control vector or the HA-tagged FOXF1 expression plasmid DNA was cotransfected with the green fluorescent protein (GFP) expression plasmid DNA into cells. Columns, mean percentages of GFP-positive cells in sub-G1 (indicating apoptosis), G1, S, and G2-M of three independent experiments; bars, SD (*, P < 0.05, statistically significant data).
To assess the effect of FOXF1 on the growth of breast cancer cells, we transfected HA-tagged FOXF1-expressing plasmids into four FOXF1-negative breast cancer cell lines, including MCF7, Hs578T, SKBR3, and T47D cell lines. FOXF1 re-expression strongly suppressed the growth of all these cell lines examined (Fig. 3B). Next, the injection of FOXF1-expressing Hs578T cells s.c. into athymic nude mice showed that tumor formation was significantly suppressed or completely blocked compared to mice injected with vector-transfected Hs578T cells (Fig. 3C). These studies, taken together, strongly suggest that FOXF1 transcription factor possesses tumor suppressor characteristics.
To begin to decipher the mechanism whereby FOXF1 suppresses cell growth, cell cycle analysis was performed. The results showed a significant increase in the percentage of FOXF1-transfected MCF7 cells in G1 phase, indicating that FOXF1 induced G1 arrest (Fig. 3D). In addition to G1 arrest, apoptosis was concurrently observed in SKBR3 cells (Fig. 3D).
FOXF1 abrogates activation of the CDK2-RB-E2F cascade
To understand FOXF1-mediated G1 arrest, protein molecules involved in G1-S transition were examined. Phosphorylation of CDK2 at Thr160 is known to activate CDK2 kinase activity, which in turn promotes G1-S transition (34). Immunofluorescence analysis showed that FOXF1 re-expression in MCF7 diminished the phosphorylation of CDK2 at Thr160 (Fig. 4A) and concurrently abrogated phosphorylation of RB, a downstream target of CDK2 at Ser807/811 (Fig. 4B). FOXF1 had no significant effect on the total levels of CDK2 protein (Fig. 4A), indicating that FOXF1-mediated abrogation of CDK2 activity was not caused by downregulating protein levels of CDK2. In the immunofluorescence analysis of total RB, a moderate decrease in total levels of RB protein was observed in FOXF1-transfected cells compared with nontransfected cells (Fig. 4B, bottom), but even so, this effect could not completely account for dramatically diminishing RB phosphorylation by FOXF1 (Fig. 4B, top). Further confirmation of immunofluorescence data was provided by flow cytometric analysis through quantitative assessment of the overall staining intensities of phosphorylated CDK2 and phosphorylated RB in FOXF1-transfected MCF7 cells compared with vector-transfected control cells (Fig. 4C). These data, taken together, indicate that FOXF1 suppresses CDK2 activation and in turn inhibits CDK2-mediated phosphorylation of RB protein.
FOXF1 re-expression in breast cancer cells inhibits the CDK2-RB-E2F cascade pathway. A, immunofluorescence analysis of MCF7 cells, transfected with the GFP or HA-FOXF1 expression plasmid DNA, using anti-HA (green) and anti-phosphorylated CDK2 (Thr160; red) antibodies (top). DNA was labeled with 4′,6-diamidino-2-phenylindole (DAPI; blue). Immunofluorescence of total CDK2 protein levels (bottom). B, immunofluorescence analysis of the same cells as in A with anti-phosphorylated RB (Ser807/811; top) and anti-RB (bottom) antibodies. C, flow cytometric analysis of the effect of FOXF1 on the levels of phosphorylated CDK2 and phosphorylated RB proteins. FACS histogram data (top). Columns, mean from three independent FACS experiments; bars, SD (*, P < 0.05, statistically significant difference; bottom). D, FOXF1 attenuates E2F1 transcriptional activity. We transfected MCF7 cells with E2F-responsive luciferase reporter DNA and different combinations of expression plasmid DNA as indicated. The luciferase data were converted to fold activation relative to the empty vector control (set as 1). Columns, mean of three independent experiments; bars, SD (*, P < 0.05).
The E2F transcription factors are downstream targets of RB action. We therefore investigated the effect of FOXF1 on E2F1 transcriptional activity by the reporter assay using a promoter reporter plasmid containing E2F-responsive DNA polymerase-α promoter upstream of the luciferase reporter gene (35). Wild-type E2F1, but not mutant E2F1 (E132; ref. 36), induced an upregulation (almost 18-fold) of the promoter activity compared with the vector control (Fig. 4D). Cotransfection of FOXF1 with wild-type E2F1 significantly suppressed E2F1-mediated transactivation (67%) of the reporter activity (Fig. 4D). Suppression was E2F1-specific because FOXF1, in the absence of E2F1, had no significant effect on promoter-linked reporter activity (Fig. 4D). Immunofluorescence data showing a consistent decrease in endogenous total RB protein levels in FOXF1-transfected cells (Fig. 4B) supported the in vitro reporter data because RB is one of the well-known E2F downstream target genes. These results, taken together, strongly suggest that FOXF1-induced inhibition of the CDK2-RB-E2F cascade contributes to a blockage of G1-S transition of the cell cycle.
FOXF1 knockdown leads to an increase in DNA re-replication and elevates the expression of E2F-induced genes
Based on the observed negative role for FOXF1 in G1-S progression, we investigated the effects of siRNA-mediated knockdown of FOXF1 on cell cycle progression of BT549 breast cancer cells. BrdUrd incorporation analysis was used to examine the effect of FOXF1 depletion on DNA replication, cell cycle progression, and apoptosis. As shown in Fig. 5A, FOXF1 siRNA depleted >80% of endogenous FOXF1 mRNA levels. Flow cytometric analysis of BrdUrd-labeled, propidium iodide–stained cells revealed that FOXF1 knockdown resulted in a significant increase in the fraction of BT549 cells with greater than 4N of BrdUrd-incorporated DNA (6.7 ± 1.1% of control cells versus 17.3 ± 2.8% of FOXF1-depleted cells, P < 0.05; Fig. 5B). The polyploid cells formed upon FOXF1 depletion could still be labeled with BrdUrd, indicating that these cells were alive and actively replicating. In fact, the appearance of more BrdUrd-labeled polyploid cells indicates an increase in cells undergoing DNA re-replication (overreplication), which occurs when the initiation of DNA replication fires more than once in each cell cycle (37). In addition to increased DNA re-replication in BT549 cells depleted of FOXF1, a 3-fold induction of apoptosis was observed (Fig. 5B). The increase in apoptosis could be attributed to the death of cells that had undergone DNA re-replication. These findings suggest that inactivation of FOXF1 results in a defect in the stringent control of DNA replication initiation, which in turn, promotes more uncontrolled DNA replication in BT549 cells.
FOXF1 knockdown in cells leads to an increase in DNA re-replication, apoptosis, and expression of E2F-induced genes. A, analysis of the knockdown efficiency of FOXF1 siRNA. Gel-based RT-PCR (top) and real-time qRT-PCR (bottom) results are shown. B, BrdU incorporation analysis of BT549 cells depleted of FOXF1. FACS two-dimensional dot plot data from BrdU incorporation experiments (top). The percentages of DNA re-replicated (right-gated cells) and apoptotic (left-gated cells) cells are shown in two-dimensional dot plots. Columns, mean from three independent experiments; bars, SD (bottom; *, P < 0.05). C, expression profiling of cell cycle genes in BT549 cells depleted of FOXF1. The common logarithms of gene expression values from control siRNA-transfected cells were plotted against those from FOXF1 siRNA-transfected cells to make the scatter plot. The upregulated (≥2-fold) and downregulated (≤2-fold) genes in FOXF1 siRNA-transfected BT549 cells are indicated by the black and gray colors, respectively. D, schematic summary of genes regulated by FOXF1.
To gain further insight into the effects of FOXF1, expression profiling of 84 cell cycle genes was performed. Of 84 genes, 9 were downregulated and 7 were upregulated by FOXF1 in BT549 cells (Fig. 5C and D; Supplementary Table S2). Surprisingly, six out of nine FOXF1-downregulated genes were E2F-induced genes (38, 39), including MCM3, RB1, CCNB1, CCNB2, KPNA2, and BIRC5 (Fig. 5D). Among these genes, MCM3 encodes a protein participating in the initiation of DNA replication, suggesting that FOXF1 might negatively regulate its expression to suppress DNA replication initiation. In addition, CDC34, required for ubiquitin-mediated degradation of cell cycle G1 regulators and for the initiation of DNA replication (40), is also a FOXF1-downregulated gene (Fig. 5C and D). Therefore, deregulation of MCM3 and CDC34 expression might contribute to DNA overreplication after FOXF1 function was silenced (Fig. 5B). Intriguingly, the expression of two genes encoding CDK5-associated factors (CDK5R1 and CDK5RAP1) that are implicated in the activation of CDK5 (41, 42) was elevated in FOXF1-depleted cells (Fig. 5C), implying that FOXF1 might abrogate CDK5 activation which is required for the survival and proliferation of invasive breast cancer cells (43). FOXF1-upregulated genes included those that encode proteins involved in DNA repair (BRCA1 and MRE11A), mitosis (ANAPC2, ANAPC4, and CUL3), transcriptional regulation (CCNT1), and G1 progression (CDK6; Fig. 5C and D), suggesting that FOXF1 might regulate multiple biological processes. Upregulation of DNA repair genes BRCA1 and MRE11A by FOXF1 further supported the role of FOXF1 in maintaining genomic stability. ANAPC2, ANAPC4, and CUL3 are known to inhibit DNA replication and G2 progression during mitosis (44, 45); FOXF1, therefore, might inhibit DNA replication and G2 progression through upregulation of these mitotic genes. Moreover, it has been reported that the expression patterns and function of CDK6 in breast cancer cells exhibited tumor-suppressive properties (46). Therefore, upregulation of CDK6 by FOXF1 might contribute to the growth-suppressive effect of FOXF1 in breast cancer cells.
Discussion
In this article, we report the identification of FOXF1 gene as a novel epigenetic target in breast cancer. DNA hypermethylation was identified as the predominant underlying mechanism for silencing FOXF1 gene expression. Our studies showed that FOXF1 re-expression in breast cancer cell lines inhibited their growth and tumorigenicity, mainly by blocking the G1-S transition of the cell cycle, whereas inactivation of FOXF1 resulted in DNA overreplication, a risk factor leading to genomic instability. Our findings support the hypothesis that aberrant epigenetic inactivation of FOXF1 in mammary epithelial cells could contribute to mammary tumorigenesis.
Gene expression analysis in silico using the microarray database revealed that, compared with normal tissues, FOXF1 is underexpressed in a variety of cancers such as lung, prostate, bladder, ovarian, and breast cancers. Hypermethylation of the FOXF1 promoter was found to be frequent (80% in cell lines, 40% in primary tumors) and is, in all likelihood, the predominant mechanism accounting for the loss or downregulation of FOXF1 expression in breast cancer. FOXF1 promoter hypermethylation positively correlated with the tumor grade and seems to be an early event across both ductal and lobular breast cancers with the same frequency (data not shown). Therefore, the FOXF1 gene is a common epigenetic target among various types of breast cancer. It will be important to decipher whether epigenetic mechanisms are also implicated in downregulating FOXF1 expression in other types of cancers.
Our studies have shed light on the mechanisms underlying FOXF1-mediated growth suppression. We found that enforced overexpression of FOXF1 led to G1 arrest concurrently with or without apoptosis dependent on the cell type examined. Inhibition of CDK2-RB-E2F signaling cascade by FOXF1 is most likely the mechanism accounting for the G1 arrest activity of FOXF1. Although mechanisms for FOXF1-mediated inactivation of the CDK2-RB-E2F cascade are still unclear, our siRNA knockdown studies provided some mechanistic clues. We found that knockdown of FOXF1 in cells led to an increase in DNA re-replication, a sign of a defect in the stringent control of DNA replication initiation. The stringency of DNA replication initiation is finely controlled by a sequential process involving DNA replication origin licensing, activation of origin to replicate DNA, and disassembly of prereplication complexes (47–49). Deregulation of these processes permits repeated firing of DNA replication origins, leading to DNA re-replication (also called overreplication; ref. 48). Expression profiling analysis of cell cycle regulatory genes suggested that FOXF1 inhibits DNA replication initiation by downregulating the expression of DNA replication initiation factor genes (e.g., MCM3 and CDC34) and by upregulating the expression of mitotic genes (e.g., ANAPC2, ANAPC4, and CUL3) involved in negatively regulating G1-S transition. In addition, other FOXF1-regulated cell cycle genes such as CDK6, CDK5R1, and CDK5RAP1 might also contribute to this effect. These findings, taken together, suggest that gain or loss of FOXF1 function alters the expression of these cell cycle regulatory genes, which in turn, leads to these observed outcomes. In addition to MCM3, the expression of five other E2F-induced genes was also elevated in cells depleted of FOXF1, suggesting that E2F transcriptional activity was activated upon silencing of FOXF1. These results are consistent with the data showing that FOXF1 suppresses E2F transcriptional activity through abrogation of the CDK2-RB-E2F cascade.
In addition to cell cycle functions, the identified FOXF1-regulated genes were also implicated in multiple biological processes such as DNA repair, transcriptional regulation, nuclear transportation, and cell survival control. For example, we found that BRCA1 is the downstream target of FOXF1 signaling. BRCA1 is a well-known tumor suppressor gene that plays multiple roles in transcription and DNA repair. BRCA1 is inactivated through genetic mutation or aberrant downregulation of its mRNA expression (50). It is tempting to speculate that epigenetic inactivation of FOXF1 is one of mechanisms leading to the downregulation of BRCA1 expression in breast cancer. Collectively, our findings link the molecular roles of FOXF1 to cell cycle regulation and other biological processes.
In summary, our studies support a novel role for FOXF1 transcription factor in cell cycle regulation; a function that is abrogated through epigenetic silencing of FOXF1 during tumorigenesis. Therefore, FOXF1 and its downstream effectors might serve as novel molecular targets for the detection and treatment of breast cancer and other human cancers.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Acknowledgments
We thank Dr. W. Kaelin, Dana-Farber Cancer Institute at Boston, for kindly providing the E2F-responsive reporter plasmid and E2F1 expression plasmids; and Dr. P. Carlsson, Göteborg University, Sweden, for kindly providing FOXF-responsive reporter plasmids.
Grant Support: NIH Specialized Programs of Research Excellence-1 P50 CA088843-08, DOD Center of Excellence W81XWH-04-1-0595, and the Avon Research Foundation (S. Sukumar).
- Received May 9, 2010.
- Accepted May 16, 2010.
- ©2010 American Association for Cancer Research.