Lawrenson, Kate

Abstract Genome-wide association studies (GWAS) have identified several susceptibility regions that confer risk of two or more cancers. A region at 19p13 is associated with both breast and ovarian cancer risk in the general population as well as in BRCA1/BRCA2 mutation carriers. Using the iCOGs custom genotyping array, we analyzed 438 SNPs in this region in 7,435 ER-negative breast cancer cases and 42,599 controls from BCAC, and 15,347 ovarian cancer cases and 30,845 controls from OCAC and 15,252 BRCA1 mutation carriers. We identified two significant peaks of association. Peak 1, consisting of 13 highly correlated SNPs (r2>0.95), spans a 19.4kb region harboring the BABAM1, ANKLE1 and ABHD8 genes. Peak 1 SNPs are associated with serous ovarian cancer (top hit rs4808075, P = 9.2×10-20), ER-negative breast cancer (rs61494113, P = 1.1×10-13) and breast cancer in BRCA1 mutation carriers (rs61494113, P = 7.7×10-16). Peak 2 spans over 100kb and contains 7 SNPs associated with breast cancer in BRCA1 mutation carriers only (rs4464206, P = 8.9×10-8). Multiple assays were used to establish the putative functional roles of risk SNPs and candidate susceptibility genes in this region. Expression quantitative trait locus (eQTL) analyses identified significant genotype-gene expression correlations between genotype at the Peak 1 SNP rs480816 and ABHD8 expression in ovarian tumors (P = 3.0×10-5) and normal breast tissues (P = 2.8×10-3). The Peak 2 SNP, rs7246243, was significantly associated with BABAM1 expression (P = 5×10-3) in normal breast tissue only. Overexpression of BABAM1 in normal ovarian and breast cells did not modulate neoplastic phenotypes, whereas overexpression of ABHD8 induced a significant reduction in migration in both cell types in vitro. Differentially expressed genes in ABHD8 models, identified by RNAseq, were associated with cellular movement (P<1.82×10-44). Chromosome conformation capture (3C) identified an interaction between the ABHD8 promoter and a region containing 4 risk SNPs (rs4808075, rs10419397, rs56069439 and rs4808076) in breast and ovarian normal and cancer cells but only rs56069439 coincides with epigenetic enhancer marks in breast and ovarian cells. CRISPR/Cas9 genome editing was used to delete a 57bp region containing rs56069439 and expression of BABAM1, ANKLE1 and ABHD8 were measured. Taken together, these data suggest a common functional mechanism underlying risk of breast and ovarian cancer at this locus, mediated by rs56069439 variants differentially regulating expression of ABHD8, a gene involved in cell motility. Citation Format: Kate Lawrenson, Siddhartha Kar, Karoline Kuckenbaeker, Stacey Edwards, Qiyuan Li, Jonathan Tyrer, Jonathan Beesley, Kyriaki Michailidou, Susan Ramus, Alison Dunning, Fergus Couch, Alvaro Monteiro, Jacques Simard, Matthew Freedman, Douglas Easton, Georgia Chenevix-Trench, Paul Pharoah, Antonis Antoniou, Simon Gayther, OCAC, BCAC and CIMBA. Common functional mechanisms underlying pleiotropy at the 19p13.1 breast and ovarian cancer cusceptibility locus. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2783. doi:10.1158/1538-7445.AM2015-2783

Abstract Background: Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations. Methods: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls). Results: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network. Conclusion: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development. Impact: Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization. Cancer Epidemiol Biomarkers Prev; 24(10); 1574–84. ©2015 AACR.

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

