Longitudinal minimal residual disease assessment in multiple myeloma patients in complete remission - results from the NMSG flow-MRD substudy within the EMN02/HO95 MM trial.

BACKGROUND: Multiple myeloma remains an incurable disease with multiple relapses due to residual myeloma cells in the bone marrow of patients after therapy. Presence of small number of cancer cells in the body after cancer treatment, called minimal residual disease, has been shown to be prognostic for progression-free and overall survival. However, for multiple myeloma, it is unclear whether patients attaining minimal residual disease negativity may be candidates for treatment discontinuation. We investigated, if longitudinal flow cytometry-based monitoring of minimal residual disease (flow-MRD) may predict disease progression earlier and with higher sensitivity compared to biochemical assessments. METHODS: Patients from the Nordic countries with newly diagnosed multiple myeloma enrolled in the European-Myeloma-Network-02/Hovon-95 (EMN02/HO95) trial and undergoing bone marrow aspiration confirmation of complete response, were eligible for this Nordic Myeloma Study Group (NMSG) substudy. Longitdudinal flow-MRD assessment of bone marrow samples was performed to identify and enumerate residual malignant plasma cells until observed clinical progression. RESULTS: Minimal residual disease dynamics were compared to biochemically assessed changes in serum free light chain and M-component. Among 20 patients, reaching complete response or stringent complete response during the observation period, and with ≥3 sequential flow-MRD assessments analysed over time, increasing levels of minimal residual disease in the bone marrow were observed in six cases, preceding biochemically assessed disease and clinical progression by 5.5 months and 12.6 months (mean values), respectively. Mean malignant plasma cells doubling time for the six patients was 1.8 months (95% CI, 1.4-2.3 months). Minimal malignant plasma cells detection limit was 4 × 10-5. CONCLUSIONS: Flow-MRD is a sensitive method for longitudinal monitoring of minimal residual disease dynamics in multiple myeloma patients in complete response. Increasing minimal residual disease levels precedes biochemically assessed changes and is an early indicator of subsequent clinical progression. TRIAL REGISTRATION: NCT01208766.

Utility of a new notation to visualize flow cytometry analysis results: first preliminary comparison with immunohistochemistry to detect CD30 expression on T-cell lymphoma cells.

BACKGROUND: It is important to confirm CD30 expression in T-cell lymphoma cases, but immunohistochemical staining for CD30 is not commonly performed and no comparison has been done between the results of flow cytometry (FCM) and immunohistochemical staining for CD30. Therefore, we devised a notation that we termed proportion of immunoreactivity/expression for FCM (PRIME-F notation), based on the cellular proportion showing different antigen-antibody reactivity. METHODS: We retrospectively compiled 211 cases of T-cell lymphoma, assessed via FCM, from major hospitals in Miyagi Prefecture from January 2012 to January 2019, and compared 52 of these cases with the immunohistochemical immunoreactive (IR) pattern of CD30 (PRIME-I notation). The PRIME-F notation was divided into five levels: notations starting with "-" followed by 3, 2, and 1 ">" correspond to level-I, level-II, or level-III; notations starting with "(dim)+" correspond to level-IV; and those starting with "+" or "(bright)+" correspond to level-V. RESULTS: The 52 cases of PRIME-F notation with "+" included 16 cases of peripheral T-cell lymphoma (PTCL/NOS), 3 of follicular T-cell lymphoma (FTL), 3 of angioimmunoblastic T-cell lymphoma (AITL), 6 of extranodal NK/T-cell lymphoma/nasal type (ENKL), 18 of adult T-cell lymphoma (ATL), and 6 cases of anaplastic large cell lymphoma (ALCL). Eight of the 52 cases were immunohistochemically CD30-negative. In the PRIME-F level-I to III group (excluding false-positive cases), 21.7% (5 out of 23 cases) were < 10% positive for CD30 upon immunohistochemistry (IHC). Contrarily, in the level-IV & -V group, no CD30 positivity rate of < 10% upon IHC was found (0%) (p = 0.0497). In level-IV, 42.9% of cases presented a CD30 negative rate > 1/3 upon IHC, while in level-V, only 7.1% (one out of 14 cases) did. The CD30 negative rate tended to be low (p = 0.0877) in level-V. CONCLUSIONS: To our knowledge, this is the first report describing the correspondence between FCM and immunohistochemistry findings for CD30 through newly proposed notations. The PRIME-F and PRIME-I notations for CD30 showed a minor positive correlation. The PRIME notation is considered universally applicable to antibodies, and notations of both FCM and IHC show great potential for big data.

