Sutphen, Rebecca

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 As pathways of the immune system play an important role in multiple malignancies, we assessed the association between single-nucleotide polymorphisms (SNPs) in inflammation-related genes and risk of ovarian cancer in a multi-site case-control study. Using a custom array, we genotyped 930 epithelial ovarian cancer cases and 1,037 controls for SNPs in ALOX12, ALOX15, ALOX5, CCL11, CCL2, CCL3, CCR3, CRP, CXCL16, IL10, IL15RA, IL18, IL1A, IL1B, IL1RN, IL4R, IL6, IL6R, IL7R, IL8RA, IL8RB, IL9, NOS3, PTGS1, PTGS2, TLR2, and TNF. After analysis by logistic regression, SNPs with p<0.10 were evaluated in an independent data set of 3,143 cases and 2,102 controls from the Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI) collaboration. Combined analysis revealed association with interleukin 1, alpha (IL1A) SNPs (rs17561 and rs4848300, r2=0.990, p=0.005) which varied by histological subtype (heterogeneity p=0.03). For example, IL1A rs17561, which results in a missense change in amino acid at position 114 and correlates with multiple inflammatory phenotypes, was inversely associated with risk of clear cell, mucinous, and endometrioid subtypes, but not with the serous subtype (Table 1). Additional data on IL1A rs17561 in an even larger collection of ovarian cancer cases and controls are pending and will be analyzed soon. We also found that genotype at rs1864414 in ALOX5 (arachidonate 5-lipoxygenase) was associated with decreased risk combined across subtypes (OR 0.86, 95% CI 0.78-0.95; p=0.002). Thus, inherited variation in IL1A and ALOX5 appears to have a role in ovarian cancer risk, though it is limited to rarer subtypes for IL1A. These data add to existing literature on the importance of inflammation in tumorigenesis and growing evidence of subtype-specific traits in ovarian cancer and may provide important clues on prevention. 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 2603. doi:1538-7445.AM2012-2603

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: Overexpression of mitotic kinases has been associated with prognosis, histologic grade, and clinical stage in ovarian cancer, but the relationship between inherited variation in these genes and ovarian cancer risk has not been well defined. Methods: We measured associations between 397 single nucleotide polymorphisms (SNPs) from 67 mitotic kinases and invasive epithelial ovarian cancer risk in two case–control studies (n = 671 cases; n = 939 controls). Thirty-six candidate SNPs (P < 0.05) were assessed in a replication analysis consisting of three additional studies (n = 1,094 cases; n = 829 controls). Results: In initial analysis, thirty-six SNPs were suggestive of association with risk of serous ovarian cancer, all subtypes of ovarian cancer, or both (P < 0.05). Replication analyses suggested an association between rs2125846 in the Nemo-like kinase (NLK) gene and ovarian cancer (serous OR = 1.36, 95% CI: 1.11–1.67, P = 1.77 × 10−3; all subtypes OR = 1.30, 95% CI: 1.08–1.56, P = 2.97 × 10−3). Furthermore, rs2125846 was associated with risk in the combined discovery and replication sets (serous OR = 1.33, 95% CI: 1.15–1.54; all subtypes OR = 1.27, 95% CI: 1.12–1.45). Conclusions: Variation in NLK may be associated with risk of invasive epithelial ovarian cancer. Further studies are needed to confirm and understand the biologic relationship between this mitotic kinase and ovarian cancer risk. Impact: An association between SNPs in NLK and ovarian cancer may provide biologic insight into the development of this disease. Cancer Epidemiol Biomarkers Prev; 21(3); 523–8. ©2012 AACR.

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 The importance of inflammation pathways to the development of many human cancers prompted us to examine the associations between single-nucleotide polymorphisms (SNP) in inflammation-related genes and risk of ovarian cancer. In a multisite case–control study, we genotyped SNPs in a large panel of inflammatory genes in 930 epithelial ovarian cancer cases and 1,037 controls using a custom array and analyzed by logistic regression. SNPs with P < 0.10 were evaluated among 3,143 cases and 2,102 controls from the Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI) post-GWAS collaboration. Combined analysis revealed association with SNPs rs17561 and rs4848300 in the interleukin gene IL1A which varied by histologic subtype (Pheterogeneity = 0.03). For example, IL1A rs17561, which correlates with numerous inflammatory phenotypes, was associated with decreased risk of clear cell, mucinous, and endometrioid subtype, but not with the most common serous subtype. Genotype at rs1864414 in the arachidonate 5-lipoxygenase ALOX5 was also associated with decreased risk. Thus, inherited variation in IL1A and ALOX5 seems to affect ovarian cancer risk which, for IL1A, is limited to rarer subtypes. Given the importance of inflammation in tumorigenesis and growing evidence of subtype-specific features in ovarian cancer, functional investigations will be important to help clarify the importance of inherited variation related to inflammation in ovarian carcinogenesis. Cancer Res; 72(5); 1064–9. ©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: Invasive ovarian cancer is a significant cause of gynecologic cancer mortality. Methods: We examined whether this mortality was associated with inherited variation in approximately 170 candidate genes/regions [993 single-nucleotide polymorphisms (SNPs)] in a multistage analysis based initially on 312 Mayo Clinic cases (172 deaths). Additional analyses used The Cancer Genome Atlas (TCGA; 127 cases, 62 deaths). For the most compelling gene, we immunostained Mayo Clinic tissue microarrays (TMA, 326 cases) and conducted consortium-based SNP replication analysis (2,560 cases, 1,046 deaths). Results: The strongest initial mortality association was in HGF (hepatocyte growth factor) at rs1800793 (HR = 1.7, 95% CI = 1.3–2.2, P = 2.0 × 10−5) and with overall variation in HGF (gene-level test, P = 3.7 × 10−4). Analysis of TCGA data revealed consistent associations [e.g., rs5745709 (r2 = 0.96 with rs1800793): TCGA HR = 2.4, CI = 1.4–4.1, P = 2.2 × 10−3; Mayo Clinic + TCGA HR = 1.6, CI = 1.3–1.9, P = 7.0 × 10−5] and suggested genotype correlation with reduced HGF mRNA levels (P = 0.01). In Mayo Clinic TMAs, protein levels of HGF, its receptor MET (C-MET), and phospho-MET were not associated with genotype and did not serve as an intermediate phenotype; however, phospho-MET was associated with reduced mortality (P = 0.01) likely due to higher expression in early-stage disease. In eight additional ovarian cancer case series, HGF rs5745709 was not associated with mortality (HR = 1.0, CI = 0.9–1.1, P = 0.87). Conclusions: We conclude that although HGF signaling is critical to migration, invasion, and apoptosis, it is unlikely that HGF genetic variation plays a major role in ovarian cancer mortality. Furthermore, any minor role is not related to genetically-determined expression. Impact: Our study shows the utility of multiple data types and multiple data sets in observational studies. Cancer Epidemiol Biomarkers Prev; 20(8); 1638–48. ©2011 AACR.