Stephens, Owen

Abstract Background: Persons of African ancestry (AA) experience a 1.5-2-fold risk of multiple myeloma (MM) compared to persons of European ancestry (EA). We assembled a set of MM patients with self-reported AA in order to evaluate the contribution of genetics to etiology in this high-risk group. Methods: Here we present the results of a meta-analysis of two GWAS in 1,813 cases and 8,871 controls of AA. We also conducted an admixture mapping scan to identify risk alleles associated with local ancestry, fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA, and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. Finally, we conducted an eQTL analysis measuring gene expression in those genes harboring a risk variant in malignant plasma cells from 292 of the patients from a single site. Results: In GWAS analysis, we identified two suggestive novel loci located at 9p24.3 and 9p13.1 at P<1 × 10-6, but no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. 20 of the 23 known EA risk variants showed directional consistency and 9 replicated at P<0.05 in AA individuals. In eight regions, we identified markers that better capture MM risk in persons of AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95%CI: 1.56, 2.11) increased MM risk compared to those with average risk (25-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P= 5.1 × 10–12). Conclusion: Our study shows that common genetic variation contributes to MM risk individuals of AA. This abstract is also being presented as Poster C040. Citation Format: Zhaohui Du, Niels Weinhold, Gregory Chi Song, Kristen A. Rand, David J. Van Den Berg, Amie E. Hwang, Xin Sheng, Victor Hom, Sikander Ailawadhi, Ajay K. Nooka, Seema Singhal, Karen Pawlish, Edward Peters, Cathryn Bock, Ann Mohrbacher, Alexander Stram, Sonja I. Berndt, William J. Blot, Graham Casey, Victoria L. Stevens, Rick Kittles, Phyllis J. Goodman, W. Ryan Diver, Anselm Hennis, Barbara Nemesure, Eric A. Klein, Benjamin A. Rybicki, Janet L. Stanford, John S. Witte, Lisa Signorello, Esther M. John, Leslie Bernstein, Antoinette Stroup, Owen W. Stephens, Maurizio Zangari, Frits Van Rhee, Andrew Olshan, Wei Zheng, Jennifer J. Hu, Regina Ziegler, Sarah J. Nyante, Sue Ann Ingles, Michael Press, John David Carpten, Stephen Chanock, Jayesh Mehta, Graham A Colditz, Jeffrey Wolf, Thomas G. Martin, Michael Tomasson, Mark A. Fiala, Howard Terebelo, Nalini Janakiraman, Laurence Kolonel, Kenneth C. Anderson, Loic Le Marchand, Daniel Auclair, Brian C.-H. Chiu, Elad Ziv, Daniel Stram, Ravi Vij, Leon Bernal-Mizrachi, Gareth J. Morgan, Jeffrey A. Zonder, Carol Ann Huff, Sagar Lonial, Robert Z. Orlowski, David V. Conti, Christopher A. Haiman, Wendy Cozen. A meta-analysis of genome-wide association study and eQTL analysis of multiple myeloma among African Americans [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr PR05.

Abstract Persons of African ancestry (AA) have a twofold higher risk for multiple myeloma (MM) compared with persons of European ancestry (EA). Genome-wide association studies (GWASs) support a genetic contribution to MM etiology in individuals of EA. Little is known about genetic risk factors for MM in individuals of AA. We performed a meta-analysis of 2 GWASs of MM in 1813 cases and 8871 controls and conducted an admixture mapping scan to identify risk alleles. We fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. In GWAS meta-analysis, we identified 2 suggestive novel loci located at 9p24.3 and 9p13.1 at P < 1 × 10−6; however, no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. Of the 23 known EA risk variants, 20 showed directional consistency, and 9 replicated at P < .05 in AA individuals. In 8 regions, we identified markers that better capture MM risk in persons with AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95% confidence interval, 1.56-2.11) increased MM risk compared with those with average risk (25%-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P = 5.1 × 10−12). Our study shows that common genetic variation contributes to MM risk in individuals with AA.

Abstract Background: Bone disease in myeloma occurs as a result of complex interactions between myeloma cells and the bone marrow microenvironment. To date, no studies have evaluated the potential impact of genetic polymorphisms upon and/or within this microenvironment. Patients and Methods: Peripheral blood DNA from 282 patients enrolled in the Total Therapy 2 (TT2) protocol was studied using the previously reported Affymetrix 3k BOAC custom chip to evaluate relevant genetic polymorphisms. DKK1 gene expression and high risk GEP gene signatures were assessed as previously reported (Blood109:4470–4477 and 109:2276–2284 2007). Patients were classified using both full skeletal x-rays and MRI findings. The lower cut-off used was absence of focal abnormalities on x-ray and/or <7 focal lesions on MRI (see JCO25:1121–1128 2007). Results: The top 200 SNPs were first evaluated based upon univariate P values linked to limited or extensive bone disease. The top 50 SNPs with the smallest P values were then selected (eliminating closely linked genes) for recursive partitioning analysis. The recursive partitioning was conducted both with and without the insertion of known biologically relevant SNPs. The pruned tree developed with recursive partitioning proved to be quite stable and incorporated 4 dominant SNPs: rs 3766934 Epoxide hydrolase (EPHX1); sr3783408 MAP kinase; sr 1062637 RNA helicase DDX18; and sr3181366 TNFSF8-TNF-α. No alternate SNP substituted for EPHX1 in the recursive model. The 4 SNP signature of EPHX1 GG, MAP4K5 AG/AA, DDX18 GG/CC plus TNFSF8 CT/TT was highly correlated with the number of MRI lesions: mean 8.66 lesions versus 3.33 for alternate SNPs (P<.001). In the univariate association of standard prognostic factors, 4 SNP signature and 70 and 17 gene models for high risk, the best correlation with bone disease was with the 4 SNPs (P<.0001) followed by the 17 gene GEP model (P=0.05). Using stepwise multivariate regression analysis the best correlation with bone disease was Log2DKK-1 expression (P=.0001) but closely followed by SNPs EPHX1 (P=.0009); MAP4K5 (P=.0071) plus TNFSF8 and serum LDH (P=.02). The various analyses training (2/3) and validation sets were used with assessment of sensitivity and specificity. The microenvironmental 4 SNP signature combined with GEP DKK1 provided the best prediction of myeloma bone disease and overall outcome. Conclusions: This first assessment of genetic polymorphisms linked to myeloma bone disease in myeloma has led to the identification of polymorphisms with both potential biologic relevance and utility in prognostic models of myeloma bone disease including risk stratification for bisphosphonate use. Validation of these findings may allow a search for potential therapeutic targets.