Bhatti, Parveen

Summary Evidence for an association between insulin‐like growth factors (IGF) and multiple myeloma (MM) is inconsistent. We examined total IGF‐I concentrations and risk of MM by combining baseline serological data among UK Biobank participants ( n = 444 187; 732 incident MM) with a two‐sample Mendelian randomisation (MR) analysis using identified genetic variants associated with circulating total IGF‐I and IGF‐binding protein 3 (IGFBP‐3) in the InterLymph consortium (2434 MM and their 2567 controls). Finally, additional lymphoid neoplasm (LN) subtypes were included for comparison with the main hypothesis. Circulating IGF‐I level was positively associated with MM risk Hazard ratio‐HR‐per one standard deviation‐SD‐increase (HR 1‐SD = 1.11, 95% confidence interval [CI]: 1.01–1.22; p ‐value = 0.03), especially closer to diagnosis. Genetically inferred IGF‐I levels were associated with increased MM risk (odds ratio [OR] = 1.27, 95% CI: 1.05–1.54) but not with any other LNs. Genetically inferred IGFBP‐3 levels showed no associations with any LN evaluated. Corroborating previous findings, in a secondary analysis, IGF‐I levels were associated with the risk of chronic lymphocytic leukaemia/small lymphocytic lymphoma (CLL/SLL) in males with higher body mass index (HR 1‐SD in obese male = 1.36, 95% CI: 1.14–1.61). Our serological and MR analyses suggest a contributing role of IGF‐I in the susceptibility of MM; the CLL/SLL findings warrant further investigation considering sex‐specific adiposity.

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.

AbstractPleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome‐wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10−8) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4‐rs34517439‐A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13‐1.32, P = 4.81 × 10−7). rs34517439‐A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA‐binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4‐rs34517439 as a potentially novel MM risk locus.

Abstract Background: Genome-wide association studies (GWAS) of multiple myeloma in populations of European ancestry (EA) identified and confirmed 24 susceptibility loci. For other cancers (e.g., colorectum and melanoma), risk loci have also been associated with patient survival. Methods: We explored the possible association of all the known risk variants and their polygenic risk score (PRS) with multiple myeloma overall survival (OS) in multiple populations of EA [the International Multiple Myeloma rESEarch (IMMEnSE) consortium, the International Lymphoma Epidemiology consortium, CoMMpass, and the German GWAS] for a total of 3,748 multiple myeloma cases. Cox proportional hazards regression was used to assess the association between each risk SNP with OS under the allelic and codominant models of inheritance. All analyses were adjusted for age, sex, country of origin (for IMMEnSE) or principal components (for the others) and disease stage (ISS). SNP associations were meta-analyzed. Results: SNP associations were meta-analyzed. From the meta-analysis, two multiple myeloma risk SNPs were associated with OS (P &lt; 0.05), specifically POT1-AS1-rs2170352 [HR = 1.37; 95% confidence interval (CI) = 1.09–1.73; P = 0.007] and TNFRSF13B-rs4273077 (HR = 1.19; 95% CI = 1.01–1.41; P = 0.04). The association between the combined 24 SNP MM-PRS and OS, however, was not significant. Conclusions: Overall, our results did not support an association between the majority of multiple myeloma risk SNPs and OS. Impact: This is the first study to investigate the association between multiple myeloma PRS and OS in multiple myeloma.

Abstract Background: A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non–Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes. Methods: In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell–mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression. Results: We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08–1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59–2.07; P-trend (Ptrend) &lt; 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell–mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend &lt; 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction. Conclusions: Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions. Impact: Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.

AbstractEvidence for the human health effects of pesticides is needed to inform risk assessment. We studied the relationship between occupational insecticide use and risk of non‐Hodgkin lymphoma (NHL) by pooling data from nine case‐control studies participating in the InterLymph Consortium, including 7909 cases and 8644 controls from North America, the European Union and Australia. Insecticide use was coded using self‐report or expert assessment, for insecticide groups (eg, organophosphates, pyrethroids) and active ingredients (eg, malathion, permethrin). Associations with insecticides were estimated using logistic regression to produce odds ratios (ORs) and 95% confidence intervals (CI) for all NHL and NHL subtypes, with adjustment for study site, demographic factors and use of other pesticides. Occupational insecticide use, overall, was not associated with risk of NHL. Use of organophosphate insecticides was associated with increased risk of all NHL and the subtype follicular lymphoma, and an association was found with diazinon, in particular (ever use: OR = 2.05, 95%CI: 1.24‐3.37). The carbamate insecticide, carbaryl, was associated with risk of all NHL, and the strongest associations were found with T‐cell NHL for ever‐use (OR = 2.44, 95%CI: 1.13‐5.28) and longer duration (>8 years vs never: OR = 2.90, 95%CI: 1.02‐8.25). There was no association of NHL with other broad groups of insecticides, including organochlorines and pyrethroids, and some inverse associations were estimated in relation to historical DDT use. Our findings contribute to the totality of evidence available to help inform risk decisions by public health and regulatory agencies of importance given continued, widespread use of organophosphate and carbamate insecticides.

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.

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 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.