Ziv, Elad

Abstract Background: Genome-wide association studies (GWAS) in European populations have identified genetic risk variants associated with multiple myeloma. Methods: We performed association testing of common variation in eight regions in 1,318 patients with multiple myeloma and 1,480 controls of European ancestry and 1,305 patients with multiple myeloma and 7,078 controls of African ancestry and conducted a meta-analysis to localize the signals, with epigenetic annotation used to predict functionality. Results: We found that variants in 7p15.3, 17p11.2, 22q13.1 were statistically significantly (P < 0.05) associated with multiple myeloma risk in persons of African ancestry and persons of European ancestry, and the variant in 3p22.1 was associated in European ancestry only. In a combined African ancestry–European ancestry meta-analysis, variation in five regions (2p23.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) was statistically significantly associated with multiple myeloma risk. In 3p22.1, the correlated variants clustered within the gene body of ULK4. Correlated variants in 7p15.3 clustered around an enhancer at the 3′ end of the CDCA7L transcription termination site. A missense variant at 17p11.2 (rs34562254, Pro251Leu, OR, 1.32; P = 2.93 × 10−7) in TNFRSF13B encodes a lymphocyte-specific protein in the TNF receptor family that interacts with the NF-κB pathway. SNPs correlated with the index signal in 22q13.1 cluster around the promoter and enhancer regions of CBX7. Conclusions: We found that reported multiple myeloma susceptibility regions contain risk variants important across populations, supporting the use of multiple racial/ethnic groups with different underlying genetic architecture to enhance the localization and identification of putatively functional alleles. Impact: A subset of reported risk loci for multiple myeloma has consistent effects across populations and is likely to be functional. Cancer Epidemiol Biomarkers Prev; 25(12); 1609–18. ©2016 AACR.

Abstract Genome-wide association studies (GWAS) of multiple myeloma (MM) in Northern Europeans have identified seven novel risk loci (2p23.3, 3p22.1, 3q26.2, 6p21.33, 7p15.3, 17p11.2, 22q13.1). We performed a multiethnic meta-analysis of these regions in 1,274 MM patients and 1,486 controls of European ancestry (EA) and 1,049 MM patients and 7,080 controls of African ancestry (AA), leveraging the differential linkage-disequilibrium of these populations in order to better localize the putative functional variants. We observed directionally consistent effects for all seven index SNPs in both populations, with four significantly associated (p<0.05) with risk in EAs (3p22.1, 7p15.3, 17p11.2, 22q13.1), and two significantly associated with risk in AAs (7p15.3 and 22q13.1). In a fixed effects meta-analysis of six regions (excluding the HLA region on chromosome 6), variation in five of the regions (2p33.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) had statistically significant associations with risk (Table 1). In one region, the index variant had the strongest association [rs4487645 at 7p15.3, (OR = 1.30, p = 8.7×10−8)]. Five of the six most significantly associated variants identified in the multiethnic analyses overlapped with biologically relevant features indicating regulatory activity based on CD20+ (B lymphocyte) cells, showing evidence of potential function; those included a missense variant in (17p11.2, rs34562254, Pro251Leu) in TNFRSF13B, which encodes a lymphocyte-specific protein in the tumor necrosis factor receptor family that interacts with the NF-kb pathway. Our study shows that these regions are important in MM risk across ethnicities and further supports the use of multiple ethnic groups in genetic studies to enhance identification of risk variants. Table 1.Most significant associations for each region in the multiethnic meta-analysis.Individuals of European AncestryIndividuals of African AncestryMulitethnic Metar2 with IndexcCHRSNPRAaFreqbORP-valueFreqbORP-valueORP-valueP-het2rs732075G0.591.222.0×10−30.621.122.0×10−21.162.6×10−40.280.09/0.283rs73069394A0.191.243.0×10−30.621.181.5×10−21.201.3×10−50.550.77/0.963rs12637184G0.761.136.0×10−20.921.192.7×10−11.151.0×10−20.640.94/1.007rs4487645C0.671.237.0×10−40.891.485.5×10−51.308.7×10−80.07-d17rs34562254A0.121.452.4×10−50.131.212.2×10−31.312.5×10−60.120.33/0.9022rs139400T0.491.224.0×10−40.531.172.1×10−31.191.2×10−60.630.63/0.96aRisk allelebFrequency of the risk allele in European and African ancestry studiescr2 metrics based on 1000 Genomes Project AFR/EUR populationsdIndex SNP Citation Format: Kristin A. Rand, Chi Song, Eric Dean, Daniel Serie, Karen Curtin, Dennis Hazelett, Amie E. Hwang, Xin Sheng, Alex Stram, David J. Van Den Berg, Carol Ann Huff, Leon Bernal-Mizrachi, Michael H. Tomasson, Sikander Ailawadhi, Anneclaire De Roos, Seema Singhal, Karen Pawlish, Edward Peters, Catherine Bock, David V. Conti, Graham Colditz, Todd Zimmerman, Scott Huntsman, John Graff, African Ancestry Prostate Cancer GWAS Consortium,African Ancestry Breast Cancer GWAS Consortium, Stephen J. Chanock, Michael Lieber, Jayesh Mehta, Eric A. Klein, Nalini Janakiraman, Richard K. Severson, Angela R. Brooks-Wilson, Vincent Rajkumar, Elizabeth E. Brown, Laurence Kolonel, Susan Slager, Brian E. Henderson, Graham G. Giles, John J. Spinelli, Brian Chiu, Kenneth C. Anderson, Jeffrey Zonder, Robert Z. Orlowski, Sagar Lonial, Nicola Camp, Celine Vachon, Elad Ziv, Dan O. Stram, Christopher A. Haiman, Wendy Cozen. Multiple myeloma susceptibility loci examined in African and European ancestry populations. [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 4629. doi:10.1158/1538-7445.AM2015-4629

Abstract Purpose: Mammographic breast density is an established and strong risk factor for breast cancer. Recently, a number of common genetic susceptibility loci have been identified as risk factors for breast cancer. We present the first report on the relative contribution of 76 validated breast cancer susceptibility loci, in the context of a polygenic risk score (PRS), to the breast density and breast cancer association. We also examine whether the PRS improves prediction of the Breast Cancer Surveillance Consortium (BCSC) 5-year risk model above and beyond breast density and clinical factors included in the model. Methods: The study population included 1643 cases and 2397 controls from three independent epidemiologic studies: the Mayo Mammography Health Study (MMHS), the Mayo Clinic Breast Cancer Study (MCBCS), and the Bavarian Breast Cancer Cases and Control Study (BBCC). Data collected on patients in each of the studies included clinical risk factor, BI-RADS breast density [a) almost entirely fat; b) scattered fibroglandular densities; c) heterogeneously dense; and d) extremely dense] and genotypes for the 76 breast cancer susceptibility loci known at the time of the study. We formed a PRS from the reported per-SNP odds ratios (OR) for 76 known breast cancer susceptibility loci, and evaluated whether BI-RADS density and the PRS were independent risk factors for breast cancer, when adjusted for age, 1/BMI and study. We also incorporated the PRS (OR) into the BCSC 5-year risk model and estimated 5-year risk for the MMHS nested case-control study of 339 invasive cases, 765 controls. Results: BI-RADS density (p<0.0001) and PRS (p<0.0001) were independent risk factors for breast cancer that together showed greater discrimination of risk (AUC=0.69) than density (AUC=0.66; ΔAUC=0.029) or PRS score alone (AUC=0.68; ΔAUC=0.013; p<0.001). Relative to those with scattered fibroglandular densities and average (2nd quartile) PRS, women with extremely dense breasts and in the highest PRS quartile, had a 2.7 fold (95%CI: 1.7-4.1) increased risk of breast cancer, while those with fatty breasts and in the lowest PRS quartile had a reduced risk (OR=0.30, 95%CI: 0.18-0.51). Incorporation of the PRS into the BCSC risk model improved discrimination (ΔAUC=0.031, p=0.001), for a net reclassification improvement of 20% (95%CI: 11%-28%), split equally among cases (9%) and controls (11%). Conclusion: BI-RADS density and the PRS are both important risk factors for breast cancer that can be included in breast cancer risk models to improve prediction of this disease. Using these models to identify high and low-risk risk groups will facilitate improved tailored screening and primary prevention interventions. Citation Format: Celine M. Vachon, V. Shane Pankratz, Christopher G. Scott, Lothar Haeberle, Elad Ziv, Matthew R. Jensen, Kathleen R. Brandt, Dana H. Whaley, Janet E. Olson, Katharina Heusinger, Carolin C. Hack, Sebastian M. Jud, Matthias W. Beckmann, Jeffrey A. Tice, Kristen S. Purrington, Thomas A. Sellers, Karla Kerlikowske, Peter A. Fasching, Fergus J. Couch. The contribution of common breast cancer susceptibility loci to the breast density and breast cancer association and the Breast Cancer Surveillance Consortium (BCSC) risk model. [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 940. doi:10.1158/1538-7445.AM2014-940

Abstract Background: Genome wide association studies (GWAS) have recently identified 8 loci identified with multiple myeloma (MM) susceptibility. We hypothesized that germ line genetic variants may also affect MM progression. We performed the first GWAS of MM survival, by conducting a meta-analysis of two existing studies at University of California San Francisco (UCSF) and the Mayo Clinic, and replicating findings in the International Multiple Myeloma rESEarch (IMMEnSE) consortium. Methods: The UCSF study included 353 MM patients genotyped on an Illumina Omni5 and Illumina 660 array. The Mayo clinic study included 239 MM patients genotyped on an Affymetrix 6.0 array. We imputed missing genotypes using the 1000 Genomes dataset and performed a GWAS for survival using proportional hazards models adjusting for age and gender and genetic ancestry using principal components analysis. We conducted a meta-analysis of the results from the two GWAS. We replicated the top SNPs in the IMMEnSE cohort which included 772 patients with survival data from 7 European countries and one North American site. We used the dataset from Grundberg et al (Nat Genetics 2012) to test for an association between SNPs associated with survival and gene expression in immortalized lymphocytes. We used gene expression data from Zhan et al (Blood 2006) to analyze the association between expression and survival. Results: In meta-analysis, we found a genome wide significant association between SNPs on a region at 16p13 and survival (rs72773978; p=2.9x10e-10). Patients with the minor allele were at increased risk for mortality (HR 2.73; 95% CI: 1.99 - 3.73). The top associated SNPs were all in strong linkage disequilibrium and were centered over the gene FOPNL. We replicated the association in the IMMEnSE cohort and found a significant association between the top SNP from the discovery dataset (rs72773978 p=0.037). We used previously published data to determine the association between gene expression and the top associated SNPs from the GWAS. The top SNP (rs72773978) was not in previously published eQTL databases, so we used rs7201759 which was also associated with survival in our data (p=3.8x10e-10) as a proxy. The minor allele for this SNP was associated with increased levels of expression of FOPNL (p=9x10e-5), but not with expression of any other genes within a 1 megabase region. Finally, we also observed an association between increased expression of FOPNL and shorter survival (HR: 2.62 per standard deviation.; 95% CI: 1.08 - 6.32). Conclusion: Germ line variants at 16p13 near the FOPNL gene are associated with survival among MM patients. FOPNL is thought to be involved in centrosomal function, and centrosomal activity has previously been associated with survival in MM patients. Our results support this previous observation. Furthermore, they suggest that germ line variants could also be useful in assessing prognosis for patients with MM. Citation Format: Elad Ziv, Eric Dean, Donglei Hu, Alessandro Martino, Daniel Serie, Daniele Campa, Blake Aftab, Paige Bracci, Gabriele Buda, Jennifer Caswell, Charles Dumontet, Marek Dudziński, Laura Fejerman, Alexandra Greenberg, Scott Hunstman, Artur Jurczyszyn, Krzysztof Jamroziak, Shaji Kumar, Herlander Marques, Thomas Martin, Joaquin Martinez-Lopez, Vincent Rajkumar, Juan Sainz, Annette Juul Vangsted, Marzena Watek, Jeffrey Wolf, Susan Slager, Federico Canzian, Celine Vachon. Genome wide association study identifies variants at 16p13 associated with survival in multiple myeloma patients. [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 5078. doi:10.1158/1538-7445.AM2014-5078