Fenstermacher, David

Background The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study. Methods We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant. Results For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10−5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10−8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants. Conclusions This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

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 DNA copy number variations (CNVs) are a significant and ubiquitous source of human genetic variation. However, the influence of CNVs on cancer susceptibility remains poorly understood. Out of 83 genome-wide association studies (GWAS) on cancer to date, only a few studies found significant associations on germline CNVs and cancer. Here we analyzed data from our ongoing two-stage GWAS of four North American case-control studies to test the hypothesis that CNVs in germline DNA from peripheral blood lymphocytes may serve as risk factors for epithelial ovarian cancer (EOC). To reduce disease heterogeneity, we focused on the subset of cases with serous histology. The analysis was therefore based on 942 serous ovarian cancer patients and 1,682 healthy controls who were genotyped using the Illumina 610K quad array. Subjects with extreme noise and genomic waviness in log R ratio (LR) were excluded prior to segmentation; principal component analysis was performed to adjust for batch effects. CNV segmentation was performed on LR data from 22 autosomes using circular binary segmentation embedded in the Copy Number Analysis Module from Golden Helix SNP Variation Suite version 7. Copy number segment covariates were discretized based on the thresholds that signify a transition between copy number states (deletion/no deletion; duplication/no duplication) after the segmentation. Unconditional logistic regression on a log-additive model was used to evaluate the association between copy number states and serous EOC risk after adjusting for study sites. By comparing single marker copy number states at 388,958 SNPs on 22 autosomes, we observed a total of 134 SNPs significantly associated with risk of serous EOC with a p value below 10−6. Associations with deletion polymorphisms were observed on chromosomes 7, 8, 14, and 18 when controlling for false discovery rate at 1%; no duplication polymorphisms were significant. We observed a large deletion at chromosome 14 that occurred in 8.9% of cases but in only 3.9% of controls, with a p value of 5.59×10−8 (Odds Ratio (OR): 2.56, 95% confidence interval (CI): 1.82-3.60). Another common deletion on chromosome 7 occurred in 8.8% of cases and only 4.3% of controls (p=3.83×10−6; OR: 2.23, 95% CI: 1.60-3.10). Two additional regions on chromosomes 8 and 18 with deletion events less than 5% were also identified. Women who harbored the deletion on chromosome 8 were at lower risk of developing serous ovarian cancer (p=4.83×10−8; OR: 0.04, 95% CI: 0.01-0.30). Women who carried the deletion on chromosome 18 had an increased risk for serous EOC (p=9.04×10−7; OR: 7.26, 95% CI: 2.94-17.97). Further validation of these four regions using independent data sets is currently underway. In summary, this is the largest reported genome-wide study of CNVs and serous EOC risk. These preliminary results suggested that germline CNVs may play an important role in ovarian cancer susceptibility. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 1649. doi:1538-7445.AM2012-1649

Abstract Mi(cro)RNAs are short non-coding RNA molecules that play a key role in carcinogenesis by regulating tumor suppressors and oncogenes. Published data has implicated miRNAs in ovarian cancer (OC) development and progression, and we hypothesize that single nucleotide polymorphisms (SNPs) in sites of miRNA: messenger RNA (mRNA) binding may influence OC risk by altering expression levels of targeted mRNAs. To evaluate associations between SNPs in miRNA binding sites and OC risk, we used bioinformatics tools and public databases to identify miRNA binding SNPs relevant to ovarian cancer. We then evaluated the frequency of approximately 1,100 SNPs in ∼20,000 epithelial OC cases and 20,000 controls represented in the international Collaborative Oncological Gene-environment Study. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI) between genotypes and case status, with adjustment for the first five principal components representing European ancestry. Log-additive genetic models were applied to each SNP, modeling the number of copies of the minor allele. Subgroup analysis was conducted for serous adenocarcinomas, the most predominant histologic subtype of epithelial OC. Preliminary analysis of 11,215 invasive epithelial OCs and 15,910 controls included four SNPs with P<10−6 and 1 SNP with P<10−5. Analysis of 6,215 invasive serous adenocarcinomas revealed two SNPs with P<10−7 and three SNPs with P<10−6. Noteworthy is rs6104808, a C>T SNP predicted to reside in a highly conserved miRNA binding site of BTBD3 (BTB POZ domain containing 3), a transcription factor suggested to be involved in cell proliferation, cancer progression, and response to platinum-based chemotherapy. rs6104808 was associated with a decreased risk of OC overall (OR (95% CI): 0.42 (0.28-0.62), P=1.26 x 10−5)), and the association was stronger among cases with serous adenocarcinoma (OR (95% CI): 0.30 (0.17-0.50), P=9.89 x 10−6)). This represents the largest, most comprehensive epidemiologic study to date to evaluate associations between miRNA binding site SNPs and OC susceptibility. Although based on a preliminary analysis, these data suggest a significant role and putative biological mechanism for common germline variants in miRNA binding sites leading to cancer development. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2927. doi:1538-7445.AM2012-2927

Abstract Epithelial ovarian cancer (EOC) is the 6th most commonly diagnosed cancer among women in the world, accounting for 4% of all female cancers. On a worldwide basis, an estimated 204,000 new cases are diagnosed and 125,000 women die of EOC annually, with 5-year survival rate at 45%. Metastasis and its complications are the usual cause of EOC-related deaths. Thus, factors related to tumor aggressiveness, response to therapy, and the underlying health state of the patient often predict survival. Candidate gene and genome-wide association studies (GWAS) for survival from EOC have discovered only a few the genetic variants, accounting for only a small proportion of the heritability. One explanation for the missing heritability is that the common analysis approach _ assessing the effect of each single nucleotide polymorphism (SNP) individually _ is not well suited to the detection of small effects of multiple SNPs. Gene set analysis (GSA) overcomes this limitation by assessing the overall evidence of association of a phenotype with all measured variation in a set of genes or pathway. Therefore, we conducted a GSA using data from an EOC GWAS, using the Illumina 610-quad SNP array, to determine novel gene sets (GSs) associated with overall survival (OS) from EOC (N cases = 1057, N deaths = 495). All subjects were non-Hispanic, non-Jewish Caucasians, with sample and SNP call rates >95%. Subjects with ambiguous gender, unresolved identical genotypes, and less than 80% European ancestry were also excluded. SNPs were mapped to genes within 20 kb, with genes then mapped to GSs defined by the following: KEGG, Gene Ontology (biological, cellular, molecular) and PharmGKB. GSA was completed in which the association of each gene with OS was assessed using a principal component gene-level test (adjusting for study site, age and population sub-structure based on the first eigenvector from EigenStrat), followed by the aggregation of the gene-level p-values using the Gamma method with a soft truncation threshold of 0.15. Empirical p-values were determined using permutation methods. Of the 2,573 GSs assessed, 30 GSs had a p-value < 0.01. The top GSs were related to: mismatch repair (p = 3.0E-4), platelet formation (p=0.001), leukocyte adhesive activation (p=0.002), antigen processing and presentation of endogenous antigens (p = 0.002) or peptide antigens (p = 0.002), platinum pathway (p = 0.003), aging and lifespan (p = 0.004, p = 0.007), immune system development (p = 0.004) and somatic diversification (p=0.007). These results illustrate the usefulness of GSA for the development of novel hypotheses that can be followed-up with additional clinical or functional studies. Moreover, application of GSA to survival from EOC allowed the detection of significant GSs (and the corresponding genes/SNPs within these GSs) that impact OS that were not detected by assessing each SNP individually, where the smallest p-value from single SNP analysis of the GWAS for OS was 1.3E-05. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 2738. doi:10.1158/1538-7445.AM2011-2738

