DeFazio, Anna

Abstract Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10−6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10−7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10−3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis. Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419–30. ©2018 AACR.

Abstract Background: Comorbidities can affect survival of ovarian cancer patients by influencing treatment efficacy. However, little evidence exists on the association between individual concurrent comorbidities and prognosis in ovarian cancer patients. Methods: Among patients diagnosed with invasive ovarian carcinoma who participated in 23 studies included in the Ovarian Cancer Association Consortium, we explored associations between histories of endometriosis; asthma; depression; osteoporosis; and autoimmune, gallbladder, kidney, liver, and neurological diseases and overall and progression-free survival. Using Cox proportional hazards regression models adjusted for age at diagnosis, stage of disease, histology, and study site, we estimated pooled HRs and 95% confidence intervals to assess associations between each comorbidity and ovarian cancer outcomes. Results: None of the comorbidities were associated with ovarian cancer outcome in the overall sample nor in strata defined by histologic subtype, weight status, age at diagnosis, or stage of disease (local/regional vs. advanced). Conclusions: Histories of endometriosis; asthma; depression; osteoporosis; and autoimmune, gallbladder, kidney, liver, or neurologic diseases were not associated with ovarian cancer overall or progression-free survival. Impact: These previously diagnosed chronic diseases do not appear to affect ovarian cancer prognosis. Cancer Epidemiol Biomarkers Prev; 26(9); 1470–3. ©2017 AACR.

Abstract Background: Many epithelial ovarian cancer (EOC) risk factors relate to hormone exposure and elevated estrogen levels are associated with obesity in postmenopausal women. Therefore, we hypothesized that gene–environment interactions related to hormone-related risk factors could differ between obese and non-obese women. Methods: We considered interactions between 11,441 SNPs within 80 candidate genes related to hormone biosynthesis and metabolism and insulin-like growth factors with six hormone-related factors (oral contraceptive use, parity, endometriosis, tubal ligation, hormone replacement therapy, and estrogen use) and assessed whether these interactions differed between obese and non-obese women. Interactions were assessed using logistic regression models and data from 14 case–control studies (6,247 cases; 10,379 controls). Histotype-specific analyses were also completed. Results: SNPs in the following candidate genes showed notable interaction: IGF1R (rs41497346, estrogen plus progesterone hormone therapy, histology = all, P = 4.9 × 10−6) and ESR1 (rs12661437, endometriosis, histology = all, P = 1.5 × 10−5). The most notable obesity–gene–hormone risk factor interaction was within INSR (rs113759408, parity, histology = endometrioid, P = 8.8 × 10−6). Conclusions: We have demonstrated the feasibility of assessing multifactor interactions in large genetic epidemiology studies. Follow-up studies are necessary to assess the robustness of our findings for ESR1, CYP11A1, IGF1R, CYP11B1, INSR, and IGFBP2. Future work is needed to develop powerful statistical methods able to detect these complex interactions. Impact: Assessment of multifactor interaction is feasible, and, here, suggests that the relationship between genetic variants within candidate genes and hormone-related risk factors may vary EOC susceptibility. Cancer Epidemiol Biomarkers Prev; 25(5); 780–90. ©2016 AACR.

Abstract Background: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC). Methods: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped). Results: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E−6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E−6; Pcorrected = 0.01). Conclusions: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed. Impact: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446–54. ©2016 AACR.

Abstract Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk. Cancer Epidemiol Biomarkers Prev; 24(11); 1680–91. ©2015 AACR.

