Olson, Sara H.

Glioblastoma (GBM) is the most common malignant brain tumor in the United States. Incidence of GBM increases with age, and younger age‐at‐diagnosis is significantly associated with improved prognosis. While the relationship between candidate GBM risk SNPs and age‐at‐diagnosis has been explored, genome‐wide association studies (GWAS) have not previously been stratified by age. Potential age‐specific genetic effects were assessed in autosomal SNPs for GBM patients using data from four previous GWAS. Using age distribution tertiles (18–53, 54–64, 65+) datasets were analyzed using age‐stratified logistic regression to generate p values, odds ratios (OR), and 95% confidence intervals (95%CI), and then combined using meta‐analysis. There were 4,512 total GBM cases, and 10,582 controls used for analysis. Significant associations were detected at two previously identified SNPs in 7p11.2 (rs723527 [p54–63 = 1.50x10−9, OR54–63 = 1.28, 95%CI54–63 = 1.18–1.39; p64+ = 2.14x10−11, OR64+ = 1.32, 95%CI64+ = 1.21–1.43] and rs11979158 [p54–63 = 6.13x10−8, OR54–63 = 1.35, 95%CI54–63 = 1.21–1.50; p64+ = 2.18x10−10, OR64+ = 1.42, 95%CI64+ = 1.27–1.58]) but only in persons >54. There was also a significant association at the previously identified lower grade glioma (LGG) risk locus at 8q24.21 (rs55705857) in persons ages 18–53 (p18–53 = 9.30 × 10−11, OR18–53 = 1.76, 95%CI18–53 = 1.49–2.10). Within The Cancer Genome Atlas (TCGA) there was higher prevalence of ‘LGG’‐like tumor characteristics in GBM samples in those 18–53, with IDH1/2 mutation frequency of 15%, as compared to 2.1% [54–63] and 0.8% [64+] (p = 0.0005). Age‐specific differences in cancer susceptibility can provide important clues to etiology. The association of a SNP known to confer risk for IDH1/2 mutant glioma and higher prevalence of IDH1/2 mutation within younger individuals 18–53 suggests that more younger individuals may present initially with ‘secondary glioblastoma.’

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3′ gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97–1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03–1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10−28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.

Abstract Glioma incidence is highest in non-Hispanic Whites, where it occurs ~2x as frequently compared with other race/ethnicity groups. Glioma GWAS to date have included European ancestry populations only, and it is unknown whether variants identified by these analyses are associated with glioma in non- European ancestry populations. African Americans and Hispanics are admixed populations with varying proportions of European ancestry. While global ancestry may be similar within admixed groups, the proportion of European ancestry at each allele can vary across the genome. As glioma is more common in European ancestry populations, the presence of increased local European ancestry in these admixed populations could be used to identify glioma risk loci. Here we assessed whether excess European ancestry at established risk loci (Melin et al, Nature Genetics, 2017) was associated with glioma risk in non-European ancestry populations. Global ancestry was estimated using fastStructure, and local ancestry was estimated using RFMix. Both methods used 1,000 genomes project reference populations (African: YRI; European: CEU; East Asian: CHB/JPT; and Native American: CLM/PEL/MXL). We evaluated differences in local European ancestry between cases and controls using logistic regression conditioned on global European ancestry within 500kb of 25 previously identified risk variants among individuals with ≥50% African ancestry, and ≥30% Native American ancestry for all gliomas, and for grade IV glioblastoma (GBM) and grade II-III non-GBM. There were 347 individuals (184 cases and 163 controls) with ≥50% global African ancestry, and 277 individuals (153 cases and 124 controls) with ≥30% global American ancestry. There was no significant difference in proportion of global European ancestry between cases and controls with ≥50% global African ancestry (cases: 18.2%, controls: 17.7%, p=0.6834), and no significant difference in proportion of global European ancestry between cases and controls with ≥30% global American ancestry (cases: 51.1%, controls: 49.0%, p=0.2123). Among individuals with >50% African ancestry, we observed a nominally significant association between all glioma and increased local European ancestry at 7p11.2 (EGFR, pmin=0.0070) and between GBM and increased local European ancestry at 22q13.1 (CSNK1E, pmin=0.0098), both near SNPs previously associated with glioblastoma in majority European-ancestry populations. The dataset used for this analysis represents the largest collection of genotyped non-European glioma cases. These results suggest that glioma risk in African Americans may be associated with an increased local European ancestry variants at glioma risk loci previously identified in majority European ancestry populations (7p11.2 and 22q13.1). Citation Format: Quinn T. Ostrom, Kathleen M. Egan, L. Burt Nabors, Travis Gerke, Reid C. Thompson, Jeffrey J. Olson, Renato LaRocca, Sajeel Chowdhary, Jeanette E. Eckel-Passow, Georgina Armstrong, John K. Wiencke, Christopher I. Amos, Jonine L. Bernstein, Elizabeth B. Claus, Dora Il'yasova, Christoffer Johansen, Daniel H. Lachance, Rose Lai, Ryan T. Merrell, Sara H. Olson, Siegal Sadetzki, Joellen Schildkraut, Sanjay Shete, Richard S. Houlston, Robert B. Jenkins, Beatrice Melin, Melissa L. Bondy, Jill S. Barnholtz-Sloan. Evaluating glioma risk associated with extent of European admixture in African-Americans and Latinos [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 233.

