The CCND1 c.870G>A polymorphism is a risk factor for t(11;14)(q13;q32) multiple myeloma

Large B-cell lymphomas with CCND1 rearrangement have different immunoglobulin gene breakpoints and genomic profile than mantle cell lymphoma

Mantle cell lymphoma (MCL) is genetically characterized by the IG::CCND1 translocation mediated by an aberrant V(D)J rearrangement. CCND1 translocations and overexpression have been identified in occasional aggressive B-cell lymphomas with unusual features for MCL. The mechanism generating CCND1 rearrangements in these tumors and their genomic profile are not known. We have reconstructed the IG::CCND1 translocations and the genomic profile of 13 SOX11-negative aggressive B-cell lymphomas using whole genome/exome and target sequencing. The mechanism behind the translocation was an aberrant V(D)J rearrangement in three tumors and by an anomalous IGH class-switch recombination (CSR) or somatic hypermutation (SHM) mechanism in ten. The tumors with a V(D)J-mediated translocation were two blastoid MCL and one high-grade B-cell lymphoma. None of them had a mutational profile suggestive of DLBCL. The ten tumors with CSR/SHM-mediated IGH::CCND1 were mainly large B-cell lymphomas, with mutated genes commonly seen in DLBCL and BCL6 rearrangements in 6. Two cases, which transformed from marginal zone lymphomas, carried mutations in KLF2, TNFAIP3 and KMT2D. These findings expand the spectrum of tumors carrying CCND1 rearrangement that may occur as a secondary event in DLBCL mediated by aberrant CSR/SHM and associated with a mutational profile different from that of MCL.

Precision medicine for multiple myeloma: The case for translocation (11;14)

The t(11;14) translocation is among the most prevalent cytogenetic abnormalities in multiple myeloma (MM), distinguished by its unique features and biology that have been thoroughly explored for decades. What further sets this MM subtype apart is its oscillating prognostic significance, from initially being considered a favorable alteration to intermediate risk and potential future reclassification as favorable risk. Despite not being inherently a high-risk alteration indicative of an aggressive phenotype, it appears that t(11;14)-MM is less responsive to novel agents like proteasome inhibitors and immunomodulatory drugs which have otherwise transformed the disease's treatment landscape, perhaps partially explained by its reduced propensity for immunoglobulin production and oligosecretory nature. However, its distinct reliance on Bcl-2 has heightened its sensitivity to venetoclax. Further subclassification based on morphological and genomic characteristics could enhance our prediction models of treatment responses and enable more tailored therapeutic strategies for patients. This review aims to encapsulate the existing research evidence in this area.

Haplotype analysis identifies functional elements in monoclonal gammopathy of unknown significance

Genome-wide association studies (GWASs) based on common single nucleotide polymorphisms (SNPs) have identified several loci associated with the risk of monoclonal gammopathy of unknown significance (MGUS), a precursor condition for multiple myeloma (MM). We hypothesized that analyzing haplotypes might be more useful than analyzing individual SNPs, as it could identify functional chromosomal units that collectively contribute to MGUS risk. To test this hypothesis, we used data from our previous GWAS on 992 MGUS cases and 2910 controls from three European populations. We identified 23 haplotypes that were associated with the risk of MGUS at the genome-wide significance level (p < 5 × 10-8) and showed consistent results among all three populations. In 10 genomic regions, strong promoter, enhancer and regulatory element-related histone marks and their connections to target genes as well as genome segmentation data supported the importance of these regions in MGUS susceptibility. Several associated haplotypes affected pathways important for MM cell survival such as ubiquitin-proteasome system (RNF186, OTUD3), PI3K/AKT/mTOR (HINT3), innate immunity (SEC14L1, ZBP1), cell death regulation (BID) and NOTCH signaling (RBPJ). These pathways are important current therapeutic targets for MM, which may highlight the advantage of the haplotype approach homing to functional units.

Deciphering the genetics and mechanisms of predisposition to multiple myeloma

Multiple myeloma (MM) is an incurable malignancy of plasma cells. Epidemiological studies indicate a substantial heritable component, but the underlying mechanisms remain unclear. Here, in a genome-wide association study totaling 10,906 cases and 366,221 controls, we identify 35 MM risk loci, 12 of which are novel. Through functional fine-mapping and Mendelian randomization, we uncover two causal mechanisms for inherited MM risk: longer telomeres; and elevated levels of B-cell maturation antigen (BCMA) and interleukin-5 receptor alpha (IL5RA) in plasma. The largest increase in BCMA and IL5RA levels is mediated by the risk variant rs34562254-A at TNFRSF13B. While individuals with loss-of-function variants in TNFRSF13B develop B-cell immunodeficiency, rs34562254-A exerts a gain-of-function effect, increasing MM risk through amplified B-cell responses. Our results represent an analysis of genetic MM predisposition, highlighting causal mechanisms contributing to MM development.

PASTRY: achieving balanced power for detecting risk and protective minor alleles in meta-analysis of association studies with overlapping subjects

BACKGROUND

Meta-analysis is a statistical method that combines the results of multiple studies to increase statistical power. When multiple studies participating in a meta-analysis utilize the same public dataset as controls, the summary statistics from these studies become correlated. To solve this challenge, Lin and Sullivan proposed a method to provide an optimal test statistic adjusted for the correlation. This method quickly became the standard practice. However, we identified an unexpected power asymmetry phenomenon in this standard framework. This can lead to unbalanced power for detecting protective minor alleles and risk minor alleles.

RESULTS

We found that the power asymmetry of the current framework is mainly due to the errors in approximating the correlation term. We then developed a meta-analysis method based on an accurate correlation estimator, called PASTRY (A method to avoid Power ASymmeTRY). PASTRY outperformed the standard method on both simulated and real datasets in terms of the power symmetry.

CONCLUSIONS

Our findings suggest that PASTRY can help to alleviate the power asymmetry problem. PASTRY is available at https://github.com/hanlab-SNU/PASTRY .

SLAMF7 as a Promising Immunotherapeutic Target in Multiple Myeloma Treatments

