Identification of multiple risk loci and regulatory mechanisms influencing susceptibility to multiple myeloma

Genome-wide association studies (GWAS) have transformed our understanding of susceptibility to multiple myeloma (MM), but much of the heritability remains unexplained. We report a new GWAS, a meta-analysis with previous GWAS and a replication series, totalling 9974 MM cases and 247,556 controls of European ancestry. Collectively, these data provide evidence for six new MM risk loci, bringing the total number to 23. Integration of information from gene expression, epigenetic profiling and in situ Hi-C data for the 23 risk loci implicate disruption of developmental transcriptional regulators as a basis of MM susceptibility, compatible with altered B-cell differentiation as a key mechanism. Dysregulation of autophagy/apoptosis and cell cycle signalling feature as recurrently perturbed pathways. Our findings provide further insight into the biological basis of MM

Application of an artificial intelligence-based tool in [18F]FDG PET/CT for the assessment of bone marrow involvement in multiple myeloma

Purpose: [18F]FDG PET/CT is an imaging modality of high performance in multiple myeloma (MM). Nevertheless, the inter-observer reproducibility in PET/CT scan interpretation may be hampered by the different patterns of bone marrow (BM) infiltration in the disease. Although many approaches have been recently developed to address the issue of standardization, none can yet be considered a standard method in the interpretation of PET/CT. We herein aim to validate a novel three-dimensional deep learning-based tool on PET/CT images for automated assessment of the intensity of BM metabolism in MM patients. Materials and methods: Whole-body [18F]FDG PET/CT scans of 35 consecutive, previously untreated MM patients were studied. All patients were investigated in the context of an open-label, multicenter, randomized, active-controlled, phase 3 trial (GMMG-HD7). Qualitative (visual) analysis classified the PET/CT scans into three groups based on the presence and number of focal [18F]FDG-avid lesions as well as the degree of diffuse [18F]FDG uptake in the BM. The proposed automated method for BM metabolism assessment is based on an initial CT-based segmentation of the skeleton, its transfer to the SUV PET images, the subsequent application of different SUV thresholds, and refinement of the resulting regions using postprocessing. In the present analysis, six different SUV thresholds (Approaches 1–6) were applied for the definition of pathological tracer uptake in the skeleton [Approach 1: liver SUVmedian 7 1.1 (axial skeleton), gluteal muscles SUVmedian 7 4 (extremities). Approach 2: liver SUVmedian 7 1.5 (axial skeleton), gluteal muscles SUVmedian 7 4 (extremities). Approach 3: liver SUVmedian 7 2 (axial skeleton), gluteal muscles SUVmedian 7 4 (extremities). Approach 4: ≥ 2.5. Approach 5: ≥ 2.5 (axial skeleton), ≥ 2.0 (extremities). Approach 6: SUVmax liver]. Using the resulting masks, subsequent calculations of the whole-body metabolic tumor volume (MTV) and total lesion glycolysis (TLG) in each patient were performed. A correlation analysis was performed between the automated PET values and the results of the visual PET/CT analysis as well as the histopathological, cytogenetical, and clinical data of the patients. Results: BM segmentation and calculation of MTV and TLG after the application of the deep learning tool were feasible in all patients. A significant positive correlation (p < 0.05) was observed between the results of the visual analysis of the PET/CT scans for the three patient groups and the MTV and TLG values after the employment of all six [18F]FDG uptake thresholds. In addition, there were significant differences between the three patient groups with regard to their MTV and TLG values for all applied thresholds of pathological tracer uptake. Furthermore, we could demonstrate a significant, moderate, positive correlation of BM plasma cell infiltration and plasma levels of β2-microglobulin with the automated quantitative PET/CT parameters MTV and TLG after utilization of Approaches 1, 2, 4, and 5. Conclusions: The automated, volumetric, whole-body PET/CT assessment of the BM metabolic activity in MM is feasible with the herein applied method and correlates with clinically relevant parameters in the disease. This methodology offers a potentially reliable tool in the direction of optimization and standardization of PET/CT interpretation in MM. Based on the present promising findings, the deep learning-based approach will be further evaluated in future prospective studies with larger patient cohorts

Genome-wide assessment of the association of rare and common copy number variations to testicular germ cell cancer.

