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

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.

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.

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.

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.

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.

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.

Human lymphoid translocation fragile zones are hypomethylated and have accessible chromatin.

Chromosomal translocations are a hallmark of hematopoietic malignancies. CG motifs within translocation fragile zones (typically 20 to 600 bp in size) are prone to chromosomal translocation in lymphomas. Here we demonstrate that the CG motifs in human translocation fragile zones are hypomethylated relative to the adjacent DNA. Using a methyltransferase footprinting assay on isolated nuclei (in vitro), we find that the chromatin at these fragile zones is accessible. We also examined in vivo accessibility using cellular expression of a prokaryotic methylase. Based on this assay, which measures accessibility over a much longer time interval than is possible with in vitro methods, these fragile zones were found to be more accessible than the adjacent DNA. Because DNA within the fragile zones can be methylated by both cellular and exogenous methyltransferases, the fragile zones are predominantly in a duplex DNA conformation. These observations permit more-refined models for why these zones are 100- to 1,000-fold more prone to undergo chromosomal translocation than the adjacent regions.

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.

Both CpG methylation and activation-induced deaminase are required for the fragility of the human bcl-2 major breakpoint region: implications for the timing of the breaks in the t(14;18) translocation.

The t(14;18) chromosomal translocation typically involves breakage at the bcl-2 major breakpoint region (MBR) to cause human follicular lymphoma. A theory to explain the striking propensity of the MBR breaks at three CpG clusters within the 175-bp MBR region invoked activation-induced deaminase (AID). In a test of that theory, we used here minichromosomal substrates in human pre-B cell lines. Consistent with the essential elements of the theory, we found that the MBR breakage process is indeed highly dependent on DNA methylation at the CpG sites and highly dependent on the AID enzyme to create lesions at peak locations within the MBR. Interestingly, breakage of the phosphodiester bonds at the AID-initiated MBR lesions is RAG dependent, but, unexpectedly, most are also dependent on Artemis. We found that Artemis is capable of nicking small heteroduplex structures and is even able to nick single-base mismatches. This raises the possibility that activated Artemis, derived from the unjoined D to J(H) DNA ends at the IgH locus on chromosome 14, nicks AID-generated TG mismatches at methyl CpG sites, and this would explain why the breaks at the chromosome 18 MBR occur within the same time window as those on chromosome 14.

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.

Heterogeneity and randomness of DNA methylation patterns in human embryonic stem cells.

DNA methylation has been proposed to be important in many biological processes and is the subject of intense study. Traditional bisulfite genomic sequencing allows detailed high-resolution methylation pattern analysis of each molecule with haplotype information across a few hundred bases at each locus, but lacks the capacity to gather voluminous data. Although recent technological developments are aimed at assessing DNA methylation patterns in a high-throughput manner across the genome, the haplotype information cannot be accurately assembled when the sequencing reads are short or when each hybridization target only includes one or two cytosine-phosphate-guanine (CpG) sites. Whether a distinct and nonrandom DNA methylation pattern is present at a given locus is difficult to discern without the haplotype information, and the DNA methylation patterns are much less apparent because the data are often obtained only as methylation frequencies at each CpG site with some of these methods. It would facilitate the interpretation of data obtained from high-throughput bisulfite sequencing if the loci with nonrandom DNA methylation patterns could be distinguished from those that are randomly methylated. In this study, we carried out traditional genomic bisulfite sequencing using the normal diploid human embryonic stem (hES) cell lines, and utilized Hamming distance analysis to evaluate the existence of a distinct and nonrandom DNA methylation pattern at each locus studied. Our findings suggest that Hamming distance is a simple, quick, and useful tool to identify loci with nonrandom DNA methylation patterns and may be utilized to discern links between biological changes and DNA methylation patterns in the high-throughput bisulfite sequencing data sets.

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.

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.

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.