Abstract Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10−6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10−7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10−3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis. Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419–30. ©2018 AACR.

Significance Much of the inherited susceptibility to pancreatic cancer remains unexplained. Germline variants that cause protein misfolding and impaired secretion of pancreatic enzymes such as CPA1 (encoding carboxypeptidase A1) can cause pancreatic acinar cell endoplasmic reticulum (ER) stress. We investigated the hypothesis that pancreatic cancer could arise from germline variants in genes encoding pancreatic secretory enzymes that induce pancreatic acinar cell stress. We find ∼1% of 1,579 patients with pancreatic cancer vs. 1 of 2,012 controls have germline variants in the genes encoding CPA1 and CPB1 (carboxypeptidase B1) that impair secretion of its protein product and induce ER stress. These findings implicate pancreatic acinar cell stress as a mechanism of pancreatic cancer susceptibility.

Abstract Objective: To elucidate the genetic architecture of gene expression in pancreatic tissues. Design: We performed expression quantitative trait locus (eQTL) and allele specific expression (ASE) analyses using RNA-sequence data and 1000 Genomes (1000G) imputed GWAS genotypes from 95 fresh frozen histologically normal pancreatic tissue samples. Data from 115 pancreatic tumor-derived tissue samples from The Cancer Genome Atlas (TCGA) was included for comparison. Results: We identified 38,615 cis-eQTLs (corresponding to 484 Genes) in histologically normal tissues and 39,713 cis-eQTL (corresponding to 237 Genes) in tumor tissues (FDR<0.1), with the strongest effects seen near transcriptional start sites (TSS). Approximately 23% and 42% of genes with significant cis-eQTLs (eGenes) appeared to be specific for tumor and normal derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in noncoding regulatory regions marked by modified histones, DNAse hypersensitivity and bound transcription factors, in particular for pancreatic tissues (1.53-3.12 fold, P≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 in the ABO gene (rs687289) was associated with ABO expression in histologically normal (P=5.8x10-8) and tumor-derived (P=8.3x10-5) pancreatic tissues. The high linkage disequilibrium (LD) between this variant and the O blood group generating deletion variant in exon 6 of ABO suggested that nonsense-mediated decay (NMD) of the “O” mRNA could explain the eQTL. However, knock-down of crucial NMD regulators did not influence decay of the ABO “O” mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL. Conclusions: We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich dataset for further studies on gene expression and regulation in pancreatic tissues. Citation Format: Laufey T. Amundadottir, Soren Lykke Andresen, Wenming Xiao, Jason Hoskins, Ashley Jermusyk, Lauren Rost, Irene Collins, Jinping Jia, Michael Mobaraki, Bin Zhu, Robert Kurtz, Hemang Parikh, Lei Song, Meredith Yeager, Torben Jensen, William Bamlet, Nilanjan Chatterjee, Brian Wolpin, Jill Smith, Sara Olson, Gloria Petersen, Jianxin Shi, Mingfeng Zhang. Analysis of cis-eQTLs in normal and tumor-derived pancreatic tissues reveals functional insights, including for the 9q34.1 ABO pancreatic cancer risk locus [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1442. doi:10.1158/1538-7445.AM2017-1442

Abstract BACKGROUND: Incidence of glioma varies significantly by sex, and most glioma histologies occur with greater incidence in males. Previous analyses have examined the impact of estrogen exposure as a risk factor for these tumors, but have found results of varying significance and low effect size. There may be differences in effect of previously discovered risk alleles that contribute to sex differences. METHODS: Using data collected for three previous glioma GWAS in European-ancestry populations (MD Anderson Cancer Center, the San Francisco Adult Glioma Study, and the Glioma International Case Control Study) we assessed sex-specific effects for 14 previously identified and 13 newly identified glioma risk SNPs (27 total) overall and for glioblastoma (GBM) and non-GBM tumors separately. There were 3,892 male cases (59% GBM), 4,522 male controls, 2,500 female cases (52% GBM) and 4,940 female controls. Sex-specific odds ratios (ORM and ORF), 95% confidence intervals (95% CI) and p values (pM and pF) were generated using stratified logistic regression models. Data from each study were analyzed separately and combined using inverse variance weighted meta-analysis. Results were considered statistically significant at p<6.2x10-4. RESULTS: In GBM, rs11979158 (7p11.2, pM=1.01x10-10, ORM=1.43 [95% CI: 1.28-1.59]; pF=3.43x10-3, ORF=1.22 [95% CI: 1.07-1.39]) and rs2562152 (16p13.3, pM=5.59x10-4, ORM=1.23 [95% CI: 1.09-1.39]; pF=1.22x10-1, ORF=1.12 [95% CI: 0.97-1.29]) had significant effect in males only. In non-GBM gliomas, rs12076373 (1q44, pM=3.56x10-7, ORM=1.41, 95% CI: 1.23-1.61; pF=1.76x10-2, ORF=1.20, 95% CI: 1.03-1.38), rs11979158 (7p11.2, pM=1.77x10-5, ORM=1.32 [95% CI: 1.16-1.49]; pF=2.73x10-1, ORF=1.08 [95% CI: 0.94-1.24]), and rs3751667 (16p13.3, pM=9.44x10-7, ORM=1.30 [95% CI: 1.17-1.44]; pF=5.18x10-2, ORF=1.13 [95% CI: 1.00-1.28]) had significant effect in males only. Effect size for rs55705857 (8q24.21) varied significantly by sex, with ORM=2.63 (95% CI: 2.24-3.09, pM=4.42x10-32), as compared to ORF=3.95 (95% CI: 3.28-4.76, pF=1.82x10-47). A sensitivity analysis was performed due to allele frequency heterogeneity by study and results did not change. CONCLUSIONS: Sex differences and other demographic differences in cancer susceptibility can provide important clues to etiology, and these differences can be leveraged for discovery in genetic association studies. Significant differences in effect size may suggest variation in genetic effect of risk alleles or in unidentified risk factors that vary in prevalence or effect by sex. There may also be differences in the distribution of molecular subtypes within each histology by sex. Further investigation using an agnostic approach may further elucidate the relationship between effect of risk alleles and sex. Citation Format: Quinn T. Ostrom, Ben Kinnersley, Margaret Wrensch, Jeanette E. Eckel-Passow, Georgina Armstrong, Terri Rice, Yanwen Chen, John Wiencke, Lucie McCoy, Helen Hansen, Christopher Amos, Jonine L. Bernstein, Elizabeth B. Claus, Dora Il'yasova, Christoffer Johansen, Daniel Lachance, Rose Lai, Ching C. Lau, Ryan T. Merrell, Sara H. Olson, Siegal Sadetzki, Joellen Schildkraut, Sanjay Shete, Richard S. Houlston, Robert B. Jenkins, Beatrice Melin, Melissa Bondy, Jill S. Barnholtz-Sloan. Estimating sex-specific effects of genetic loci associated with glioma risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1315. doi:10.1158/1538-7445.AM2017-1315

Abstract BACKGROUND: Glioma accounts for ~27% of all primary brain tumors and is responsible for ~13,000 cancer-related deaths in the US each year. Glioma tumors can be broadly classified into glioblastoma (GBM) and lower-grade non-GBM. Typically gliomas have a poor prognosis irrespective of medical care, with the most common form, GBM, having a five-year survival rate of only 5%. While genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, individual studies have had limited power to identify risk loci. METHODS: We performed the largest glioma GWAS to date, comprising a meta-analysis of six existing GWAS (6,405 cases, 14,100 controls) as well as new GWAS from the Glioma International Case Control Consortium (GICC; 4,572 cases and 3,286 controls) and University of California, San Francisco (UCSF)-Mayo (1,519 cases, 804 controls), totaling 12,496 cases (6,191 classified as GBM, 5,819 as non-GBM) and 18,190 controls. RESULTS: We identified five new risk loci for GBM at 1p31.3 (rs12752552; near JAK1, P=2.04×10-9, odds ratio (OR)=1.22), 11q14.1 (rs11233250; P=9.95×10-10, OR=1.24), 16p13.3 (rs2562152; near MPG, P=1.93x10-8, OR=1.21), 16q12.1 (rs10852606; HEATR3, P=1.29×10-11, OR=1.18), 22q13.1 (rs2235573; P=1.76×10-10, OR=1.15) and eight for non-GBM at 1q32.1 (rs4252707; MDM4, P=3.34×10-9, OR=1.19), 1q44 (rs12076373; AKT3, P=2.63×10-10, OR=1.23), 2q33.3 (rs7572263; near IDH1, P=2.18×10-10, OR=1.20), 3p14.1 (rs11706832; LRIG1, P=7.66×10-9, OR=1.15), 10q24.33 (rs11598018; OBFC1, P=3.39×10-8, OR=1.14), 11q21 (rs7107785; P=3.87×10-10, OR=1.16), 14q12 (rs10131032; P=5.07x10-11, OR=1.33) and 16p13.3 (rs3751667; P=2.61×10-9, OR=1.18). Case-only analyses confirmed the specificity of 11q14.1, 16p13.3 and 22q13.1 associations for GBM and 1q44, 2q33.3, 3p14.1, 11q21 and 14q12 for non-GBM tumors. In the combined meta-analysis, among previously published glioma risk SNPs, those for all glioma at 17p13.1 (TP53), GBM at 5p15.33 (TERT), 7p11.2 (EGFR), 9p21.3 (CDKN2B-AS1) and 20q13.33 (RTEL1) and for non-GBM at 8q24.21 (CCDC26), 11q23.2, 11q23.3 (PHLDB1) and 15q24.2 (ETFA) showed even greater evidence for association. SNPs at 10q25.2 and 12q12.1 for non-GBM tumors retained genome-wide significance (i.e. P<5.0x10-8). Associations at the previously reported loci for GBM at 3q26.2 (near TERC) and 12q23.33 (POLR3B) did not retain statistical significance. CONCLUSIONS: Our findings substantiate genetic susceptibility to GBM and non-GBM glioma being highly distinct, consistent with their distinctive molecular profiles presumably resulting from different etiological pathways. Functional analyses should lead to further insights into the biological basis of the different glioma histologies. Such information can inform gene discovery initiatives and therefore have a measurable impact on the successful development of new therapeutic agents. Citation Format: Ben Kinnersley, Beatrice S. Melin, Jill S. Barnholtz-Sloan, Margaret R. Wrensch, Christoffer Johansen, Dora Il’yasova, Quinn Ostrom, and members of GICC, Karim Labreche, Jeanette E. Eckel-Passow, Paul A. Decker, Marianne Labussière, Ahmed Idbaih, Khe Hoang-Xuan, Anna-Luisa Di Stefano, Karima Mokhtari, Jean-Yves Delattre, Pilar Galan, Konstantinos Gousias, Johannes Schramm, Minouk J. Schoemaker, Sarah J. Fleming, Stefan Herms, Stefanie Heilmann, Marcus M. Nöthen, Heinz-Erich Wichmann, Stefan Schreiber, Anthony Swerdlow, Mark Lathrop, Matthias Simon, Marc Sanson, Preetha Rajaraman, Stephen Chanock, Martha Linet, Zhaoming Wang, Meredith Yeager, Rose K. Lai, Elizabeth B. Claus, Sara H. Olson, Robert B. Jenkins, Richard S. Houlston, Melissa L. Bondy. Genome-wide association study of glioma reveals specific differences in genetic susceptibility to glioblastoma and non-glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1302. doi:10.1158/1538-7445.AM2017-1302