Abstract Since most GWAS hits fall in non-coding regions of the genome and long non-coding RNAs (lncRNAs) are emerging as drivers of carcinogenesis, we hypothesized that SNPs in lncRNAs influence epithelial ovarian cancer (EOC) risk. A comprehensive lncRNA database (lncRNA db) that contained 104 human lncRNAs when downloaded in November 2011 was used to identify 1,737 variants in 63 unique lncRNAs that were represented in three genome-wide association studies from North America, the United Kingdom, and Poland (3,995 EOC cases and 3,277 controls of European background). SNPs with MAF <5% were excluded, and missing genotypes were inferred with Mach using the HapMap CEU population. Unconditional logistic regression treating alternate alleles as an ordinal variable was used to evaluate individual SNP-EOC risk associations. Subgroup analysis was conducted for serous adenocarcinomas (N= 2,066), the main histologic subtype of EOC. Models were adjusted for study and European ancestry using the first two principal components (PCs). We reduced correlated signals (r2 ≥ 0.8) within each lncRNA and kept the most statistically significant variants based on Wald tests, resulting in 651 SNPs for further analyses. Fisher's product method was used to aggregate evidence of multiple variants within each lncRNA using 10,000 permutations. To adjust for multiple comparisons, a false discovery rate (FDR) of 10% was used. A replication phase was carried out for top-ranked SNPs from an additional 9,854 EOC cases (5,802 serous) and 17,633 controls genotyped through an international effort known as the Collaborative Oncological Gene-Environment study (COGS). Four lncRNAs, WT1-AS (p = 0.005, n = 2 SNPs), Hoxa11as (p = 0.013, n = 2), AK082072 (p = 0.016, n = 4) and H19 (p = 0.020, n = 2), were significantly associated with EOC risk among all histologies (FDR ≤ 10%). The SNP in WT1-AS most significantly associated with EOC risk was rs3809061 (T>C) (minor allele frequency= 36%, odds ratio (OR) = 0.91, CI: 0.84-0.97, p = 0.008). Among serous cases, rs3809061 was also protective (OR = 0.92, p = 0.04). These findings were replicated by a correlated SNP rs2301250 (r2 = 0.84) evaluated in COGS among serous cases (OR = 0.94, CI: 0.90-0.99, p = 0.009) and less significant among all histologies (OR = 0.96, CI: 0.93-1.00, p = 0.054). We further interrogated its functional role by performing an expression quantitative trait locus analysis on WTI-AS mRNA expression using the Cancer Genome Atlas (TCGA) data (N = 462 serous cases). WT1-AS expression was significantly higher among carriers of the variant allele (compared to TT homozygotes) using the Mann-Whitney U test (p = 0.001). WT1-AS is the antisense and possible regulator of WT1, a tumor suppressor gene reported to be prognostic in advanced EOC. These findings implicate germline lncRNA variants associated with EOC risk that merit further investigation. Citation Format: Yian Ann Chen, Zhihua Chen, Jennifer Permuth-Wey, Ya-Yu Tsai, Hui-Yi Lin, Xiaotao Qu, Kate Lawrenson, David Fenstermacher, Catherine M. Phelan, Alvaro Monteiro, Simon A. Gayther, Steven A. Narod, Rebecca Sutphen, Michael J. Birrer, Nicolas Wentzensen, Joellen M. Schildkraut, Ellen L. Goode, Paul Pharoah, Thomas Sellers. Variants in long non-coding RNAs are associated with epithelial ovarian cancer risk in a pooled analysis of three genome-wide association studies. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 4579. doi:10.1158/1538-7445.AM2013-4579

Abstract The presence of regulatory T cells (Treg) in solid tumors is known to play a role in patient survival in ovarian cancer and other malignancies. We assessed inherited genetic variations via 749 tag single-nucleotide polymorphisms (SNP) in 25 Treg-associated genes (CD28, CTLA4, FOXP3, IDO1, IL10, IL10RA, IL15, 1L17RA, IL23A, IL23R, IL2RA, IL6, IL6R, IL8, LGALS1, LGALS9, MAP3K8, STAT5A, STAT5B, TGFB1, TGFB2, TGFB3, TGFBR1, TGRBR2, and TGFBR3) in relation to ovarian cancer survival. We analyzed genotype and overall survival in 10,084 women with invasive epithelial ovarian cancer, including 5,248 high-grade serous, 1,452 endometrioid, 795 clear cell, and 661 mucinous carcinoma cases of European descent across 28 studies from the Ovarian Cancer Association Consortium (OCAC). The strongest associations were found for endometrioid carcinoma and IL2RA SNPs rs11256497 [HR, 1.42; 95% confidence interval (CI), 1.22–1.64; P = 5.7 × 10−6], rs791587 (HR, 1.36; 95% CI, 1.17–1.57; P = 6.2 × 10−5), rs2476491 (HR, = 1.40; 95% CI, 1.19–1.64; P = 5.6 × 10−5), and rs10795763 (HR, 1.35; 95% CI, 1.17–1.57; P = 7.9 × 10−5), and for clear cell carcinoma and CTLA4 SNP rs231775 (HR, 0.67; 95% CI, 0.54–0.82; P = 9.3 × 10−5) after adjustment for age, study site, population stratification, stage, grade, and oral contraceptive use. The rs231775 allele associated with improved survival in our study also results in an amino acid change in CTLA4 and previously has been reported to be associated with autoimmune conditions. Thus, we found evidence that SNPs in genes related to Tregs seem to play a role in ovarian cancer survival, particularly in patients with clear cell and endometrioid epithelial ovarian cancer. Cancer Immunol Res; 2(4); 332–40. ©2014 AACR.