Deep phenotyping of T cell populations under long-term treatment of tacrolimus and rapamycin in patients receiving renal transplantations by mass cytometry.

Tacrolimus (FK506) and rapamycin (RAPA) are widely used to maintain long-term immunosuppression after organ transplantation. However, the impact of accumulative drug administration on the recipients' immune systems remains unclear. We investigated the impact of 3-year FK506 or RAPA treatment after renal transplantation on the human immune systems. A discovery cohort of 30 patients was first recruited, and we discovered two distinctive T lineage suppressive regulatory patterns induced by chronic treatment of FK506 and RAPA. The increased percentage of senescent CD8+ CD57+ T lineages and less responsive T cell receptor (TCR) pathway in the FK506 group indicate better graft acceptance. Meanwhile, percentages of regulatory T cells (Tregs) and expression of CTLA-4 were both up to two-fold higher in the RAPA group, suggesting the inconsistent reactivation potential of the FK506 and RAPA groups when an anti-tumour or anti-infection immune response is concerned. Additionally, up-regulation of phosphorylated signaling proteins in T lineages after in vitro CD3/CD28 stimulation suggested more sensitive TCR-signaling pathways reserved in the RAPA group. An independent validation cohort of 100 renal transplantation patients was further investigated for the hypothesis that long-term RAPA administration mitigates the development of tumours and infections during long-term intake of immunosuppressants. Our results indicate that RAPA administration indeed results in less clinical oncogenesis and infection. The deep phenotyping of T-cell lineages, as educated by the long-term treatment of different immunosuppressants, provides new evidence for personalized precision medicine after renal transplantations.

T-cell engager antibodies enable T cells to control HBV infection and to target HBsAg-positive hepatoma in mice.

BACKGROUND & AIMS: Current antiviral therapies control but rarely eliminate HBV, leaving chronic HBV carriers at risk of developing hepatocellular carcinoma (HCC). Lacking or dysfunctional virus-specific adaptive immunity prevents control of HBV and allows the virus to persist. Restoring antiviral T-cell immunity could lead to HBV elimination and cure of chronically infected patients. METHODS: We constructed bispecific T-cell engager antibodies that are designed to induce antiviral immunity through simultaneous binding of HBV envelope proteins (HBVenv) on infected hepatocytes and CD3 or CD28 on T cells. T-cell engager antibodies were employed in co-cultures with healthy donor lymphocytes and HBV-infected target cells. Activation of the T-cell response was determined by detection of pro-inflammatory cytokines, effector function (by cytotoxicity) and antiviral effects. To study in vivo efficacy, immune-deficient mice were transplanted with HBVenv-positive and -negative hepatoma cells. RESULTS: The 2 T-cell engager antibodies synergistically activated T cells to become polyfunctional effectors that in turn elicited potent antiviral effects by killing infected cells and in addition controlled HBV via non-cytolytic, cytokine-mediated antiviral mechanisms. In vivo in mice, the antibodies attracted T cells specifically to the tumors expressing HBVenv resulting in T-cell activation, tumor infiltration and reduction of tumor burden. CONCLUSION: This study demonstrates that the administration of HBVenv-targeting T-cell engager antibodies facilitates a robust T-cell redirection towards HBV-positive target cells and provides a feasible and promising approach for the treatment of chronic viral hepatitis and HBV-associated HCC. LAY SUMMARY: T-cell engager antibodies are an interesting, novel therapeutic tool to restore immunity in patients with chronic hepatitis B. As bispecific antibodies, they bind envelope proteins on the surface of the hepatitis B virus (HBV) and CD3 or CD28 on T cells. This way, they induce a potent antiviral and cytotoxic T-cell response that leads to the elimination of HBV-positive cells. These bispecific T-cell engager antibodies are exciting therapeutic candidates for chronic hepatitis B and HBV-associated hepatocellular carcinoma.

Flow cytometry multiplexed method for the detection of neutralizing human antibodies to the native SARS-CoV-2 spike protein.