Testicular germ cell cancer (TGCC) is one of the most heritable forms of cancer. Previous genome-wide association studies have focused on single nucleotide polymorphisms (SNPs), largely ignoring the influence of copy number variants (CNVs). Here we present a genome-wide study of copy number variation on a cohort of 212 cases and 437 controls from Denmark, which was genotyped at ~1.8 million markers, half of which were non-polymorphic copy number markers. No association of common variants were found, whereas analysis of rare variants (present in less than 1% of the samples) initially indicated a single gene with significantly higher accumulation of rare CNVs in cases as compared to controls, at the gene PTPN1 (P=3.8*10-2, 0.9% of cases and 0% of controls). However, the CNV could not be verified by qPCR in the affected samples. The CNV calling of the array data was validated by sequencing of the GSTM1 gene, which showed that the CNV frequency was in complete agreement between the two platforms. This study therefore disconfirms the hypothesis that there exists a single CNV locus with a major effect size that predisposes to TGCC. Genome-wide pathway association analysis indicated a weak association of rare CNVs related to cell migration (FDR=0.021, 1.8% of cases and 1.1% of controls). Dysregulation during migration of primordial germ cells has previously been suspected to be a part of TGCC development and this set of multiple rare variants may thereby have a minor contribution to an increased susceptibility of TGCCs

Enrichment of B cell receptor signaling and epidermal growth factor receptor pathways in monoclonal gammopathy of undetermined significance: a genome-wide genetic interaction study

Background: Recent identification of 10 germline variants predisposing to monoclonal gammopathy of undetermined significance (MGUS) explicates genetic dependency of this asymptomatic precursor condition with multiple myeloma (MM). Yet much of genetic burden as well as functional links remain unexplained. We propose a workflow to expand the search for susceptibility loci with genome-wide interaction and for subsequent identification of genetic clusters and pathways. Methods: Polygenic interaction analysis on 243 cases/1285 controls identified 14 paired risk loci belonging to unique chromosomal bands which were then replicated in two independent sets (case only study, 82 individuals; case/control study 236 cases/ 2484 controls). Further investigation on gene-set enrichment, regulatory pathway and genetic network was carried out with stand-alone in silico tools separately for both interaction and genome-wide association study-detected risk loci. Results: Intronic-PREX1 (20q13.13), a reported locus predisposing to MM was confirmed to have contribution to excess MGUS risk in interaction with SETBP1, a well-established candidate predisposing to myeloid malignancies. Pathway enrichment showed B cell receptor signaling pathway (P < 5.3 × 10− 3) downstream to allograft rejection pathway (P < 5.6 × 10− 4) and autoimmune thyroid disease pathway (P < 9.3 × 10− 4) as well as epidermal growth factor receptor regulation pathway (P < 2.4 × 10− 2) to be differentially regulated. Oncogene ALK and CDH2 were also identified to be moderately interacting with rs10251201 and rs16966921, two previously reported risk loci for MGUS. Conclusions: We described novel pathways and variants potentially causal for MGUS. The methodology thus proposed to facilitate our search streamlines risk locus-based interaction, genetic network and pathway enrichment analyses

Analyzing Longitudinal wb-MRI Data and Clinical Course in a Cohort of Former Smoldering Multiple Myeloma Patients: Connections between MRI Findings and Clinical Progression Patterns

The purpose of this study was to analyze size and growth dynamics of focal lesions (FL) as well as to quantify diffuse infiltration (DI) in untreated smoldering multiple myeloma (SMM) patients and correlate those MRI features with timepoint and cause of progression. We investigated 199 whole-body magnetic resonance imaging (wb-MRI) scans originating from longitudinal imaging of 60 SMM patients and 39 computed tomography (CT) scans for corresponding osteolytic lesions (OL) in 17 patients. All FLs >5 mm were manually segmented to quantify volume and growth dynamics, and DI was scored, rating four compartments separately in T1- and fat-saturated T2-weighted images. The majority of patients with at least two FLs showed substantial spatial heterogeneity in growth dynamics. The volume of the largest FL (p = 0.001, c-index 0.72), the speed of growth of the fastest growing FL (p = 0.003, c-index 0.75), the DI score (DIS, p = 0.014, c-index 0.67), and its dynamic over time (DIS dynamic, p < 0.001, c-index 0.67) all significantly correlated with the time to progression. Size and growth dynamics of FLs correlated significantly with presence/appearance of OL in CT within 2 years after the respective MRI assessment (p = 0.016 and p = 0.022). DIS correlated with decrease of hemoglobin (p < 0.001). In conclusion, size and growth dynamics of FLs correlate with prognosis and local bone destruction. Connections between MRI findings and progression patterns (fast growing FL—OL; DIS—hemoglobin decrease) might enable more precise diagnostic and therapeutic approaches for SMM patients in the future

