Huntsman, Scott

Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.

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