Mosaic chromosomal alterations in hematopoietic cells and clinical outcomes in patients with multiple myeloma.

Mosaic chromosomal alterations (mCAs) in hematopoietic cells increase mortality and risk of hematological cancers and infections. We investigated the landscape of mCAs and their clinical consequences in 976 patients with multiple myeloma undergoing high-dose chemotherapy and autologous stem cell support (ASCT) with median 6.4 years of follow-up. mCAs were detected in the stem cell harvest product of 158 patients (16.2%). Autosomal aberrations were found in 60 patients (6.1%) and affected all chromosomes. Loss of chromosome X was found in 51 females (12.7%) and loss of chromosome Y in 55 males (9.6%). Overall survival and progression were similar between carriers of autosomal mCAs and non-carriers. In contrast, female patients with loss of the X chromosome had longer overall survival (age-adjusted[a.a.] HR 0.54, 95% CI 0.32-0.93, p = 0.02), lower risk of progression (a.a. HR 0.55, 95% CI 0.35-0.87; p = 0.01), and better post-transplant response (higher degree of complete response (CR) or very good partial response (VGPR)). The reason for this substantial effect is unknown. Additionally, myeloma clones in the stem cell product was confirmed by mCA analysis in the few patients with multiple mCAs (n = 12 patients). Multiple mCAs conferred inferior overall survival (a.a. HR 2.0, 95% CI 1.02-3.84; p = 0.04) and higher risk of myeloma progression (a.a. HR 3.36, 95% CI 1.67-6.81; p < 0.001), which is presumed to be driven by suspected myeloma contaminants.

Sequence variants influencing the regulation of serum IgG subclass levels.

Immunoglobulin G (IgG) is the main isotype of antibody in human blood. IgG consists of four subclasses (IgG1 to IgG4), encoded by separate constant region genes within the Ig heavy chain locus (IGH). Here, we report a genome-wide association study on blood IgG subclass levels. Across 4334 adults and 4571 individuals under 18 years, we discover ten new and identify four known variants at five loci influencing IgG subclass levels. These variants also affect the risk of asthma, autoimmune diseases, and blood traits. Seven variants map to the IGH locus, three to the Fcγ receptor (FCGR) locus, and two to the human leukocyte antigen (HLA) region, affecting the levels of all IgG subclasses. The most significant associations are observed between the G1m (f), G2m(n) and G3m(b*) allotypes, and IgG1, IgG2 and IgG3, respectively. Additionally, we describe selective associations with IgG4 at 16p11.2 (ITGAX) and 17q21.1 (IKZF3, ZPBP2, GSDMB, ORMDL3). Interestingly, the latter coincides with a highly pleiotropic signal where the allele associated with lower IgG4 levels protects against childhood asthma but predisposes to inflammatory bowel disease. Our results provide insight into the regulation of antibody-mediated immunity that can potentially be useful in the development of antibody based therapeutics.

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.

An exploratory study of the damage markers NfL, GFAP, and t-Tau, in cerebrospinal fluid and other findings from a patient cohort enriched for suspected autoimmune psychiatric disease.

There is growing evidence suggesting that immunological mechanisms play a significant role in the development of psychiatric symptoms in certain patient subgroups. However, the relationship between clinical red flags for suspected autoimmune psychiatric disease and signs of central nervous system (CNS) pathology (e.g., routine cerebrospinal fluid (CSF) alterations, CNS damage markers, neurophysiological or neuroimaging findings) has received limited attention. Here, we aimed to describe the prevalence and distribution of potential CNS pathologies in psychiatric patients in relation to clinical red flags for autoimmune psychiatric disease and psychiatric symptoms. CSF routine findings and CNS damage markers; neurofilament light chain protein (NfL), glial fibrillary acidic protein (GFAP) and total Tau (t-Tau), in CSF from 127 patients with psychiatric disease preselected for suspected immunological involvement were related to recently proposed clinical red flags, psychiatric features, and MRI and EEG findings. Twenty-one percent had abnormal routine CSF findings and 27% had elevated levels of CNS damage markers. Six percent had anti-neuronal antibodies in serum and 2% had these antibodies in the CSF. Sixty-six percent of patients examined with MRI (n = 88) had alterations, mostly atrophy or nonspecific white matter lesions. Twenty-seven percent of patients with EEG recordings (n = 70) had abnormal findings. Elevated NfL levels were associated with comorbid autoimmunity and affective dysregulation symptoms. Elevated t-Tau was associated with catatonia and higher ratings of agitation/hyperactivity. Elevated GFAP was associated with acute onset, atypical presentation, infectious prodrome, tics, depressive/anxiety symptom ratings and overall greater psychiatric symptom burden. In conclusion, preselection based on suspected autoimmune psychiatric disease identifies a population with a high prevalence of CSF alterations suggesting CNS pathology. Future studies should examine the value of these markers in predicting treatment responses.

