Pharoah, Paul D.P.

Abstract Background: Large-scale genome-wide association studies (GWAS) have identified ~35 loci associated with epithelial ovarian cancer risk. However, these genetic loci together explain only a small portion of the heritability of this malignancy. The large majority of the GWAS-identified variants are located in non-coding regions, thus possible causal genes in these loci remain largely unknown. We performed a transcriptome-wide association study (TWAS) to search for novel genetic loci for ovarian cancer risk and plausible causal genes at GWAS-identified loci. Method: Transcriptome data from normal ovarian tissue samples (n=68) and all tissue samples (n = 369), along with their high-density genotyping data, obtained from only European descendants included in the Genotype-Tissue Expression Project (GTEx), were used to build ovarian linear prediction models and cross-tissue models (to increase statistical power) using the elastic net method. Based on model performance, we evaluated 17,121 genes for their cis-predicted gene expressions in relation to ovarian cancer risk using summary statistics data generated in GWAS of ovarian cancer from 25,509 cases and 40,941 controls. MetaXcan was used to integrate gene expression prediction model with summary statistics. Results: We identified 35 genes with predicted expression levels associated with ovarian cancer risk at P value < 2.2 × 10-6, the Bonferroni corrected significance level for multiple comparisons. Of these, 12 genes at 4 genetic loci are located ≥500kb away from risk SNPs previously reported in GWAS, representing potential novel genetic loci for ovarian cancer risk. The remaining 23 genes at 12 loci are located within known ovarian cancer risk loci. Fifteen of these genes at the 12 loci have not been reported in previous studies. Analyses by ovarian cancer histological subtypes showed that the majority of these 35 genes are associated with serous invasive carcinoma. Several new associations were identified in histological subtype analyses. Conclusion: In this TWAS we identified multiple genes with predicted expressions related to ovarian cancer risk and provide substantial new information to enhance the understanding of ovarian cancer biology and genetics. Citation Format: Yingchang Lu, Joellen M. Schildkraut, Thomas A. Sellers, Lang Wu, Xingyi Guo, Bingshan Li, Y. Ann Chen, Jennifer B. Doherty, Simon Gayther, Ellen L. Goode, Hae Kyung Im, Siddhartha Kar, Kate Lawrenson, Ani W. Manichaikul, Jennifer B. Permuth, Brett M. Reid, Jamie K. Teer, Paul Pharoah, Wei Zheng, Jirong Long. Transcriptome-wide association study among 66,450 women to identify candidate susceptible genes for ovarian cancer risk [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 1308. doi:10.1158/1538-7445.AM2017-1308

Abstract Epithelial ovarian cancer (EOC) is a rare but deadly disease for which there is notably poorer survival in women of African Ancestry (AA) compared to women of European Ancestry (EA). Previous EA-based genome-wide association studies (GWAS) have identified 30 common, low penetrant EOC susceptibility alleles. Using the custom-designed 533,631 SNP Illumina OncoArray and imputation to ~12 million genetic variants in the 1000 Genomes Phase 3, we conducted a GWAS in 755 AA EOC cases, including 537 high-grade serous ovarian cancer (HGSOC) cases, and 1,235 AA controls from the Ovarian Cancer Association Consortium (OCAC). We identified novel susceptibility loci with suggestive evidence of association with EOC, based on a threshold of P<1x10-6, at three distinct loci including 10p15.1 (lead SNP rs4525119, intronic to AKR1C3, P=4.9 x 10-7, minor allele frequency [MAF]=0.33), 3p25.3 (lead SNP rs7643459, intronic to LOC101927394, P=8.4 x 10-7, MAF=0.36), and 4q13.3 (lead SNP rs4286604, 15 kb 3’ of UGT2A1, P=8.5x10-7, MAF=0.27). In analysis of HGSOC, we identified suggestive evidence of association at six distinct loci including 5q11.2 (lead SNP rs37792, 132 kb 5’ of FST, P=6.0x10-8, MAF=0.34), Xq27.2 (lead SNP rs57403204, 81 kb 3' of MAGEC1, P=1.7x10-7, MAF=0.06), 10p15.1 (lead SNP rs79079890, LOC105376360 intronic, P=3.0x10-7, MAF=0.03), 17p25.1 (lead SNP rs66459581, 5 kb 5’ of PRPSAP1, P=5.1x10-7, MAF=0.23), 15p12 (lead SNP rs116046250, GABRG3 intronic, P=8.7x10-7, MAF=0.05), and 4q21.21 (lead SNP rs192876988, 30 kb 3’ of GK2, P=9.2x10-7, MAF=0.05). These SNPs showed no evidence of association with risk in the OCAC EA GWAS including up to 23,543 EOC cases and 29,444 controls. Of the SNPs identified with suggestive evidence of association in our GWAS, two are located near genes known to regulate hormones and previously reported in relation to diseases of the ovary (AKR1C3 and FST), two have been linked to cancer (AKR1C3 and MAGEC1) and one represents testis-specific expression (GK2). Of the 30 SNPs identified previously in a EA EOC GWAS, we observed nominally significant associations (P<0.05) with consistent direction of effect for eight SNPs in AA women. Our investigation presents evidence of (1) variants for EOC shared among EA and AA women and (2) novel EOC risk loci of importance for AA women. Citation Format: Ani Manichaikul, Lauren C. Peres, Xin-Qun Wang, Kate Lawrenson, Sarah Abbott, Ann G. Schwartz, Anna H. Wu, Edward Peters, Patricia G. Moorman, Michele L. Cote, Melissa Bondy, Linda Kelemen, Ellen L. Goode, Jill Barnholtz-Sloan, Simon A. Gather, Andrew Berchuck, Jennifer A. Doherty, Paul Pharoah, Joellen Schildkraut. Identification of novel epithelial ovarian cancer loci in women of African ancestry from the Ovarian Cancer Association Consortium [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 1307. doi:10.1158/1538-7445.AM2017-1307

Investigating the most likely causal variants identified by fine‐mapping analyses may improve the power to detect gene–environment interactions. We assessed the interplay between 70 single nucleotide polymorphisms identified by genetic fine‐scale mapping of susceptibility loci and 11 epidemiological breast cancer risk factors in relation to breast cancer. Analyses were conducted on up to 58,573 subjects (26,968 cases and 31,605 controls) from the Breast Cancer Association Consortium, in one of the largest studies of its kind. Analyses were carried out separately for estrogen receptor (ER) positive (ER+) and ER negative (ER–) disease. The Bayesian False Discovery Probability (BFDP) was computed to assess the noteworthiness of the results. Four potential gene–environment interactions were identified as noteworthy (BFDP < 0.80) when assuming a true prior interaction probability of 0.01. The strongest interaction result in relation to overall breast cancer risk was found between CFLAR‐rs7558475 and current smoking (ORint = 0.77, 95% CI: 0.67–0.88, pint = 1.8 × 10−4). The interaction with the strongest statistical evidence was found between 5q14‐rs7707921 and alcohol consumption (ORint =1.36, 95% CI: 1.16–1.59, pint = 1.9 × 10−5) in relation to ER– disease risk. The remaining two gene–environment interactions were also identified in relation to ER– breast cancer risk and were found between 3p21‐rs6796502 and age at menarche (ORint = 1.26, 95% CI: 1.12–1.43, pint =1.8 × 10−4) and between 8q23‐rs13267382 and age at first full‐term pregnancy (ORint = 0.89, 95% CI: 0.83–0.95, pint = 5.2 × 10−4). While these results do not suggest any strong gene–environment interactions, our results may still be useful to inform experimental studies. These may in turn, shed light on the potential interactions observed.

Abstract In this review, we summarize current progress in the genetic epidemiology of epithelial ovarian cancer (EOC), focusing exclusively on elucidating the role of common germline genetic variation in conferring susceptibility to EOC. We provide an overview of the more than 30 EOC risk loci identified to date by genome-wide association studies (GWAS) and describe the contribution of large-scale, cross-cancer type, custom genotyping projects, such as the OncoArray and the Collaborative Oncological Gene-Environment Study, to locus discovery and replication. We discuss the histotype-specific nature of these EOC risk loci, pleiotropy, or overlapping genetic effects between EOC and other hormone-related cancer types, and the application of findings to polygenic risk prediction for EOC. The second part of the article offers a concise review of primarily laboratory-based studies that have led to the identification of several putative EOC susceptibility genes using common variants at the known EOC risk loci as starting points. More global biological insights emerging from network- and pathway-based analyses of GWAS for EOC susceptibility are also highlighted. Finally, we delve into potential future directions, including the need to identify EOC risk loci in non-European populations and the next generation of GWAS functional studies that are likely to involve genome editing to establish the cell type–specific carcinogenic effects of EOC risk variants Cancer Epidemiol Biomarkers Prev; 27(4); 395–404. ©2018 AACR. See all articles in this CEBP Focus section, “Genome-Wide Association Studies in Cancer.”

Abstract Breast cancer risks conferred by many germline missense variants in the BRCA1 and BRCA2 genes, often referred to as variants of uncertain significance (VUS), have not been established. In this study, associations between 19 BRCA1 and 33 BRCA2 missense substitution variants and breast cancer risk were investigated through a breast cancer case–control study using genotyping data from 38 studies of predominantly European ancestry (41,890 cases and 41,607 controls) and nine studies of Asian ancestry (6,269 cases and 6,624 controls). The BRCA2 c.9104A&gt;C, p.Tyr3035Ser (OR = 2.52; P = 0.04), and BRCA1 c.5096G&gt;A, p.Arg1699Gln (OR = 4.29; P = 0.009) variant were associated with moderately increased risks of breast cancer among Europeans, whereas BRCA2 c.7522G&gt;A, p.Gly2508Ser (OR = 2.68; P = 0.004), and c.8187G&gt;T, p.Lys2729Asn (OR = 1.4; P = 0.004) were associated with moderate and low risks of breast cancer among Asians. Functional characterization of the BRCA2 variants using four quantitative assays showed reduced BRCA2 activity for p.Tyr3035Ser compared with wild-type. Overall, our results show how BRCA2 missense variants that influence protein function can confer clinically relevant, moderately increased risks of breast cancer, with potential implications for risk management guidelines in women with these specific variants. Cancer Res; 77(11); 2789–99. ©2017 AACR.

Abstract Background: Genome-wide association studies (GWAS) have identified multiple loci associated with epithelial ovarian cancer (EOC) susceptibility, but further progress requires integration of epidemiology and biology to illuminate true risk loci below genome-wide significance levels (P &lt; 5 × 10−8). Most risk SNPs lie within non–protein-encoding regions, and we hypothesize that long noncoding RNA (lncRNA) genes are enriched at EOC risk regions and represent biologically relevant functional targets. Methods: Using imputed GWAS data from about 18,000 invasive EOC cases and 34,000 controls of European ancestry, the GENCODE (v19) lncRNA database was used to annotate SNPs from 13,442 lncRNAs for permutation-based enrichment analysis. Tumor expression quantitative trait locus (eQTL) analysis was performed for sub-genome-wide regions (1 × 10−5 &gt; P &gt; 5 × 10−8) overlapping lncRNAs. Results: Of 5,294 EOC-associated SNPs (P &lt; 1.0 × 10−5), 1,464 (28%) mapped within 53 unique lncRNAs and an additional 3,484 (66%) SNPs were correlated (r2 &gt; 0.2) with SNPs within 115 lncRNAs. EOC-associated SNPs comprised 130 independent regions, of which 72 (55%) overlapped with lncRNAs, representing a significant enrichment (P = 5.0 × 10−4) that was more pronounced among a subset of 5,401 lncRNAs with active epigenetic regulation in normal ovarian tissue. EOC-associated lncRNAs and their putative promoters and transcription factors were enriched for biologically relevant pathways and eQTL analysis identified five novel putative risk regions with allele-specific effects on lncRNA gene expression. Conclusions: lncRNAs are significantly enriched at EOC risk regions, suggesting a mechanistic role for lncRNAs in driving predisposition to EOC. Impact: lncRNAs represent key candidates for integrative epidemiologic and functional studies. Further research on their biologic role in ovarian cancer is indicated. Cancer Epidemiol Biomarkers Prev; 26(1); 116–25. ©2016 AACR.