Multiple myeloma (MM) is a common hematological malignancy that has fostered several new therapeutic approaches to combat newly diagnosed or relapsed MM. While the field has advanced over the past 2 decades, the majority of patients will develop resistance to these treatments, causing the need for new therapeutic targets. SLAMF7 is an attractive therapeutic target in multiple myeloma, and a monoclonal antibody that targets SLAMF7 has shown consistent beneficial outcomes in clinical trials to date. In this review, we will focus on the structure and regulation of SLAMF7 and its mechanism of action. The most recent clinical trials will be reviewed to further understand the clinical implications and improve the prognosis of MM. Furthermore, the efficacy of anti-SLAMF7 monoclonal antibodies combined with standard therapies and possible resistance mechanisms will be discussed. This review aimed to provide a detailed summary of the role of SLAMF7 in the pathogenesis of patients with MM and the rationale for further investigation into SLAMF7-mediated molecular pathways associated with MM development.

D-Type Cyclins in Development and Disease

D-type cyclins encode G1/S cell cycle checkpoint proteins, which play a crucial role in defining cell cycle exit and progression. Precise control of cell cycle exit is vital during embryonic development, with defects in the pathways regulating intracellular D-type cyclins resulting in abnormal initiation of stem cell differentiation in a variety of different organ systems. Furthermore, stabilisation of D-type cyclins is observed in a wide range of disorders characterized by cellular over-proliferation, including cancers and overgrowth disorders. In this review, we will summarize and compare the roles played by each D-type cyclin during development and provide examples of how their intracellular dysregulation can be an underlying cause of disease.

Investigation of Rare Non-Coding Variants in Familial Multiple Myeloma

Multiple myeloma (MM) is a plasma cell malignancy whereby a single clone of plasma cells over-propagates in the bone marrow, resulting in the increased production of monoclonal immunoglobulin. While the complex genetic architecture of MM is well characterized, much less is known about germline variants predisposing to MM. Genome-wide sequencing approaches in MM families have started to identify rare high-penetrance coding risk alleles. In addition, genome-wide association studies have discovered several common low-penetrance risk alleles, which are mainly located in the non-coding genome. Here, we further explored the genetic basis in familial MM within the non-coding genome in whole-genome sequencing data. We prioritized and characterized 150 upstream, 5' untranslated region (UTR) and 3' UTR variants from 14 MM families, including 20 top-scoring variants. These variants confirmed previously implicated biological pathways in MM development. Most importantly, protein network and pathway enrichment analyses also identified 10 genes involved in mitogen-activated protein kinase (MAPK) signaling pathways, which have previously been established as important MM pathways.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

BACKGROUND

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

CASE PRESENTATION

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

CONCLUSIONS

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

A polygenic risk score for multiple myeloma risk prediction

There is overwhelming epidemiologic evidence that the risk of multiple myeloma (MM) has a solid genetic background. Genome-wide association studies (GWAS) have identified 23 risk loci that contribute to the genetic susceptibility of MM, but have low individual penetrance. Combining the SNPs in a polygenic risk score (PRS) is a possible approach to improve their usefulness. Using 2361 MM cases and 1415 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium, we computed a weighted and an unweighted PRS. We observed associations with MM risk with OR = 3.44, 95% CI 2.53-4.69, p = 3.55 × 10-15 for the highest vs. lowest quintile of the weighted score, and OR = 3.18, 95% CI 2.1 = 34-4.33, p = 1.62 × 10-13 for the highest vs. lowest quintile of the unweighted score. We found a convincing association of a PRS generated with 23 SNPs and risk of MM. Our work provides additional validation of previously discovered MM risk variants and of their combination into a PRS, which is a first step towards the use of genetics for risk stratification in the general population.

Functional dissection of inherited non-coding variation influencing multiple myeloma risk

Thousands of non-coding variants have been associated with increased risk of human diseases, yet the causal variants and their mechanisms-of-action remain obscure. In an integrative study combining massively parallel reporter assays (MPRA), expression analyses (eQTL, meQTL, PCHiC) and chromatin accessibility analyses in primary cells (caQTL), we investigate 1,039 variants associated with multiple myeloma (MM). We demonstrate that MM susceptibility is mediated by gene-regulatory changes in plasma cells and B-cells, and identify putative causal variants at six risk loci (SMARCD3, WAC, ELL2, CDCA7L, CEP120, and PREX1). Notably, three of these variants co-localize with significant plasma cell caQTLs, signaling the presence of causal activity at these precise genomic positions in an endogenous chromosomal context in vivo. Our results provide a systematic functional dissection of risk loci for a hematologic malignancy.

Genome Instability in Multiple Myeloma: Facts and Factors

Multiple myeloma (MM) is a malignant neoplasm of terminally differentiated immunoglobulin-producing B lymphocytes called plasma cells. MM is the second most common hematologic malignancy, and it poses a heavy economic and social burden because it remains incurable and confers a profound disability to patients. Despite current progress in MM treatment, the disease invariably recurs, even after the transplantation of autologous hematopoietic stem cells (ASCT). Biological processes leading to a pathological myeloma clone and the mechanisms of further evolution of the disease are far from complete understanding. Genetically, MM is a complex disease that demonstrates a high level of heterogeneity. Myeloma genomes carry numerous genetic changes, including structural genome variations and chromosomal gains and losses, and these changes occur in combinations with point mutations affecting various cellular pathways, including genome maintenance. MM genome instability in its extreme is manifested in mutation kataegis and complex genomic rearrangements: chromothripsis, templated insertions, and chromoplexy. Chemotherapeutic agents used to treat MM add another level of complexity because many of them exacerbate genome instability. Genome abnormalities are driver events and deciphering their mechanisms will help understand the causes of MM and play a pivotal role in developing new therapies.

Increased complexity of t(11;14) rearrangements in plasma cell neoplasms compared with mantle cell lymphoma

Plasma cell neoplasms (PCN) and mantle cell lymphoma (MCL) can both harbor t(11;14)(q13;q32) (CCND1/IGH), usually resulting in cyclin D1 overexpression. In some cases, particularly at low levels of disease, it can be morphologically challenging to distinguish between these entities in the bone marrow (BM) since PCN with t(11;14) are often CD20-positive with lymphoplasmacytic cytology, while MCL can rarely have plasmacytic differentiation. We compared the difference in CCND1/IGH by fluorescence in situ hybridization (FISH) in PCN and MCL to evaluate for possible differentiating characteristics. We identified 326 cases of MCL with t(11;14) and 279 cases of PCN with t(11;14) from either formalin-fixed, paraffin-embedded tissue or fresh BM specimens. The "typical," balanced CCND1/IGH FISH signal pattern was defined as three total CCND1 signals, three total IGH signals, and two total fusion signals. Any deviation from the "typical" pattern was defined as an "atypical" pattern, which was further stratified into "gain of fusion" vs "complex" patterns. There was a significantly higher proportion of cases that showed an atypical FISH pattern in PCN compared with MCL (53% vs 27%, P < .0001). There was also a significantly higher proportion of cases that showed a complex FISH pattern in PCN compared with MCL (47% vs 17%, P < .0001). We confirmed these findings using mate-pair sequencing of 25 PCN and MCL samples. PCN more often have a complex CCND1/IGH FISH pattern compared with MCL, suggesting possible differences in the genomic mechanisms underlying these rearrangements in plasma cells compared with B cells.

