Dumontet, Charles

Background/Objectives: Family history is a known risk factor for multiple myeloma (MM) and its precursor condition, monoclonal gammopathy of undetermined significance (MGUS). Previous genome-wide association studies (GWASs) have identified 35 common loci associated with MM risk and 21 associated with MGUS. The objective of this study was to identify less common and rare genetic loci predisposing to MM/MGUS through whole-exome sequencing (WES)-based linkage analysis. Methods:Multipoint linkage analysis was conducted using the Multipoint Engine for Rapid Likelihood Inference (MERLIN) with the Lander–Green algorithm on germline WES data from 79 pedigrees with 2 or more affected relatives (120 MM, 86 MGUS, and 21 unaffected). Genome-wide linkage was evaluated using 12,946 independent single-nucleotide variants (linkage disequilibrium r2 < 0.05). Results: Significant linkage was observed at chromosome 6q22.33–q24.2 by the non-parametric model (logarithm-of-odds (LOD) = 3.3) and suggestive linkage by the dominant parametric model (heterogeneity LOD (HLOD) = 2.5). Fourteen rare variants within this region were prioritized using family-specific partial LOD scores and in silico functional prediction tools. Nine of these variants, REPS1, THEMIS, TAAR6, AHI1, VNN1, VNN3, MTFR2/FAM54A, LAMA2, and PHACTR2, overlapped immune-regulatory regions in blood cell lines and were not previously identified in GWASs. Conclusions: This study demonstrates the utility of applying a linkage analysis framework to familial WES data for identifying genomic regions and candidate genes that may contribute to MM/MGUS predisposition. These findings provide new insight into the inherited risk and etiology of familial MM and MGUS.

Introduction: WM is a rare subtype of lymphoma, identified by the presence of IgM paraprotein, infiltration of clonal, small lymphoplasmacytic B cells in the bone marrow, and the somatic mutation, MYD88 L265P, which is observed in over 90% of WM patients. We previously identified two rare variants at 6p25.3 (EXOC2) and 14q32.13 (TCL1) that were significantly associated with an increased risk of WM/LPL (PMID: 30305637). However, the genetic basis for WM/LPL predisposition remains largely unknown. We conducted the largest WM/LPL GWAS to discover further germline susceptibility loci for WM/LPL. Methods: Our meta-analysis discovery set included 1,707 cases and 58,328 controls of European ancestry from the InterLymph Consortium, genotyped on a high-density SNP array and imputed to TOPMed, providing over 10 million common variants for association testing. Germline genetic loci with suggestive significance (P&lt;5x10-6) were selected for replication in four independent sets (570 cases and 771,770 controls; UK Biobank, Copenhagen Hospital Biobank, FinnGen, Mayo Clinic), giving us a total of 2,277 WM/LPL cases and 831,098 controls. Conditional and sensitivity analyses were performed considering family history and sex. A transcriptome-wide association study (TWAS) was performed using FUSION with expression quantitative trait loci (eQTL) data from GTEx and Young Finns Study in whole blood, lymphocytes, and other tissues. Functional annotations were done using FUMA, F-MAGMA, FORGE, and OpenTarget. Colocalization was measured using CoLoC with whole blood gene expression data. Results & Discussion: We identified 11 germline genetic loci associated with WM/LPL risk that replicated and reached genome-wide significance (P&lt; 5x10-8), including the two known loci on 6p25.3 and 14q32.13 and nine novel loci in the following regions: 6p25.3 (n=2 independent loci), 2q37.1, 9p13.2, 6p21.32, 10q23.31, 16q21, 1p31.1, and 15q15.1. Four of the discovered variants showed strong associations with WM/LPL [odds ratios (ORs) &gt;10]. Enrichment was observed for enhancer and promoter elements in B-cells (FDR q &lt; 0.01), and gene-based and colocalization analyses identified several potential biological pathways. Colocalization analyses show the intronic PAX5 risk variant on 9p (OR=1.5, P= 2.3x10-25) associated with whole blood eQTL expression of PAX5, which plays an essential role in regulating lymphoid progenitors to the B-lymphocyte lineage, as well as RP11-297B17.3 (RNA gene, antisense to PAX5). In addition, GWAS have shown this variant is in linkage disequilibrium (LD) with GWAS loci for CD20 on IgD+ CD38- B cell that represents a specific immunophenotype that helps identify and characterize B cells involved with WM and other similar conditions. The 2q37.1 (SP140) locus (OR=1.6, P=1.2x10-33) contains variants linked to chronic lymphocytic leukemia (CLL) and autoimmune diseases, including multiple sclerosis and Crohn's disease. This locus is important in chromatin-mediated gene regulation, suggesting a significant role in the immune system. The specific variant influences exon-skipping and potentially protein expression, altering the gene's function. The variants on 10q23.31 in FAS (OR=1.3, P=4.4x10-13) and 15q15.1 (OR=1.2, P=8.5x10-11) in BCL2 modifying factor (BMF), are in strong LD with known CLL susceptibility variants. The 15q variant is also linked with multiple lymphoid leukemia GW significant SNPs. Sex-specific analyses identified one novel locus at 20q13.33 associated with the expression of ZGPAT and UCKL1 in males (OR=1.9, P=9.5x10-10). TWAS identified associations with CRHR1 and CDKN2A; overexpression of both are shown to be associated with an increased risk of WM/LPL. Expression of CDKN2A has been linked to childhood acute lymphoblastic leukemia, while CRHR1 has been associated with breast cancer risk. Conclusions: Our study, the largest GWAS of WM to date, has identified 11 susceptibility loci, including two previously known and nine novel regions. The findings underscore the genetic complexity of WM, with the PAX5 and SP140 loci implicated in both chromatin-mediated gene regulation and immune system functions. These insights not only enhance our understanding of WM predisposition but also pave the way for future research focused on functional understanding of these loci, potentially leading to new therapeutic strategies.

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.

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 Inherited genetic risk factors play a role in multiple myeloma (MM), yet considerable missing heritability exists. Rare risk variants at genome-wide association study (GWAS) loci are a new avenue to explore. Pleiotropy between lymphoid neoplasms (LNs) has been suggested in family history and genetic studies, but no studies have interrogated sequencing for pleiotropic genes or rare risk variants. Sequencing genetically enriched cases can help discover rarer variants. We analyzed exome sequencing in familial or early-onset MM cases to identify rare, functionally relevant variants near GWAS loci for a range of LNs. A total of 149 distinct and significant LN GWAS loci have been published. We identified six recurrent, rare, potentially deleterious variants within 5 kb of significant GWAS single nucleotide polymorphisms in 75 MM cases. Mutations were observed in BTNL2, EOMES, TNFRSF13B, IRF8, ACOXL and TSPAN32. All six genes replicated in an independent set of 255 early-onset MM or familial MM or precursor cases. Expansion of our analyses to the full length of these six genes resulted in a list of 39 rare and deleterious variants, seven of which segregated in MM families. Three genes also had significant rare variant burden in 733 sporadic MM cases compared with 935 control individuals: IRF8 (P = 1.0 × 10−6), EOMES (P = 6.0 × 10−6) and BTNL2 (P = 2.1 × 10−3). Together, our results implicate six genes in MM risk, provide support for genetic pleiotropy between LN subtypes and demonstrate the utility of sequencing genetically enriched cases to identify functionally relevant variants near GWAS loci.

