Variants in ELL2 influencing immunoglobulin levels associate with multiple myeloma

FaMMily Affairs: Dissecting inherited contributions to multiple myeloma risk

Etiological links to multiple myeloma (MM) remain poorly understood, though emerging evidence suggests a significant hereditary component. This review integrates current literature on inherited factors contributing to MM risk, synthesizing both epidemiologic and genomic data. We examine familial clustering patterns, assess genome-wide association studies (GWAS) that reveal common genetic variants linked to MM, and explore rare, high-penetrance variants in key susceptibility genes. Additionally, we advocate for routine germline screening in high-risk MM populations, particularly those with a strong family history of cancer, a personal history of cancer, or early-onset disease. By elucidating the inherited influences on MM predisposition, this review seeks to inform future research and refine risk assessment strategies in this population.

Germline predisposition in multiple myeloma

We present a study of rare germline predisposition variants in 954 unrelated individuals with multiple myeloma (MM) and 82 MM families. Using a candidate gene approach, we identified such variants across all age groups in 9.1% of sporadic and 18% of familial cases. Implicated genes included genes suggested in other MM risk studies as potential risk genes (DIS3, EP300, KDM1A, and USP45); genes involved in predisposition to other cancers (ATM, BRCA1/2, CHEK2, PMS2, POT1, PRF1, and TP53); and BRIP1, EP300, and FANCM in individuals of African ancestry. Variants were characterized using loss of heterozygosity (LOH), biallelic events, and gene expression analyses, revealing 31 variants in 3.25% of sporadic cases for which pathogenicity was supported by multiple lines of evidence. Our results suggest that the disruption of DNA damage repair pathways may play a role in MM susceptibility. These results will inform improved surveillance in high-risk groups and potential therapeutic strategies.

Developing and Verifying an Effective Diagnostic Model Linked to Immune Infiltration in Stanford Type A Aortic Dissection

BACKGROUND

The deadly cardiovascular condition known as Stanford type A aortic dissection (TAAD) carries a high risk of morbidity and mortality. One important step in the pathophysiology of the condition is the influx of immune cells into the aorta media, which causes medial degeneration. The purpose of this work was to investigate the potential pathogenic significance of immune cell infiltration in TAAD and to test for associated biomarkers.

METHODS

The National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database provided the RNA sequencing microarray data (GSE153434, GPL20795, GSE52093). Immune cell infiltration abundance was predicted using ImmuCellAI. GEO2R was used to select differentially expressed genes (DEGs), which were then processed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Additionally, hub genes linked to immune infiltration were found using functional and pathway enrichment, least absolute shrinkage and selection operator (LASSO), weighted gene co-expression network analysis (WGCNA), and differential expression analysis. Lastly, hub genes were validated and assessed using receiver operating characteristic (ROC) curves in the microarray dataset GSE52093. The hub gene expression and its connection to immune infiltration in TAAD were confirmed using both animal models and clinic data.

RESULTS

We identified the most important connections between macrophages, T helper cell 17 (Th17), iTreg cells, B cells, natural killer cells and TAAD. And screened seven hub genes associated with immune cell infiltration: ABCG2, FAM20C, ELL2, MTHFD2, ANKRD6, GLRX, and CDCP1. The diagnostic model in TAAD diagnosis with the area under ROC (AUC) was 0.996, and the sensitivity was 99.21%, the specificity was 98.67%, which demonstrated a surprisingly strong diagnostic power of TAAD in the validation datasets. The expression pattern of four hub DEGs (ABCG2, FAM20C, MTHFD2, CDCP1) in clinic samples and animal models matched bioinformatics analysis, and ABCG2, FAM20C, MTHFD2 up-regulation, and the of CDCP1 down-regulation were also linked to poor cardiovascular function.

CONCLUSIONS

This study developed and verified an effective diagnostic model linked to immune infiltration in TAAD, providing new approaches to studying the potential pathogenesis of TAAD and discovering new medication intervention targets.