Epidemiology, genetics and treatment of multiple myeloma and precursor diseases

Multiple myeloma (MM) is a hematological malignancy caused by the clonal expansion of plasma cells. The incidence of MM worldwide is increasing with greater than 140 000 people being diagnosed with MM per year. Whereas 5-year survival after a diagnosis of MM has improved from 28% in 1975 to 56% in 2012, the disease remains essentially incurable. In this review, we summarize our current understanding of MM including its epidemiology, genetics and biology. We will also provide an overview of MM management that has led to improvements in survival, including recent changes to diagnosis and therapies. Areas of unmet need include the management of patients with high-risk MM, those with reduced performance status and those refractory to standard therapies. Ongoing research into the biology and early detection of MM as well as the development of novel therapies, such as immunotherapies, has the potential to influence MM practice in the future.

Genetic determinants of multiple myeloma risk within the Wnt/beta-catenin signaling pathway

BACKGROUND

Aberrant Wnt/beta-catenin pathway activation is implicated in Multiple Myeloma (MM) development, but little is known if genetic variants within this pathway contribute to MM susceptibility.

METHODS

We performed a discovery candidate pathway analysis in 269 non-Hispanic white MM cases and 272 controls focusing on 171 variants selected from 26 core genes within the Wnt/beta-catenin pathway. Significant candidate variants (P < 0.05) were selected for validation in internal and external non-Hispanic white populations totaling 818 cases and 1209 controls. We also examined significant variants in non-Hispanic black and Hispanic case/control study populations to identify potential differences by race/ethnicity. Possible biological functions of candidate variants were predicted in silico.

RESULTS

Seven variants were significantly associated with MM risk in non-Hispanic whites in the discovery population, of which LRP6:rs7966410 (OR: 0.57; 95 % CI: 0.38-0.88; P = 9.90 × 10^-3) and LRP6:rs7956971 (OR: 0.64; 95 % CI: 0.44-0.95; P = 0.027) remained significant in the internal and external populations. CSNK1D:rs9901910 replicated among all three racial/ethnic groups, with 2-6 fold increased risk of MM (OR: 2.40; 95 % CI: 1.67-3.45; P = 2.43 × 10^-6 - non-Hispanic white; OR: 6.42; 95 % CI: 2.47-16.7; P = 3.14 × 10^-4 - non-Hispanic black; OR: 4.31; 95 % CI: 1.83-10.1; P = 8.10 × 10^-4 - Hispanic). BTRC:rs7916830 was associated with a significant 37 % and 24 % reduced risk of MM in the non-Hispanic white (95 % CI: 0.49-0.82; P = 5.60 × 10^-4) and non-Hispanic Black (95 % CI: 0.60-0.97; P = 0.028) population, respectively. In silico tools predicted that these loci altered function through via gene regulation.

CONCLUSION

We identified several variants within the Wnt/beta-catenin pathway associated with MM susceptibility. Findings of this study highlight the potential genetic role of Wnt/beta-catenin signaling in MM etiology among a diverse patient population.

Sequencing at lymphoid neoplasm susceptibility loci maps six myeloma risk genes

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.

Duo Shared Genomic Segment analysis identifies a genome-wide significant risk locus at 18q21.33 in myeloma pedigrees

AIM

High-risk pedigrees (HRPs) are a powerful design to map highly penetrant risk genes. We previously described Shared Genomic Segment (SGS) analysis, a mapping method for single large extended pedigrees that also addresses genetic heterogeneity inherent in complex diseases. SGS identifies shared segregating chromosomal regions that may inherit in only a subset of cases. However, single large pedigrees that are individually powerful (at least 15 meioses between studied cases) are scarce. Here, we expand the SGS strategy to incorporate evidence from two extended HRPs by identifying the same segregating risk locus in both pedigrees and allowing for some relaxation in the size of each HRP.

METHODS

Duo-SGS is a procedure to combine single-pedigree SGS evidence. It implements statistically rigorous duo-pedigree thresholding to determine genome-wide significance levels that account for optimization across pedigree pairs. Single-pedigree SGS identifies optimal segments shared by case subsets at each locus across the genome, with nominal significance assessed empirically. Duo-SGS combines the statistical evidence for SGS segments at the same genomic location in two pedigrees using Fisher's method. One pedigree is paired with all others and the best duo-SGS evidence at each locus across the genome is established. Genome-wide significance thresholds are determined through distribution-fitting and the Theory of Large Deviations. We applied the duoSGS strategy to eleven extended, myeloma HRPs.

RESULTS

We identified one genome-wide significant region at 18q21.33 (0.85 Mb, P = 7.3 × 10-9) which contains one gene, CDH20. Thirteen regions were genome-wide suggestive: 1q42.2, 2p16.1, 3p25.2, 5q21.3, 5q31.1, 6q16.1, 6q26, 7q11.23, 12q24.31, 13q13.3, 18p11.22, 18q22.3 and 19p13.12.

CONCLUSION

Our results provide novel risk loci with segregating evidence from multiple HRPs and offer compelling targets and specific segment carriers to focus a future search for functional variants involved in inherited risk formyeloma.

A meta-analysis of genome-wide association studies of multiple myeloma among men and women of African ancestry

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.

Genetic Predictors of Mortality in Patients with Multiple Myeloma

Multiple myeloma (MM) is a heterogeneous disease featured by clonal plasma cell proliferation and genomic instability. The advent of next-generation sequencing allowed unraveling the complex genomic landscape of the disease. Several recurrent genomic aberrations including immunoglobulin genes translocations, copy number abnormalities, complex chromosomal events, transcriptomic and epigenomic deregulation, and mutations define various molecular subgroups with distinct outcomes. In this review, we describe the recurrent genomic events identified in MM impacting patients’ outcome and survival. These genomic aberrations constitute new markers that could be incorporated into a prognostication model to eventually guide therapy at every stage of the disease.

Genetically determined telomere length and multiple myeloma risk and outcome