Abstract Introduction Multiple myeloma (MM) is a result of a malignant transformation of plasma cells that is preceded by the presence of an asymptomatic clonal plasma cell expansion, a condition referred to as monoclonal gammopathy of undetermined significance (MGUS). We and others have shown familial aggregation of MM and MGUS. Evidence from epidemiologic, family and genome-wide association studies (GWAS) suggests a genetic component underlying MM etiology. GWAS have successfully established 17 common genetic risk loci for MM to date and recently, rare inherited susceptibility variants in the LSD1 / KDM1A and USP45 genes were identified in familial MM / MGUS kindreds. Family-based approaches may be used to elucidate genetic variation contributing to familial MM. Genetic linkage analysis has historically been used to detect the chromosomal location of disease genes. The objective of this study was to conduct a linkage analysis of MM / MGUS families to identify genomic regions for MM / MGUS. Methods Linkage analysis was performed on whole exome sequencing (WES) data generated from germline DNA extracted from peripheral blood from MM / MGUS families ascertained at four sites, of which 79 were selected with 2 or more affected members with any combination of MM or MGUS. Whole exome capture was performed using Agilent SureSelect 38 Mb paired end sequencing and ran on Illumina HiSeq 2000s/2500s; standard analyses for alignment (GRCh37) and quality control were conducted. The Genome Analysis Toolkit's (GATK) HaplotypeCaller in per-sample mode was used to jointly call all the samples together. Quality control included removing variants that had <75% call rate, <8x coverage, or minor allele frequency <0.01 (1KGenomes). The jointly called WES data from the 79 families consisted of 230 individuals, including 119 MM, 93 MGUS and 18 unaffected relatives. Additional early onset (defined as age <60) MM cases (n=400), sporadic MM cases (n=879) and controls (n=2389) were also included in the joint calling to be used for follow-up of linkage regions. We included independent variants defined as linkage disequilibrium (LD) r2 <0.05 using both PLINK 1.07's LD based variant pruner and also MERLIN's pairwise r2 cluster marker-marker modeling to ensure variant independence. After quality control and LD filtering: 12,946 variants remained. We conducted multipoint linkage analysis using MERLIN (Multipoint Engine for Rapid Likelihood Inference using the Lander-Green approach) and considered both non-parametric and parametric (dominant and recessive) models to test for co-segregation of chromosomal regions with MM / MGUS. A logarithm of the odds (LOD) score greater than 3.3 was significant evidence for linkage. Within the same families, we interrogated any suggestive regions by performing gene-level association tests using the PEDGENE R package to calculate burden and kernel statistics, accounting for family relationships. Follow-up of these genes in the early onset and sporadic MM cases and controls and interrogation of single variants driving the associations of these genes is ongoing. Results Among the 79 families, 28 consisted of two or more MM, 41 had two or more MM and MGUS cases and 10 had two or more MGUS. Analyses of the 79 MM / MGUS families identified chromosome 6 at 123420001-149070000 BP region, with evidence for linkage by the non-parametric model (LOD score=3.3) and supportive evidence for linkage by the dominant parametric model (LOD=2.5) (Figure). The chromosome 6 region, which is outside the HLA-region, contains a total of 72 genes. Using gene-level association tests, we identified four genes within the region (MOXD1, TRDN, ADAT2, LTV1) that were associated within MM / MGUS in the family cases (p<0.05). Discussion We found evidence for a locus on chromosome 6 linked to MM / MGUS. Four genes in this region may be associated with MM / MGUS in these families. Follow-up of these genes/variants in high-risk individuals (early-onset / age <60), sporadic cases and controls will help elucidate genetic mechanisms contributing to familial MM / MGUS. In addition, we demonstrate the benefit of applying a linkage analysis framework to familial WES data for discovery of genomic regions influencing MM / MGUS. Disclosures Dumontet: Roche: Research Funding; Merck: Consultancy, Membership on an entity's Board of Directors or advisory committees; Janssen: Honoraria; Sanofi: Honoraria. Kumar:AbbVie: Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees, Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; AbbVie: Membership on an entity's Board of Directors or advisory committees, Research Funding; Oncopeptides: Membership on an entity's Board of Directors or advisory committees; Takeda: Membership on an entity's Board of Directors or advisory committees; Merck: Membership on an entity's Board of Directors or advisory committees, Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; KITE: Membership on an entity's Board of Directors or advisory committees, Research Funding; Novartis: Research Funding; Roche: Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees, Research Funding; KITE: Membership on an entity's Board of Directors or advisory committees, Research Funding.

