Birrer, Michael J.

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