Iversen, Edwin S.

Abstract Background: Variation in mitochondrial genes in cancer cells alters the mitochondrial bioenergetic and biosynthetic state, which mediates cross-talk between mitochondria and the nucleus to modulate many cellular processes. The objective of this study was to comprehensively evaluate mitochondria-nuclear cross-talk in terms of SNP-SNP interactions associated with invasive epithelial ovarian cancer (EOC) risk, which may play an important role in unveiling the underlying mechanism of this complex lethal disease. Methods: Our analyses were based on 1,938 invasive EOC cases and 2,009 controls with greater than 80% European ancestry from five case-control studies genotyped using the Illumina 610quad chip. We evaluated 2,488 SNPs in over 1,500 genes in eight mitochondria-related pathways (apoptosis, oxidation-reduction, epithelial-mesenchymal transition, immunomoduation, transmembrane transport, small GTPases, BRCA1-interactors and circadian rhythm). SNPs with strong linkage disequilibrium of r2 > 0.8 were excluded. We examined 2-way SNP-SNP interactions associated with invasive EOC risk using the two-stage Random Forests plus Multivariate Adaptive Regression Splines (TRM) approach. Bootstrapping was further performed to select important variables and reduce false positive findings. Based on the scree plot using 1,000 bootstrap samples, the factors with more than 4% of bootstrap frequency were defined as significant factors. We then included each identified factors separately in model adjusting for study site and first principal component representing European ancestry. Results: Five significant SNP-SNP interactions were observed: (rs2270799 in IQSEC3 (12p13.33) + rs933305 in CUX2 (12q24.3), (p= 6.3x10−5); rs6063251 in PREX1 (20q13.3) + rs9705 in SLC39A8 (4q22-q24), (p= 3.5x10−5); rs933305 in CUX2 (12q24.3) + rs1712957 in ARHGAP15 (2q22.2-q22.3), (p=4.2x10−7), rs2861828 in SLIT1 (10q23.3-q24) + rs6063251 in PREX1 (20q13.3), (p=2.5x10−7); rs10954593 in SEMA3C (7q21-q31) + rs933305 in CUX2 (12q24.3), (p=2.8x10−3). Four of the five interactions involved SNPs from genes in the small GTPases pathway (IQSEC3; PREX1 (twice) and ARHGAP15). In addition, three of the five interaction pairs involved the CUX2 gene. This gene encodes a protein which contains three CUT domains and a homeodomain, which are DNA-binding motifs. The main mitochondrial SNP-SNP interaction involved the nuclear mitochondria gene SLC39A8. Conclusions: Our results suggest that genetic variation in mitochondrial-related pathways may be associated with invasive EOC risk. These findings highlight important genes/pathways that may interact as a network to influence invasive EOC risk and may generate novel therapeutic targets. External validation will be done using other independent studies for future analyses. Citation Format: Hui-Yi Lin, Ya-Yu Tsai, Y. Ann Chen, Zhihua Chen, Xiaotao Qu, Ellen L. Goode, Joellen Schildkraut, Edwin Iversen, Steven A. Narod, Thomas A. Sellers, Catherine M. Phelan. SNP-SNP interactions in mitochondria-related pathways are associated with invasive Epithelial Ovarian Cancer (EOC) risk . [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 4839. doi:10.1158/1538-7445.AM2013-4839

ScopeWe reevaluated previously reported associations between variants in pathways of one‐carbon (1‐C) (folate) transfer genes and ovarian carcinoma (OC) risk, and in related pathways of purine and pyrimidine metabolism, and assessed interactions with folate intake.Methods and resultsOdds ratios (OR) for 446 genetic variants were estimated among 13 410 OC cases and 22 635 controls, and among 2281 cases and 3444 controls with folate information. Following multiple testing correction, the most significant main effect associations were for dihydropyrimidine dehydrogenase (DPYD) variants rs11587873 (OR = 0.92; p = 6 × 10−5) and rs828054 (OR = 1.06; p = 1 × 10−4). Thirteen variants in the pyrimidine metabolism genes, DPYD, DPYS, PPAT, and TYMS, also interacted significantly with folate in a multivariant analysis (corrected p = 9.9 × 10−6) but collectively explained only 0.2% of OC risk. Although no other associations were significant after multiple testing correction, variants in SHMT1 in 1‐C transfer, previously reported with OC, suggested lower risk at higher folate (pinteraction = 0.03–0.006).ConclusionVariation in pyrimidine metabolism genes, particularly DPYD, which was previously reported to be associated with OC, may influence risk; however, stratification by folate intake is unlikely to modify disease risk appreciably in these women. SHMT1 SNP‐by‐folate interactions are plausible but require further validation. Polymorphisms in selected genes in purine metabolism were not associated with OC.

Abstract Background: There are several well-established environmental risk factors for ovarian cancer, and recent genome-wide association studies have also identified six variants that influence disease risk. However, the interplay between such risk factors and susceptibility loci has not been studied. Methods: Data from 14 ovarian cancer case–control studies were pooled, and stratified analyses by each environmental risk factor with tests for heterogeneity were conducted to determine the presence of interactions for all histologic subtypes. A genetic “risk score” was created to consider the effects of all six variants simultaneously. A multivariate model was fit to examine the association between all environmental risk factors and genetic risk score on ovarian cancer risk. Results: Among 7,374 controls and 5,566 cases, there was no statistical evidence of interaction between the six SNPs or genetic risk score and the environmental risk factors on ovarian cancer risk. In a main effects model, women in the highest genetic risk score quartile had a 65% increased risk of ovarian cancer compared with women in the lowest [95% confidence interval (CI), 1.48–1.84]. Analyses by histologic subtype yielded risk differences across subtype for endometriosis (Phet < 0.001), parity (Phet < 0.01), and tubal ligation (Phet = 0.041). Conclusions: The lack of interactions suggests that a multiplicative model is the best fit for these data. Under such a model, we provide a robust estimate of the effect of each risk factor that sets the stage for absolute risk prediction modeling that considers both environmental and genetic risk factors. Further research into the observed differences in risk across histologic subtype is warranted. Cancer Epidemiol Biomarkers Prev; 22(5); 880–90. ©2013 AACR.

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