An Integrated Multi-omic Single-Cell Atlas of Human B Cell Identity.

B cells are capable of a wide range of effector functions including antibody secretion, antigen presentation, cytokine production, and generation of immunological memory. A consistent strategy for classifying human B cells by using surface molecules is essential to harness this functional diversity for clinical translation. We developed a highly multiplexed screen to quantify the co-expression of 351 surface molecules on millions of human B cells. We identified differentially expressed molecules and aligned their variance with isotype usage, VDJ sequence, metabolic profile, biosynthesis activity, and signaling response. Based on these analyses, we propose a classification scheme to segregate B cells from four lymphoid tissues into twelve unique subsets, including a CD45RB+CD27- early memory population, a class-switched CD39+ tonsil-resident population, and a CD19hiCD11c+ memory population that potently responds to immune activation. This classification framework and underlying datasets provide a resource for further investigations of human B cell identity and function.

Human chromosomal translocations at CpG sites and a theoretical basis for their lineage and stage specificity.

We have assembled, annotated, and analyzed a database of over 1700 breakpoints from the most common chromosomal rearrangements in human leukemias and lymphomas. Using this database, we show that although the CpG dinucleotide constitutes only 1% of the human genome, it accounts for 40%-70% of breakpoints at pro-B/pre-B stage translocation regions-specifically, those near the bcl-2, bcl-1, and E2A genes. We do not observe CpG hotspots in rearrangements involving lymphoid-myeloid progenitors, mature B cells, or T cells. The stage specificity, lineage specificity, CpG targeting, and unique breakpoint distributions at these cluster regions may be explained by a lesion-specific double-strand breakage mechanism involving the RAG complex acting at AID-deaminated methyl-CpGs.

Bone marrow morphology is a strong discriminator between chronic eosinophilic leukemia, not otherwise specified and reactive idiopathic hypereosinophilic syndrome.

Chronic eosinophilic leukemia, not otherwise specified can be difficult to distinguish from idiopathic hypereosinophilic syndrome according to the current World Health Organization guideline. To examine whether the morphological features of bone marrow might aid in the differential diagnosis of these two entities, we studied a total of 139 patients with a diagnosis of chronic eosinophilic leukemia, not otherwise specified (n=17) or idiopathic hypereosinophilic syndrome (n=122). As a group, abnormal bone marrow morphological features, resembling myelodysplastic syndromes, myeloproliferative neoplasm or myelodysplastic/myeloproliferative neoplasm, were identified in 40/139 (27%) patients: 16 (94%) of those with chronic eosinophilic leukemia and 24 (20%) of those with hypereosinophilic syndrome. Abnormal bone marrow correlated with older age (P<0.001), constitutional symptoms (P<0.001), anemia (P=0.041), abnormal platelet count (P=0.002), organomegaly (P=0.008), elevated lactate dehydrogenase concentration (P=0.005), abnormal karyotype (P<0.001), as well as the presence of myeloid neoplasm-related mutations (P<0.001). Patients with abnormal bone marrow had shorter survival (48.1 months versus not reached, P<0.001), a finding which was independent of other confounding factors (P<0.001). The association between abnormal bone marrow and shorter survival was also observed in hypereosinophilic syndrome patients alone. In summary, most patients with chronic eosinophilic leukemia, not otherwise specified and a proportion of those with idiopathic hypereosinophilic syndrome show abnormal bone marrow features similar to the ones encountered in patients with myelodysplastic syndromes, myelodysplastic/myeloproliferative neoplasm or BCR-ABL1-negative myeloproliferative neoplasm. Among patients who are currently considered to have idiopathic hypereosinophilic syndrome, abnormal bone marrow is a strong indicator of clonal hematopoiesis. Similar to other myeloid neoplasms, bone marrow morphology should be one of the major criteria to distinguish patients with chronic eosinophilic leukemia, not otherwise specified or clonal hypereosinophilic syndrome from those with truly reactive idiopathic hypereosinophilic syndrome.

H3K4me3 stimulates the V(D)J RAG complex for both nicking and hairpinning in trans in addition to tethering in cis: implications for translocations.

The PHD finger of the RAG2 polypeptide of the RAG1/RAG2 complex binds to the histone H3 modification, trimethylated lysine 4 (H3K4me3), and in some manner increases V(D)J recombination. In the absence of biochemical studies of H3K4me3 on purified RAG enzyme activity, the precise role of H3K4me3 remains unclear. Here, we find that H3K4me3 stimulates purified RAG enzymatic activity at both the nicking (2- to 5-fold) and hairpinning (3- to 11-fold) steps of V(D)J recombination. Remarkably, this stimulation can be achieved with free H3K4me3 peptide (in trans), indicating that H3K4me3 functions via two distinct mechanisms. It not only tethers the RAG enzyme complex to a region of DNA, but it also induces a substantial increase in the catalytic turnover number (k(cat)) of the RAG complex. The H3K4me3 catalytic stimulation applies to suboptimal cryptic RSS sites located at H3K4me3 peaks that are critical in the inception of human T cell acute lymphoblastic lymphomas.