The genetics and epidemiology of N- and O-immunoglobulin A glycomics

BACKGROUND

Immunoglobulin (Ig) glycosylation modulates the immune response and plays a critical role in ageing and diseases. Studies have mainly focused on IgG glycosylation, and little is known about the genetics and epidemiology of IgA glycosylation.

METHODS

We generated, using a novel liquid chromatography-mass spectrometry method, the first large-scale IgA glycomics dataset in serum from 2423 twins, encompassing 71 N- and O-glycan species.

RESULTS

We showed that, despite the lack of a direct genetic template, glycosylation is highly heritable, and that glycopeptide structures are sex-specific, and undergo substantial changes with ageing. We observe extensive correlations between the IgA and IgG glycomes, and, exploiting the twin design, show that they are predominantly influenced by shared genetic factors. A genome-wide association study identified eight loci associated with both the IgA and IgG glycomes (ST6GAL1, ELL2, B4GALT1, ABCF2, TMEM121, SLC38A10, SMARCB1, and MGAT3) and two novel loci specifically modulating IgA O-glycosylation (C1GALT1 and ST3GAL1). Validation of our findings in an independent cohort of 320 individuals from Qatar showed that the underlying genetic architecture is conserved across ancestries.

CONCLUSIONS

Our study delineates the genetic landscape of IgA glycosylation and provides novel potential functional links with the aetiology of complex immune diseases, including genetic factors involved in IgA nephropathy risk.

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.

Screening in Multiple Myeloma and Its Precursors: Are We There Yet?

BACKGROUND

Multiple myeloma (MM) is a hematological malignancy that develops over years from the asymptomatic precursors, monoclonal gammopathy of undetermined significance, and smoldering multiple myeloma. Recent evidence shows that by initiating treatment at an asymptomatic stage, outcomes in MM can be significantly improved. However, a vast majority of MM patients are diagnosed after the development of symptomatic end-organ damage and cannot reap the benefits of early treatment. The precursors of MM are easily detected by serum protein electrophoresis and free light chain assay of the serum, raising the question of whether population-based screening could detect MM at an asymptomatic stage and significantly expand the availability of early treatment in MM. Screening is a hallmark of care in many malignancies, and there are accepted criteria for when screening is appropriate.

CONTENT

Here we review the available relevant evidence for the introduction of screening and discuss whether screening for MM and its precursors fulfills these criteria. We also highlight gaps in our current knowledge, most notably a lack of data on the benefits and harms of screening and the lack of a defined target population. There are ongoing studies that may fill these critical gaps in the literature, but their results are still pending.

SUMMARY

Screening could lead to a paradigm shift in the care of patients with MM, but critical scientific questions need to be answered before screening of healthy individuals can be recommended. In short, we should not screen for MM and its precursors-yet.

Design and quality control of large-scale two-sample Mendelian randomization studies

Background

Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors.

Methods

We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set.

Results

We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls.

Conclusions

In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk

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

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.

Genetic regulation of serum IgA levels and susceptibility to common immune, infectious, kidney, and cardio-metabolic traits

Immunoglobulin A (IgA) mediates mucosal responses to food antigens and the intestinal microbiome and is involved in susceptibility to mucosal pathogens, celiac disease, inflammatory bowel disease, and IgA nephropathy. We performed a genome-wide association study of serum IgA levels in 41,263 individuals of diverse ancestries and identified 20 genome-wide significant loci, including 9 known and 11 novel loci. Co-localization analyses with expression QTLs prioritized candidate genes for 14 of 20 significant loci. Most loci encoded genes that produced immune defects and IgA abnormalities when genetically manipulated in mice. We also observed positive genetic correlations of serum IgA levels with IgA nephropathy, type 2 diabetes, and body mass index, and negative correlations with celiac disease, inflammatory bowel disease, and several infections. Mendelian randomization supported elevated serum IgA as a causal factor in IgA nephropathy. African ancestry was consistently associated with higher serum IgA levels and greater frequency of IgA-increasing alleles compared to other ancestries. Our findings provide novel insights into the genetic regulation of IgA levels and its potential role in human disease.

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.