A correct identification of seropositive individuals for the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection is of paramount relevance to assess the degree of protection of a human population to present and future outbreaks of the COVID-19 pandemic. We describe here a sensitive and quantitative flow cytometry method using the cytometer-friendly non-adherent Jurkat T-cell line that stably expresses the full-length native spike "S" protein of SARS-CoV-2 and a truncated form of the human EGFR that serves a normalizing role. S protein and huEGFRt coding sequences are separated by a T2A self-cleaving sequence, allowing to accurately quantify the presence of anti-S immunoglobulins by calculating a score based on the ratio of fluorescence intensities obtained by double-staining with the test sera and anti-EGFR. The method allows to detect immune individuals regardless of the result of other serological tests or even repeated PCR monitoring. As examples of its use, we show that as much as 28% of the personnel working at the CBMSO in Madrid is already immune. Additionally, we show that anti-S antibodies with protective neutralizing activity are long-lasting and can be detected in sera 8 months after infection.

Understanding the Activation of Platelets in Diabetes and Its Modulation by Allyl Methyl Sulfide, an Active Metabolite of Garlic.

BACKGROUND: Diabetes mellitus (DM) is a chronic metabolic disorder associated with higher risk of having cardiovascular disease. Platelets play a promising role in the pathogenesis of cardiovascular complications in diabetes. Since last several decades, garlic and its bioactive components are extensively studied in diabetes and its complications. Our aim was to explore the antiplatelet property of allyl methyl sulfide (AMS) focusing on ameliorating platelet activation in diabetes. METHOD: We used streptozotocin- (STZ-) induced diabetic rats as model for type 1 diabetes. We have evaluated the effect of allyl methyl sulfide on platelet activation by administrating AMS to diabetic rats for 10 weeks. Flow cytometry-based analysis was used to evaluate the platelet activation, platelet aggregation, platelet macrophage interaction, and endogenous ROS generation in the platelets obtained from control, diabetes, and AMS- and aspirin-treated diabetic rats. RESULTS: AMS treatment for 10 weeks effectively reduced the blood glucose levels in diabetic rats. Three weeks of AMS (50 mg/kg/day) treatment did not reduce the activation of platelets but a significant (p < 0.05) decrease was observed after 10 weeks of treatment. Oral administration of AMS significantly (p < 0.05) reduced the baseline and also reduced ADP-induced aggregation of platelets after 3 and 10 weeks of treatment. Furthermore, 10 weeks of AMS treatment in diabetic rats attenuated the endogenous ROS content (p < 0.05) of platelets and platelet macrophage interactions. The inhibition of platelet activation in diabetic rats after AMS treatment was comparable with aspirin treatment (30 mg/kg/day). CONCLUSION: We observed an inhibitory effect of allyl methyl sulfide on platelet aggregation, platelet activation, platelet macrophage interaction, and increased ROS levels in type 1 diabetes. Our data suggests that AMS can be useful to control cardiovascular complication in diabetes via inhibition of platelet activation.

Flow Cytometry for Detection and Quantification of Micrometastases in Sentinel Lymph Nodes from Patients with Primary Melanoma.