Characterization of rare germline variants in familial multiple myeloma

Multiple myeloma (MM) is a malignancy of plasma cells, characterized by the presence of monoclonal immunoglobulin, known as M protein1. MM is preceded by monoclonal gammopathy of undetermined significance (MGUS) which is also a precursor of immunoglobulin light chain (AL) amyloidosis. Previous studies have reported a 2- to 4-fold increased risk of MGUS or MM in first-degree relatives of MM or MGUS patients, suggesting the existence of inherited susceptibility. For many years, high-risk germline predisposing genes have been lacking for MM. However, recent sequencing efforts have proposed a few novel candidates, most notably loss-offunction (LoF) variants in the tumor suppressor gene DIS3 and in the histone demethylase gene KDM1A, and others as recently reviewed in detail in Pertesi et al. In addition to the suspected rare, high-penetrance variants, genome-wide association studies have identified over 20 common, low-penetrance variants associated with the risk of MM; these were estimated to account for about 15% of the familial MM risk

A genome-wide association study of men with symptoms of testicular dysgenesis syndrome and its network biology interpretation

Background Testicular dysgenesis syndrome (TDS) is a common disease that links testicular germ cell cancer, cryptorchidism and some cases of hypospadias and male infertility with impaired development of the testis. The incidence of these disorders has increased over the last few decades, and testicular cancer now affects 1% of the Danish and Norwegian male population. Methods To identify genetic variants that span the four TDS phenotypes, the authors performed a genome-wide association study (GWAS) using Affymetrix Human SNP Array 6.0 to screen 488 patients with symptoms of TDS and 439 selected controls with excellent reproductive health. Furthermore, they developed a novel integrative method that combines GWAS data with other TDS-relevant data types and identified additional TDS markers. The most significant findings were replicated in an independent cohort of 671 Nordic men. Results Markers located in the region of TGFBR3 and BMP7 showed association with all TDS phenotypes in both the discovery and replication cohorts. An immunohistochemistry investigation confirmed the presence of transforming growth factor beta receptor type III (TGFBR3) in peritubular and Leydig cells, in both fetal and adult testis. Single-nucleotide polymorphisms in the KITLG gene showed significant associations, but only with testicular cancer. Conclusions The association of single-nucleotide polymorphisms in the TGFBR3 and BMP7 genes, which belong to the transforming growth factor b signalling pathway, suggests a role for this pathway in the pathogenesis of TDS. Integrating data from multiple layers can highlight findings in GWAS that are biologically relevant despite having border significance at currently accepted statistical levels

Comprehensive genomic analysis of refractory multiple myeloma reveals a complex mutational landscape associated with drug resistance and novel therapeutic vulnerabilities

The outcomes of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) remain poor. In this study, we performed whole genome and transcriptome sequencing of 39 heavily pretreated relapsed/refractory MM (RRMM) patients to identify mechanisms of resistance and potential therapeutic targets. We observed a high mutational load and indications of increased genomic instability. Recurrently mutated genes in RRMM, which had not been previously reported or only observed at a lower frequency in newly diagnosed MM, included NRAS, BRAF, TP53, SLC4A7, MLLT4, EWSR1, HCFC2, and COPS3. We found multiple genomic regions with bi-allelic events affecting tumor suppressor genes and demonstrated a significant adverse impact of bi-allelic TP53 alterations on survival. With regard to potentially resistance conferring mutations, recurrently mutated gene networks included genes with relevance for PI and IMiD activity; the latter particularly affecting members of the Cereblon and the COP9 signalosome complex. We observed a major impact of signatures associated with exposure to melphalan or impaired DNA double-strand break homologous recombination repair in RRMM. The latter coincided with mutations in genes associated with PARP inhibitor sensitivity in 49% of RRMM patients; a finding with potential therapeutic implications. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies