A functional link between FOXA1 and breast cancer SNPs

Genome-wide association studies have revealed a multitude of breast cancer-associated SNPs. The majority of these SNPs are located in noncoding regions of the genome. Yet how they contribute to breast cancer development is unknown. Recently, a groundbreaking study by the Lupien group has shown that risk-associated SNPs of breast cancer are enriched for FOXA1 binding sites, which influences the function of this transcription factor.


Background
Breast cancer is a heterogeneous disease and the most frequent tumors are positive for the expression of estrogen receptor (ER) [1] or the epidermal growth factor receptor 2 (HER2) [2]. Th e transcription factor FOXA1 is a common feature for the majority of these tumors [3,4]. Genomic analyses have revealed that the binding of FOXA1 shows a tissue-specifi c occupancy and that it is associated with the binding of ER or androgen receptor [5][6][7][8]. FOXA1 recruitment to chromatin is mediated by the epigenetic signature consisting of monomethylated and dimethylated histone H3 on lysine 4 (H3K4me1,2) [6].
In the last years, genomic studies have uncovered a multitude of SNPs associated with the risk of breast cancer. However, given that these SNPs are often found in non-coding regions [9], it was unclear how they could be contributing to the development of the disease [10].

Article
In a seminal study published in Nature Genetics, Cowper-Sal Lari and colleagues have integrated genomic methodologies to identify the functional relationship between risk-associated SNPs with high linkage disequilibrium and the transcription factor binding or histone modifi cations in ER-positive breast cancer sub type [11]. To check the enrichment of risk-associated SNPs in chromatin regions, the authors employed variant-set enrichment analysis. Using this statistical analysis tool, they identifi ed that the breast cancer-associated variant set is strongly enriched for FOXA1, ER binding regions and histone modifi cation H3K4me1 marks. Importantly, the associated variant set from other diseases (for example, prostate cancer, colorectal cancer and bone cancer) did not show the same correlation. Based on these fi ndings, the authors state that the enrich ment of the breast cancer-associated SNPs for these transcription factors and histone mark are both cancer type specifi c and tissue specifi c.
Having shown that breast cancer risk-associated SNPs are enriched with FOXA1 binding sites, the authors tested how these SNPs modulate the binding of the transcription factor to chromatin. Th ey used the intragenomic replicates method, which establishes the eff ect of individual SNPs on the affi nity of a transcription factor binding to chromatin. Th e authors focused on one of the most studied risk-associated SNPs: rs4784227, which maps within the forkhead motif for FOXA1 genomic interaction. Interestingly, FOXA1 binding affi nity was increased for the [T] variant allele compared with the [C] reference allele. Extensive intragenomic replicates analysis and variant-set enrichment analysis on all SNP clusters indicated that more than one-half of them are associated with altered affi nity of FOXA1 to chromatin.
To verify the eff ect of the SNP on breast cancer progression the authors focused on the TOX3 gene, which is localized in the rs4784227 risk-associated region. Th ey used chromatin conformation capture assays and they confi rmed a physical interaction between the rs4784227 region and the promoter of the TOX3 gene. Th is result confi rms that TOX3 is a gene target of the regulatory region harboring the SNP. Furthermore, the authors analyzed the impact of the diff erent alleles on TOX3 gene expression performing allele-specifi c expression assay. With this method they showed that the [T] variant allele has a repressive eff ect on TOX3 expression. Finally, the Abstract Genome-wide association studies have revealed a multitude of breast cancer-associated SNPs. The majority of these SNPs are located in noncoding regions of the genome. Yet how they contribute to breast cancer development is unknown. Recently, a groundbreaking study by the Lupien group has shown that risk-associated SNPs of breast cancer are enriched for FOXA1 binding sites, which infl uences the function of this transcription factor.

Viewpoint
Th e published study opens a new way to understand the impact of risk-associated SNPs on breast cancer disease. Importantly, the authors have identifi ed a subset of FOXA1 binding regions enriched for risk-associated SNPs, which immediately raises the question of whether the remaining FOXA1 regions also play a critical role in the progression of the disease. On the contrary, not all breast cancer-associated SNPs were associated with FOXA1 binding regions. In future studies it will therefore be of interest to investigate the eff ect of these SNPs on other transcription factors.
Furthermore, the study concludes that more than onehalf of the risk-associated SNPs infl uence FOXA1 binding [11]. FOXA1 has been described to function as a pioneer factor that allows the binding of ER and Androgen Receptor [5]. Recently, the expression of FOXA1 both in ER and HER2 breast cancer cell lines has been also reported to repress transcription [12]. Whether the SNPs are infl uencing both or any of these functions needs to be determined.
Finally, the study also reports that the depletion of TOX3 increased signifi cantly the proliferation of breast cancer cells, suggesting that the rs4784227 riskassociated region might be used as a clinical marker for cancer prognosis. Investigating how other SNPs could impact current breast cancer treatments will be important. Th e use of SNPs as biomarkers to predict anticancer drug effi cacy and risk levels in breast cancer patients is a promising path to explore.