Wang, Zhaoming

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

Abstract Background: Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. Methods: The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background. Results: The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis. Conclusions: Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures. Impact: Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126–35. ©2016 AACR.

Abstract Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103–14. ©2016 AACR.

Abstract Histiocytic neoplasms are clonal, hematopoietic disorders characterized by an accumulation of abnormal, monocyte-derived dendritic cells or macrophages in Langerhans cell histiocytosis (LCH) and non-Langerhans cell histiocytosis (non-LCH), respectively. The discovery of BRAFV600E mutations in approximately 50% of these patients provided the first molecular therapeutic target in histiocytosis. However, recurrent driving mutations in the majority of patients with BRAFV600E–wild-type non-LCH are unknown, and recurrent cooperating mutations in non-MAP kinase pathways are undefined for the histiocytic neoplasms. Through combined whole-exome and transcriptome sequencing, we identified recurrent kinase fusions involving BRAF, ALK, and NTRK1, as well as recurrent, activating MAP2K1 and ARAF mutations in patients with BRAFV600E–wild-type non-LCH. In addition to MAP kinase pathway lesions, recurrently altered genes involving diverse cellular pathways were identified. Treatment of patients with MAP2K1- and ARAF-mutated non-LCH using MEK and RAF inhibitors, respectively, resulted in clinical efficacy, demonstrating the importance of detecting and targeting diverse kinase alterations in these disorders. Significance: We provide the first description of kinase fusions in systemic histiocytic neoplasms and activating ARAF and MAP2K1 mutations in non-Langerhans histiocytic neoplasms. Refractory patients with MAP2K1- and ARAF-mutant histiocytoses had clinical responses to MEK inhibition and sorafenib, respectively, highlighting the importance of comprehensive genomic analysis of these disorders. Cancer Discov; 6(2); 154–65. ©2015 AACR. This article is highlighted in the In This Issue feature, p. 109

Abstract Genome wide association studies (GWAS) have mapped multiple independent cancer risk loci (n = 6) to a small region on chr5p15.33 for at least ten distinct cancers, including bladder, breast, glioma, lung, melanoma, non-melanoma skin, ovarian, pancreas, prostate, and testicular germ cell cancer. This region harbors two plausible target genes, TERT which encodes the catalytic subunit of telomerase reverse transcriptase which maintains chromosome ends by adding telomeres repeats, and CLPTM1L which encodes the cleft lip and palate transmembrane protein 1-like protein which promotes cancer cell growth, protect cells from apoptosis, and can induce abnormal cytokinesis. The most significant pancreatic cancer GWAS SNP on chr5p15.33 was rs401681 (P = 3.7×10−7, OR = 1.19), located in the 13th intron of the CLPTM1L gene and within the independent multi-cancer risk locus in this region most distant from TERT. Imputation and fine mapping in pancreatic cancer resolved this signal by three orders of magnitude (P = 1.4×10−10, OR = 1.30) to a set of seven highly correlated SNPs. Electrophoretic mobility shift assays (EMSA) and luciferase assays across these SNPs led us to a single SNP, rs36115365, that showed allele-specific effects in both assays and overlapped with prominent ENCODE marks indicating regulatory potential. Examining this SNP across cell lines derived from multiple types of cancer linked to this locus, specifically pancreatic, testicular, lung, and melanoma, established a consistent pattern of a specific protein binding and enhanced regulatory activity for the risk allele across eight cancer cell lines (two cell lines per cancer type). siRNA-based targeting of the putative gene regulatory-region harboring rs36115365 resulted in a strong inhibition of TERT expression (average 60%, range 45-85%) but no effect on CLPTM1L expression. This inhibition of TERT expression was stronger from the pancreatic cancer risk as compared to the protective allele. Quantitative mass spectrometry identified Zinc finger protein 148 (ZNF148) as the protein specifically binding the risk allele, further supported by direct and allele-specific binding of purified recombinant ZNF148 in EMSA assays. Lastly, siRNA-mediated knockdown of ZNF148 resulted in significant reduction of TERT but not CLPTM1L expression across multiple cancer cell lines. Our results indicate that the pancreatic cancer GWAS signal at the single independent chr5p15.33 multi-cancer risk locus most distant from TERT may be explained by a single SNP that influences TERT expression via ZNF148. Citation Format: Kevin M. Brown, Jun Fang, Jinping Jia, Zhaoming Wang, Matthew Makowski, Tongwu Zhang, Jason Hoskins, Jiyeon Choi, Younghun Han, Mingfeng Zhang, Mai Xu, Peter Kanetsky, Andresson Thorkell, Gloria M. Petersen, Katherine L. Nathanson, Christopher I. Amos, Maria T. Landi, Stephen J. Chanock, Michiel Vermeulen, Laufey T. Amundadottir. Functional characterization of a multicancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4610. doi:10.1158/1538-7445.AM2015-4610