Functional and molecular profiling of hematopoietic stem cells during regeneration.

Hematopoietic stem cells (HSCs) enable hematopoietic stem cell transplantation (HCT) through their ability to replenish the entire blood system. Proliferation of HSCs is linked to decreased reconstitution potential, and a precise regulation of actively dividing HSCs is thus essential to ensure long-term functionality. This regulation becomes important in the transplantation setting where HSCs undergo proliferation followed by a gradual transition to quiescence and homeostasis. Although mouse HSCs have been well studied under homeostatic conditions, the mechanisms regulating HSC activation under stress remain unclear. Here, we analyzed the different phases of regeneration after transplantation. We isolated bone marrow from mice at 8 time points after transplantation and examined the reconstitution dynamics and transcriptional profiles of stem and progenitor populations. We found that regenerating HSCs initially produced rapidly expanding progenitors and displayed distinct changes in fatty acid metabolism and glycolysis. Moreover, we observed molecular changes in cell cycle, MYC and mTOR signaling in both HSCs, and progenitor subsets. We used a decay rate model to fit the temporal transcription profiles of regenerating HSCs and identified genes with progressively decreased or increased expression after transplantation. These genes overlapped to a large extent with published gene sets associated with key aspects of HSC function, demonstrating the potential of this data set as a resource for identification of novel HSC regulators. Taken together, our study provides a detailed functional and molecular characterization of HSCs at different phases of regeneration and identifies a gene set associated with the transition from proliferation to quiescence.

AliGater: a framework for the development of bioinformatic pipelines for large-scale, high-dimensional cytometry data.

Motivation: AliGater is an open-source framework to accelerate the development of bioinformatic pipelines for the analysis of large-scale, high-dimensional flow cytometry data. AliGater provides a Python package for automatic feature extraction workflows, as well as building blocks to construct analysis pipelines. Results: We illustrate the use of AliGater in a high-resolution flow cytometry-based genome-wide association study on 46 immune cell populations in 14 288 individuals. Availability and implementation: Source code and documentation at https://github.com/LudvigEk/aligater and https://aligater.readthedocs.io.

Sequence enrichment profiles enable target-agnostic antibody generation for a broad range of antigens.

Phenotypic drug discovery (PDD) enables the target-agnostic generation of therapeutic drugs with novel mechanisms of action. However, realizing its full potential for biologics discovery requires new technologies to produce antibodies to all, a priori unknown, disease-associated biomolecules. We present a methodology that helps achieve this by integrating computational modeling, differential antibody display selection, and massive parallel sequencing. The method uses the law of mass action-based computational modeling to optimize antibody display selection and, by matching computationally modeled and experimentally selected sequence enrichment profiles, predict which antibody sequences encode specificity for disease-associated biomolecules. Applied to a phage display antibody library and cell-based antibody selection, ∼105 antibody sequences encoding specificity for tumor cell surface receptors expressed at 103-106 receptors/cell were discovered. We anticipate that this approach will be broadly applicable to molecular libraries coupling genotype to phenotype and to the screening of complex antigen populations for identification of antibodies to unknown disease-associated targets.

Genetic regulation of fetal hemoglobin across global populations.

Human genetic variation has enabled the identification of several key regulators of fetal-to-adult hemoglobin switching, including BCL11A, resulting in therapeutic advances. However, despite the progress made, limited further insights have been obtained to provide a fuller accounting of how genetic variation contributes to the global mechanisms of fetal hemoglobin (HbF) gene regulation. Here, we have conducted a multi-ancestry genome-wide association study of 28,279 individuals from several cohorts spanning 5 continents to define the architecture of human genetic variation impacting HbF. We have identified a total of 178 conditionally independent genome-wide significant or suggestive variants across 14 genomic windows. Importantly, these new data enable us to better define the mechanisms by which HbF switching occurs in vivo. We conduct targeted perturbations to define BACH2 as a new genetically-nominated regulator of hemoglobin switching. We define putative causal variants and underlying mechanisms at the well-studied BCL11A and HBS1L-MYB loci, illuminating the complex variant-driven regulation present at these loci. We additionally show how rare large-effect deletions in the HBB locus can interact with polygenic variation to influence HbF levels. Our study paves the way for the next generation of therapies to more effectively induce HbF in sickle cell disease and ß-thalassemia.

Genome-wide association study on 13 167 individuals identifies regulators of blood CD34+cell levels.