Telomeres are involved in processes like cellular growth, chromosomal stability, and proper segregation to daughter cells. Telomere length measured in leukocytes (LTL) has been investigated in different cancer types, including multiple myeloma (MM). However, LTL measurement is prone to heterogeneity due to sample handling and study design (retrospective vs. prospective). LTL is genetically determined; genome-wide association studies identified 11 SNPs that, combined in a score, can be used as a genetic instrument to measure LTL and evaluate its association with MM risk. This approach has been already successfully attempted in various cancer types but never in MM. We tested the "teloscore" in 2407 MM patients and 1741 controls from the International Multiple Myeloma rESEarch (IMMeNSE) consortium. We observed an increased risk for longer genetically determined telomere length (gdTL) (OR = 1.69; 95% CI 1.36-2.11; P = 2.97 × 10^-6 for highest vs. lowest quintile of the score). Furthermore, in a subset of 1376 MM patients we tested the relationship between the teloscore and MM patients survival, observing a better prognosis for longer gdTL compared with shorter gdTL (HR = 0.93; 95% CI 0.86-0.99; P = 0.049). In conclusion, we report convincing evidence that longer gdTL is a risk marker for MM risk, and that it is potentially involved in increasing MM survival.

Discovering Genotype Variants in an Infant with VACTERL through Clinical Exome Sequencing: A Support for Personalized Risk Assessment and Disease Prevention

Congenital anomalies may have an increased risk of noncommunicable diseases (NCDs) We performed a clinical exome analysis in an infant affected by "Vertebral, Anorectal, Cardiac, Tracheoesophageal, Genitourinary, and Limb" (VACTERL) malformation association to identify potential biomarkers that may be helpful for preventing malignancy risk or other chronic processes. Among the variants, six variants that may be linked with VACTERL were identified in the exome analysis. The variants c.501G>C on OLR1 and c.-8C>G on PSMA6 were previously associated with myocardial infarction. The variants c.1936A>G on AKAP10 and c.575A>G on PON1 are linked to defects in cardiac conduction and artery disease, respectively. Alterations in metabolism were also suggested by the variants c.860G>A on EPHX2 and c.214C>A on GHRL. In addition, three variants associated with colon cancer were discovered. Specifically, the reported variants were c.723G>A on CCND1 and c.91T>A on AURKA proto-oncogenes as well as c.827A>C in the tumor suppressor PTPRJ. A further inspection identified 15 rare variants carried by cancer genes. Specifically, these mutations are located on five tumor suppressors (SDHA, RB1CC1, PTCH1, DMBT1, BCR) and eight proto-oncogenes (MERTK, CSF1R, MYB, ROS1, PCM1, FGFR2, MYH11, BRCC3) and have an allele frequency lower than 0.01 in the Genome Aggregation Database (GnomAD). We observed that the cardiac and metabolic phenotypic traits are linked with the genotype of the patient. In addition, the risk of developing neoplasia cannot be excluded a priori. Long-term surgical issues of patients with VATER syndrome could benefit from the clinical exome sequencing of a personalized risk assessment for the appearance of further disease in pubertal timing and adult age.

A Journey Through Myeloma Evolution: From the Normal Plasma Cell to Disease Complexity

The knowledge of cancer origin and the subsequent tracking of disease evolution represent unmet needs that will soon be within clinical reach. This will provide the opportunity to improve patient's stratification and to personalize treatments based on cancer biology along its life history. In this review, we focus on the molecular pathogenesis of multiple myeloma (MM), a hematologic malignancy with a well-known multi-stage disease course, where such approach can sooner translate into a clinical benefit. We describe novel insights into modes and timing of disease initiation. We dissect the biology of the preclinical and pre-malignant phases, elucidating how knowledge of the genomics of the disease and the composition of the microenvironment allow stratification of patients based on risk of disease progression. Then, we explore cell-intrinsic and cell-extrinsic drivers of MM evolution to symptomatic disease. Finally, we discuss how this may relate to the development of refractory disease after treatment. By integrating an evolutionary view of myeloma biology with the recent acquisitions on its clonal heterogeneity, we envision a way to drive the clinical management of the disease based on its detailed biological features more than surrogates of disease burden.

Mouse tumor susceptibility genes identify drug combinations for multiple myeloma

Long-term genetic studies utilizing backcross and congenic strain analyses coupled with positional cloning strategies and functional studies identified Cdkn2a, Mtor, and Mndal as mouse plasmacytoma susceptibility/resistance genes. Tumor incidence data in congenic strains carrying the resistance alleles of Cdkn2a and Mtor led us to hypothesize that drug combinations affecting these pathways are likely to have an additive, if not synergistic effect in inhibiting tumor cell growth. Traditional and novel systems-level genomic approaches were used to assess combination activity, disease specificity, and clinical potential of a drug combination involving rapamycin/everolimus, an Mtor inhibitor, with entinostat, an histone deacetylase inhibitor. The combination synergistically repressed oncogenic MYC and activated the Cdkn2a tumor suppressor. The identification of MYC as a primary upstream regulator led to the identification of small molecule binders of the G-quadruplex structure that forms in the NHEIII region of the MYC promoter. These studies highlight the importance of identifying drug combinations which simultaneously upregulate tumor suppressors and downregulate oncogenes.

Coinherited genetics of multiple myeloma and its precursor, monoclonal gammopathy of undetermined significance

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.

Genome-wide study on uveal melanoma patients finds association to DNA repair gene TDP1

Uveal melanoma is a life-threatening disease for which data on germline predisposition are essentially limited to mutations in the BAP1 gene. Many risk factors are shared between uveal melanoma and cutaneous melanoma, and these include fair skin color and light eye color. We carried out a genome-wide association study on 590 uveal melanoma patients and 5199 controls. Using a P-value limit of 10 we identified 11 loci with related odds ratios for the risk alleles ranging from 1.32 to 1.78. The smallest P-value in the overall analysis reached 1.07 × 10 for rs3759710 at 14q32.11, which is intronic to TDP1 (tyrosyl-DNA phosphodiesterase 1). This locus emerged as a genome-wide significant association for uveal melanoma clinical subtypes with any chromosomal aberrations (P = 10) and presence of epithelioid cells (P = 10). TDP1 is a DNA repair enzyme capable of repairing many types of DNA damage, including oxidative DNA lesions which may be relevant for uveal melanoma. We additionally wanted to replicate the previous candidate locus for uveal melanoma at chromosome 5p15.33 intronic to the CLPTM1L gene. Our analysis gave an odds ratio of 1.23 (95% confidence interval: 1.09-1.38; P = 0.0008) for the C allele of rs421284 and 1.21 (95% confidence interval: 1.07-1.36; P = 0.002) for the C allele of rs452932. Our data thus replicated the association of uveal melanoma with the CLPTM1L locus. Our data on TDP1 offer an attractive model positing that oxidative damage in pigmented tissue may be an initiation event in uveal melanoma and the level of damage may be regulated by the degree and type of iris pigmentation.