Abstract The etiology of ovarian cancer is poorly understood but there is clearly a heritable component. Efforts to identify susceptibility alleles rarely consider interactions among alleles or their joint effects because it quickly becomes computationally intractable. In this study, we sought to identify multi-SNP effects jointly with a pathway-based analysis (GSEA-SNP) of 1,952 EOC cases and 2,042 frequency-matched controls genotyped with the Illumina 610K array. All subjects were self-reported non-Hispanic non-Jewish Caucasians, with SNP and sample call rates > 95%. Subjects with ambiguous gender, unresolved identical genotypes and < 80% European ancestry were excluded. SNPs with MAF < 1% were excluded. Missing genotypes were inferred using Mach based on the HapMap CEU population. SNPs not within introns were annotated to genes within 100 bp. We retrieved the following databases of gene sets (GSs) from a compiled database, MsigDB: human chromosome and cytogenetic band, chemical and genetic perturbations, canonical pathways, microRNA binding targets, transcription factor targets (TFT), cancer gene neighborhood (CGN), Cancer modules and Gene Ontology. The GSEA-SNP approach calculates and rank orders the trend statistic for association between each SNP and EOC risk. The Enrichment Score (ES) estimates the overrepresentation of top-ranked SNPs for each GS; the statistical significance was estimated using 10,000 permutations. The ES was normalized according to the size of GS to yield Normalized Enrichment Score (NES). The false discovery (FDR) rate was estimated using NES to adjust for multiple hypothesis testing. A total of 5181 gene-sets were included in the analysis. When controlling the FDR at 15%, 14 of the GSs were highly enriched with association signals, including two chromosomal regions, 8q13 (p = 0.0055, FDR = 0.15) and 6p24 (p = 0.0012, FDR = 0.08). BIOCARTA_LYM_PATHWAY, a pathway related to cell adhesion and diapedesis of lymphocytes, was also significant (p = 0.0002, FDR = 0.12). A TFT GS, composed of genes with promoter regions around transcription start sites containing the motif GGCNRNWCTTYS, was associated with risk (p = 0.0004, FDR = 0.08). Currently, no known transcription factors bind to this computationally predicted motif. Molecular function of double stranded RNA binding was also enriched (p = 0.005, FDR = 0.10). The remaining enriched sets were CGN sets, for which the neighborhoods around the cancer genes were originally defined using correlation of gene expression from 4 large data sets of various cancer types. The enriched CGN included: CD48 (p = 0.011), CD53 (p = 0.009), CD97 (p = 0.011), INPP5D (p = 0.010), PTPN6 (p = 0.008), VAV1 (p = 0.011), ITGAL (p = 0.012), PTPRC (p = 0.012), and STAT6 (p = 0.010). In summary, these analyses detected biologically plausible GSs related to etiology of EOC, highlighting SNPs in core enrichment groups that were not identified using individual SNP tests. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4729. doi:10.1158/1538-7445.AM2011-4729

Abstract Background: Genome-wide association studies (GWAS) for epithelial ovarian cancer (EOC), the most lethal gynecologic malignancy, have identified novel susceptibility loci. GWAS for survival after EOC have had more limited success. The association of each single-nucleotide polymorphism (SNP) individually may not be well suited to detect small effects of multiple SNPs, such as those operating within the same biologic pathway. Gene set analysis (GSA) overcomes this limitation by assessing overall evidence for association of a phenotype with all measured variation in a set of genes. Methods: To determine gene sets associated with EOC overall survival, we conducted GSA using data from two large GWAS (N cases = 2,813, N deaths = 1,116), with a novel Principal Component-Gamma GSA method. Analysis was completed for all cases and then separately for high-grade serous histologic subtype. Results: Analysis of the high-grade serous subjects resulted in 43 gene sets with P < 0.005 (1.7%); of these, 21 gene sets had P < 0.10 in both GWAS, including intracellular signaling pathway (P = 7.3 × 10−5) and macrolide binding (P = 6.2 × 10−4) gene sets. The top gene sets in analysis of all cases were meiotic mismatch repair (P = 6.3 × 10−4) and macrolide binding (P = 1.0 × 10−3). Of 18 gene sets with P < 0.005 (0.7%), eight had P < 0.10 in both GWAS. Conclusion: This research detected novel gene sets associated with EOC survival. Impact: Novel gene sets associated with EOC survival might lead to new insights and avenues for development of novel therapies for EOC and pharmacogenomic studies. Cancer Epidemiol Biomarkers Prev; 21(3); 529–36. ©2012 AACR.