Stem cell transplantation is a cornerstone in the treatment of blood malignancies. The most common method to harvest stem cells for transplantation is by leukapheresis, requiring mobilization of CD34+ hematopoietic stem and progenitor cells (HSPCs) from the bone marrow into the blood. Identifying the genetic factors that control blood CD34+ cell levels could reveal new drug targets for HSPC mobilization. Here we report the first large-scale, genome-wide association study on blood CD34+ cell levels. Across 13 167 individuals, we identify 9 significant and 2 suggestive associations, accounted for by 8 loci (PPM1H, CXCR4, ENO1-RERE, ITGA9, ARHGAP45, CEBPA, TERT, and MYC). Notably, 4 of the identified associations map to CXCR4, showing that bona fide regulators of blood CD34+ cell levels can be identified through genetic variation. Further, the most significant association maps to PPM1H, encoding a serine/threonine phosphatase never previously implicated in HSPC biology. PPM1H is expressed in HSPCs, and the allele that confers higher blood CD34+ cell levels downregulates PPM1H. Through functional fine-mapping, we find that this downregulation is caused by the variant rs772557-A, which abrogates an MYB transcription factor-binding site in PPM1H intron 1 that is active in specific HSPC subpopulations, including hematopoietic stem cells, and interacts with the promoter by chromatin looping. Furthermore, PPM1H knockdown increases the proportion of CD34+ and CD34+90+ cells in cord blood assays. Our results provide the first large-scale analysis of the genetic architecture of blood CD34+ cell levels and warrant further investigation of PPM1H as a potential inhibition target for stem cell mobilization.

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.

ZBTB33 is mutated in clonal hematopoiesis and myelodysplastic syndromes and impacts RNA splicing.

Clonal hematopoiesis results from somatic mutations in cancer driver genes in hematopoietic stem cells. We sought to identify novel drivers of clonal expansion using an unbiased analysis of sequencing data from 84,683 persons and identified common mutations in the 5-methylcytosine reader, ZBTB33, as well as in YLPM1, SRCAP, and ZNF318. We also identified these mutations at low frequency in myelodysplastic syndrome patients. Zbtb33 edited mouse hematopoietic stem and progenitor cells exhibited a competitive advantage in vivo and increased genome-wide intron retention. ZBTB33 mutations potentially link DNA methylation and RNA splicing, the two most commonly mutated pathways in clonal hematopoiesis and MDS.

Germline variants at SOHLH2 influence multiple myeloma risk.

Multiple myeloma (MM) is caused by the uncontrolled, clonal expansion of plasma cells. While there is epidemiological evidence for inherited susceptibility, the molecular basis remains incompletely understood. We report a genome-wide association study totalling 5,320 cases and 422,289 controls from four Nordic populations, and find a novel MM risk variant at SOHLH2 at 13q13.3 (risk allele frequency = 3.5%; odds ratio = 1.38; P = 2.2 × 10-14). This gene encodes a transcription factor involved in gametogenesis that is normally only weakly expressed in plasma cells. The association is represented by 14 variants in linkage disequilibrium. Among these, rs75712673 maps to a genomic region with open chromatin in plasma cells, and upregulates SOHLH2 in this cell type. Moreover, rs75712673 influences transcriptional activity in luciferase assays, and shows a chromatin looping interaction with the SOHLH2 promoter. Our work provides novel insight into MM susceptibility.

Accelerating target deconvolution for therapeutic antibody candidates using highly parallelized genome editing.

Therapeutic antibodies are transforming the treatment of cancer and autoimmune diseases. Today, a key challenge is finding antibodies against new targets. Phenotypic discovery promises to achieve this by enabling discovery of antibodies with therapeutic potential without specifying the molecular target a priori. Yet, deconvoluting the targets of phenotypically discovered antibodies remains a bottleneck; efficient deconvolution methods are needed for phenotypic discovery to reach its full potential. Here, we report a comprehensive investigation of a target deconvolution approach based on pooled CRISPR/Cas9. Applying this approach within three real-world phenotypic discovery programs, we rapidly deconvolute the targets of 38 of 39 test antibodies (97%), a success rate far higher than with existing approaches. Moreover, the approach scales well, requires much less work, and robustly identifies antibodies against the major histocompatibility complex. Our data establish CRISPR/Cas9 as a highly efficient target deconvolution approach, with immediate implications for the development of antibody-based drugs.

Inherited myeloproliferative neoplasm risk affects haematopoietic stem cells.

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P < 5.0 × 10-8), 7 of which have not been previously reported. We find that there is a shared genetic architecture between MPN risk and several haematopoietic traits from distinct lineages; that there is an enrichment for MPN risk variants within accessible chromatin of HSCs; and that increased MPN risk is associated with longer telomere length in leukocytes and other clonal haematopoietic states-collectively suggesting that MPN risk is associated with the function and self-renewal of HSCs. We use gene mapping to identify modulators of HSC biology linked to MPN risk, and show through targeted variant-to-function assays that CHEK2 and GFI1B have roles in altering the function of HSCs to confer disease risk. Overall, our results reveal a previously unappreciated mechanism for inherited MPN risk through the modulation of HSC function.

Correction: S100A6 is a critical regulator of hematopoietic stem cells.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.