ShinySOM: graphical SOM-based analysis of single-cell cytometry data.

SUMMARY: ShinySOM offers a user-friendly interface for reproducible, high-throughput analysis of high-dimensional flow and mass cytometry data guided by self-organizing maps. The software implements a FlowSOM-style workflow, with improvements in performance, visualizations and data dissection possibilities. The outputs of the analysis include precise statistical information about the dissected samples, and R-compatible metadata useful for the batch processing of large sample volumes. AVAILABILITY AND IMPLEMENTATION: ShinySOM is free and open-source, available online at gitlab.com/exaexa/ShinySOM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DUX4r, ZNF384r and PAX5-P80R mutated B-cell precursor acute lymphoblastic leukemia frequently undergo monocytic switch.

Recently, we described B-cell precursor acute lymphoblastic leukemia (BCP-ALL) subtype with early switch to the monocytic lineage and loss of the B-cell immunophenotype, including CD19 expression. Thus far, the genetic background has remained unknown. Among 726 children consecutively diagnosed with BCP-ALL, 8% patients experienced switch detectable by flow cytometry (FC). Using exome and RNA sequencing, switch was found to positively correlate with three different genetic subtypes: PAX5-P80R mutation (5 cases with switch out of 5), rearranged DUX4 (DUX4r; 30 cases of 41) and rearranged ZNF384 (ZNF384r; 4 cases of 10). Expression profiles or phenotypic patterns correlated with genotypes, but within each genotype they could not identify cases who subsequently switched. If switching was not taken into account, the B-cell-oriented FC assessment underestimated the minimal residual disease level. For patients with PAX5-P80R, a discordance between FC-determined and PCR-determined MRD was found on day 15, resulting from a rapid loss of the B-cell phenotype. Discordance on day 33 was observed in all the DUX4r, PAX5-P80R and ZNF384r subtypes. Importantly, despite the substantial phenotypic changes, possibly even challenging the appropriateness of BCP-ALL therapy, the monocytic switch was not associated with a higher incidence of relapse and poorer prognosis in patients undergoing standard ALL treatment.

Nuclear pore protein TPR associates with lamin B1 and affects nuclear lamina organization and nuclear pore distribution.

The organization of the nuclear periphery is crucial for many nuclear functions. Nuclear lamins form dense network at the nuclear periphery and play a substantial role in chromatin organization, transcription regulation and in organization of nuclear pore complexes (NPCs). Here, we show that TPR, the protein located preferentially within the nuclear baskets of NPCs, associates with lamin B1. The depletion of TPR affects the organization of lamin B1 but not lamin A/C within the nuclear lamina as shown by stimulated emission depletion microscopy. Finally, reduction of TPR affects the distribution of NPCs within the nuclear envelope and the effect can be reversed by simultaneous knock-down of lamin A/C or the overexpression of lamin B1. Our work suggests a novel role for the TPR at the nuclear periphery: the TPR contributes to the organization of the nuclear lamina and in cooperation with lamins guards the interphase assembly of nuclear pore complexes.

What matters in chronic Burkholderia cenocepacia infection in cystic fibrosis: Insights from comparative genomics.