Targeted next-generation sequencing identifies a subset of idiopathic hypereosinophilic syndrome with features similar to chronic eosinophilic leukemia, not otherwise specified.

The distinction between chronic eosinophilic leukemia, not otherwise specified and idiopathic hypereosinophilic syndrome largely relies on clonality assessment. Prior to the advent of next-generation sequencing, clonality was usually determined by cytogenetic analysis. We applied targeted next-generation sequencing panels designed for myeloid neoplasms to bone marrow specimens from a cohort of idiopathic hypereosinophilic syndrome patients (n=51), and assessed the significance of mutations in conjunction with clinicopathological features. The findings were further compared with those of 17 chronic eosinophilic leukemia, not otherwise specified patients defined by their abnormal cytogenetics and/or increased blasts. Mutations were detected in 14/51 idiopathic hypereosinophilic syndrome patients (idiopathic hypereosinophilic syndrome/next-generation sequencing-positive) (28%), involving single gene in 7 and ≥2 in 7 patients. The more frequently mutated genes included ASXL1 (43%), TET2 (36%), EZH2 (29%), SETBP1 (22%), CBL (14%), and NOTCH1 (14%). Idiopathic hypereosinophilic syndrome/next-generation sequencing-positive patients showed a number of clinical features and bone marrow findings resembling chronic eosinophilic leukemia, not otherwise specified. Chronic eosinophilic leukemia, not otherwise specified patients showed a disease-specific survival of 14.4 months, markedly inferior to idiopathic hypereosinophilic syndrome/next-generation sequencing-negative (P<0.001), but not significantly different from idiopathic hypereosinophilic syndrome/next-generation sequencing-positive (P=0.117). These data suggest that targeted next-generation sequencing helps to establish clonality in a subset of patients with hypereosinophilia that would otherwise be classified as idiopathic hypereosinophilic syndrome. In conjunction with other diagnostic features, mutation data can be used to establish a diagnosis of chronic eosinophilic leukemia, not otherwise specified in patients presenting with hypereosinophilia.

BCL6 breaks occur at different AID sequence motifs in Ig-BCL6 and non-Ig-BCL6 rearrangements.

BCL6 translocations are common in B-cell lymphomas and frequently have chromosomal breaks in immunoglobulin heavy chain (IgH) switch regions, suggesting that they occur during class-switch recombination. We analyze 120 BCL6 translocation breakpoints clustered in a 2156-bp segment of BCL6 intron 1, including 62 breakpoints (52%) joined to IgH, 12 (10%) joined to Ig light chains, and 46 (38%) joined to non-Ig partners. The BCL6 breaks in Ig-BCL6 translocations prefer known activation-induced cytosine deaminase (AID) hotspots such as WGCW and WRC (W = A/T, R = A/G), whereas BCL6 breaks in non-Ig rearrangements occur at CpG/CGC sites in addition to WGCW. Unlike previously identified CpG breaks in pro-B/pre-B-cell translocations, the BCL6 breaks do not show evidence of recombination activating gene or terminal deoxynucleotidyl transferase activity. Both WGCW/WRC and CpG/CGC breaks at BCL6 are most likely initiated by AID in germinal center B-cells, and their differential use suggests subtle mechanistic differences between Ig-BCL6 and non-Ig-BCL6 rearrangements.

IgH partner breakpoint sequences provide evidence that AID initiates t(11;14) and t(8;14) chromosomal breaks in mantle cell and Burkitt lymphomas.

Previous studies have implicated activation-induced cytidine deaminase (AID) in B-cell translocations but have failed to identify any association between their chromosomal breakpoints and known AID target sequences. Analysis of 56 unclustered IgH-CCND1 translocations in mantle cell lymphoma across the ~ 344-kb bcl-1 breakpoint locus demonstrates that half of the CCND1 breaks are near CpG dinucleotides. Most of these CpG breaks are at CGC motifs, and half of the remaining breaks are near WGCW, both known AID targets. These findings provide the strongest evidence to date that AID initiates chromosomal breaks in translocations that occur in human bone marrow B-cell progenitors. We also identify WGCW breaks at the MYC locus in Burkitt lymphoma translocations and murine IgH-MYC translocations, both of which arise in mature germinal center B cells. Finally, we propose a developmental model to explain the transition from CpG breaks in early human B-cell progenitors to WGCW breaks in later stage B cells.

