Mendoza-Fandiño, Gustavo

Abstract Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. Significance: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439

Abstract Genome Wide Association Studies (GWAS) are designed to identify single nucleotide polymorphisms (SNPs) mapping to genomic regions associated with a particular trait or disease. Currently, 19 genomic loci associated with risk for estrogen receptor negative breast cancer have been identified. However, for most of these loci the mechanisms driving susceptibility have remained elusive. The 2p23.2 locus, defined by SNP rs67073037 (p = 4.76 x 10-9), contains 79 SNPs exhibiting genome-wide significant associations as well as four candidate genes: SPDYA, TRMT61B, WDR43, and FAM179A. Post-GWAS functional analysis revealed two candidate causal regulatory elements: an enhancer-like element defined by SNP rs44017214 located in the first intron of the WDR43 gene; and a segment in the 3’-UTR region of WDR43 defined by SNPs rs11680458 and rs11331880 which are predicted to have allele specific binding of non-coding microRNAs miR-548 and miR-376, respectively. We validated allele specific transcription factor binding to rs44017214 using competitive binding of USF1/USF2 to the risk allele in both MCF10A and CAL51 cell lines in Electrophoretic Mobility Shift Assays (EMSA). Additionally, rs44017214 was removed using CRISPR-cas9, and expression of neighboring genes was evaluated by qRT-PCR. No change in expression was found for any gene in a 2 Mb region surrounding the SNP, with the exception of PLB1. This suggests that the enhancer element in the intron of WDR43 has PLB1 as its target. To assess the miRNA binding site impact on WDR43 expression the pMIR-REPORT Luciferase allele-specific vector was co-transfected with miR-141/200, miR548, or miR-376 in MCF10A cells. The rs1131880 allele C in the WDR43 3’-UTR lead to the suppression of luciferase expression by miR-376 compared to the G allele. This study reveals two independent regulatory mechanisms within the 2p23.2 locus; (i) modulation of enhancer activity by allele specific binding of transcription factors USF1/USF2 to rs4401714, and (ii) decreasing the expression of WDR43 through the formation of a miR-376 binding site at the 3’-UTR of the WDR43 gene. Elucidating the functional significance of these common variants in the novel susceptibility loci 2p23.2 furthers our knowledge of the etiology of estrogen receptor negative breast cancer. Citation Format: Carly M. Harro, Gustavo Mendoza-Fandino, Nicholas T. Woods, Xueli Li, Fergus J. Couch, Alvaro N. Monteiro. Two distinct regulatory mechanisms underlie estrogen receptor negative breast cancer susceptibility at the 2p23.2 locus [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1306. doi:10.1158/1538-7445.AM2017-1306

Abstract An ovarian cancer Genome Wide Association Study (GWAS) identified 9p22.2 as a novel susceptibility locus with the most statistically significant single nucleotide polymorphisms (SNPs) located in an intergenic region near the Basonuclin 2 (BNC2) gene, which codes for a putative transcription regulator containing three pairs of zinc finger (ZF) domains. The minor alleles are protective in terms of ovarian cancer susceptibility. However, the molecular mechanisms by which these SNPs modify susceptibility remain largely unknown. The significant SNPs in the 9p22.2 locus lie in non-coding regions and therefore are hypothesized to affect the activity of regulatory elements and modify the expression of a target gene(s). In order to test this hypothesis we conducted fine mapping for the locus which delimited a ∼74kb region containing SNPs with a p-value less than 10-8. FAIRE-Seq and ChIP-Seq experiments for histone markers conducted in immortalized ovarian surface epithelial cells (IOSE) and fallopian tube epithelial cells (IFTE) were used to prioritize functional SNP candidates which overlap with regulatory elements. Luciferase assays tested the ability of these regulatory elements to activate transcription. Chromosome conformation capture (3C) experiments demonstrate a physical interaction between the BNC2 promoter and candidate regulatory elements containing risk-associated SNPs. Methylation Quantification at Trait Loci (mQTL) revealed an association between decreased methylation at the BNC2 promoter and the protective minor alleles. Expression analysis shows decreased expression of BNC2 in cancer versus normal tissue implicating tumor suppressor function of BNC2. Therefore the protective minor allele likely increases expression of BNC2 which in turn contains tumor suppressing properties that decrease ovarian cancer risk. BNC2 protein is in complex with transcriptional regulatory proteins, in particular the NURD complex components. Moreover, BNC2 acts as a transcriptional repressor in in vitro transfection assays indicating that it functions in transcriptional repression. We then used protein binding microarrays and CHIP-Seq experiments to identify its putative DNA binding motifs and its downstream target genes. Analysis of this data suggests that BNC2 functions in a regulatory transcription network that impacts on genes implicated in ovarian cancer with enrichment for genes in the TGF-beta response pathway. Citation Format: Melissa A. Buckley, Howard C. Shen, Gustavo A. Mendoza-Fandino, Nicholas T. Woods, Anxhela Gjyshi, Juliet French, Kate Lawrenson, Honglin Song, Jonathan Tyrer, Renato S. Carvalho, Alexandra Valle, Ann Chen, Sean Yoder, Gregory Bloom, Ya-Yu Tsai, Ally Yang, Timothy R. Hughes, Xiaotao Qu, Mine Cicek, Melissa Larson, Ellen Goode, Brooke Fridley, Susan Ramus, Georgia Chenevix-Trench, Paul Pharoah, Thomas A. Sellers, Simon Gayther, Alvaro N.A. Monteiro, Ovarian Cancer Association Consortium. Functional analysis of the 9p22 locus implicates the transcriptional regulation of BNC2 as a mechanism in ovarian cancer predisposition. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3285. doi:10.1158/1538-7445.AM2014-3285