Burkholderia cenocepacia causes severe pulmonary infections in cystic fibrosis (CF) patients. Since the bacterium is virtually untreatable by antibiotics, chronic infections persist for years and might develop into fatal septic pneumonia (cepacia syndrome, CS). To devise new strategies to combat chronic B. cenocepacia infections, it is essential to obtain comprehensive knowledge about their pathogenesis. We conducted a comparative genomic analysis of 32 Czech isolates of epidemic clone B. cenocepacia ST32 isolated from various stages of chronic infection in 8 CF patients. High numbers of large-scale deletions were found to occur during chronic infection, affecting preferentially genomic islands and nonessential replicons. Recombination between insertion sequences (IS) was inferred as the mechanism behind deletion formation; the most numerous IS group was specific for the ST32 clone and has undergone transposition burst since its divergence. Genes functionally related to transition metal metabolism were identified as hotspots for deletions and IS insertions. This functional category was also represented among genes where nonsynonymous point mutations and indels occurred parallelly among patients. Another category exhibiting parallel mutations was oxidative stress protection; mutations in catalase KatG resulted in impaired detoxification of hydrogen peroxide. Deep sequencing revealed substantial polymorphism in genes of both categories within the sputum B. cenocepacia ST32 populations, indicating extensive adaptive evolution. Neither oxidative stress response nor transition metal metabolism genes were previously reported to undergo parallel evolution during chronic CF infection. Mutations in katG and copper metabolism genes were overrepresented in patients where chronic infection developed into CS. Among professional phagocytes, macrophages use both hydrogen peroxide and copper for their bactericidal activity; our results thus tentatively point to macrophages as suspects in pathogenesis towards the fatal CS.

ERG deletions in childhood acute lymphoblastic leukemia with DUX4 rearrangements are mostly polyclonal, prognostically relevant and their detection rate strongly depends on screening method sensitivity.

ERG-deletions occur recurrently in acute lymphoblastic leukemia, especially in the DUX4-rearranged subtype. The ERG-deletion was shown to positively impact prognosis of patients with IKZF1-deletion and its presence precludes assignment into IKZF1 plus group, a novel high-risk category on AIEOP-BFM ALL trials. We analyzed the impact of different methods on ERG-deletion detection rate, evaluated ERG-deletion as a potential marker for DUX4-rearranged leukemia, studied its associations with molecular and clinical characteristics within this leukemia subtype, and analyzed its clonality. Using single-nucleotide-polymorphism array, genomic polymerase chain reaction (PCR) and amplicon-sequencing we found ERG-deletion in 34% (16 of 47), 66% (33 of 50) and 78% (39 of 50) of DUX4-rearranged leukemia, respectively. False negativity of ERG-deletion by single-nucleotide-polymorphism array caused IKZF1 plus misclassification in 5 patients. No ERG-deletion was found outside the DUX4-rearranged cases. Within DUX4-rearranged leukemia, the ERG-deletion was associated with higher total number of copy-number aberrations, and, importantly, the ERG-deletion positivity by PCR was associated with better outcome [5-year event-free survival (EFS), ERG-deletion-positive 93% vs. ERG-deletion-negative 68%, P=0.022; 5-year overall survival (OS), ERG-deletion-positive 97% vs. ERG-deletion-negative 75%, P=0.029]. Ultra-deep amplicon-sequencing revealed distinct co-existing ERG-deletions in 22 of 24 patients. In conclusion, our data demonstrate inadequate sensitivity of single-nucleotide-polymorphism array for ERG-deletion detection, unacceptable for proper IKZF1 plus classification. Even using more sensitive methods (PCR/amplicon-sequencing) for its detection, ERG-deletion is absent in 22-34% of DUX4-rearranged leukemia and does not represent an adequately sensitive marker of this leukemia subtype. Importantly, the ERG-deletion potentially stratifies the DUX4-rearranged leukemia into biologically/clinically distinct subsets. Frequent polyclonal pattern of ERG-deletions shows that late origin of this lesion is more common than has been previously described.

Genomic landscape of pediatric B-other acute lymphoblastic leukemia in a consecutive European cohort.