Abstract A genome-wide association study (GWAS) of prostate cancer in Kaiser Permanente health plan members (7,783 cases, 38,595 controls; 80.3% non-Hispanic white, 4.9% African-American, 7.0% East Asian, and 7.8% Latino) revealed a new independent risk indel rs4646284 at the previously identified locus 6q25.3 that replicated in PEGASUS (N = 7,539) and the Multiethnic Cohort (N = 4,679) with an overall P = 1.0 × 10−19 (OR, 1.18). Across the 6q25.3 locus, rs4646284 exhibited the strongest association with expression of SLC22A1 (P = 1.3 × 10−23) and SLC22A3 (P = 3.2 × 10−52). At the known 19q13.33 locus, rs2659124 (P = 1.3 × 10−13; OR, 1.18) nominally replicated in PEGASUS. A risk score of 105 known risk SNPs was strongly associated with prostate cancer (P < 1.0 × 10−8). Comparing the highest to lowest risk score deciles, the OR was 6.22 for non-Hispanic whites, 5.82 for Latinos, 3.77 for African-Americans, and 3.38 for East Asians. In non-Hispanic whites, the 105 risk SNPs explained approximately 7.6% of disease heritability. The entire GWAS array explained approximately 33.4% of heritability, with a 4.3-fold enrichment within DNaseI hypersensitivity sites (P = 0.004). Significance: Taken together, our findings of independent risk variants, ethnic variation in existing SNP replication, and remaining unexplained heritability have important implications for further clarifying the genetic risk of prostate cancer. Our findings also suggest that there may be much promise in evaluating understudied variation, such as indels and ethnically diverse populations. Cancer Discov; 5(8); 878–91. ©2015 AACR. This article is highlighted in the In This Issue feature, p. 783

A small number of common susceptibility loci have been identified for pancreatic cancer, one of which is marked by rs401681 in the TERT–CLPTM1L gene region on chromosome 5p15.33. Because this region is characterized by low linkage disequilibrium, we sought to identify whether additional single nucleotide polymorphisms (SNPs) could be related to pancreatic cancer risk, independently of rs401681. We performed an in‐depth analysis of genetic variability of the telomerase reverse transcriptase (TERT) and the telomerase RNA component (TERC) genes, in 5,550 subjects with pancreatic cancer and 7,585 controls from the PANcreatic Disease ReseArch (PANDoRA) and the PanScan consortia. We identified a significant association between a variant in TERT and pancreatic cancer risk (rs2853677, odds ratio = 0.85; 95% confidence interval = 0.80–0.90, p = 8.3 × 10−8). Additional analysis adjusting rs2853677 for rs401681 indicated that the two SNPs are independently associated with pancreatic cancer risk, as suggested by the low linkage disequilibrium between them (r2 = 0.07, D′ = 0.28). Three additional SNPs in TERT reached statistical significance after correction for multiple testing: rs2736100 (p = 3.0 × 10−5), rs4583925 (p = 4.0 × 10−5) and rs2735948 (p = 5.0 × 10−5). In conclusion, we confirmed that the TERT locus is associated with pancreatic cancer risk, possibly through several independent variants.