Terminal deoxynucleotidyl transferase and CD84 identify human multi-potent lymphoid progenitors.

Lymphoid specification in human hematopoietic progenitors is not fully understood. To better associate lymphoid identity with protein-level cell features, we conduct a highly multiplexed single-cell proteomic screen on human bone marrow progenitors. This screen identifies terminal deoxynucleotidyl transferase (TdT), a specialized DNA polymerase intrinsic to VDJ recombination, broadly expressed within CD34+ progenitors prior to B/T cell emergence. While these TdT+ cells coincide with granulocyte-monocyte progenitor (GMP) immunophenotype, their accessible chromatin regions show enrichment for lymphoid-associated transcription factor (TF) motifs. TdT expression on GMPs is inversely related to the SLAM family member CD84. Prospective isolation of CD84lo GMPs demonstrates robust lymphoid potentials ex vivo, while still retaining significant myeloid differentiation capacity, akin to LMPPs. This multi-omic study identifies human bone marrow lymphoid-primed progenitors, further defining the lympho-myeloid axis in human hematopoiesis.

Rapid Setup of Tissue Microarray and Tiled Area Imaging on the Multiplexed Ion Beam Imaging Microscope using the Tile/SED/Array Interface.

Multiplexed ion beam imaging (MIBI) is a next-generation mass spectrometry-based microscopy technique that generates 40+ plex images of protein expression in histologic tissues, enabling detailed dissection of cellular phenotypes and histoarchitectural organization. A key bottleneck in operation occurs when users select the physical locations on the tissue for imaging. As the scale and complexity of MIBI experiments have increased, the manufacturer-provided interface and third-party tools have become increasingly unwieldy for imaging large tissue microarrays and tiled tissue areas. Thus, a web-based, interactive, what-you-see-is-what-you-get (WYSIWYG) graphical interface layer - the tile/SED/array Interface (TSAI) - was developed for users to set imaging locations using familiar and intuitive mouse gestures such as drag-and-drop, click-and-drag, and polygon drawing. Written according to web standards already built into modern web browsers, it requires no installation of external programs, extensions, or compilers. Of interest to the hundreds of current MIBI users, this interface dramatically simplifies and accelerates the setup of large, complex MIBI runs.

Longitudinal study of 2 patients with cyclic thrombocytopenia, STAT3 and MPL mutations.

Cyclic thrombocytopenia (CTP) is a rare disease of periodic platelet count oscillations. The pathogenesis of CTP remains elusive. To study the underlying pathophysiology and genetic and cellular associations with CTP, we applied systems biology approaches to 2 patients with stable platelet cycling and reciprocal thrombopoietin (TPO) cycling at multiple time points through 2 cycles. Blood transcriptome analysis revealed cycling of platelet-specific genes, which are in parallel with and precede platelet count oscillation, indicating that cyclical platelet production leads platelet count cycling in both patients. Additionally, neutrophil and erythrocyte-specific genes also showed fluctuations correlating with platelet count changes, consistent with TPO effects on hematopoietic progenitors. Moreover, we found novel genetic associations with CTP. One patient had a novel germline heterozygous loss-of-function (LOF) thrombopoietin receptor (MPL) c.1210G>A mutation, and both had pathogenic somatic gain-of-function (GOF) variants in signal transducer and activator of transcription 3 (STAT3). In addition, both patients had clonal T-cell populations that remained stable throughout platelet count cycles. These mutations and clonal T cells may potentially involve in the pathogenic baseline in these patients, rendering exaggerated persistent thrombopoiesis oscillations of their intrinsic rhythm upon homeostatic perturbations. This work provides new insights into the pathophysiology of CTP and possible therapies.

Supervised dimensionality reduction for exploration of single-cell data by HSS-LDA.

Single-cell technologies generate large, high-dimensional datasets encompassing a diversity of omics. Dimensionality reduction captures the structure and heterogeneity of the original dataset, creating low-dimensional visualizations that contribute to the human understanding of data. Existing algorithms are typically unsupervised, using measured features to generate manifolds, disregarding known biological labels such as cell type or experimental time point. We repurpose the classification algorithm, linear discriminant analysis (LDA), for supervised dimensionality reduction of single-cell data. LDA identifies linear combinations of predictors that optimally separate a priori classes, enabling the study of specific aspects of cellular heterogeneity. We implement feature selection by hybrid subset selection (HSS) and demonstrate that this computationally efficient approach generates non-stochastic, interpretable axes amenable to diverse biological processes such as differentiation over time and cell cycle. We benchmark HSS-LDA against several popular dimensionality-reduction algorithms and illustrate its utility and versatility for the exploration of single-cell mass cytometry, transcriptomics, and chromatin accessibility data.