Novel biological subtypes and clinically important genetic aberrations (druggable lesions, prognostic factors) have been described in B-other acute lymphoblastic leukemia (ALL) during the last decade; however, due to a lack of studies on unselected cohorts, their population frequency and mutual associations still have to be established. We studied 110 consecutively diagnosed and uniformly treated childhood B-other patients using single nucleotide polymorphism arrays and whole exome/transcriptome sequencing. The frequency of DUX4-rearranged, BCR-ABL1-like, ZNF384-rearranged, ETV6-RUNX1-like, iAMP21 and MEF2D-rearranged subtypes was 27%, 15%, 5%, 5%, 4%, and 2%, respectively; 43% of cases were not classified into any of these subtypes (B-rest). We found worse early response to treatment in DUX4-rearranged leukemia and a strong association of ZNF384-rearranged leukemia with B-myeloid immunophenotype. Of the druggable lesions, JAK/STAT-class and RAS/RAF/MAPK-class aberrations were found in 21% and 43% of patients, respectively; an ABL-class aberration was found in one patient. A recently described negative prognostic factor, IKZF1plus , was found in 14% of patients and was enriched in (but not exclusive for) BCR-ABL1-like subtype. PAX5 fusions (including 4 novel), intragenic amplifications and P80R mutations were mutually exclusive and only occurred in the B-rest subset, altogether accounting for 20% of the B-other group. PAX5 P80R was associated with a specific gene expression signature, potentially defining a novel leukemia subtype. Our study shows unbiased European population-based frequencies of novel ALL subtypes, recurrent (cyto)genetic aberrations and their mutual associations. This study also strengthens and widens the current knowledge of B-other ALL and provides an objective basis for optimization of current genetic diagnostics.

flowIO: Flow cytometry standard conformance testing, editing, and export tool.