Next-Generation Sequencing for Clinical Management of Multiple Myeloma: Ready for Prime Time?

Personalized treatment is an attractive strategy that promises increased efficacy with reduced side effects in cancer. The feasibility of such an approach has been greatly boosted by next-generation sequencing (NGS) techniques, which can return detailed information on the genome and on the transcriptome of each patient's tumor, thus highlighting biomarkers of response or druggable targets that may differ from case to case. However, while the number of cancers sequenced is growing exponentially, much fewer cases are amenable to a molecularly-guided treatment outside of clinical trials to date. In multiple myeloma, genomic analysis shows a variety of gene mutations, aneuploidies, segmental copy-number changes, translocations that are extremely heterogeneous, and more numerous than other hematological malignancies. Currently, in routine clinical practice we employ reduced FISH panels that only capture three high-risk features as part of the R-ISS. On the contrary, recent advances have suggested that extending genomic analysis to the full spectrum of recurrent mutations and structural abnormalities in multiple myeloma may have biological and clinical implications. Furthermore, increased efficacy of novel treatments can now produce deeper responses, and standard methods do not have enough sensitivity to stratify patients in complete biochemical remission. Consequently, NGS techniques have been developed to monitor the size of the clone to a sensitivity of up to a cell in a million after treatment. However, even these techniques are not within reach of standard laboratories. In this review we will recapitulate recent advances in multiple myeloma genomics, with special focus on the ones that may have immediate translational impact. We will analyze the benefits and pitfalls of NGS-based diagnostics, highlighting crucial aspects that will need to be taken into account before this can be implemented in most laboratories. We will make the point that a new era in myeloma diagnostics and minimal residual disease monitoring is close and conventional genetic testing will not be able to return the required information. This will mandate that even in routine practice NGS should soon be adopted owing to a higher informative potential with increasing clinical benefits.

Genetic predisposition for multiple myeloma

Multiple myeloma (MM) is the second most common blood malignancy. Epidemiological family studies going back to the 1920s have provided evidence for familial aggregation, suggesting a subset of cases have an inherited genetic background. Recently, studies aimed at explaining this phenomenon have begun to provide direct evidence for genetic predisposition to MM. Genome-wide association studies have identified common risk alleles at 24 independent loci. Sequencing studies of familial cases and kindreds have begun to identify promising candidate genes where variants with strong effects on MM risk might reside. Finally, functional studies are starting to give insight into how identified risk alleles promote the development of MM. Here, we review recent findings in MM predisposition field, and highlight open questions and future directions.

Meta-Analysis of Genome-Wide Association Studies Identifies Three Loci Associated With Stiffness Index of the Calcaneus

The stiffness index (SI) from quantitative ultrasound measurements is a good indicator of BMD and may be used to predict the risk of osteoporotic fracture. We conducted a genomewide association study (GWAS) for SI using 7742 individuals from the Taiwan Biobank, followed by a replication study in a Korean population (n = 2955). Approximately 6.1 million SNPs were subjected to association analysis, and SI-associated variants were identified. We further conducted a meta-analysis of Taiwan Biobank significant SNPs with a Korean population-based cohort. Candidate genes were prioritized according to epigenetic annotations, gene ontology, protein-protein interaction, GWAS catalog, and expression quantitative trait loci analyses. Our results revealed seven significant single-nucleotide polymorphisms (SNPs) within three loci: 7q31.31, 17p13.3, and 11q14.2. Conditional analysis showed that three SNPs, rs2536195 (CPED1/WNT16), rs1231207 (SMG6), and rs4944661 (LOC10050636/TMEM135), were the most important signals within these regions. The associations for the three SNPs were confirmed in a UK Biobank estimated BMD GWAS; these three cytobands were replicated successfully after a meta-analysis with a Korean population cohort as well. However, two SNPs were not replicated. After prioritization, we identified two novel genes, RAB15 and FNTB, as strong candidates for association with SI. Our study identified three SI-associated SNPs and two novel SI-related genes. Overall, these results provide further insight into the genetic architecture of osteoporosis. Further studies in larger East Asian populations are needed. © 2019 American Society for Bone and Mineral Research.

Germline Risk Contribution to Genomic Instability in Multiple Myeloma

Genomic instability, a well-established hallmark of human cancer, is also a driving force in the natural history of multiple myeloma (MM) - a difficult to treat and in most cases fatal neoplasm of immunoglobulin producing plasma cells that reside in the hematopoietic bone marrow. Long recognized manifestations of genomic instability in myeloma at the cytogenetic level include abnormal chromosome numbers (aneuploidy) caused by trisomy of odd-numbered chromosomes; recurrent oncogene-activating chromosomal translocations that involve immunoglobulin loci; and large-scale amplifications, inversions, and insertions/deletions (indels) of genetic material. Catastrophic genetic rearrangements that either shatter and illegitimately reassemble a single chromosome (chromotripsis) or lead to disordered segmental rearrangements of multiple chromosomes (chromoplexy) also occur. Genomic instability at the nucleotide level results in base substitution mutations and small indels that affect both the coding and non-coding genome. Sometimes this generates a distinctive signature of somatic mutations that can be attributed to defects in DNA repair pathways, the DNA damage response (DDR) or aberrant activity of mutator genes including members of the APOBEC family. In addition to myeloma development and progression, genomic instability promotes acquisition of drug resistance in patients with myeloma. Here we review recent findings on the genetic predisposition to myeloma, including newly identified candidate genes suggesting linkage of germline risk and compromised genomic stability control. The role of ethnic and familial risk factors for myeloma is highlighted. We address current research gaps that concern the lack of studies on the mechanism by which germline risk alleles promote genomic instability in myeloma, including the open question whether genetic modifiers of myeloma development act in tumor cells, the tumor microenvironment (TME), or in both. We conclude with a brief proposition for future research directions, which concentrate on the biological function of myeloma risk and genetic instability alleles, the potential links between the germline genome and somatic changes in myeloma, and the need to elucidate genetic modifiers in the TME.