An optimized protocol for phenotyping human granulocytes by mass cytometry.

Granulocytes encompass diverse roles, from fighting off pathogens to regulating inflammatory processes in allergies. These roles are represented by distinct cellular phenotypes that we captured with mass cytometry (CyTOF). Our protocol enables simultaneous evaluation of human basophils, eosinophils, and neutrophils under homeostasis and upon immune activation by anti-Immunoglobulin E (anti-IgE) or interleukin-3 (IL-3). Granulocyte integrity and detection of protein markers were optimized so that rare granulocyte populations could be deeply characterized by single cell mass cytometry. For complete details on the use and execution of this protocol, please refer to Vivanco Gonzalez et al. (2020).

Nonhomologous DNA end joining (NHEJ) and chromosomal translocations in humans.

Double-strand breaks (DSBs) arise in dividing cells about ten times per cell per day. Causes include replication across a nick, free radicals of oxidative metabolism, ionizing radiation, and inadvertent action by enzymes of DNA metabolism (such as failures of type II topoisomerases or cleavage by recombinases at off-target sites). There are two major double-strand break repair pathways. Homologous recombination (HR) can repair double-strand breaks, but only during S phase and typically only if there are hundreds of base pairs of homology. The more commonly used pathway is nonhomologous DNA end joining, abbreviated NHEJ. NHEJ can repair a DSB at any time during the cell cycle and does not require any homology, although a few nucleotides of terminal microhomology are often utilized by the NHEJ enzymes, if present. The proteins and enzymes of NHEJ include Ku, DNA-PKcs, Artemis, DNA polymerase mu (Pol micro), DNA polymerase lambda (Pol lambda), XLF (also called Cernunnos), XRCC4, and DNA ligase IV. These enzymes constitute what some call the classical NHEJ pathway, and in wild type cells, the vast majority of joining events appear to proceed using these components. NHEJ is present in many prokaryotes, as well as all eukaryotes, and very similar mechanistic flexibility evolved both convergently and divergently. When two double-strand breaks occur on different chromosomes, then the rejoining is almost always done by NHEJ. The causes of DSBs in lymphomas most often involve the RAG or AID enzymes that function in the specialized processes of antigen receptor gene rearrangement.

Mechanisms of chromosomal rearrangement in the human genome.

Many human cancers are associated with characteristic chromosomal rearrangements, especially hematopoietic cancers such as leukemias and lymphomas. The first and most critical step in the rearrangement process is the induction of two DNA double-strand breaks (DSB). In all cases, at least one of the two DSBs is generated by a pathologic process, such as (1) randomly-positioned breaks due to ionizing radiation, free radical oxidative damage, or spontaneous hydrolysis; (2) breaks associated with topoisomerase inhibitor treatment; or (3) breaks at direct or inverted repeat sequences, mediated by unidentified strand breakage mechanisms. In lymphoid cells, one of the two requisite DSBs is often physiologic, the result of V(D)J recombination or class switch recombination (CSR) at the lymphoid antigen receptor loci. The RAG complex, which causes the DSBs in V(D)J recombination, can cause (4) sequence-specific, pathologic DSBs at sites that fit the consensus of their normal V(D)J recombination signal targets; or (5) structure-specific, pathologic DSBs at regions of single- to double-strand transition. CSR occurs specifically in the B-cell lineage, and requires (6) activation-induced cytidine deaminase (AID) action at sites of single-stranded DNA, which may occur pathologically outside of the normal target loci of class switch recombination regions and somatic hypermutation (SHM) zones. Recent work proposes a seventh mechanism: the sequential action of AID and the RAG complex at CpG sites provides a coherent model for the pathologic DSBs at some of the most common sites of translocation in human lymphoma - the bcl-2 gene in follicular lymphoma and diffuse large B-cell lymphoma, and the bcl-1 gene in mantle cell lymphoma.

Mass Cytometry Phenotyping of Human Granulocytes Reveals Novel Basophil Functional Heterogeneity.