Abstract Genome-wide association analyses have identified an ovarian cancer susceptibility locus near BNC2 at 9p22.2 (OR 0.82, p 5.1 × 10−19, Nat Genet 41:996-1000, 2009), based on 1,772 and 2,354 Stage I cases and controls, 4,803 and 5,237 Stage 2 cases and controls, and extensive Stage 3 analysis in independent samples. Here, we report on Stage 3 analysis of nine additional regions with Stage 2 results at p < 10−4. Stage 1 samples were genotyped using the Illumina 610Quad and imputed at up to 2 million SNPs; Stage 2 samples were genotyped on an Illumina iSelect of 22,790 SNPs. Stage 3 samples from 4,335 cases (including 2,197 serous subtype) and 5,951 controls (14 studies) were genotyped using Illumina or Sequenom at up to 32 SNPs on chromosomes 1, 2, 3, 7, 8, 11, and 17. No evidence of heterogeneity by study site was observed. Logistic regression revealed four regions with ordinal p ≤ 10−7 in all cases or serous subtype only. The most significant finding for all cases was for rs2072590 at 2q31.1 (all cases OR 1.15, p=8.1 × 10−13; serous only OR 1.19 p=9.0 × 10−13) which is located 3.8 kb 3-prime of HOXD3, a homeobox gene which also had an intronic SNP (rs2301301) with p=3.1 × 10−5 (all cases). Second, we observed a locus at 8q24.21 highlighted by rs10088218 (all cases OR 0.84, p=1.3 × 10−9; serous only OR 0.75 p=2.2 × 10−14). This locus resides on the opposite side of MYC and 1.3 Mb from the 8q24 regions harboring loci for other cancers, and is independent of the other cancer-associated variants (r2 < 0.02). Third, at 3q25.31, promising results were observed at rs2665390, a SNP intronic to TIPARP (TiPARP, TCDD-inducible poly(ADP-ribose) polymerase; all cases OR 1.19, p=8.7 × 10−7; serous only OR 1.26, p=3.2 × 10−8). Finally, at 17q21.32, rs9303542 was associated with risk particularly among serous cases (all cases OR 1.13, p=1.5 × 10−6; serous only OR 1.16 p=4.3 × 10−7); this SNP is intronic to SKAP1 (src kinase associated phosphoprotein 1). We evaluated the functional role of genes at these loci in ovarian cancer and normal ovarian epithelial cells and in an in vitro model of neoplastic initiation and progression from normal ovarian epithelial cells after CYMC and KRAS activation. Both BNC2 and TiPARP were down regulated in cancer versus normal cells (p=5 × 10−9 and p=3.4 × 10−5, respectively) suggesting that loss of function is required for cancer development. Over-expressing CMYC (at 8q24) induced neoplastic transformation of normal ovarian cells. Both BNC2 and TiPARP were down-regulated following CMYC over-expression and showed significant loss of function with further transformation after mutant KRAS activation (p=7.8 × 10−7 and p=6.3 × 10−4, respectively). These results emphasize the importance of genome-wide association studies and consortia for the identification of moderate-risk alleles in rare cancers. Importantly, our findings unravel the genetic components of this disease and suggest functionally-significant novel targets. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 3860.