The Impact of Tumor Heterogeneity on Diagnostics and Novel Therapeutic Strategies in Multiple Myeloma

Myeloma is characterized by extensive inter-patient genomic heterogeneity due to multiple different initiating events. A recent multi-region sequencing study demonstrated spatial differences, with progression events, such as TP53 mutations, frequently being restricted to focal lesions. In this review article, we describe the clinical impact of these two types of tumor heterogeneity. Target mutations are often dominant at one site but absent at other sites, which poses a significant challenge to personalized therapy in myeloma. The same holds true for high-risk subclones, which can be locally restricted, and as such not detectable at the iliac crest, which is the usual sampling site. Imaging can improve current risk classifiers and monitoring of residual disease, but does not allow for deciphering the molecular characteristics of tumor clones. In the era of novel immunotherapies, the clinical impact of heterogeneity certainly needs to be re-defined. Yet, preliminary observations indicate an ongoing impact of spatial heterogeneity on the efficacy of monoclonal antibodies. In conclusion, we recommend combining molecular tests with imaging to improve risk prediction and monitoring of residual disease. Overcoming intra-tumor heterogeneity is the prerequisite for curing myeloma. Novel immunotherapies are promising but research addressing their impact on the spatial clonal architecture is highly warranted.

Genetic polymorphisms in genes of class switch recombination and multiple myeloma risk and survival: an IMMEnSE study

Genetic variants in genes acting during the maturation process of immature B-cell to differentiated plasma cell could influence the risk of developing multiple myeloma (MM). During B-cell maturation, several programmed genetic rearrangements occur to increase the variation of the immunoglobulin chains. Class switch recombination (CSR) is one of the most important among these mechanisms. Germline polymorphisms altering even subtly this process could play a role in the etiology and outcome of MM. We performed an association study of 30 genetic variants in the key CSR genes, using 2632 MM patients and 2848 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium, the Heidelberg MM Group and the ESTHER cohort. We found an association between LIG4-rs1555902 and decreased MM risk, which approached statistical significance, as well as significant associations between AICDA-rs3794318 and better outcome. Our results add to our knowledge on the genetic component of MM risk and survival.

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.

Risk of MGUS in relatives of multiple myeloma cases by clinical and tumor characteristics

We and others have shown increased risk of monoclonal gammopathy of undetermined significance (MGUS) in first-degree relatives of patients with multiple myeloma (MM). Whether familial risk of MGUS differs by the MM proband's age at onset, tumor or clinical characteristics is unknown. MM and smoldering MM (SMM) cases (N = 430) were recruited from the Mayo Clinic in Rochester, Minnesota between 2005-2015. First-degree relatives over age 40 provided serum samples for evaluation of MGUS (N = 1179). Age and sex specific rates of MGUS among first-degree relatives were compared to a population-based sample. Cytogenetic subtypes were classified by Fluorescence in situ hybridization. MGUS was detected in 75 first-degree relatives for an age- and sex- adjusted prevalence of 5.8% (95% CI: 4.5-7.2). Prevalence of MGUS in first-degree relatives was 2.4 fold (95% CI: 1.9-2.9) greater than expected rates. Familial risk did not differ by proband's age at diagnosis, gender, isotype, IgH translocation, or trisomy. This study confirms first-degree relatives of MM cases have a significantly higher risk of MGUS compared to the general population, regardless of age, gender, or tumor characteristics. In selected situations, such as multiple affected first-degree relatives, screening of first-degree relatives of MM cases could be considered for follow-up and prevention strategies.

Genetic variation of acquired structural chromosomal aberrations

Human malignancies are often hallmarked with genomic instability, which itself is also considered a causative event in malignant transformation. Genomic instability may manifest itself as genetic changes in the nucleotide sequence of DNA, or as structural or numerical changes of chromosomes. Unrepaired or insufficiently repaired DNA double-strand breaks, as well as telomere shortening, are important contributors in the formation of structural chromosomal aberrations (CAs). In the present review, we discuss potential mechanisms behind the formation of CAs and their relation to cancer. Based on our own studies, we also illustrate how inherited genetic variation may modify the frequency and types of CAs occurring in humans. Recently, we published a series of studies on variations in genes relevant to maintaining genomic integrity, such as those encoding xenobiotic-metabolising enzymes, DNA repair, the tumour suppressor TP53, the spindle assembly checkpoint, and cyclin D1 (CCND1). While individually genetic variation in these genes exerted small modulating effects, in interactions they were associated with CA frequencies in peripheral blood lymphocytes of healthy volunteers. Moreover, we observed opposite associations between the CCND1 splice site polymorphism rs9344 G870A and the frequency of CAs compared to their association with translocation t(11,14). We discuss the functional consequences of the CCND1 gene in interplay with DNA damage response and DNA repair during malignant transformation. Our review summarizes existing evidence that gene variations in relevant cellular pathways modulate the frequency of CAs, predominantly in a complex interaction. More functional/mechanistic studies elucidating these observations are required. Several questions emerge, such as the role of CAs in malignancies with respect to a particular phenotype and heterogeneity, the formation of CAs during the process of malignant transformation, and the formation of CAs in individual types of lymphocytes in relation to the immune response.

Genetic Predisposition to Multiple Myeloma at 5q15 Is Mediated by an ELL2 Enhancer Polymorphism

Multiple myeloma (MM) is a malignancy of plasma cells. Genome-wide association studies have shown that variation at 5q15 influences MM risk. Here, we have sought to decipher the causal variant at 5q15 and the mechanism by which it influences tumorigenesis. We show that rs6877329 G > C resides in a predicted enhancer element that physically interacts with the transcription start site of ELL2. The rs6877329-C risk allele is associated with reduced enhancer activity and lowered ELL2 expression. Since ELL2 is critical to the B cell differentiation process, reduced ELL2 expression is consistent with inherited genetic variation contributing to arrest of plasma cell development, facilitating MM clonal expansion. These data provide evidence for a biological mechanism underlying a hereditary risk of MM at 5q15.

Search for rare protein altering variants influencing susceptibility to multiple myeloma