Basophils, the rarest granulocyte, play critical roles in parasite- and allergen-induced inflammation. We applied mass cytometry (CyTOF) to simultaneously asses 44 proteins to phenotype and functionally characterize neutrophils, eosinophils, and basophils from 19 healthy donors. There was minimal heterogeneity seen in eosinophils and neutrophils, but data-driven analyses revealed four unique subpopulations within phenotypically basophilic granulocytes (PBG; CD45+HLA-DR-CD123+). Through CyTOF and fluorescence-activated cell sorting (FACS), we classified these four PBG subpopulations as (I) CD16lowFcεRIhighCD244high (88.5 ± 1.2%), (II) CD16highFcεRIhighCD244high (9.1 ± 0.4%), (III) CD16lowFcεRIlowCD244low (2.3 ± 1.3), and (IV) CD16highFcεRIlowCD244low (0.4 ± 0.1%). Prospective isolation confirmed basophilic-morphology of PBG I-III, but neutrophilic-morphology of PBG IV. Functional interrogation via IgE-crosslinking or IL-3 stimulation demonstrated that PBG I-II had significant increases in CD203c expression, whereas PBG III-IV remained unchanged compared with media-alone conditions. Thus, PBG III-IV could serve roles in non-IgE-mediated immunity. Our findings offer new perspectives in human basophil heterogeneity and the varying functional potential of these new subsets in health and disease.

Multiplexed single-cell morphometry for hematopathology diagnostics.

The diagnosis of lymphomas and leukemias requires hematopathologists to integrate microscopically visible cellular morphology with antibody-identified cell surface molecule expression. To merge these into one high-throughput, highly multiplexed, single-cell assay, we quantify cell morphological features by their underlying, antibody-measurable molecular components, which empowers mass cytometers to 'see' like pathologists. When applied to 71 diverse clinical samples, single-cell morphometric profiling reveals robust and distinct patterns of 'morphometric' markers for each major cell type. Individually, lamin B1 highlights acute leukemias, lamin A/C helps distinguish normal from neoplastic mature T cells, and VAMP-7 recapitulates light-cytometric side scatter. Combined with machine learning, morphometric markers form intuitive visualizations of normal and neoplastic cellular distribution and differentiation. When recalibrated for myelomonocytic blast enumeration, this approach is superior to flow cytometry and comparable to expert microscopy, bypassing years of specialized training. The contextualization of traditional surface markers on independent morphometric frameworks permits more sensitive and automated diagnosis of complex hematopoietic diseases.

Single-stranded DNA ligation and XLF-stimulated incompatible DNA end ligation by the XRCC4-DNA ligase IV complex: influence of terminal DNA sequence.

The double-strand DNA break repair pathway, non-homologous DNA end joining (NHEJ), is distinctive for the flexibility of its nuclease, polymerase and ligase activities. Here we find that the joining of ends by XRCC4-ligase IV is markedly influenced by the terminal sequence, and a steric hindrance model can account for this. XLF (Cernunnos) stimulates the joining of both incompatible DNA ends and compatible DNA ends at physiologic concentrations of Mg2+, but only of incompatible DNA ends at higher concentrations of Mg2+, suggesting charge neutralization between the two DNA ends within the ligase complex. XRCC4-DNA ligase IV has the distinctive ability to ligate poly-dT single-stranded DNA and long dT overhangs in a Ku- and XLF-independent manner, but not other homopolymeric DNA. The dT preference of the ligase is interesting given the sequence bias of the NHEJ polymerase. These distinctive properties of the XRCC4-DNA ligase IV complex explain important aspects of its in vivo roles.

Impact of somatic and germline mutations on the outcome of systemic mastocytosis.

Systemic mastocytosis (SM) is a highly heterogeneous disease with indolent and aggressive forms, with the mechanisms leading to malignant transformation still remaining to be elucidated. Here, we investigated the presence and frequency of genetic variants in 34 SM patients with multilineal KIT D816V mutations. Initial screening was performed by targeted sequencing of 410 genes in DNA extracted from purified bone marrow cells and hair from 12 patients with nonadvanced SM and 8 patients with advanced SM, followed by whole-genome sequencing (WGS) in 4 cases. Somatic mutations were further investigated in another 14 patients with advanced SM. Despite the fact that no common mutation other than KIT D816V was found in WGS analyses, targeted next-generation sequencing identified 67 nonsynonymous genetic variants involving 39 genes. Half of the mutations were somatic (mostly multilineal), whereas the other half were germline variants. The presence of ≥1 multilineal somatic mutation involving genes other than KIT D816V, ≥3 germline variants, and ≥1 multilineal mutation in the SRSF2, ASXL1, RUNX1, and/or EZH2 genes (S/A/R/E genes), in addition to skin lesions, splenomegaly, thrombocytopenia, low hemoglobin levels, and increased alkaline phosphatase and ß2-microglobulin serum levels, were associated with a poorer patient outcome. However, the presence of ≥1 multilineal mutation, particularly involving S/A/R/E genes, was the only independent predictor for progression-free survival and overall survival in our cohort.