Hu, Donglei

Background: Multiple myeloma (MM) treatment has advanced considerably with proteasome inhibitors, immunomodulatory drugs (IMiDs), and, most recently, monoclonal antibodies. However, treatment response is quite heterogeneous, and many patients still progress even with novel combination therapies. Thus, understanding the factors that underlie response to treatment is a priority. Several somatic genomic aberrations and mutations predict poor response to therapy. However, germline variation has not previously been investigated as a predictor of treatment response in MM. We used genome-wide association and transcriptome-wide association approaches to identify germline genetic predictors of treatment response in MM. Methods: We included 510 MM patients from Mayo Clinic and 324 MM patients from UCSF diagnosed from 1999-2015. Basic demographics, laboratory values and pathology at diagnosis, type of initial therapy, duration of therapy, and follow-up were ascertained by chart review. Patients were grouped into categories of treatment based on their first line therapy: proteasome inhibitor-based regimen, IMiD based regimen, combination proteasome and IMiD based regimen, or other. Response was assessed using International Myeloma Working Group (IMWG) response criteria after 4-6 cycles of induction. As such, response was categorized into progressive disease (PD), minimal response (MR), partial response (PR), very good partial response (VGPR), or complete response (CR). Germline samples were genotyped using Illumina or Affymetrix arrays and were imputed based on the Haplotype Reference Consortium (HRC). For the genome-wide association study (GWAS), we included all SNPs with minor allele frequency >0.01 and imputation r2 of >0.5. We tested each SNP for association with treatment response using a linear regression model that adjusted for age, gender, and genetic ancestry (from principal components analysis (PCA)). To perform the transcriptome-wide association study (TWAS), we calculated predicted gene expression data using PREDIXCAN software on a reference cohort of 922 individuals with genotype and RNA expression data from peripheral blood. We then tested the association between predicted gene expression and response to therapy using linear regression models. Both the GWAS and TWAS were performed on subgroups of patients who received either proteasome inhibitors or IMiD therapies. The analyses of patients on proteasome inhibitors were adjusted for IMiD use as a covariate and analyses of patients on IMiDs were adjusted for proteasome inhibitor use. The threshold for genome-wide significance loci was set at P = 5 x 10-8 and the threshold for suggestive loci was set at P = 10-6. The threshold for significance for TWAS was set at P = 4 x 10-6 by using a Bonferroni correction for the number of genes for which genetic models of expression could be developed by PREDIXCAN. Results: Overall, 42.7% (59 of 138) of patients on proteasome inhibitors alone, 32.5% (66 of 203) of patients on IMiDs alone, and 58.1% (50 of 86) of patients on combination achieved at least a VGPR. There were no significant differences in response across centers in analyses that adjusted for age, sex and types of therapy. There were no genome wide significant loci to predict for response. We identified 8 suggestive SNPs associated with proteasome inhibitor response and 4 suggestive SNPs associated with IMiD response. TWAS identified ZNF622 as a candidate genetic modifier of proteasome inhibitor effect that was significant after correction for multiple hypothesis testing (P = 1.6 x 10-6). Higher genetically predicted expression was associated with improved response to proteasome inhibitor therapy. Among patients above the median of predicted expression of ZNF622 on proteasome inhibitors, 62.6% achieved at least VGPR; among patients at or below the median of expression, only 37.6% achieved at least VGPR. Conclusions: We identified an association between predicted expression of ZNF622 and clinical response to proteasome inhibitor therapy among MM patients. ZNF622 is a zinc finger binding protein which is known to be co-activator of B Myb activity and to affect apoptosis in response to oxidative stress. Our work highlights the potential importance of pharmacogenetic modifiers of treatment response in MM. Disclosures Shah: Nkarta: Consultancy, Membership on an entity's Board of Directors or advisory committees; Indapta Therapeutics: Equity Ownership; University of California, San Francisco: Employment; Celgene, Janssen, Bluebird Bio, Sutro Biopharma: Research Funding; Genentech, Seattle Genetics, Oncopeptides, Karoypharm, Surface Oncology, Precision biosciences GSK, Nektar, Amgen, Indapta Therapeutics, Sanofi: Membership on an entity's Board of Directors or advisory committees; Amgen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees; Poseida: Research Funding; Bristol-Myers Squibb: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees; Kite: Consultancy, Membership on an entity's Board of Directors or advisory committees; Teneobio: Consultancy, Membership on an entity's Board of Directors or advisory committees. Wong:Celgene: Research Funding; Janssen: Research Funding; Roche: Research Funding; Fortis: Research Funding. Martin:Roche and Juno: Consultancy; Amgen, Sanofi, Seattle Genetics: Research Funding. Kumar:Celgene: Consultancy, Research Funding; Takeda: Research Funding; Janssen: Consultancy, Research Funding.

Abstract Introduction The 2-3 fold excess risk of multiple myeloma (MM) among family members of cases suggests a heritable contribution to risk. Recently, a genome-wide association study (GWAS) identified two genome-wide significant and one promising novel loci associated with multiple myeloma risk. To confirm these associations and identify additional novel risk loci, we performed a four-center, genome-wide association meta-analysis. Methods A fixed effects model was used for the meta-analysis which included a total of 1248 cases and 1485 controls, all of European descent, genotyped and analyzed at four separate centers with samples contributed by 10 studies. After quality control and imputation using the 1000 Genomes Project, the analysis included ∼9.5 million variants (λ=1.024). Associations between (single nucleotide polymorphisms) SNP genotypes and MM risk were evaluated under a log-additive model of inheritance, with each study adjusting for age, sex, and up to 10 principal components to control for population stratification. Promising results were replicated in an independent set of 1587 cases and 1770 controls using TaqMan, for a total of 2835 and 3255 cases and controls, respectively, in a combined meta-analysis. Results The discovery meta-analysis did not reveal any genome-wide significant associations (defined as p<5 x 10-8). We used a novel pruning algorithm to identify the top 35 most promising single nucleotide polymorphisms (SNPs) to advance to replication. We successfully genotyped 22 SNPs in the replication set. In the combined discovery and replication meta-analysis, rs1345359 at 12q23.1 was the most strongly associated SNP (P=9 x 10-8, Table 1). The variant allele C was associated with reduced risk (odds ratio discovery set [OR]= 0.69, OR replication set = 0.78, OR combined = 0.74). A second locus at 20q13.2 (rs150220835), was associated with a two-fold increased risk (P=1.22 x 10-6), a borderline increased risk (P=0.0900) and 45% increased risk (P=2.44 X 10-5) in the discovery, replication, and combined analysis sets respectively (Table 1). We also confirmed the association between MM risk and two previously published SNPs (rs4487645, p=0.0007and rs105251, p=0.0044) (Broderick et al., Nat. Genet., 2011). The third previously suggested SNP (rs6746082) was of nominal significance (p=0.0517) in the meta-analysis. Discussion We confirmed the association between MM risk and two previously published SNPs and identified a possible association with a novel SNP in chromosome 12q23.1 (rs1345359). This SNP is not located in a gene nor associated with biofeatures in ENCODE, thus further examination of correlated SNPs is necessary to identify a functional SNP linked to this locus. We also found suggestive evidence for a second locus at 20q13.2 requiring additional replication. Larger studies would improve risk variant discovery for this rare hematologic malignancy. Disclosures: Wolf: Celgene: Honoraria, Research Funding; Millenium: Honoraria; Onyx: Honoraria. Anderson:Celgene, Millennium, BMS, Onyx: Membership on an entity’s Board of Directors or advisory committees; Acetylon, Oncopep: Scientific Founder , Scientific Founder Other.

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 &lt; 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 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