The genetic basis underlying the inherited risk of developing multiple myeloma (MM) is largely unknown. To examine the impact of rare protein altering variants on the risk of developing MM we analyzed high-coverage exome sequencing data on 513 MM cases and 1,569 healthy controls, performing both single variant and gene burden tests. We did not identify any recurrent coding low-frequency alleles (1–5%) with moderate effect that were statistically associated with MM. In a gene burden analysis we did however identify a promising relationship between variation in the marrow kinetochore microtubule stromal gene KIF18A, which plays a role in control mitotic chromosome positioning dynamics, and risk of MM (P =3.6×10⁻⁶). Further analysis showed KIF18A displays a distinct pattern of expression across molecular subgroups of MM as well as being associated with patient survival. Our results inform future study design and provide a resource for contextualizing the impact of candidate MM susceptibility genes.

Novel pedigree analysis implicates DNA repair and chromatin remodeling in multiple myeloma risk

The high-risk pedigree (HRP) design is an established strategy to discover rare, highly-penetrant, Mendelian-like causal variants. Its success, however, in complex traits has been modest, largely due to challenges of genetic heterogeneity and complex inheritance models. We describe a HRP strategy that addresses intra-familial heterogeneity, and identifies inherited segments important for mapping regulatory risk. We apply this new Shared Genomic Segment (SGS) method in 11 extended, Utah, multiple myeloma (MM) HRPs, and subsequent exome sequencing in SGS regions of interest in 1063 MM / MGUS (monoclonal gammopathy of undetermined significance-a precursor to MM) cases and 964 controls from a jointly-called collaborative resource, including cases from the initial 11 HRPs. One genome-wide significant 1.8 Mb shared segment was found at 6q16. Exome sequencing in this region revealed predicted deleterious variants in USP45 (p.Gln691* and p.Gln621Glu), a gene known to influence DNA repair through endonuclease regulation. Additionally, a 1.2 Mb segment at 1p36.11 is inherited in two Utah HRPs, with coding variants identified in ARID1A (p.Ser90Gly and p.Met890Val), a key gene in the SWI/SNF chromatin remodeling complex. Our results provide compelling statistical and genetic evidence for segregating risk variants for MM. In addition, we demonstrate a novel strategy to use large HRPs for risk-variant discovery more generally in complex traits.

Genome-Wide Association Studies in Glioma

Since the first reports in 2009, genome-wide association studies (GWAS) have been successful in identifying germline variants associated with glioma susceptibility. In this review, we describe a chronological history of glioma GWAS, culminating in the most recent study comprising 12,496 cases and 18,190 controls. We additionally summarize associations at the 27 glioma-risk SNPs that have been reported so far. Future efforts are likely to be principally focused on assessing association of germline-risk SNPs with particular molecular subgroups of glioma, as well as investigating the functional basis of the risk loci in tumor formation. These ongoing studies will be important to maximize the impact of research into glioma susceptibility, both in terms of insight into tumor etiology as well as opportunities for clinical translation. Cancer Epidemiol Biomarkers Prev; 27(4); 418-28. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

A common variant within the HNF1B gene is associated with overall survival of multiple myeloma patients: results from the IMMEnSE consortium and meta-analysis

Diabetogenic single nucleotide polymorphisms (SNPs) have recently been associated with multiple myeloma (MM) risk but their impact on overall survival (OS) of MM patients has not been analysed yet. In order to investigate the impact of 58 GWAS-identified variants for type 2 diabetes (T2D) on OS of patients with MM, we analysed genotyping data of 936 MM patients collected by the International Multiple Myeloma rESEarch (IMMENSE) consortium and an independent set of 700 MM patients recruited by the University Clinic of Heidelberg. A meta-analysis of the cox regression results of the two sets showed that rs7501939 located in the HNF1B gene negatively impacted OS (HRRec= 1.44, 95% CI = 1.18-1.76, P = 0.0001). The meta-analysis also showed a noteworthy gender-specific association of the SLC30A8rs13266634 SNP with OS. The presence of each additional copy of the minor allele at rs13266634 was associated with poor OS in men whereas no association was seen in women (HRMen-Add = 1.32, 95% CI 1.13-1.54, P = 0.0003). In conclusion, these data suggest that the HNF1Brs7501939 SNP confers poor OS in patients with MM and that a SNP in SLC30A8 affect OS in men.

Profile of common prostate cancer risk variants in an unscreened Romanian population

To find sequence variants affecting prostate cancer (PCA) susceptibility in an unscreened Romanian population we use a genome‐wide association study (GWAS). The study population included 990 unrelated pathologically confirmed PCA cases and 1034 male controls. DNA was genotyped using Illumina SNP arrays, and 24.295.558 variants were imputed using the 1000 Genomes data set. An association test was performed between the imputed markers and PCA. A systematic literature review for variants associated with PCA risk identified 115 unique variants that were tested in the Romanian sample set. Thirty of the previously reported SNPs replicated (P‐value < 0.05), with the strongest associations observed at: 8q24.21, 11q13.3, 6q25.3, 5p15.33, 22q13.2, 17q12 and 3q13.2. The replicated variants showing the most significant association in Romania are rs1016343 at 8q24.21 (P = 2.2 × 10⁻⁴), rs7929962 at 11q13.3 (P = 2.7 × 10⁻⁴) and rs9364554 at 6q25.2 (P = 4.7 × 10⁻⁴). None of the variants tested in the Romanian GWAS reached genome‐wide significance (P‐value <5 × 10⁻⁸) but 807 markers had P‐values <1 × 10⁻⁴. Here, we report the results of the first GWAS of PCA performed in a Romanian population. Our study provides evidence that a substantial fraction of previously validated PCA variants associate with risk in this unscreened Romanian population.

FOLD: a method to optimize power in meta-analysis of genetic association studies with overlapping subjects

Motivation

In genetic association studies, meta-analyses are widely used to increase the statistical power by aggregating information from multiple studies. In meta-analyses, participating studies often share the same individuals due to the shared use of publicly available control data or accidental recruiting of the same subjects. As such overlapping can inflate false positive rate, overlapping subjects are traditionally split in the studies prior to meta-analysis, which requires access to genotype data and is not always possible. Fortunately, recently developed meta-analysis methods can systematically account for overlapping subjects at the summary statistics level.

Results

We identify and report a phenomenon that these methods for overlapping subjects can yield low power. For instance, in our simulation involving a meta-analysis of five studies that share 20% of individuals, whereas the traditional splitting method achieved 80% power, none of the new methods exceeded 32% power. We found that this low power resulted from the unaccounted differences between shared and unshared individuals in terms of their contributions towards the final statistic. Here, we propose an optimal summary-statistic-based method termed as FOLD that increases the power of meta-analysis involving studies with overlapping subjects.

Availability and implementation

