Spinelli, John J.

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10−9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10−4−5.79 × 10−14). Mechanistically, we found that the ULK4rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10−4), whereas the IKBKErs17433804 SNP correlated with the number of transitional CD24+CD38+ B cells (p = 4.8 × 10−4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10−4). We also found that the CD46rs1142469 SNP correlated with numbers of CD19+ B cells, CD19+CD3− B cells, CD5+IgD− cells, IgM− cells, IgD−IgM− cells, and CD4−CD8− PBMCs (p = 4.9 × 10−4−8.6 × 10−4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2Ars2811710 SNP correlated with levels of CD4+EMCD45RO+CD27− cells (p = 9.3 × 10−4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3−, MCP-2−, and IL20-dependent pathways.

AbstractGene expression profiling can be used for predicting survival in multiple myeloma (MM) and identifying patients who will benefit from particular types of therapy. Some germline single nucleotide polymorphisms (SNPs) act as expression quantitative trait loci (eQTLs) showing strong associations with gene expression levels. We performed an association study to test whether eQTLs of genes reported to be associated with prognosis of MM patients are directly associated with measures of adverse outcome. Using the genotype‐tissue expression portal, we identified a total of 16 candidate genes with at least one eQTL SNP associated with their expression with P < 10−7 either in EBV‐transformed B‐lymphocytes or whole blood. We genotyped the resulting 22 SNPs in 1327 MM cases from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and examined their association with overall survival (OS) and progression‐free survival (PFS), adjusting for age, sex, country of origin and disease stage. Three polymorphisms in two genes (TBRG4‐rs1992292, TBRG4‐rs2287535 and ENTPD1‐rs2153913) showed associations with OS at P < .05, with the former two also associated with PFS. The associations of two polymorphisms in TBRG4 with OS were replicated in 1277 MM cases from the International Lymphoma Epidemiology (InterLymph) Consortium. A meta‐analysis of the data from IMMEnSE and InterLymph (2579 cases) showed that TBRG4‐rs1992292 is associated with OS (hazard ratio = 1.14, 95% confidence interval 1.04‐1.26, P = .007). In conclusion, we found biologically a plausible association between a SNP in TBRG4 and OS of MM patients.

Abstract So far, 23 germline susceptibility loci have been associated with multiple myeloma (MM) risk. It is unclear whether the genetic variation associated with MM susceptibility also predisposes to its precursor, monoclonal gammopathy of undetermined significance (MGUS). Leveraging 2434 MM cases, 754 MGUS cases, and 2 independent sets of controls (2567/879), we investigated potential shared genetic susceptibility of MM and MGUS by (1) performing MM and MGUS genome-wide association studies (GWAS); (2) validating the association of a polygenic risk score (PRS) based on 23 established MM loci (MM-PRS) with risk of MM, and for the first time with MGUS; and (3) examining genetic correlation of MM and MGUS. Heritability and genetic estimates yielded 17% (standard error [SE] ±0.04) and 15% (SE ±0.11) for MM and MGUS risk, respectively, and a 55% (SE ±0.30) genetic correlation. The MM-PRS was associated with risk of MM when assessed continuously (odds ratio [OR], 1.17 per SD; 95% confidence interval [CI], 1.13-1.21) or categorically (OR, 1.70; 95% CI, 1.38-2.09 for highest; OR, 0.71; 95% CI, 0.55-0.90 for lowest compared with middle quintile). The MM-PRS was similarly associated with MGUS (OR, 1.19 per SD; 95% CI, 1.14-1.26 as a continuous measure, OR, 1.77, 95%CI: 1.29-2.43 for highest and OR, 0.70, 95%CI: 0.50-0.98 for lowest compared with middle quintile). MM and MGUS associations did not differ by age, sex, or MM immunoglobulin isotype. We validated a 23-SNP MM-PRS in an independent series of MM cases and provide evidence for its association with MGUS. Our results suggest shared common genetic susceptibility to MM and MGUS.

Abstract We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.

Background Genome-wide association studies (GWAS) conducted in populations of European ancestry (EA) have identified and confirmed 23 germline susceptibility loci for multiple myeloma (MM). The effect sizes of single nucleotide polymorphisms (SNPs) at these loci are small, therefore combining them into a single summary measure, known as a polygenic risk score (PRS), may provide a more meaningful risk factor. We have previously shown a PRS comprised of the 23 SNPs for MM contributes to increased risk of MM, with a 2.7-fold increase for highest vs. lowest PRS quintiles. Whether the MM-PRS is also associated with overall survival (OS) in MM cases has not been evaluated. We examined the association between MM-PRS and OS in two EA studies. Methods The first study consisted of 2,179 EA MM cases from ten studies included in the Multiple Myeloma Working Group within the International Lymphoma Consortium (InterLymph). Cases were diagnosed between 1970 and 2015 and genotyped using multiple platforms (Oncoarray, Affymetrix, Human660W-quad Beadchip, and Illumina arrays); 885 cases also had stage [based on International Staging System (ISS)] available. Each of the GWAS was subjected to rigorous standard quality control independently (prior to imputation via the Michigan imputation server based on the Haplotype Reference Consortium (HRC). The second study consisted of 515 newly diagnosed EA MM cases from CoMMpass (Relating Clinical Outcomes in Multiple Myeloma to Personal Assessment of Genetic Profile), diagnosed from 2011-2013, who had whole genome sequencing (WGS) performed on germline DNA. The WGS data was used to call common germline genetic variants through the Mayo Clinic bioinformatics pipeline. Briefly, genetic variants were detected with GenomeGPS, aligned to the hg19 reference genome, called using the GATK (V3.6) Haplotype Caller, and merged for multiple-sample joint calling. To reduce the false positive variants, variant quality score recalibration (VQSR) was applied for both SNPs and INDELs. After quality control, 458 EA samples remained. Follow-up was available for both studies and consisted of time from MM diagnosis date until death or date of last known follow-up. The PRS was constructed from the 23 MM SNPs using the published per allele odds ratio associated with MM risk. The published log odds ratios for each SNP were multiplied by the number of risk alleles (0, 1, 2) for the corresponding SNP, and summed, resulting in a unique score per person. Kaplan-Meier curves and Cox proportional hazard models were used to assess the association between PRS with MM OS considering two models: 1) adjusted for age, sex, study and 2) additional adjustment by stage (ISS). Hazard ratios (OR) and 95% confidence intervals (CI) were estimated. The PRS was evaluated both as a continuous variable, per standard deviation (SD), and as a categorical variable (quintiles). Results MM cases (N=2,179) in the InterLymph study were 59% male and 41% female and the median age was 61.0 years (26-90 years). Median follow-up time was 57.2 months (1.0-509.0 months) with 868 reported deaths. MM cases with stage information available consisted of 20% stage I (n=178), 53% stage II (n=466), and 27% stage III (n=241). No association was observed between PRS and OS in MM patients regardless of adjustment for stage (continuous PRS (HR: 1.03, 95% CI: 0.83-1.28, P=0.80) or by quintile PRS (p>0.05)) (Table). The CoMMpass EA MM cases (n=458) had similar distributions for sex (61% male and 39% females) but were slightly older 65 years (27-93 years) and had shorter follow-up time (median=39.75 months (0.13-77.2)) with 117 deaths. Stage was available for 96% of CoMMpass cases including 36% stage I (n=159), 33% stage II (n=146), and 31% stage III (n=134). We also observed no association of PRS and OS in the CoMMpass study (HR=1.02, 95% CI: 0.72 -1.46, P= 0.89), adjusted for age, sex, and stage (Table). Discussion A PRS score for MM risk is not associated with OS for MM cases in two EA populations. Given that prior studies have shown association of genetic variation with MM survival, efforts to identify additional loci associated with OS or MM specific survival are warranted. Future studies should also consider germline variants impact on molecular subtypes, specific therapies, and outcomes. Disclosures Kumar: Celgene: Consultancy, Research Funding; Janssen: Consultancy, Research Funding; Takeda: 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 Genome-wide association studies (GWAS) conducted among populations of European ancestry (EA) have identified 23 common single nucleotide polymorphisms (SNPs) associated with multiple myeloma (MM) risk. We hypothesize that the combination of these SNPs in a polygenic risk score (PRS) is likely to be a strong risk factor for MM. However, it is unclear whether the genetic variation associated with MM susceptibility also predisposes to monoclonal gammopathy of undetermined significance (MGUS). Thus, we calculated a PRS and evaluated the association with risk of MM and its precursor, MGUS. We pooled genotype data for 2434 MM and 3446 controls from ten MM GWAS of individuals of EA within the Interlymph Consortium, for 23 MM risk SNPs identified by prior GWAS. An additional 754 MGUS cases were ascertained from Mayo Clinic and MD Anderson clinical practices. To calculate the PRS, we used the risk estimates corresponding to the 23 SNP associations from the largest published MM GWAS. The log of the odds ratio (OR) for each SNP was multiplied by the respective number of risk alleles and summed to generate a PRS for each individual. The PRS was examined continuously, per one standard deviation (SD), and as quintiles, based on the PRS distribution in the controls. Associations of PRS with MM and MGUS risk were examined separately, using multivariable logistic regression assuming an additive model to assess ORs and 95% confidence intervals adjusted for age, sex, and site. We also evaluated age and sex stratified models. The distribution of sex within MM cases, MGUS cases and controls were each ~60% male and ~40% female. The median age was 61, 66, and 66 years for MM cases, MGUS cases and controls, respectively. PRS ranged from 1.52-4.91, with a median PRS of 3.21 for MM cases, 3.19 for MGUS cases, and 3.05 for controls. PRS was significantly associated with MM risk when assessed continuously (OR=1.19 per SD, p=2.2x10-16) and categorically; compared with the middle quintile (Q3), individuals in the highest quintile (Q5) had a 66% increased MM risk (OR=1.66, p=2.3x10-9) and those in the lowest quintile (Q1) had a 38% decreased MM risk (OR=0.62, p=1.3x10-6). PRS was also significantly associated with MGUS risk (OR=1.19 per SD, p=1.7x10-11); individuals with the highest PRS (Q5) had a 77% increased risk (OR=1.77, p=4.0x10-4) and those with lowest PRS (Q1) had 30% decreased risk (OR=0.70, p=0.04), compared with Q3. When stratified by age and sex, similar associations and trends were found. Using an independent sample of MM / MGUS cases and controls, we showed that a PRS constructed from 23 common genetic variants for MM risk is associated with risk of both MM and MGUS, regardless of age or sex. A future direction of this work is testing associations with PRS and clinical characteristics of the MM cases, as well as differences between MGUS cases that progress and those that do not. Our results suggest that common genetic variation may predispose to MGUS as the precursor to MM. Citation Format: Alyssa I. Clay-Gilmour, Michelle A. Hildebrandt, Nicola J. Camp, Elad Ziv, Elizabeth E. Brown, Jonathan N. Hofmann, John J. Spinelli, Graham G. Giles, Parveen Bhatti, Wendy Cozen, Xifeng Wu, Dennis P. Robinson, Aaron D. Norman, Jason P. Sinnwell, Shaji K. Kumar, S Vincent Rajkumar, Susan L. Slager, Celine M. Vachon. Associations between a polygenic risk score and risk of multiple myeloma and its precursor [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2686.

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.