The Flow Cytometry Standard (FCS) format is a widely accepted norm for storing Flow Cytometry (FCM) data. Its goal as a standard is to allow FCM data sharing and re-analysis. Over more than three decades of its existence FCS has evolved into a well-defined, flexible file format reflecting technical changes in the FCM field. Its flexibility as well as rising numbers of instrument vendors leads to suboptimal implementations of FCS in some cases. Such situations compromise the primary goal of the standard and hinder the ability to reproduce FCM analyses. It is further underlined by rapid rise of advanced FCM analyses, often carried out outside traditional software tools and heavily relying on standard data storage and presentation. We have developed flowIO, an R package which tests FCS file conformance with the standard as defined by International Society for Advancement of Cytometry (ISAC) normative. Along with the package we provide a web based application (also at http://bioinformin.cesnet.cz/flowIO/) allowing user friendly access to the conformance testing as well as FCS file editing and export for further analysis.

High-resolution Antibody Array Analysis of Childhood Acute Leukemia Cells.

Acute leukemia is a disease pathologically manifested at both genomic and proteomic levels. Molecular genetic technologies are currently widely used in clinical research. In contrast, sensitive and high-throughput proteomic techniques for performing protein analyses in patient samples are still lacking. Here, we used a technology based on size exclusion chromatography followed by immunoprecipitation of target proteins with an antibody bead array (Size Exclusion Chromatography-Microsphere-based Affinity Proteomics, SEC-MAP) to detect hundreds of proteins from a single sample. In addition, we developed semi-automatic bioinformatics tools to adapt this technology for high-content proteomic screening of pediatric acute leukemia patients.To confirm the utility of SEC-MAP in leukemia immunophenotyping, we tested 31 leukemia diagnostic markers in parallel by SEC-MAP and flow cytometry. We identified 28 antibodies suitable for both techniques. Eighteen of them provided excellent quantitative correlation between SEC-MAP and flow cytometry (p< 0.05). Next, SEC-MAP was applied to examine 57 diagnostic samples from patients with acute leukemia. In this assay, we used 632 different antibodies and detected 501 targets. Of those, 47 targets were differentially expressed between at least two of the three acute leukemia subgroups. The CD markers correlated with immunophenotypic categories as expected. From non-CD markers, we found DBN1, PAX5, or PTK2 overexpressed in B-cell precursor acute lymphoblastic leukemias, LAT, SH2D1A, or STAT5A overexpressed in T-cell acute lymphoblastic leukemias, and HCK, GLUD1, or SYK overexpressed in acute myeloid leukemias. In addition, OPAL1 overexpression corresponded to ETV6-RUNX1 chromosomal translocation.In summary, we demonstrated that SEC-MAP technology is a powerful tool for detecting hundreds of proteins in clinical samples obtained from pediatric acute leukemia patients. It provides information about protein size and reveals differences in protein expression between particular leukemia subgroups. Forty-seven of SEC-MAP identified targets were validated by other conventional method in this study.

Myocardial ischemia and reperfusion leads to transient CD8 immune deficiency and accelerated immunosenescence in CMV-seropositive patients.

RATIONALE: There is mounting evidence of a higher incidence of coronary heart disease in cytomegalovirus-seropositive individuals. OBJECTIVE: The aim of this study was to investigate whether acute myocardial infarction triggers an inflammatory T-cell response that might lead to accelerated immunosenescence in cytomegalovirus-seropositive patients. METHODS AND RESULTS: Thirty-four patients with acute myocardial infarction undergoing primary percutaneous coronary intervention were longitudinally studied within 3 months after reperfusion (Cohort A). In addition, 54 patients with acute myocardial infarction and chronic myocardial infarction were analyzed in a cross-sectional study (Cohort B). Cytomegalovirus-seropositive patients demonstrated a greater fall in the concentration of terminally differentiated CD8 effector memory T cells (TEMRA) in peripheral blood during the first 30 minutes of reperfusion compared with cytomegalovirus-seronegative patients (-192 versus -63 cells/µL; P=0.008), correlating with the expression of programmed cell death-1 before primary percutaneous coronary intervention (r=0.8; P=0.0002). A significant proportion of TEMRA cells remained depleted for ≥3 months in cytomegalovirus-seropositive patients. Using high-throughput 13-parameter flow cytometry and human leukocyte antigen class I cytomegalovirus-specific dextramers, we confirmed an acute and persistent depletion of terminally differentiated TEMRA and cytomegalovirus-specific CD8(+) cells in cytomegalovirus-seropositive patients. Long-term reconstitution of the TEMRA pool in chronic cytomegalovirus-seropositive postmyocardial infarction patients was associated with signs of terminal differentiation including an increase in killer cell lectin-like receptor subfamily G member 1 and shorter telomere length in CD8(+) T cells (2225 versus 3397 bp; P<0.001). CONCLUSIONS: Myocardial ischemia and reperfusion in cytomegalovirus-seropositive patients undergoing primary percutaneous coronary intervention leads to acute loss of antigen-specific, terminally differentiated CD8 T cells, possibly through programmed cell death-1-dependent programmed cell death. Our results suggest that acute myocardial infarction and reperfusion accelerate immunosenescence in cytomegalovirus-seropositive patients.

Mantle cell lymphoma-variant Richter syndrome: Detailed molecular-cytogenetic and backtracking analysis reveals slow evolution of a pre-MCL clone in parallel with CLL over several years.

Richter syndrome represents the transformation of the chronic lymphocytic leukemia (CLL) into an aggressive lymphoma, most frequently the diffuse large B-cell lymphoma (DLBCL). In this report we describe a patient with CLL, who developed a clonally-related pleomorphic highly-aggressive mantle cell lymphoma (MCL) after five cycles of a fludarabine-based second-line therapy for the first relapse of CLL. Molecular cytogenetic methods together with whole-exome sequencing revealed numerous gene alterations restricted to the MCL clone (apart from the canonical t(11;14)(q13;q32) translocation) including gain of one copy of ATM gene or emergence of TP53, CREBBP, NUP214, FUBP1 and SF3B1 gene mutations. Similarly, gene expression analysis revealed vast differences between the MCL and CLL transcriptome, including overexpression of cyclin D1, downregulation of cyclins D2 and D3, or downregulation of IL4R in the MCL clone. Backtracking analysis using quantitative PCR specifically detecting an MCL-restricted focal deletion of TP53 revealed that the pre-MCL clone appeared in the bone marrow and peripheral blood of the patient approximately 4 years before the clinical manifestation of MCL. Both molecular cytogenetic and sequencing data support the hypothesis of a slow development of the pre-MCL clone in parallel to CLL over several years, and thereby exclude the possibility that the transformation event occurred at the stage of the CLL relapse clone by mere t(11;14)(q13;q32) acquisition.

Detection and monitoring of normal and leukemic cell populations with hierarchical clustering of flow cytometry data.

Flow cytometry is a valuable tool in research and diagnostics including minimal residual disease (MRD) monitoring of hematologic malignancies. However, its gradual advancement toward increasing numbers of fluorescent parameters leads to information rich datasets, which are challenging to analyze by standard gating and do not reflect the multidimensionality of the data. We have developed a novel method to analyze complex flow cytometry data, based on hierarchical clustering analysis (HCA) but with a new underlying algorithm, using Mahalanobis distance measure. HCA is scalable to analyze complex multiparameter datasets (here demonstrated on up to 12 color flow cytometry and on a 20-parameter synthetic dataset). We have validated this method by comparison with standard gating approaches when performed independently by expert cytometrists. Acute lymphoblastic leukemia blast populations were analyzed in diagnostic and follow-up datasets (n = 123) from three centers. HCA results correlated very well (Passing-Bablok correlation coefficient = 0.992, slope = 1, intercept = -0.01) with standard gating data obtained by the I-BFM FLOW-MRD study group. To further improve the performance in follow-up samples with low MRD levels and to automate MRD detection, we combined HCA with support vector machine (SVM) learning. HCA in combination with SVM provides a novel diagnostic tool that not only allows analysis of increasingly complex flow cytometry data but also is less observer-dependent compared with classical gating and has potential for automation.

An automated analysis of highly complex flow cytometry-based proteomic data.

The combination of color-coded microspheres as carriers and flow cytometry as a detection platform provides new opportunities for multiplexed measurement of biomolecules. Here, we developed a software tool capable of automated gating of color-coded microspheres, automatic extraction of statistics from all subsets and validation, normalization, and cross-sample analysis. The approach presented in this article enabled us to harness the power of high-content cellular proteomics. In size exclusion chromatography-resolved microsphere-based affinity proteomics (Size-MAP), antibody-coupled microspheres are used to measure biotinylated proteins that have been separated by size exclusion chromatography. The captured proteins are labeled with streptavidin phycoerythrin and detected by multicolor flow cytometry. When the results from multiple size exclusion chromatography fractions are combined, binding is detected as discrete reactivity peaks (entities). The information obtained might be approximated to a multiplexed western blot. We used a microsphere set with >1,000 subsets, presenting an approach to extract biologically relevant information. The R-project environment was used to sequentially recognize subsets in two-dimensional space and gate them. The aim was to extract the median streptavidin phycoerythrin fluorescence intensity for all 1,000+ microsphere subsets from a series of 96 measured samples. The resulting text files were subjected to algorithms that identified entities across the 24 fractions. Thus, the original 24 data points for each antibody were compressed to 1-4 integrated values representing the areas of individual antibody reactivity peaks. Finally, we provide experimental data on cellular protein changes induced by treatment of leukemia cells with imatinib mesylate. The approach presented here exemplifies how large-scale flow cytometry data analysis can be efficiently processed to employ flow cytometry as a high-content proteomics method.

Standardization of Workflow and Flow Cytometry Panels for Quantitative Expression Profiling of Surface Antigens on Blood Leukocyte Subsets: An HCDM CDMaps Initiative.

Background: The Human Cell Differentiation Molecules (HCDM) organizes Human Leukocyte Differentiation Antigen (HLDA) workshops to test and name clusters of antibodies that react with a specific antigen. These cluster of differentiation (CD) markers have provided the scientific community with validated antibody clones, consistent naming of targets and reproducible identification of leukocyte subsets. Still, quantitative CD marker expression profiles and benchmarking of reagents at the single-cell level are currently lacking. Objective: To develop a flow cytometric procedure for quantitative expression profiling of surface antigens on blood leukocyte subsets that is standardized across multiple research laboratories. Methods: A high content framework to evaluate the titration and reactivity of Phycoerythrin (PE)-conjugated monoclonal antibodies (mAbs) was created. Two flow cytometry panels were designed: an innate cell tube for granulocytes, dendritic cells, monocytes, NK cells and innate lymphoid cells (12-color) and an adaptive lymphocyte tube for naive and memory B and T cells, including TCRγδ+, regulatory-T and follicular helper T cells (11-color). The potential of these 2 panels was demonstrated via expression profiling of selected CD markers detected by PE-conjugated antibodies and evaluated using 561 nm excitation. Results: Using automated data annotation and dried backbone reagents, we reached a robust workflow amenable to processing hundreds of measurements in each experiment in a 96-well plate format. The immunophenotyping panels enabled discrimination of 27 leukocyte subsets and quantitative detection of the expression of PE-conjugated CD markers of interest that could quantify protein expression above 400 units of antibody binding capacity. Expression profiling of 4 selected CD markers (CD11b, CD31, CD38, CD40) showed high reproducibility across centers, as well as the capacity to benchmark unique clones directed toward the same CD3 antigen. Conclusion: We optimized a procedure for quantitative expression profiling of surface antigens on blood leukocyte subsets. The workflow, bioinformatics pipeline and optimized flow panels enable the following: 1) mapping the expression patterns of HLDA-approved mAb clones to CD markers; 2) benchmarking new antibody clones to established CD markers; 3) defining new clusters of differentiation in future HLDA workshops.