Multiple Myeloma: Bioinformatic Analysis for Identification of Key Genes and Pathways

Multiple myeloma (MM) is a hematological malignancy in which monoclonal plasma cells multiply in the bone marrow and monoclonal immunoglobulins are overproduced in older people. Several molecular and cytogenetic advances allow scientists to identify several genetic and chromosomal abnormalities that cause the disease. The comprehension of the pathophysiology of MM requires an understanding of the characteristics of malignant clones and the changes in the bone marrow microenvironment. This study aims to identify the central genes and to determine the key signaling pathways in MM by in silico approaches. A list of 114 differentially expressed genes (DEGs) is important in the prognosis of MM. The DEGs are collected from scientific publications and databases (https://www.ncbi.nlm.nih.gov/). These data are analyzed by Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) software (https://string-db.org/) through the construction of protein-protein interaction (PPI) networks and enrichment analysis of the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, by CytoHubba, AutoAnnotate, Bingo Apps plugins in Cytoscape software (https://cytoscape.org/) and by DAVID database (https://david.ncifcrf.gov/). The analysis of the results shows that there are 7 core genes, including TP53; MYC; CDND1; IL6; UBA52; EZH2, and MDM2. These top genes appear to play a role in the promotion and progression of MM. According to functional enrichment analysis, these genes are mainly involved in the following signaling pathways: Epstein-Barr virus infection, microRNA pathway, PI3K-Akt signaling pathway, and p53 signaling pathway. Several crucial genes, including TP53, MYC, CDND1, IL6, UBA52, EZH2, and MDM2, are significantly correlated with MM, which may exert their role in the onset and evolution of MM.

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.

A genetic risk score of alleles related to MGUS interacts with socioeconomic position in a population-based cohort

Environmental, genetic, and social factors are suggested to jointly influence monoclonal gammopathy of undetermined significance (MGUS), a precursor of multiple myeloma. Aim of this study was to investigate interactions between MGUS-related genetic variants and socioeconomic position (SEP) indicators education and income on MGUS in a population-based study. Two different MGUS-related genetic risk allele sum scores (GRS) were calculated based on recent genome-wide meta-analyses. Odds Ratios (OR) were estimated in 4329 participants including 238 MGUS cases to assess associations and multiplicative interaction. The relative excess risk due to interaction (RERI) was calculated to assess additive interaction. Both GRSs were associated with MGUS. A multiplicative interaction between one GRS and education was observed with genetic effects of OR 1.34 (95% CI 1.11–1.62) per risk allele in the highest and OR 1.06 (95% CI 0.86–1.31) in the lowest education group. A RERI of 0.10 (95% CI 0.05–0.14) also indicated additive interaction. Further, additive GRS by income interaction (RERI 0.07; 95% CI 0.01–0.13) for the same GRS was also indicated. Results indicate interaction between MGUS-related genetic risk and SEP. Non-genetic MGUS risk factors more common in higher education groups may influence the expression of MGUS-related genetic variants.

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.

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.

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.

The Significance of mRNA in the Biology of Multiple Myeloma and Its Clinical Implications

Multiple myeloma (MM) is a genetically complex disease that results from a multistep transformation of normal to malignant plasma cells in the bone marrow. However, the molecular mechanisms responsible for the initiation and heterogeneous evolution of MM remain largely unknown. A fundamental step needed to understand the oncogenesis of MM and its response to therapy is the identification of driver mutations. The introduction of gene expression profiling (GEP) in MM is an important step in elucidating the molecular heterogeneity of MM and its clinical relevance. Since some mutations in myeloma occur in non-coding regions, studies based on the analysis of mRNA provide more comprehensive information on the oncogenic pathways and mechanisms relevant to MM biology. In this review, we discuss the role of gene expression profiling in understanding the biology of multiple myeloma together with the clinical manifestation of the disease, as well as its impact on treatment decisions and future directions.

The Norwegian childhood cancer biobank

Background

The rapidly expanding era of “omics” research is highly dependent on the availability of quality‐proven biological material, especially for rare conditions such as pediatric malignancies. Professional biobanks provide such material, focusing on standardized collection and handling procedures, distinctive quality measurements, traceability of storage conditions, and accessibility. For pediatric malignancies, traditional tumor biobanking is challenging due to the rareness and limited amount of tissue and blood samples. The higher molecular heterogeneity, lower mutation rates, and unique genomic landscapes, however, renders biobanking of this tissue even more crucial.

Aim

The aim of this study was to test and establish methods for a prospective and centralized biobank for infants, children, and adolescents up to 18 years of age diagnosed with cancer in Norway.

Methods

Obtain judicial and ethical approvals and administration through a consortium, steering committee, and advisory board. Develop pipelines including SOPs for all aspects in the biobank process, including collection, processing and storing of samples and data, as well of quality controlling, safeguarding, distributing, and transport.

Results

The childhood cancer biobanking started at Oslo University Hospital in March 2017 and was from 2019 run as a national Norwegian Childhood Cancer Biobank. Informed consent and biological samples are collected regionally and stored centrally. Approximately 12 000 samples from 510 patients and have been included by January 1, 2021, representing a 96% consent and participation rate among our newly diagnosed patients.

Conclusion

A well‐functioning nationwide collection and centralized biobank with standardized procedures and national storage for pediatric malignancies has been established with a high acceptance among families.

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.

Search for multiple myeloma risk factors using Mendelian randomization

The etiology of multiple myeloma (MM) is poorly understood. Summary data from genome-wide association studies (GWASs) of multiple phenotypes can be exploited in a Mendelian randomization (MR) phenome-wide association study (PheWAS) to search for factors influencing MM risk. We performed an MR-PheWAS analyzing 249 phenotypes, proxied by 10 225 genetic variants, and summary genetic data from a GWAS of 7717 MM cases and 29 304 controls. Odds ratios (ORs) per 1 standard deviation increase in each phenotype were estimated under an inverse variance weighted random effects model. A Bonferroni-corrected threshold of P = 2 × 10-4 was considered significant, whereas P < .05 was considered suggestive of an association. Although no significant associations with MM risk were observed among the 249 phenotypes, 28 phenotypes showed evidence suggestive of association, including increased levels of serum vitamin B6 and blood carnitine (P = 1.1 × 10-3) with greater MM risk and ω-3 fatty acids (P = 5.4 × 10-4) with reduced MM risk. A suggestive association between increased telomere length and reduced MM risk was also noted; however, this association was primarily driven by the previously identified risk variant rs10936599 at 3q26 (TERC). Although not statistically significant, increased body mass index was associated with increased risk (OR, 1.10; 95% confidence interval, 0.99-1.22), supporting findings from a previous meta-analysis of prospective observational studies. Our study did not provide evidence supporting any modifiable factors examined as having a major influence on MM risk; however, it provides insight into factors for which the evidence has previously been mixed.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Multi-omics tumor profiling technologies to develop precision medicine in multiple myeloma

Multiple myeloma (MM), the second most common hematologic cancer, is caused by accumulation of aberrant plasma cells in the bone marrow. Its molecular causes are not fully understood and its great heterogeneity among patients complicates therapeutic decision-making. In the past decades, development of new therapies and drugs have significantly improved survival of MM patients. However, resistance to drugs and relapse remain the most common causes of mortality and are the major challenges to overcome. The advent of high throughput omics technologies capable of analyzing big amount of clinical and biological data has changed the way to diagnose and treat MM. Integration of omics data (gene mutations, gene expression, epigenetic information, and protein and metabolite levels) with clinical histories of thousands of patients allows to build scores to stratify the risk at diagnosis and predict the response to treatment, helping clinicians to make better educated decisions for each particular case. There is no doubt that the future of MM treatment relies on personalized therapies based on predictive models built from omics studies. This review summarizes the current treatments and the use of omics technologies in MM, and their importance in the implementation of personalized medicine.

Multi-omics tumor profiling technologies to develop precision medicine in multiple myeloma

Multiple myeloma (MM), the second most common hematologic cancer, is caused by accumulation of aberrant plasma cells in the bone marrow. Its molecular causes are not fully understood and its great heterogeneity among patients complicates therapeutic decision-making. In the past decades, development of new therapies and drugs have significantly improved survival of MM patients. However, resistance to drugs and relapse remain the most common causes of mortality and are the major challenges to overcome. The advent of high throughput omics technologies capable of analyzing big amount of clinical and biological data has changed the way to diagnose and treat MM. Integration of omics data (gene mutations, gene expression, epigenetic information, and protein and metabolite levels) with clinical histories of thousands of patients allows to build scores to stratify the risk at diagnosis and predict the response to treatment, helping clinicians to make better educated decisions for each particular case. There is no doubt that the future of MM treatment relies on personalized therapies based on predictive models built from omics studies. This review summarizes the current treatments and the use of omics technologies in MM, and their importance in the implementation of personalized medicine.

Do myeloproliferative neoplasms and multiple myeloma share the same genetic susceptibility loci?

Myeloproliferative neoplasms (MPNs) are a group of diseases that cause myeloid hematopoietic cells to overproliferate. Epidemiological and familial studies suggest that genetic factors contribute to the risk of developing MPN, but the genetic susceptibility of MPN is still not well known. Indeed, only few loci are known to have a clear role in the predisposition to this disease. Some studies reported a diagnosis of MPNs and multiple myeloma (MM) in the same patients, but the biological causes are still unclear. We tested the hypothesis that the two diseases share at least partly the same genetic risk loci. In the context of a European multicenter study with 460 cases and 880 controls, we analyzed the effect of the known MM risk loci, individually and in a polygenic risk score (PRS). The most significant result was obtained among patients with chronic myeloid leukemia (CML) for PS0RS1C1-rs2285803, which showed to be associated with an increased risk (OR = 3.28, 95% CI 1.79-6.02, P = .00012, P = .00276 when taking into account multiple testing). Additionally, the PRS showed an association with MPN risk when comparing the last with the first quartile of the PRS (OR = 2.39, 95% CI 1.64-3.48, P = 5.98 × 10^-6 ). In conclusion, our results suggest a potential common genetic background between MPN and MM, which needs to be further investigated.

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.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

MPRAscore: robust and non-parametric analysis of massively parallel reporter assays

MOTIVATION

Massively parallel reporter assays (MPRA) enable systematic screening of DNA sequence variants for effects on transcriptional activity. However, convenient analysis tools are still needed.

RESULTS

We introduce MPRAscore, a novel tool to infer allele-specific effects on transcription from MPRA data. MPRAscore uses a weighted, variance-regularized method to calculate variant effect sizes robustly, and a permutation approach to test for significance without assuming normality or independence.

AVAILABILITY AND IMPLEMENTATION

Source code (C++), precompiled binaries and data used in the paper at https://github.com/abhisheknrl/MPRAscore and https://www.ncbi.nlm.nih.gov/bioproject/PRJNA554195.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alternative polyadenylation mediates genetic regulation of gene expression

Little is known about co-transcriptional or post-transcriptional regulatory mechanisms linking noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3’ Seq to study the impact of genetic variation on alternative polyadenylation (APA) in the nuclear and total mRNA fractions of 52 HapMap Yoruba human lymphoblastoid cell lines. We mapped 602 APA quantitative trait loci (apaQTLs) at 10% FDR, of which 152 were nuclear specific. Effect sizes at intronic apaQTLs are negatively correlated with eQTL effect sizes. These observations suggest genetic variants can decrease mRNA expression levels by increasing usage of intronic PAS. We also identified 24 apaQTLs associated with protein levels, but not mRNA expression. Finally, we found that 19% of apaQTLs can be associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression, protein expression, and disease risk, and reveals uncharted modes of genetic regulation.

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.

Assessing thyroid cancer risk using polygenic risk scores

Genome-wide association studies (GWASs) have identified at least 10 single-nucleotide polymorphisms (SNPs) associated with papillary thyroid cancer (PTC) risk. Most of these SNPs are common variants with small to moderate effect sizes. Here we assessed the combined genetic effects of these variants on PTC risk by using summarized GWAS results to build polygenic risk score (PRS) models in three PTC study groups from Ohio (1,544 patients and 1,593 controls), Iceland (723 patients and 129,556 controls), and the United Kingdom (534 patients and 407,945 controls). A PRS based on the 10 established PTC SNPs showed a stronger predictive power compared with the clinical factors model, with a minimum increase of area under the receiver-operating curve of 5.4 percentage points (P ≤ 1.0 × 10^-9). Adding an extended PRS based on 592,475 common variants did not significantly improve the prediction power compared with the 10-SNP model, suggesting that most of the remaining undiscovered genetic risk in thyroid cancer is due to rare, moderate- to high-penetrance variants rather than to common low-penetrance variants. Based on the 10-SNP PRS, individuals in the top decile group of PRSs have a close to sevenfold greater risk (95% CI, 5.4-8.8) compared with the bottom decile group. In conclusion, PRSs based on a small number of common germline variants emphasize the importance of heritable low-penetrance markers in PTC.

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

Transcriptome-wide association study of multiple myeloma identifies candidate susceptibility genes

BACKGROUND

While genome-wide association studies (GWAS) of multiple myeloma (MM) have identified variants at 23 regions influencing risk, the genes underlying these associations are largely unknown. To identify candidate causal genes at these regions and search for novel risk regions, we performed a multi-tissue transcriptome-wide association study (TWAS).

RESULTS

GWAS data on 7319 MM cases and 234,385 controls was integrated with Genotype-Tissue Expression Project (GTEx) data assayed in 48 tissues (sample sizes, N = 80-491), including lymphocyte cell lines and whole blood, to predict gene expression. We identified 108 genes at 13 independent regions associated with MM risk, all of which were in 1 Mb of known MM GWAS risk variants. Of these, 94 genes, located in eight regions, had not previously been considered as a candidate gene for that locus.

CONCLUSIONS

Our findings highlight the value of leveraging expression data from multiple tissues to identify candidate genes responsible for GWAS associations which provide insight into MM tumorigenesis. Among the genes identified, a number have plausible roles in MM biology, notably APOBEC3C, APOBEC3H, APOBEC3D, APOBEC3F, APOBEC3G, or have been previously implicated in other malignancies. The genes identified in this TWAS can be explored for follow-up and validation to further understand their role in MM biology.

ELL2 Is Downregulated and Associated with Galactose-Deficient IgA1 in IgA Nephropathy

Background

Galactose-deficient IgA1 (Gd-IgA1) is an important causal factor in IgA nephropathy; however, the underlying mechanism for the production of Gd-IgA1 is unknown. The elongation factor for RNA polymerase II (ELL2), which encoded a key component of the superelongation complex (SEC), drives secretory-specific Ig mRNA production.

Methods

We enrolled 21 patients with IgAN, 18 healthy controls, and 20 patients with non-IgAN glomerulonephritis. The differential expression of ELL2 was compared using publically available data from Gene Expression Omnibus (GEO) datasets. The relationship between ELL2 expressions and galactose-deficient IgA1 (Gd-IgA1) levels in serum were also studied. At last, the results were validated by shELL2 treatment experiment.

Results

We found that the number of CD19+ B cells was increased in IgAN patients compared to healthy controls. The expression level of ELL2 in patients with IgAN was significantly lower than that of healthy control and disease control. Consistent with present results, the lower ELL2 expression in IgAN patients was observed in microarray expression profiles from GEO datasets. Pearson correlation analysis showed that ELL2 expression negatively correlated with Gd-IgA1 levels. Furthermore, in an in vitro experiment, we found that loss of ELL2 function in human B lymphoma DAKIKI cells, an IgA1-producing cell line, increased the levels of Gd-IgA1, which confirmed that ELL2 modulated the levels of Gd-IgA1.

Conclusion

Our findings implied that decreased ELL2 expression was negatively correlated with the numbers of B cells and aberrant glycosylation of IgA1 in IgAN.

Regions of homozygosity as risk factors for multiple myeloma

Genomic regions of homozygosity (ROH), detectable in outbred populations, have been implicated as determinants of inherited risk. To examine whether ROH is associated with risk of multiple myeloma (MM), we performed whole-genome homozygosity analysis using single-nucleotide polymorphism genotype data from 2,282 MM cases and 5,197 controls, with replication in an additional 878 MM cases and 7,083 controls. Globally, the distribution of ROH between cases and controls was not significantly different. However, one ROH at chromosome 9q21, harboring the B-cell transcription factor gene KLF9, showed evidence of a consistent association and may therefore warrant further investigation as a candidate risk factor for MM. Overall, our analysis provides little support for a homozygosity signature being a significant factor in MM risk.

Identification of Target Genes at Juvenile Idiopathic Arthritis GWAS Loci in Human Neutrophils

Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic disease among children which could cause severe disability. Genomic studies have discovered substantial number of risk loci for JIA, however, the mechanism of how these loci affect JIA development is not fully understood. Neutrophil is an important cell type involved in autoimmune diseases. To better understand the biological function of genetic loci in neutrophils during JIA development, we took an integrated multi-omics approach to identify target genes at JIA risk loci in neutrophils and constructed a protein-protein interaction network via a machine learning approach. We identified genes likely to be JIA risk loci targeted genes in neutrophils which could contribute to JIA development.

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.

ELL2 Influences Transcription Elongation, Splicing, Ig Secretion and Growth

ELL2 was previously discovered as a member of the Super Elongation Complex. It is involved in driving the maturation of B cells to plasma cells through shifting patterns of RNA processing, favoring generation of the secretory form of heavy chain immunoglobulin (IgH) associated with plasma cells. ELL2 influences the expression and splicing patterns of more than 4,000 genes in antibody secreting cells. The ELL2 gene has been implicated in cancers such as multiple myeloma and salivary gland carcinoma. A member of the ELL family (ELL1) was recently proven to act as an E3 ubiquitin ligase to known proto-oncogene, c-Myc, through a highly conserved cysteine residue in the C-terminal CEYLH region. Comparison of sequence homology shows this region is conserved between the three members of the ELL family, leading us to hypothesize that the other two ELLs (2 and 3) could serve the same role. In this review, we summarize what is known about ELL2 with respect to its role in driving B cell to plasma cell differentiation as well as its potential role in tumor suppression.

Genetic correlation between multiple myeloma and chronic lymphocytic leukaemia provides evidence for shared aetiology

The clustering of different types of B-cell malignancies in families raises the possibility of shared aetiology. To examine this, we performed cross-trait linkage disequilibrium (LD)-score regression of multiple myeloma (MM) and chronic lymphocytic leukaemia (CLL) genome-wide association study (GWAS) data sets, totalling 11,734 cases and 29,468 controls. A significant genetic correlation between these two B-cell malignancies was shown (Rg = 0.4, P = 0.0046). Furthermore, four of the 45 known CLL risk loci were shown to associate with MM risk and five of the 23 known MM risk loci associate with CLL risk. By integrating eQTL, Hi-C and ChIP-seq data, we show that these pleiotropic risk loci are enriched for B-cell regulatory elements and implicate B-cell developmental genes. These data identify shared biological pathways influencing the development of CLL and, MM and further our understanding of the aetiological basis of these B-cell malignancies.

RNA Splicing in the Transition from B Cells to Antibody-Secreting Cells: The Influences of ELL2, Small Nuclear RNA, and Endoplasmic Reticulum Stress

In the transition from B cells to Ab-secreting cells (ASCs) many genes are induced, such as ELL2, Irf4, Prdm1, Xbp1, whereas other mRNAs do not change in abundance. Nonetheless, using splicing array technology and mouse splenic B cells plus or minus LPS, we found that induced and "uninduced" genes can show large differences in splicing patterns between the cell stages, which could influence ASC development. We found that ∼55% of these splicing changes depend on ELL2, a transcription elongation factor that influences expression levels and splicing patterns of ASC signature genes, genes in the cell-cycle and N-glycan biosynthesis and processing pathways, and the secretory versus membrane forms of the IgH mRNA. Some of these changes occur when ELL2 binds directly to the genes encoding those mRNAs, whereas some of the changes are indirect. To attempt to account for the changes that occur in RNA splicing before or without ELL2 induction, we examined the amount of the small nuclear RNA molecules and found that they were significantly decreased within 18 h of LPS stimulation and stayed low until 72 h. Correlating with this, at 18 h after LPS, endoplasmic reticulum stress and Ire1 phosphorylation are induced. Inhibiting the regulated Ire1-dependent mRNA decay with 4u8C correlates with the reduction in small nuclear RNA and changes in the normal splicing patterns at 18 h. Thus, we conclude that the RNA splicing patterns in ASCs are shaped early by endoplasmic reticulum stress and Ire1 phosphorylation and later by ELL2 induction.

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.

Human genetic variants and age are the strongest predictors of humoral immune responses to common pathogens and vaccines

BACKGROUND

Humoral immune responses to infectious agents or vaccination vary substantially among individuals, and many of the factors responsible for this variability remain to be defined. Current evidence suggests that human genetic variation influences (i) serum immunoglobulin levels, (ii) seroconversion rates, and (iii) intensity of antigen-specific immune responses. Here, we evaluated the impact of intrinsic (age and sex), environmental, and genetic factors on the variability of humoral response to common pathogens and vaccines.

METHODS

We characterized the serological response to 15 antigens from common human pathogens or vaccines, in an age- and sex-stratified cohort of 1000 healthy individuals (Milieu Intérieur cohort). Using clinical-grade serological assays, we measured total IgA, IgE, IgG, and IgM levels, as well as qualitative (serostatus) and quantitative IgG responses to cytomegalovirus, Epstein-Barr virus, herpes simplex virus 1 and 2, varicella zoster virus, Helicobacter pylori, Toxoplasma gondii, influenza A virus, measles, mumps, rubella, and hepatitis B virus. Following genome-wide genotyping of single nucleotide polymorphisms and imputation, we examined associations between ~ 5 million genetic variants and antibody responses using single marker and gene burden tests.

RESULTS

We identified age and sex as important determinants of humoral immunity, with older individuals and women having higher rates of seropositivity for most antigens. Genome-wide association studies revealed significant associations between variants in the human leukocyte antigen (HLA) class II region on chromosome 6 and anti-EBV and anti-rubella IgG levels. We used HLA imputation to fine map these associations to amino acid variants in the peptide-binding groove of HLA-DRβ1 and HLA-DPβ1, respectively. We also observed significant associations for total IgA levels with two loci on chromosome 2 and with specific KIR-HLA combinations.

CONCLUSIONS

Using extensive serological testing and genome-wide association analyses in a well-characterized cohort of healthy individuals, we demonstrated that age, sex, and specific human genetic variants contribute to inter-individual variability in humoral immunity. By highlighting genes and pathways implicated in the normal antibody response to frequently encountered antigens, these findings provide a basis to better understand disease pathogenesis.

TRIALS REGISTRATION

ClinicalTrials.gov , NCT01699893.

Germline Lysine-Specific Demethylase 1 (LSD1/KDM1A) Mutations Confer Susceptibility to Multiple Myeloma

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.