Abstract Purpose: Chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that germline polymorphisms might be associated with clinical outcome. Experimental Design: We analyzed approximately 2.8 million genotyped and imputed SNPs from the iCOGS experiment for progression-free survival (PFS) and overall survival (OS) in 2,901 European epithelial ovarian cancer (EOC) patients who underwent first-line treatment of cytoreductive surgery and chemotherapy regardless of regimen, and in a subset of 1,098 patients treated with ≥4 cycles of paclitaxel and carboplatin at standard doses. We evaluated the top SNPs in 4,434 EOC patients, including patients from The Cancer Genome Atlas. In addition, we conducted pathway analysis of all intragenic SNPs and tested their association with PFS and OS using gene set enrichment analysis. Results: Five SNPs were significantly associated (P ≤ 1.0 × 10−5) with poorer outcomes in at least one of the four analyses, three of which, rs4910232 (11p15.3), rs2549714 (16q23), and rs6674079 (1q22), were located in long noncoding RNAs (lncRNAs) RP11-179A10.1, RP11-314O13.1, and RP11-284F21.8, respectively (P ≤ 7.1 × 10−6). ENCODE ChIP-seq data at 1q22 for normal ovary show evidence of histone modification around RP11-284F21.8, and rs6674079 is perfectly correlated with another SNP within the super-enhancer MEF2D, expression levels of which were reportedly associated with prognosis in another solid tumor. YAP1- and WWTR1 (TAZ)-stimulated gene expression and high-density lipoprotein (HDL)-mediated lipid transport pathways were associated with PFS and OS, respectively, in the cohort who had standard chemotherapy (pGSEA ≤6 × 10−3). Conclusions: We have identified SNPs in three lncRNAs that might be important targets for novel EOC therapies. Clin Cancer Res; 21(23); 5264–76. ©2015 AACR.

Abstract Primary and acquired chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that heritable polymorphisms might be associated with clinical outcome. Through the Ovarian Cancer Association Consortium, we have collected detailed clinical data from a large number of ovarian cancer patients. We have conducted association analyses of progression-free survival (PFS) and overall survival (OS) for 2,901 European women with ovarian carcinoma for ∼8 million imputed SNPs with minor allele frequency ≥ 0.02 and good imputation quality (imputation r2 ≥ 0.3) from the iCOGS experiment (Pharoah et al Nat Genet. 2013). Cox regression models adjusted for known prognostic factors were used to assess the log-additive effects of SNPs on PFS and OS in 2,825 patients who underwent cytoreductive surgery and received any chemotherapy (‘all chemo’) and in the subset of 1,098 patients who underwent cytoreductive surgery followed by ≥ 4 cycles of carboplatin AUC 5 or 6 and paclitaxel 175 or 135 mg/m2 intravenously (‘standard chemo’). We identified five SNPs at four loci that were significantly associated with outcome (p ≤ 10-7) in at least one of four analyses (Table). Analysis of expression data for genes within a 1Mb region flanking each lead SNP indicates significant associations between nearby genes and ovarian cancer survival. Of particular interest is the heat shock protein crystalline αB (CRYAB) at 11q23.1. We have also identified significant associations with KIAA1598 and RAVER2. Further analysis of these loci in a larger genotyped sample is necessary to definitively identify SNP(s) associated with outcome. Our study adds to the growing evidence that genetic variants are associated with ovarian cancer prognosis. Table: SNPs associated with ovarian cancer outcome and gene expression Overall Survival (P)*Progression-Free Survival (P)*Lead SNPChromosome LocusImputation R²Expression AnalysisAll chemoStandard chemoAll chemoStandard chemors1429768671p31.30.742.14E-0512.69E-044.68E-055.11E-073.64E-09rs559375217q32.30.600.198.58E-082.32E-052.80E-031.85E-03rs7885346210q26.110.602.23E-0724.57E-036.57E-081.65E-016.85E-03rs1089121911q23.10.841.53E-0531.51E-059.50E-081.50E-042.32E-04rs3579103211q23.10.811.53E-0531.66E-055.88E-086.46E-051.87E-04*P-values are adjusted for known confounders of ovarian cancer outcomes; FIGO stage, grade, histology, residual disease, and additionally for age at diagnosis for OS.1 RAVER2; 2 KIAA1598; 3 CRYAB Citation Format: Sharon E. Johnatty, Jonathan Tyrer, Jonathan Beesley, Bo Gao, Yi Lu, Stuart MacGregor, Anna deFazio, Paul Pharoah, Ellen Goode, Georgia Chenevix-Trench. Identification of genetic loci associated with ovarian cancer prognosis. [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 3286. doi:10.1158/1538-7445.AM2014-3286