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 4729: Pathway and gene set analyses for epithelial ovarian cancer (EOC) genome-wide association study (GWAS)
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 5616: Assessment of differential DNA methylation by genotype in ovarian cancer cases and controls
Abstract Ovarian carcinoma is a significant cause of cancer mortality in women in the United States. Increasing evidence has demonstrated that differential methylation status in lymphocyte DNA may be linked to risk of carcinogenesis. We are investigating the differential DNA methylation in germline DNA obtained from pre-treatment ovarian cancer cases (n=256) and individually-age matched controls (n=274) using the Illumina Infinium Human Methylation 27K BeadChip. We are first studying methylation and the common moderate-risk alleles for ovarian cancer on chromosomes 2q31, 3q25, 8q24, 9p22, 17q21, and 19p13. Within each locus, we correlated single nucleotide polymorphism (SNP) genotypes with methylation values at relevant nearby CpG sites. Preliminary results are available on the first set of 84 cases and 90 controls. At 2q31, cases with TG or TT genotype at rs2072590 were hypomethylated at HOXD1 (a homeobox transcription factor) compared to cases with GG genotype (mean beta TG/TT=0.08, GG=0.09, p=0.034); and at 17q21, cases with AG or GG genotype at rs9303542 were hypomethylated at SKAP1 (a T-cell adaptor protein) compared to cases with AA genotype (mean beta AG/GG=0.14, AA=0.16, p=0.045). At these two loci, similar trends were seen among controls. Logistic regression adjusted for age and methylation chip was then used to assess the association between SNP genotype and ovarian cancer risk, with and without inclusion of the relevant methylation value. At HOXD1 rs2072590, the per-allele ovarian cancer risk of 1.34 attenuated to 1.25 with inclusion of the HOXD1 methylation value as a covariate, suggesting that the disease association may be mediated in part through methylation. Though based on preliminary data, results are consistent with a role for inherited variation on germline methylation at HOXD1 and SKAP1. Analysis of methylation on additional genotyped cases and controls is underway. Finally, because, trans effects may also exist, we are assessing relationships between genotype and methylation at other loci throughout the genome. Such correlative molecular studies have the potential to help elucidate mechanisms underlying associations from genome-wide association studies. 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 5616. doi:10.1158/1538-7445.AM2011-5616
MicroRNA Processing and Binding Site Polymorphisms Are Not Replicated in the Ovarian Cancer Association Consortium
Abstract Mitochondria have been implicated in carcinogenesis because of their central role in apoptosis, free radical production, and cellular energy metabolism. Although studies have identified somatic mitochondrial DNA mutations in ovarian cancers, little attention has been directed to the possible role of germline mitochondrial variants in mitochondrial dysfunction and subsequent ovarian cancer development. Through an ongoing epithelial ovarian cancer genome-wide association study of 1,856 invasive cases (59% with serous carcinomas) and 1,912 controls (frequency-matched on age, race, and study site), we genotyped 138 SNPs tagging variation in the maternally-inherited mitochondrial genome (mtDNA) to evaluate the hypothesis that single nucleotide polymorphisms (SNPs) in mtDNA are associated with ovarian cancer risk. All subjects were non-Hispanic, non-Jewish Caucasians, and were genotyped with the Illumina 610-quad Beadchip. Samples and SNPs with call rates less than 95% were excluded. Subjects with ambiguous gender, unresolved identical genotypes, and less than 80% European ancestry were excluded. For each SNP, logistic regression was used to estimate odds ratios (OR) and corresponding 95% confidence intervals (CI) between carriers of the minor versus major maternally-inherited allele and case status. Subgroup analysis was conducted to estimate allele-specific risks between serous-only cases and controls. Haplotype analysis was performed for regions containing statistically significant SNPs (p<0.05). Two rare SNPs, C16329T in the displacement loop gene (MT-D-loop) (OR: 6.72, 95%CI: 1.5-29.8, p=0.004; minor allele frequency (MAF)=0.001) and C15905T in the tRNA threonine gene (MT-TT) (OR: 1.60, 95%CI: 1.1-2.3, p=0.009; MAF=0.027), were associated with increased ovarian cancer risk; whereas T6777C in the cytochrome c oxidase subunit 1 gene (MT-CO1) (OR: 0.68, 95%CI: 0.51-0.91, p=0.009; MAF=0.063) was associated with decreased risk. All three SNPs remained statistically significant after adjustment for multiple comparisons within the mitochondrial SNPs using permutation testing. C16329T and T6777C were also significantly associated among cases with serous histology, with ORs (95% CIs) of 6.13 (1.2-29.6) and 0.61 (0.43-0.87), respectively. Haplotype analysis for each of the three regions investigated revealed significant haplotype-disease associations; however, findings suggest that the respective associations were driven by the most significant SNP(s) listed above. This is the first large-scale epidemiologic study to provide evidence that SNPs in mitochondrial genes may be novel risk factors for ovarian cancer. Future investigation of additional mitochondrial variants (i.e. nuclear-encoded mitochondrial proteins) is necessary to more comprehensively investigate this association. Furthermore, replication of the most promising SNPs in a larger population is warranted to validate these findings. 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 2835.
Abstract 897: Association of inherited variants in HGF and ovarian cancer survival time
Abstract Using data from patients with invasive epithelial ovarian cancer, we sought to examine whether overall survival was associated with common inherited variation in candidate genes and genomic regions including those involved in angiogenesis, mitosis, and regions of interest from genome-wide association and expression studies. We conducted a two-stage analysis based on Mayo Clinic patients with replication of peak results in the Cancer Genome Atlas (TCGA). Ovarian cancer patients (N=312) seen at the Mayo Clinic diagnosed from 1999 to 2006 with median follow-up of 3.7 years (range, 0.1 − 9.5 years) and 172 observed deaths were analyzed at 993 single nucleotide polymorphisms (SNPs) in 179 genes. For all patients and for serous patients only (N=192 with 129 observed deaths), Cox regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for each SNP, adjusted for other prognostic factors. Twenty-nine SNPs in ten genes revealed p < 0.01 among all cases or among serous cases only, including 21 independent SNP at r2=0.9; analyses within four genes (HGF, PLG, PRKACB, and DCTN5) yielded p < 0.001. The strongest result was in HGF at rs1800793 (all cases: HR 1.7, 95% CI 1.3 − 2.2, p-trend 1.9 E-5; serous cases: HR 1.6, 95% CI 1.2 − 2.2, p-trend 1.8 E-3). A correlated SNP rs2214825 also showed association (r2=0.76; all cases: HR 1.4, 95% CI 1.1 − 1.8, p-trend 3.3 E-3; serous cases: HR 1.4, 95% CI 1.0 − 1.9, p-trend 0.03), and PCA gene-level tests of the HGF locus were also significantly associated with survival (all cases, p=3.7 E-4; serous cases, p=0.01). Among all cases, differential survival was also associated with DCTN5 rs12447304 (HR 1.8, 95% CI 1.3-2.6; p=7.1 E-4), and among serous cases, with PLG rs783173 (HR 1.6, 95% CI 1.2 − 2.0, p = 8.1 E-4) and PRKACB rs1402694 (HR 1.7, 95% CI 1.3 − 2.2, p=3.2 E-4). We then examined SNPs in the most suggestive gene, HGF, using data from TCGA which had genotyped four overlapping HGF SNPs in 350 invasive cases (188 deaths). Covariate-adjusted TCGA analysis at rs2214825 revealed consistent results (HR 1.7, 95% CI 1.3-2.2, p=1.4 E-4) as did analyses of three other HGF SNPs (HR range 1.3 to 1.5, p range 3.4 E-3 to 0.049); correlation between SNP genotypes and HGF expression was also observed (p=1.4 E-3, rs2214825, probe 210997_at). While this is the first report of SNPs in HGF, it has long been known that the HGF signaling pathway, activated by HGF ligation of c-MET, plays a key role ovarian cancer cell growth, migration, and invasion. c-MET is overexpressed in 40%-60% of ovarian tumors, and high overexpression is associated with lower overall survival. In addition, HGF has been shown to enhance the response of ovarian cancer cells to low doses of paclitaxel and cisplatin by activating the intrinsic apoptotic pathway, leading to cancer cell death. These results now suggest that inherited differences in HGF among women with invasive disease may play a role in survival, holding potential clues for tailored avenues of treatment. 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 897.
Abstract 4724: A novel region on 8q24.21 is associated with ovarian cancer susceptibility
Abstract Recent GWAS have identified and confirmed associations between several SNPs in the gene desert located on chromosome 8q24 with risk of cancers of the breast, prostate, glioma, bladder and colorectum. Evidence supporting an association with ovarian cancer is limited and inconsistent. Here we evaluate preliminary results in the 8q24 region from stage I of an ongoing invasive epithelial ovarian cancer GWAS among 1864 cases (1096 serous cases) and 1912 controls. All subjects were self-reported non-Hispanic non-Jewish Caucasian, and were genotyped with the Illumina 610quad Chip. Samples and SNPs with call rates less than 95% were excluded (9.6% and 9.9%, respectively) as well as 128 subjects with ambiguous gender, unresolved identical genotypes and less than 80% European ancestry were also excluded. Our most significant hit in the region was rs6651252 located at 8q24.21 with p = 0.000125 for all ovarian cancer, and 1.02*10−5 for serous cases only. This region is located approximately 849 kb downstream of the previous 8q24.21 findings on other cancers, and was not in linkage disequilibrium with any of the haplotype blocks reported previously. We further examined the block structure of this region with eight genotyped SNPs located within 129.61 to 129.64 MB (rs10088218, rs1516973, rs1516975, rs16903097, rs10098821, rs6651252, rs1561928, and rs1561927) and performed haplotype analysis to evaluate the association with ovarian cancer risk. One haplotype (AAGCAGGG, estimated frequency=0.11) was significantly associated with decreased ovarian cancer risk with an odds ratio of 0.67 (all cases, 95% CI: 0.78-0.91, p=0.0026, empirical p= 0.0020). This region contains no known genes, but it is located upstream and overlaps with a known copy number variation (CNV) reported by Redon et al. (2006). Previous published literature on array CGH and CNV results from TCGA (The Cancer Genome Atlas) also indicate that the 8q24.21 region is associated with frequent gains/amplification in ovarian tumors. Combined analysis of SNPs and CNV identified in this region is underway to elucidate the joint effect of these two types of genetic variation on ovarian cancer risk. Our result indicates a novel locus on 8q24.21 is associated with epithelial ovarian cancer susceptibility. Haplotype analysis reveals one specific haplotype that is associated with decreased ovarian cancer risk compared to the most common haplotype. Given previous 8q24 findings from other cancers and evidence from publicly available resources, fine-mapping and resequencing of this region is warranted. 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 4724.
Abstract 4718: Inherited variants in relapse-associated genes and risk of ovarian cancer
Abstract Background: We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk. Methods and Findings: Women with and without invasive ovarian cancer (749 cases, 1,041 controls) were genotyped at 136 single nucleotide polymorphisms (SNPs) within 13 candidate genes. Risk was estimated for each SNP and for overall variation within each gene. At the gene-level, variation within MSL1 (male-specific lethal-1 homolog) was associated with risk of serous cancer (p=0.03); haplotypes within PRPF31 (PRP31 pre-mRNA processing factor 31 homolog) were associated with risk of invasive disease (p=0.03). MSL1 rs7211770 was associated with decreased risk of serous disease (OR 0.81, 95% CI 0.66-0.98; p=0.03). SNPs in MFSD7, BTN3A3, ZNF200, PTPRS, and CCND1A were inversely associated with risk (p<0.05), and there was increased risk at HEXIM1 rs1053578 (p=0.04, OR 1.40, 95% CI 1.02-1.91). Conclusions: Tumor studies can reveal novel genes worthy of follow-up for cancer susceptibility. Here, we found that inherited markers in the gene encoding MSL1, part of a complex that modifies the histone H4, may decrease risk of invasive serous ovarian cancer. 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 4718.
Abstract 3860: New ovarian cancer susceptibility loci identified
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.