Our method is available at http://software.buhmhan.com/FOLD.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-wide association studies of cancer: current insights and future perspectives

Genome-wide association studies (GWAS) provide an agnostic approach for investigating the genetic basis of complex diseases. In oncology, GWAS of nearly all common malignancies have been performed, and over 450 genetic variants associated with increased risks have been identified. As well as revealing novel pathways important in carcinogenesis, these studies have shown that common genetic variation contributes substantially to the heritable risk of many common cancers. The clinical application of GWAS is starting to provide opportunities for drug discovery and repositioning as well as for cancer prevention. However, deciphering the functional and biological basis of associations is challenging and is in part a barrier to fully unlocking the potential of GWAS.

Chromosomal instability and acquired drug resistance in multiple myeloma

Chromosomal instability (CIN) is an important hallmark of human cancer. CIN not only contributes to all stages of tumor development (initiation, promotion and progression) but also drives, in large measure, the acquisition of drug resistance by cancer cells. Although CIN is a cornerstone of the complex mutational architecture that underlies neoplastic cell development and tumor heterogeneity and has been tightly associated with treatment responses and survival of cancer patients, it may be one of the least understood features of the malignant phenotype in terms of genetic pathways and molecular mechanisms. Here we review new insights into the type of CIN seen in multiple myeloma (MM), a blood cancer of terminally differentiated, immunoglobulin-producing B-lymphocytes called plasma cells that remains incurable in the great majority of cases. We will consider bona fide myeloma CIN genes, methods for measuring CIN in myeloma cells, and novel approaches to CIN-targeted treatments of patients with myeloma. The new findings generate optimism that enhanced understanding of CIN will lead to the design and testing of new therapeutic strategies to overcome drug resistance in MM in the not-so-distant future.

Evolutionary biology of high-risk multiple myeloma

The outcomes for the majority of patients with myeloma have improved over recent decades, driven by treatment advances. However, there is a subset of patients considered to have high-risk disease who have not benefited. Understanding how high-risk disease evolves from more therapeutically tractable stages is crucial if we are to improve outcomes. This can be accomplished by identifying the genetic mechanisms and mutations driving the transition of a normal plasma cell to one with the features of the following disease stages: monoclonal gammopathy of undetermined significance, smouldering myeloma, myeloma and plasma cell leukaemia. Although myeloma initiating events are clonal, subsequent driver lesions often occur in a subclone of cells, facilitating progression by Darwinian selection processes. Understanding the co-evolution of the clones within their microenvironment will be crucial for therapeutically manipulating the process. The end stage of progression is the generation of a state associated with treatment resistance, increased proliferation, evasion of apoptosis and an ability to grow independently of the bone marrow microenvironment. In this Review, we discuss these end-stage high-risk disease states and how new information is improving our understanding of their evolutionary trajectories, how they may be diagnosed and the biological behaviour that must be addressed if they are to be treated effectively.

Direct evidence for a polygenic etiology in familial multiple myeloma

Although common risk alleles for multiple myeloma (MM) were recently identified, their contribution to familial MM is unknown. Analyzing 38 familial cases identified primarily by linking Swedish nationwide registries, we demonstrate an enrichment of common MM risk alleles in familial compared with 1530 sporadic cases (P = 4.8 × 10^-2 and 6.0 × 10^-2, respectively, for 2 different polygenic risk scores) and 10 171 population-based controls (P = 1.5 × 10^-4 and 1.3 × 10^-4, respectively). Using mixture modeling, we estimate that about one-third of familial cases result from such enrichments. Our results provide the first direct evidence for a polygenic etiology in a familial hematologic malignancy.

Bortezomib-induced peripheral neuropathy: A genome-wide association study on multiple myeloma patients

The proteasome-inhibitor bortezomib was introduced into the treatment of multiple myeloma more than a decade ago. It is clinically beneficial, but peripheral neuropathy (PNP) is a side effect that may limit its use in some patients. To examine the possible genetic predisposing factors to PNP, we performed a genome-wide association study on 646 bortezomib-treated German multiple myeloma patients. Our aim was to identify genetic risk variants associated with the development of PNP as a serious side effect of the treatment. We identified 4 new promising loci for bortezomib-induced PNP at 4q34.3 (rs6552496), 5q14.1 (rs12521798), 16q23.3 (rs8060632), and 18q21.2 (rs17748074). Even though the results did not reach genome-wide significance level, they support the idea of previous studies, suggesting a genetic basis for neurotoxicity. The identified single nucleotide polymorphisms map to genes or next to genes involved in the development and function of the nervous system (CDH13, DCC, and TENM3). As possible functional clues, 2 of the variants, rs12521798 and rs17748074, affect enhancer histone marks in the brain. The rs12521798 may also impact expression of THBS4, which affects specific signal trasduction pathways in the nervous system. Further research is needed to clarify the mechanism of action of the identified single nucleotide polymorphisms in the development of drug-induced PNP and to functionally validate our in silico predictions.

NCOA1 is a novel susceptibility gene for multiple myeloma in the Chinese population: A case-control study

Multiple myeloma (MM) is an incurable malignancy of mature B-lymphoid cells, and its pathogenesis is only partially understood. Previous studies have demonstrated that a number of Non-Hodgkin Lymphoma (NHL) associated genes also show susceptibility to MM, suggesting malignancies originating from B cells may share similar genetic susceptibility. Several recent large-scale genome-wide association studies (GWAS) have identified HLA-I, HLA-II, CXCR5, ETS1, LPP and NCOA1 genes as genetic risk factors associated with NHL, and this study aimed to investigate whether these genes polymorphisms confer susceptibility with MM in the Chinese Han population. In 827 MM cases and 709 healthy controls of Chinese Han, seven single nucleotide polymorphisms (SNPs) in the HLA–I region (rs6457327), the HLA–II region (rs2647012 and rs7755224), the CXCR5 gene (rs4938573), the ETS1 gene (rs4937362), the LPP gene (rs6444305), and the NCOA1 region (rs79480871) were genotyped using the Sequenom platform. Our study indicated that genotype and allele frequencies of rs79480871 showed strong associations with MM patients (p_a = 3.5×10⁻⁴ and p_a = 1.5×10⁻⁴), and the rs6457327 genotype was more readily associated with MM patients than with controls (p_a = 4.9×10⁻³). This study was the first to reveal the correlation between NCOA1 gene polymorphisms and MM patients, indicating that NCOA1 might be a novel susceptibility gene for MM patients in the Chinese Han population.