Abstract Background: Single nucleotide polymorphisms (SNP) in microRNA-related genes have been associated with epithelial ovarian cancer (EOC) risk in two reports, yet associated alleles may be inconsistent across studies. Methods: We conducted a pooled analysis of previously identified SNPs by combining genotype data from 3,973 invasive EOC cases and 3,276 controls from the Ovarian Cancer Association Consortium. We also conducted imputation to obtain dense coverage of genes and comparable genotype data for all studies. In total, 226 SNPs within 15 kb of 4 miRNA biogenesis genes (DDX20, DROSHA, GEMIN4, and XPO5) and 23 SNPs located within putative miRNA binding sites of 6 genes (CAV1, COL18A1, E2F2, IL1R1, KRAS, and UGT2A3) were genotyped or imputed and analyzed in the entire dataset. Results: After adjustment for European ancestry, no overall association was observed between any of the analyzed SNPs and EOC risk. Conclusions: Common variants in these evaluated genes do not seem to be strongly associated with EOC risk. Impact: This analysis suggests earlier associations between EOC risk and SNPs in these genes may have been chance findings, possibly confounded by population admixture. To more adequately evaluate the relationship between genetic variants and cancer risk, large sample sizes are needed, adjustment for population stratification should be carried out, and use of imputed SNP data should be considered. Cancer Epidemiol Biomarkers Prev; 20(8); 1793–7. ©2011 AACR.

Abstract Background: Mitochondria contribute to oxidative stress, a phenomenon implicated in ovarian carcinogenesis. We hypothesized that inherited variants in mitochondrial-related genes influence epithelial ovarian cancer (EOC) susceptibility. Methods: Through a multicenter study of 1,815 Caucasian EOC cases and 1,900 controls, we investigated associations between EOC risk and 128 single nucleotide polymorphisms (SNPs) from 22 genes/regions within the mitochondrial genome (mtDNA) and 2,839 nuclear-encoded SNPs localized to 138 genes involved in mitochondrial biogenesis (BIO, n = 35), steroid hormone metabolism (HOR, n = 13), and oxidative phosphorylation (OXP, n = 90) pathways. Unconditional logistic regression was used to estimate OR and 95% CI between genotype and case status. Overall significance of each gene and pathway was evaluated by using Fisher's method to combine SNP-level evidence. At the SNP level, we investigated whether lifetime ovulation, hormone replacement therapy (HRT), and cigarette smoking were confounders or modifiers of associations. Results: Interindividual variation involving BIO was most strongly associated with EOC risk (empirical P = 0.050), especially for NRF1, MTERF, PPARGC1A, ESRRA, and CAMK2D. Several SNP-level associations strengthened after adjustment for nongenetic factors, particularly for MTERF. Statistical interactions with cigarette smoking and HRT use were observed with MTERF and CAMK2D SNPs, respectively. Overall variation within mtDNA, HOR, and OXP was not statistically significant (empirical P > 0.10). Conclusion: We provide novel evidence to suggest that variants in mitochondrial biogenesis genes may influence EOC susceptibility. Impact: A deeper understanding of the complex mechanisms implicated in mitochondrial biogenesis and oxidative stress may aid in developing strategies to reduce morbidity and mortality from EOC. Cancer Epidemiol Biomarkers Prev; 20(6); 1131–45. ©2011 AACR.

Abstract Defective microRNA (miRNA) biogenesis contributes to the development and progression of epithelial ovarian cancer (EOC). In this study, we examined the hypothesis that single nucleotide polymorphisms (SNP) in miRNA biogenesis genes may influence EOC risk. In an initial investigation, 318 SNPs in 18 genes were evaluated among 1,815 EOC cases and 1,900 controls, followed up by a replicative joint meta-analysis of data from an additional 2,172 cases and 3,052 controls. Of 23 SNPs from 9 genes associated with risk (empirical P < 0.05) in the initial investigation, the meta-analysis replicated 6 SNPs from the DROSHA, FMR1, LIN28, and LIN28B genes, including rs12194974 (G>A), an SNP in a putative transcription factor binding site in the LIN28B promoter region (summary OR = 0.90, 95% CI: 0.82–0.98; P = 0.015) which has been recently implicated in age of menarche and other phenotypes. Consistent with reports that LIN28B overexpression in EOC contributes to tumorigenesis by repressing tumor suppressor let-7 expression, we provide data from luciferase reporter assays and quantitative RT-PCR to suggest that the inverse association among rs12194974 A allele carriers may be because of reduced LIN28B expression. Our findings suggest that variants in LIN28B and possibly other miRNA biogenesis genes may influence EOC susceptibility. Cancer Res; 71(11); 3896–903. ©2011 AACR.