Abstract Background: Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. Relatives of DLBCL patients are at elevated risk for DLBCL. In small genome-wide association studies (GWAS) in which all NHL subtypes were combined, no conclusive loci for DLBCL have been previously identified in individuals of European background, whereas a recent study conducted in East Asia identified a locus at 3q27. Methods: To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of three new GWAS and one prior scan, totaling 3,857 DLBCL cases and 7,666 controls of European ancestry. Imputation was used to combine all the data using the 1000 Genomes Project release v3 and IMPUTE2. Based on these findings, we selected 6 promising SNPs for genotyping in an additional 202 cases and 3,431 controls, which were then included in a second meta-analysis. SNPTEST was used to estimate odds ratios (ORs), using the additive model and adjusting for age, gender and significant eigenvectors. Results: In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; OR=2.33, P=1.60x10-20), rs13255292 and rs4733601 at 8q24.21 (PVT1; OR=1.23, P=4.86x10-12 and OR=1.18, P=1.63x10-8, respectively), rs2523607 at 6p21.33 (HLA-B; OR=1.34, 3.35x10-9) and rs2681416 at 3q13.33 (CD86; OR=1.20, P=6.04x10-9). The two 8q24.21 SNPs displayed minimal linkage disequilibrium (r2=0.03) and conditional analysis supported their independence. We did not observe a notable signal for a locus previously reported for DLBCL on 3q27 in East Asia, rs6773854 (OR=1.06, P=0.81). We estimated that common SNPs, including but not limited to the loci discovered in this study, explain approximately 16% of the variance for DLBCL risk overall. Discussion: The susceptibility locus at 6p25.3 maps near a plausible DLBCL candidate gene, EXOC2 (exocyst complex component 2), which interacts with Ral proteins, and disruption of the Ral-exocyst regulatory node impacts cellular and developmental processes associated with malignant transformation and progression. The two 8q24.21 variants are approximately 1Mb telomeric to the 8q24 region linked with multiple cancers including CLL. Both variants are positioned in close proximity to PVT1, which is a non-coding RNA implicated in the MYC activation, and MYC is known to be deregulated in some DLBCLs. The locus at 3q13.33 maps to CD86 (B7-2), a B-7 family member that along with CD80 (B7-1) is essential in cognate interactions between T- and B-lymphocytes; both play key roles in immunodeficiency, autoimmune diseases, and in anti-tumor immunity. Finally, we had multiple signals in the HLA region, the strongest of which was at 6p21.33 near HLA-B. HLA-B encodes the HLA class I heavy chain paralogue, which heterodimerizes with a light chain (β2 microglobulin) and plays a central role in presenting intracellularly processed self or foreign antigens to CD8+ cytotoxic T lymphocytes. Class I molecules have been extensively linked to the risk of a variety of immune-mediated diseases and cancers including DLBCL. Conclusions: These data provide substantial new evidence for genetic susceptibility to this B-cell malignancy, and point towards pathways involved in immune recognition and immune function in the pathogenesis of DLBCL. Citation Format: James R. Cerhan, Sonja I. Berndt, Joseph Vijai, Hervé Ghesquières, James McKay, Sophia S. Wang, Zhaoming Wang, Meredith Yeager, Alexandra Nieters, David Cox, Alain Monnereau, Christopher Flowers, Anneclaire J. De Roos, Angela R. Brooks-Wilson, Qing Lan, Gianluca Severi, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark Purdue, Claire Vajdic, Demetrius Albanes, Kimberly A. Bertrand, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Paolo Vineis, Susan L. Slager, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, on behalf of the NHL GWAS Project. Genome-wide association study identifies multiple susceptibility loci for diffuse large B-cell lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr LB-272. doi:10.1158/1538-7445.AM2014-LB-272

Abstract Despite limited discovery stages (<1,125 cases), genome-wide association studies (GWAS) have successfully identified 13 loci associated with risk of chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL); however, these loci only explain a portion of the heritability. To identify additional genetic loci associated with CLL susceptibility, we conducted the largest meta-analysis of GWAS, to date, including 3,100 CLL cases and 7,667 controls of European ancestry from 25 independent studies genotyped on four different SNP microarrays. Association testing was conducted assuming a log-additive genetic model, adjusting for age, sex, and significant principal components. Meta-analyses were performed using the fixed effects inverse variance method based on the beta estimates and standard errors from each study. Promising loci were taken forward for replication in three additional studies totaling 788 cases and 4,872 controls. In the meta-analysis of both the discovery and replication studies, we conclusively identified ten independent SNPs in nine novel loci that reached genome-wide significance (P<5 x 10−8): 10q23.31 (FAS; P=1.22x10−14), 18q21.33 (BCL2; P=2.66x10−12 and P=7.76x10−11), 11p15.5 (C11orf21; P=2.15x10−10), 4q25 (LEF1; P=4.24x10−10), 2q33.1 (CASP8; P=2.50x10−9), 9p21.3 (CDKN2B-AS1; P=1.27x10−8), 18q21.32 (PMAIP1; P=2.51x10−8), 15q15.1 (BMF; P=2.71x10−10), and 2p22.2 (QPCT; P=1.68x10−8) as well as an independent signal at an established locus (2q13, ACOXL, P=2.08x10−18). We also found evidence for two additional promising loci that reached marginal genome-wide significance (P<2.0x10−7) at 8q22.3 (ODF1; P=5.40x10−8) and 5p15.33 (TERT; P=1.92x10−7). Although further studies are required, the proximity of five of the novel loci to genes involved in apoptosis suggests a possible unifying underlying biological mechanism. Together with the previously established loci, these genetic loci explain ~17% of the familial risk of CLL. Citation Format: Sonja I. Berndt, Christine F. Skibola, Vijai Joseph, Nicola J. Camp, Alexandra Nieters, Zhaoming Wang, Wendy Cozen, Alain Monnereau, Rachel S. Kelly, Silvia de Sanjose, Martha S. Linet, Elio Riboli, Paolo Vineis, Henrik Hjalgrim, James R. Cerhan, Stephen J. Chanock, Nathaniel Rothman, Susan L. Slager. Meta-analysis of genome-wide association studies identifies multiple loci associated with chronic lymphocytic leukemia. [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 LB-23. doi:10.1158/1538-7445.AM2013-LB-23