Phelan, Catherine M.

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 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. The hepatocyte growth factor (HGF) signaling pathway regulates key cellular functions related to tumor invasion and metastasis including cell motility, cell proliferation, and morphogenesis. We previously found that minor alleles at a single-nucleotide polymorphism (SNP; rs5745709) in the gene encoding the cytokine HGF was associated with poorer overall survival in women with invasive epithelial ovarian cancer (EOC) enrolled at Mayo Clinic (MAY) and The Cancer Genome Atlas (TCGA) (HR 1.6, 95% CI 1.3-1.9, p=8.9×10−5). In addition, we found that that genotype correlated with HGF expression in ovarian cancer tissue from TCGA women, suggesting that genetically-regulated expression of HGF may play a role in post-diagnosis survival. To further examine this possibility, we evaluated intensity of antibody immunostaining for HGF, MET (HGF's transmembrane tyrosine kinase receptor), and phospho-MET (pMET; MET's activated form) in tumor cells using tissue microarrays (TMAs) created from genotyped primary peritoneal and EOC cases, and correlated these expression levels with overall survival and HGF SNP genotypes. Methods. Participants were 326 white non-Hispanic women over 20 years of age with histologically-confirmed invasive EOC diagnosed between 1998 and 2009. All cases were enrolled at MAY within one year of diagnosis. Five μm sections of TMAs including three 0.6 mm cores per case were immunostained with HGF, MET, and pMET primary antibodies. Staining levels were classified according to combinations of staining extent and intensity as negative (< 10% cells stained), weak, moderate, or strong. The maximal expression value recorded over the multiple cores for each case was used. Cox proportional hazards regression was used to examine associations of overall survival with HGF, MET and pMET expression levels, adjusted for relevant clinical variables. Associations between expression levels and HGF SNP rs5745709 genotypes were examined using linear regression models. Results. HGF and MET expression were not associated with overall survival in our series of EOC cases (p=0.48 and 0.90, respectively). pMET expression was associated with survival, but in the opposite direction of that hypothesized (HR 0.77, 95% CI 0.39-1.52 for moderate expression; HR 0.55, 95% CI 0.27-1.10 for strong expression; each referent to negative/weak, p=0.01). Expression levels were not significantly associated with HGF rs5745709 genotypes (p=0.80 for HGF expression, p=0.89 for MET, p=0.41 for pMET). Conclusion. HGF expression levels were not associated with overall survival or with HGF SNP genotypes that correlated with survival in a previous study. Functional studies such as this one are crucial as a follow-up to discovery-based studies in order to assess biologic plausibility of observed results. 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 893. doi:10.1158/1538-7445.AM2011-893

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 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.

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.