Nucleoporin TPR Affects C2C12 Myogenic Differentiation via Regulation of Myh4 Expression.

The nuclear pore complex (NPC) has emerged as a hub for the transcriptional regulation of a subset of genes, and this type of regulation plays an important role during differentiation. Nucleoporin TPR forms the nuclear basket of the NPC and is crucial for the enrichment of open chromatin around NPCs. TPR has been implicated in the regulation of transcription; however, the role of TPR in gene expression and cell differentiation has not been described. Here we show that depletion of TPR results in an aberrant morphology of murine proliferating C2C12 myoblasts (MBs) and differentiated C2C12 myotubes (MTs). The ChIP-Seq data revealed that TPR binds to genes linked to muscle formation and function, such as myosin heavy chain (Myh4), myocyte enhancer factor 2C (Mef2C) and a majority of olfactory receptor (Olfr) genes. We further show that TPR, possibly via lysine-specific demethylase 1 (LSD1), promotes the expression of Myh4 and Olfr376, but not Mef2C. This provides a novel insight into the mechanism of myogenesis; however, more evidence is needed to fully elucidate the mechanism by which TPR affects specific myogenic genes.

A novel class of ZNF384 aberrations in acute leukemia.

Fusion of the ZNF384 gene as the 3' partner to several different 5' partner genes occurs recurrently in B-cell precursor acute lymphoblastic and mixed phenotype B/myeloid leukemia. These canonical fusions (ZNF384r) contain the complete ZNF384 coding sequence and are associated with a specific gene expression signature. Cases with this signature, but without canonical ZNF384 fusions (ZNF384r-like cases), have been described previously. Although some have been shown to harbor ZNF362 fusions, the primary aberrations remain unknown in a major proportion. We studied 3 patients with the ZNF384r signature and unknown primary genetic background and identified a previously unknown class of genetic aberration affecting the last exon of ZNF384 and resulting in disruption of the C-terminal portion of the ZNF384 protein. Importantly, in 2 cases, the ZNF384 aberration, indel, was missed during the bioinformatic analysis but revealed by the manual, targeted reanalysis. Two cases with the novel aberrations had a mixed (B/myeloid) immunophenotype commonly associated with canonical ZNF384 fusions. In conclusion, we present leukemia cases with a novel class of ZNF384 aberrations that phenocopy leukemia with ZNF384r. Therefore, we show that part of the so-called ZNF384r-like cases represent the same genetic subtype as leukemia with canonical ZNF384 fusions.