Abstract Given the frequent and largely incurable occurrence of multiple myeloma, identification of germline genetic mutations that predispose cells to multiple myeloma may provide insight into disease etiology and the developmental mechanisms of its cell of origin, the plasma cell (PC). Here, we identified familial and early-onset multiple myeloma kindreds with truncating mutations in lysine-specific demethylase 1 (LSD1/KDM1A), an epigenetic transcriptional repressor that primarily demethylates histone H3 on lysine 4 and regulates hematopoietic stem cell self-renewal. In addition, we found higher rates of germline truncating and predicted deleterious missense KDM1A mutations in patients with multiple myeloma unselected for family history compared with controls. Both monoclonal gammopathy of undetermined significance (MGUS) and multiple myeloma cells have significantly lower KDM1A transcript levels compared with normal PCs. Transcriptome analysis of multiple myeloma cells from KDM1A mutation carriers shows enrichment of pathways and MYC target genes previously associated with myeloma pathogenesis. In mice, antigen challenge followed by pharmacologic inhibition of KDM1A promoted PC expansion, enhanced secondary immune response, elicited appearance of serum paraprotein, and mediated upregulation of MYC transcriptional targets. These changes are consistent with the development of MGUS. Collectively, our findings show that KDM1A is the first autosomal-dominant multiple myeloma germline predisposition gene providing new insights into its mechanistic roles as a tumor suppressor during post-germinal center B-cell differentiation. Significance: KDM1A is the first germline autosomal dominant predisposition gene identified in multiple myeloma and provides new insights into multiple myeloma etiology and the mechanistic role of KDM1A as a tumor suppressor during post-germinal center B-cell differentiation. Cancer Res; 78(10); 2747–59. ©2018 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