BACKGROUND: Detection of micrometastases in the regional lymph nodes is one of the most important prognostic factors for melanoma patients. Our aim was to evaluate the suitability of flow cytometry for rapid quantification of disseminated melanoma cells in sentinel lymph nodes (SLN). METHODS: 132 SLNs from 104 patients diagnosed with melanoma were analyzed by flow cytometry, utilizing the extracellular marker melanoma-associated chondroitin sulfate proteoglycan, in addition to quantitative immunocytology and conventional histopathology, including immunohistochemistry. For quantification, the number of melanoma-positive cells per million lymph node cells (disseminated cancer cell density, DCCD) detected by flow cytometry was compared to the DCCD obtained by immunocytology. RESULTS: Compared to histopathology and immunocytology, flow cytometry exhibited a sensitivity of 50% and a specificity of 85%. DCCDs of immunocytology and flow cytometry of the 37 immunocytologically positive SLNs showed a positive correlation (Spearman's ρ = 0.7, P < 0.0001). In 10 SLNs from 9 patients with high tumor load, the flow cytometric DCCD was 8-fold higher on average than the immunocytologic DCCD. CONCLUSIONS: Although flow cytometry is not yet suitable for early detection of metastatic melanoma, it promises to become a valuable tool for rapidly quantifying tumor load in high-risk patients.

Neutrophilic HGF-MET Signalling Exacerbates Intestinal Inflammation.

BACKGROUND AND AIMS: Ulcerative colitis [UC] is associated with excessive neutrophil infiltration and collateral tissue damage, but the link is not yet completely understood. Since c-MET receptor tyrosine kinase [MET] is required for neutrophil chemoattraction and cytotoxicity in response to its ligand hepatocyte growth factor [HGF], we aimed to identify the function of HGF-MET signalling in neutrophils in UC patients and in mice during intestinal inflammation. METHODS: Serum and colonic biopsies from healthy controls and UC patients with active [Mayo endoscopic subscore 2-3] and inactive [Mayo endoscopic subscore 0-1] disease were collected to assess the level of serum and colonic HGF. Disease progression and immune cell infiltration were assessed during dextran sodium sulphate [DSS] colitis in wild-type and MRP8-Cre MET-LoxP mice. RESULTS: Increased mucosal HGF expression was detected in patients with active UC, and in mice during the inflammatory phase of DSS colitis. Similarly, serum HGF was significantly increased in active UC patients and positively correlated with C-reactive protein and blood neutrophil counts. Flow cytometric analysis also demonstrated an upregulation of colonic MET+ neutrophils during DSS colitis. Genetic ablation of MET in neutrophils reduced the severity of DSS-induced colitis. Concomitantly, there was a decreased number of TH17 cells, which could be due to a decreased production of IL-1ß by MET-deficient neutrophils. CONCLUSIONS: These data highlight the central role of neutrophilic HGF-MET signalling in exacerbating damage during intestinal inflammation. Hence, selective blockade of this pathway in neutrophils could be considered as a novel therapeutic approach in UC.

SwiftReg cluster registration automatically reduces flow cytometry data variability including batch effects.

Biological differences of interest in large, high-dimensional flow cytometry datasets are often obscured by undesired variations caused by differences in cytometers, reagents, or operators. Each variation type requires a different correction strategy, and their unknown contributions to overall variability hinder automated correction. We now describe swiftReg, an automated method that reduces undesired sources of variability between samples and particularly between batches. A high-resolution cluster map representing the multidimensional data is generated using the SWIFT algorithm, and shifts in cluster positions between samples are measured. Subpopulations are aligned between samples by displacing cell parameter values according to registration vectors derived from independent or locally-averaged cluster shifts. Batch variation is addressed by registering batch control or consensus samples, and applying the resulting shifts to individual samples. swiftReg selectively reduces batch variation, enhancing detection of biological differences. swiftReg outputs registered datasets as standard .FCS files to facilitate further analysis by other tools.

Mincle/Syk Signalling Promotes Intestinal Mucosal Inflammation Through Induction of Macrophage Pyroptosis in Crohn's Disease.

BACKGROUND: Macrophage-inducible C-type lectin [Mincle] signalling plays a proinflammatory role in different organs such as the brain and liver, but its role in intestinal inflammation, including Crohn's disease [CD], remains unknown. METHODS: The characteristics of Mincle signalling expression in CD patients and experimental colitis were examined. The functional role of Mincle signalling in the intestine was addressed in experimental colitis models in vivo by using Mincle knock-out [Mincle-/-] mice. In addition, neutralising anti-Mincle antibody, downstream spleen tyrosine kinase [Syk] inhibitor, and Mincle pharmacological agonist were used to study the Mincle signalling in intestine. Bone marrow-derived macrophages were collected from mice and used to further verify the effect of Mincle signalling in macrophages. RESULTS: This study has shown that Mincle signalling was significantly elevated in active human CD and experimental colitis, and macrophages were the principal leukocyte subset that upregulate Mincle signalling. Mincle deficiency and Syk pharmacological inhibition ameliorated the colitis by reducing induced macrophage pyroptosis, and activation of Mincle with the agonist aggravated the intestinal inflammation. The ex vivo studies demonstrated that activation of Mincle signalling promoted the release of proinflammatory cytokines, whereas its absence restricted release of proinflammatory cytokines from pyroptosis of macrophages. In addition, Mincle/Syk signalling in macrophages could promote the production of chemokines to recruit neutrophils by activating mitogen-activated protein kinase [MAPK] during intestinal inflammation. CONCLUSIONS: Mincle signalling promotes intestinal mucosal inflammation by inducing macrophage pyroptosis. Modulation of the Mincle/Syk axis emerges as a potential therapeutic strategy to target inflammation and treat CD.

Effects of the Positive Threshold and Data Analysis on Human MOG Antibody Detection by Live Flow Cytometry.

Human autoantibodies targeting myelin oligodendrocyte glycoprotein (MOG Ab) have become a useful clinical biomarker for the diagnosis of a spectrum of inflammatory demyelinating disorders. Live cell-based assays that detect MOG Ab against conformational MOG are currently the gold standard. Flow cytometry, in which serum binding to MOG-expressing cells and control cells are quantitively evaluated, is a widely used observer-independent, precise, and reliable detection method. However, there is currently no consensus on data analysis; for example, seropositive thresholds have been reported using varying standard deviations above a control cohort. Herein, we used a large cohort of 482 sera including samples from patients with monophasic or relapsing demyelination phenotypes consistent with MOG antibody-associated demyelination and other neurological diseases, as well as healthy controls, and applied a series of published analyses involving a background subtraction (delta) or a division (ratio). Loss of seropositivity and reduced detection sensitivity were observed when MOG ratio analyses or when 10 standard deviation (SD) or an arbitrary number was used to establish the threshold. Background binding and MOG ratio value were negatively correlated, in which patients seronegative by MOG ratio had high non-specific binding, a characteristic of serum that must be acknowledged. Most MOG Ab serostatuses were similar across analyses when optimal thresholds obtained by ROC analyses were used, demonstrating the robust nature and high discriminatory power of flow cytometry cell-based assays. With increased demand to identify MOG Ab-positive patients, a consensus on analysis is vital to improve patient diagnosis and for cross-study comparisons to ultimately define MOG Ab-associated disorders.

A Virtual Crossmatch-based Strategy Facilitates Sharing of Deceased Donor Kidneys for Highly Sensitized Recipients.

BACKGROUND: It is estimated that 19.2% of kidneys exported for candidates with >98% calculated panel reactive antibodies are transplanted into unintended recipients, most commonly due to positive physical crossmatch (PXM). We describe the application of a virtual crossmatch (VXM) that has resulted in a very low rate of transplantation into unintended recipients. METHODS: We performed a retrospective review of kidneys imported to our center to assess the reasons driving late reallocation based on the type of pretransplant crossmatch used for the intended recipient. RESULTS: From December 2014 to October 2017, 254 kidneys were imported based on our assessment of a VXM. Of these, 215 (84.6%) were transplanted without a pretransplant PXM. The remaining 39 (15.4%) recipients required a PXM on admission using a new sample because they did not have an HLA antibody test within the preceding 3 months or because they had a recent blood transfusion. A total of 93% of the imported kidneys were transplanted into intended recipients. There were 18 late reallocations: 9 (3.5%) due to identification of a new recipient medical problem upon admission, 5 (2%) due to suboptimal organ quality on arrival, and only 4 (1.6%) due to a positive PXM or HLA antibody concern. A total of 42% of the recipients of imported kidneys had a 100% calculated panel reactive antibodies. There were no hyperacute rejections and very infrequent acute rejection in the first year suggesting no evidence for immunologic memory response. CONCLUSIONS: Seamless sharing is within reach, even when kidneys are shipped long distances for highly sensitized recipients. Late reallocations can be almost entirely avoided with a strategy that relies heavily on VXM.

Inferring reaction network structure from single-cell, multiplex data, using toric systems theory.

The goal of many single-cell studies on eukaryotic cells is to gain insight into the biochemical reactions that control cell fate and state. In this paper we introduce the concept of Effective Stoichiometric Spaces (ESS) to guide the reconstruction of biochemical networks from multiplexed, fixed time-point, single-cell data. In contrast to methods based solely on statistical models of data, the ESS method leverages the power of the geometric theory of toric varieties to begin unraveling the structure of chemical reaction networks (CRN). This application of toric theory enables a data-driven mapping of covariance relationships in single-cell measurements into stoichiometric information, one in which each cell subpopulation has its associated ESS interpreted in terms of CRN theory. In the development of ESS we reframe certain aspects of the theory of CRN to better match data analysis. As an application of our approach we process cytomery- and image-based single-cell datasets and identify differences in cells treated with kinase inhibitors. Our approach is directly applicable to data acquired using readily accessible experimental methods such as Fluorescence Activated Cell Sorting (FACS) and multiplex immunofluorescence.

Nonspecific probe binding and automatic gating in flow cytometry and fluorescence activated cell sorting (FACS).

Flow cytometry is extensively used in cell biology to differentiate cells of interest (mutants) from control cells (wild-types). For mutant cells characterized by expression of a distinct membrane surface structure, fluorescent marker probes can be designed to bind specifically to these structures while the cells are in suspension, resulting in a sufficiently high fluorescence intensity measurement by the cytometer to identify a mutant cell. However, cell membranes may have relatively weak, nonspecific binding affinity to the probes, resulting in false positive results. Furthermore, the same effect would be present on mutant cells, allowing both specific and nonspecific binding to a single cell. We derive and analyze a kinetic model of fluorescent probe binding dynamics by tracking populations of mutant and wild-type cells with differing numbers of probes bound specifically and nonspecifically. By assuming the suspension is in chemical equilibrium prior to cytometry, we use a two-species Langmuir adsorption model to analyze the confounding effects of non-specific binding on the assay. Furthermore, we analytically derive an expectation maximization method to infer an appropriate estimate of the total number of mutant cells as an alternative to existing, heuristic methods. Lastly, using our model, we propose a new method to infer physical and experimental parameters from existing protocols. Our results provide improved ways to quantitatively analyze flow cytometry data.

Clinical and Molecular Features of Chronic Granulomatous Disease in Mainland China and a XL-CGD Female Infant Patient After Prenatal Diagnosis.

PURPOSE: Chronic granulomatous disease (CGD) is the most common phagocyte defect disease. Here, we describe 114 CGD patients in our center and report a rare female infant with XL-CGD to provide a better understanding of diagnosis, treatment, and prenatal diagnosis of CGD. METHOD: Patients were diagnosed by DHR-1,2,3 flow cytometry assays and gene analysis. X chromosome inactivation analysis and gp91phox protein test were used for a female infant with XL-CGD. RESULTS: XL-CGD accounts for the majority of cases in China and results in higher susceptibility to some infections than AR-CGD. The DHR assay can help diagnose CGD quickly, and atypical results should be combined with clinical manifestations, genetic analysis, and regular follow-up. For prenatal diagnosis, both gDNA and cDNA genotypes of amniotic fluid cells should be identified, and cord blood DHR assays should be performed to identify female XL-CGD patients.

diffcyt: Differential discovery in high-dimensional cytometry via high-resolution clustering.

High-dimensional flow and mass cytometry allow cell types and states to be characterized in great detail by measuring expression levels of more than 40 targeted protein markers per cell at the single-cell level. However, data analysis can be difficult, due to the large size and dimensionality of datasets as well as limitations of existing computational methods. Here, we present diffcyt, a new computational framework for differential discovery analyses in high-dimensional cytometry data, based on a combination of high-resolution clustering and empirical Bayes moderated tests adapted from transcriptomics. Our approach provides improved statistical performance, including for rare cell populations, along with flexible experimental designs and fast runtimes in an open-source framework.

Determination of essential phenotypic elements of clusters in high-dimensional entities-DEPECHE.

Technological advances have facilitated an exponential increase in the amount of information that can be derived from single cells, necessitating new computational tools that can make such highly complex data interpretable. Here, we introduce DEPECHE, a rapid, parameter free, sparse k-means-based algorithm for clustering of multi- and megavariate single-cell data. In a number of computational benchmarks aimed at evaluating the capacity to form biologically relevant clusters, including flow/mass-cytometry and single cell RNA sequencing data sets with manually curated gold standard solutions, DEPECHE clusters as well or better than the currently available best performing clustering algorithms. However, the main advantage of DEPECHE, compared to the state-of-the-art, is its unique ability to enhance interpretability of the formed clusters, in that it only retains variables relevant for cluster separation, thereby facilitating computational efficient analyses as well as understanding of complex datasets. DEPECHE is implemented in the open source R package DepecheR currently available at github.com/Theorell/DepecheR.

Computational approaches for high-throughput single-cell data analysis.

During the past decade, the number of novel technologies to interrogate biological systems at the single-cell level has skyrocketed. Numerous approaches for measuring the proteome, genome, transcriptome and epigenome at the single-cell level have been pioneered, using a variety of technologies. All these methods have one thing in common: they generate large and high-dimensional datasets that require advanced computational modelling tools to highlight and interpret interesting patterns in these data, potentially leading to novel biological insights and hypotheses. In this work, we provide an overview of the computational approaches used to interpret various types of single-cell data in an automated and unbiased way.