VoPo leverages cellular heterogeneity for predictive modeling of single-cell data

Solving the Puzzle of Preterm Birth

Solving the puzzle of preterm birth has been challenging and will require novel integrative solutions as preterm birth likely arises from many etiologies. It has been demonstrated that many sociodemographic and psychological determinants of preterm birth relate to its complex biology. It is this understanding that has enabled the development of a novel preventative strategy, which integrates the omics profile (genome, epigenome, transcriptome, proteome, metabolome, microbiome) with sociodemographic, environmental, and psychological determinants of individual pregnant people to solve the puzzle of preterm birth.

Morphological Signal Processing for Phenotype Recognition of Human Pluripotent Stem Cells Using Machine Learning Methods

Human pluripotent stem cells have the potential for unlimited proliferation and controlled differentiation into various somatic cells, making them a unique tool for regenerative and personalized medicine. Determining the best clone selection is a challenging problem in this field and requires new sensing instruments and methods able to automatically assess the state of a growing colony ('phenotype') and make decisions about its destiny. One possible solution for such label-free, non-invasive assessment is to make phase-contrast images and/or videos of growing stem cell colonies, process the morphological parameters ('morphological portrait', or signal), link this information to the colony phenotype, and initiate an automated protocol for the colony selection. As a step in implementing this strategy, we used machine learning methods to find an effective model for classifying the human pluripotent stem cell colonies of three lines according to their morphological phenotype ('good' or 'bad'), using morphological parameters from the previously published data as predictors. We found that the model using cellular morphological parameters as predictors and artificial neural networks as the classification method produced the best average accuracy of phenotype prediction (67%). When morphological parameters of colonies were used as predictors, logistic regression was the most effective classification method (75% average accuracy). Combining the morphological parameters of cells and colonies resulted in the most effective model, with a 99% average accuracy of phenotype prediction. Random forest was the most efficient classification method for the combined data. We applied feature selection methods and showed that different morphological parameters were important for phenotype recognition via either cellular or colonial parameters. Our results indicate a necessity for retaining both cellular and colonial morphological information for predicting the phenotype and provide an optimal choice for the machine learning method. The classification models reported in this study could be used as a basis for developing and/or improving automated solutions to control the quality of human pluripotent stem cells for medical purposes.

Cross-species comparative analysis of single presynapses

Comparing brain structure across species and regions enables key functional insights. Leveraging publicly available data from a novel mass cytometry-based method, synaptometry by time of flight (SynTOF), we applied an unsupervised machine learning approach to conduct a comparative study of presynapse molecular abundance across three species and three brain regions. We used neural networks and their attractive properties to model complex relationships among high dimensional data to develop a unified, unsupervised framework for comparing the profile of more than 4.5 million single presynapses among normal human, macaque, and mouse samples. An extensive validation showed the feasibility of performing cross-species comparison using SynTOF profiling. Integrative analysis of the abundance of 20 presynaptic proteins revealed near-complete separation between primates and mice involving synaptic pruning, cellular energy, lipid metabolism, and neurotransmission. In addition, our analysis revealed a strong overlap between the presynaptic composition of human and macaque in the cerebral cortex and neostriatum. Our unique approach illuminates species- and region-specific variation in presynapse molecular composition.

Multi-omic analysis of the tumor microenvironment shows clinical correlations in Ph1 study of atezolizumab +/- SoC in MM

Multiple myeloma (MM) remains incurable, and treatment of relapsed/refractory (R/R) disease is challenging. There is an unmet need for more targeted therapies in this setting; deep cellular and molecular phenotyping of the tumor and microenvironment in MM could help guide such therapies. This phase 1b study (NCT02431208) evaluated the safety and efficacy of the anti-programmed death-ligand 1 monoclonal antibody atezolizumab (Atezo) alone or in combination with the standard of care (SoC) treatments lenalidomide (Len) or pomalidomide (Pom) and/or daratumumab (Dara) in patients with R/R MM. Study endpoints included incidence of adverse events (AEs) and overall response rate (ORR). A novel unsupervised integrative multi-omic analysis was performed using RNA sequencing, mass cytometry immunophenotyping, and proteomic profiling of baseline and on-treatment bone marrow samples from patients receiving Atezo monotherapy or Atezo+Dara. A similarity network fusion (SNF) algorithm was applied to preprocessed data. Eighty-five patients were enrolled. Treatment-emergent deaths occurred in 2 patients; both deaths were considered unrelated to study treatment. ORRs ranged from 11.1% (Atezo+Len cohorts, n=18) to 83.3% (Atezo+Dara+Pom cohort, n=6). High-dimensional multi-omic profiling of the tumor microenvironment and integrative SNF analysis revealed novel correlations between cellular and molecular features of the tumor and immune microenvironment, patient selection criteria, and clinical outcome. Atezo monotherapy and SoC combinations were safe in this patient population and demonstrated some evidence of clinical efficacy. Integrative analysis of high dimensional genomics and immune data identified novel clinical correlations that may inform patient selection criteria and outcome assessment in future immunotherapy studies for myeloma.

In-silico generation of high-dimensional immune response data in patients using a deep neural network

Technologies for single-cell profiling of the immune system have enabled researchers to extract rich interconnected networks of cellular abundance, phenotypical and functional cellular parameters. These studies can power machine learning approaches to understand the role of the immune system in various diseases. However, the performance of these approaches and the generalizability of the findings have been hindered by limited cohort sizes in translational studies, partially due to logistical demands and costs associated with longitudinal data collection in sufficiently large patient cohorts. An evolving challenge is the requirement for ever-increasing cohort sizes as the dimensionality of datasets grows. We propose a deep learning model derived from a novel pipeline of optimal temporal cell matching and overcomplete autoencoders that uses data from a small subset of patients to learn to forecast an entire patient's immune response in a high dimensional space from one timepoint to another. In our analysis of 1.08 million cells from patients pre- and post-surgical intervention, we demonstrate that the generated patient-specific data are qualitatively and quantitatively similar to real patient data by demonstrating fidelity, diversity, and usefulness.

Maternal vaccination to prevent adverse pregnancy outcomes: An underutilized molecular immunological intervention?

Adverse pregnancy outcomes including maternal mortality, stillbirth, preterm birth, intrauterine growth restriction cause millions of deaths each year. More effective interventions are urgently needed. Maternal immunization could be one such intervention protecting the mother and newborn from infection through its pathogen-specific effects. However, many adverse pregnancy outcomes are not directly linked to the infectious pathogens targeted by existing maternal vaccines but rather are linked to pathological inflammation unfolding during pregnancy. The underlying pathogenesis driving such unfavourable outcomes have only partially been elucidated but appear to relate to altered immune regulation by innate as well as adaptive immune responses, ultimately leading to aberrant maternal immune activation. Maternal immunization, like all immunization, impacts the immune system beyond pathogen-specific immunity. This raises the possibility that maternal vaccination could potentially be utilised as a pathogen-agnostic immune modulatory intervention to redirect abnormal immune trajectories towards a more favourable phenotype providing pregnancy protection. In this review we describe the epidemiological evidence surrounding this hypothesis, along with the mechanistic plausibility and present a possible path forward to accelerate addressing the urgent need of adverse pregnancy outcomes.

Comparative analysis between high-grade serous ovarian cancer and healthy ovarian tissues using single-cell RNA sequencing

Introduction

High-grade serous ovarian cancer (HGSOC) is the most common histological subtype of ovarian cancer, and is associated with high mortality rates.

Methods

In this study, we analyzed specific cell subpopulations and compared different gene functions between healthy ovarian and ovarian cancer cells using single-cell RNA sequencing (ScRNA-seq). We delved deeper into the differences between healthy ovarian and ovarian cancer cells at different levels, and performed specific analysis on endothelial cells.

Results

We obtained scRNA-seq data of 6867 and 17056 cells from healthy ovarian samples and ovarian cancer samples, respectively. The transcriptional profiles of the groups differed at various stages of ovarian cell development. A detailed comparison of the cell cycle, and cell communication of different groups, revealed significant differences between healthy ovarian and ovarian cancer cells. We also found that apoptosis-related genes, URI1, PAK2, PARP1, CLU and TIMP3, were highly expressed, while immune-related genes, UBB, RPL11, CAV1, NUPR1 and Hsp90ab1, were lowly expressed in ovarian cancer cells. The results of the ScRNA-seq were verified using qPCR.

Discussion

Our findings revealed differences in function, gene expression and cell interaction patterns between ovarian cancer and healthy ovarian cell populations. These findings provide key insights on further research into the treatment of ovarian cancer.

Unbiased immune profiling reveals a natural killer cell-peripheral nerve axis in fibromyalgia

ABSTRACT

The pathophysiology of fibromyalgia syndrome (FMS) remains elusive, leading to a lack of objective diagnostic criteria and targeted treatment. We globally evaluated immune system changes in FMS by conducting multiparametric flow cytometry analyses of peripheral blood mononuclear cells and identified a natural killer (NK) cell decrease in patients with FMS. Circulating NK cells in FMS were exhausted yet activated, evidenced by lower surface expression of CD16, CD96, and CD226 and more CD107a and TIGIT. These NK cells were hyperresponsive, with increased CCL4 production and expression of CD107a when co-cultured with human leukocyte antigen null target cells. Genetic and transcriptomic pathway analyses identified significant enrichment of cell activation pathways in FMS driven by NK cells. Skin biopsies showed increased expression of NK activation ligand, unique long 16-binding protein, on subepidermal nerves of patients FMS and the presence of NK cells near peripheral nerves. Collectively, our results suggest that chronic activation and redistribution of circulating NK cells to the peripheral nerves contribute to the immunopathology associated with FMS.

Signature for Pain Recovery IN Teens (SPRINT): protocol for a multisite prospective signature study in chronic musculoskeletal pain

INTRODUCTION

Current treatments for chronic musculoskeletal (MSK) pain are suboptimal. Discovery of robust prognostic markers separating patients who recover from patients with persistent pain and disability is critical for developing patient-specific treatment strategies and conceiving novel approaches that benefit all patients. Given that chronic pain is a biopsychosocial process, this study aims to discover and validate a robust prognostic signature that measures across multiple dimensions in the same adolescent patient cohort with a computational analysis pipeline. This will facilitate risk stratification in adolescent patients with chronic MSK pain and more resourceful allocation of patients to costly and potentially burdensome multidisciplinary pain treatment approaches.

METHODS AND ANALYSIS

Here we describe a multi-institutional effort to collect, curate and analyse a high dimensional data set including epidemiological, psychometric, quantitative sensory, brain imaging and biological information collected over the course of 12 months. The aim of this effort is to derive a multivariate model with strong prognostic power regarding the clinical course of adolescent MSK pain and function.

ETHICS AND DISSEMINATION

The study complies with the National Institutes of Health policy on the use of a single internal review board (sIRB) for multisite research, with Cincinnati Children's Hospital Medical Center Review Board as the reviewing IRB. Stanford's IRB is a relying IRB within the sIRB. As foreign institutions, the University of Toronto and The Hospital for Sick Children (SickKids) are overseen by their respective ethics boards. All participants provide signed informed consent. We are committed to open-access publication, so that patients, clinicians and scientists have access to the study data and the signature(s) derived. After findings are published, we will upload a limited data set for sharing with other investigators on applicable repositories.

TRIAL REGISTRATION NUMBER

NCT04285112.

Integrated Single-cell and Plasma Proteomic Modeling to Predict Surgical Site Complications: A Prospective Cohort Study

OBJECTIVE

The aim of this study was to determine whether single-cell and plasma proteomic elements of the host's immune response to surgery accurately identify patients who develop a surgical site complication (SSC) after major abdominal surgery.

SUMMARY BACKGROUND DATA

SSCs may occur in up to 25% of patients undergoing bowel resection, resulting in significant morbidity and economic burden. However, the accurate prediction of SSCs remains clinically challenging. Leveraging high-content proteomic technologies to comprehensively profile patients' immune response to surgery is a promising approach to identify predictive biological factors of SSCs.

METHODS

Forty-one patients undergoing non-cancer bowel resection were prospectively enrolled. Blood samples collected before surgery and on postoperative day one (POD1) were analyzed using a combination of single-cell mass cytometry and plasma proteomics. The primary outcome was the occurrence of an SSC, including surgical site infection, anastomotic leak, or wound dehiscence within 30 days of surgery.

RESULTS

A multiomic model integrating the single-cell and plasma proteomic data collected on POD1 accurately differentiated patients with (n = 11) and without (n = 30) an SSC [area under the curve (AUC) = 0.86]. Model features included coregulated proinflammatory (eg, IL-6- and MyD88- signaling responses in myeloid cells) and immunosuppressive (eg, JAK/STAT signaling responses in M-MDSCs and Tregs) events preceding an SSC. Importantly, analysis of the immunological data obtained before surgery also yielded a model accurately predicting SSCs (AUC = 0.82).

CONCLUSIONS

The multiomic analysis of patients' immune response after surgery and immune state before surgery revealed systemic immune signatures preceding the development of SSCs. Our results suggest that integrating immunological data in perioperative risk assessment paradigms is a plausible strategy to guide individualized clinical care.

Single-synapse analyses of Alzheimer's disease implicate pathologic tau, DJ1, CD47, and ApoE

Synaptic molecular characterization is limited for Alzheimer’s disease (AD). Our newly invented mass cytometry–based method, synaptometry by time of flight (SynTOF), was used to measure 38 antibody probes in approximately 17 million single-synapse events from human brains without pathologic change or with pure AD or Lewy body disease (LBD), nonhuman primates (NHPs), and PS/APP mice. Synaptic molecular integrity in humans and NHP was similar. Although not detected in human synapses, Aβ was in PS/APP mice single-synapse events. Clustering and pattern identification of human synapses showed expected disease-specific differences, like increased hippocampal pathologic tau in AD and reduced caudate dopamine transporter in LBD, and revealed previously unidentified findings including increased hippocampal CD47 and lowered DJ1 in AD and higher ApoE in AD with dementia. Our results were independently supported by multiplex ion beam imaging of intact tissue. This highlights the higher depth and breadth of insight on neurodegenerative diseases obtainable through SynTOF.

Human influenza virus challenge identifies cellular correlates of protection for oral vaccination

Developing new influenza vaccines with improved performance and easier administration routes hinges on defining correlates of protection. Vaccine-elicited cellular correlates of protection for influenza in humans have not yet been demonstrated. A phase-2 double-blind randomized placebo and active (inactivated influenza vaccine) controlled study provides evidence that a human-adenovirus-5-based oral influenza vaccine tablet (VXA-A1.1) can protect from H1N1 virus challenge in humans. Mass cytometry characterization of vaccine-elicited cellular immune responses identified shared and vaccine-type-specific responses across B and T cells. For VXA-A1.1, the abundance of hemagglutinin-specific plasmablasts and plasmablasts positive for integrin α4β7, phosphorylated STAT5, or lacking expression of CD62L at day 8 were significantly correlated with protection from developing viral shedding following virus challenge at day 90 and contributed to an effective machine learning model of protection. These findings reveal the characteristics of vaccine-elicited cellular correlates of protection for an oral influenza vaccine.

Measuring the human immune response to surgery: multiomics for the prediction of postoperative outcomes

PURPOSE OF REVIEW

Postoperative complications including infections, cognitive impairment, and protracted recovery occur in one-third of the 300 million surgeries performed annually worldwide. Complications cause personal suffering along with a significant economic burden on our healthcare system. However, the accurate prediction of postoperative complications and patient-targeted interventions for their prevention remain as major clinical challenges.

RECENT FINDINGS

Although multifactorial in origin, the dysregulation of immunological mechanisms that occur in response to surgical trauma is a key determinant of postoperative complications. Prior research, primarily focusing on inflammatory plasma markers, has provided important clues regarding their pathogenesis. However, the recent advent of high-content, single-cell transcriptomic, and proteomic technologies has considerably improved our ability to characterize the immune response to surgery, thereby providing new means to understand the immunological basis of postoperative complications and to identify prognostic biological signatures.

SUMMARY

The comprehensive and single-cell characterization of the human immune response to surgery has significantly advanced our ability to predict the risk of postoperative complications. Multiomic modeling of patients' immune states holds promise for the discovery of preoperative predictive biomarkers, ultimately providing patients and surgeons with actionable information to improve surgical outcomes. Although recent studies have generated a wealth of knowledge, laying the foundation for a single-cell atlas of the human immune response to surgery, larger-scale multiomic studies are required to derive robust, scalable, and sufficiently powerful models to accurately predict the risk of postoperative complications in individual patients.

Analyzing high-dimensional cytometry data using FlowSOM

The dimensionality of cytometry data has strongly increased in the last decade, and in many situations the traditional manual downstream analysis becomes insufficient. The field is therefore slowly moving toward more automated approaches, and in this paper we describe the protocol for analyzing high-dimensional cytometry data using FlowSOM, a clustering and visualization algorithm based on a self-organizing map. FlowSOM is used to distinguish cell populations from cytometry data in an unsupervised way and can help to gain deeper insights in fields such as immunology and oncology. Since the original FlowSOM publication (2015), we have validated the tool on a wide variety of datasets, and to write this protocol, we made use of this experience to improve the user-friendliness of the package (e.g., comprehensive functions replacing commonly required scripts). Where the original paper focused mainly on the algorithm description, this protocol offers user guidelines on how to implement the procedure, detailed parameter descriptions and troubleshooting recommendations. The protocol provides clearly annotated R code, and is therefore relevant for all scientists interested in computational high-dimensional analyses without requiring a strong bioinformatics background. We demonstrate the complete workflow, starting from data preparation (such as compensation, transformation and quality control), including detailed discussion of the different FlowSOM parameters and visualization options, and concluding with how the results can be further used to answer biological questions, such as statistical comparison between groups of interest. An average FlowSOM analysis takes 1-3 h to complete, though quality issues can increase this time considerably.

CyAnno: A semi-automated approach for cell type annotation of mass cytometry datasets

MOTIVATION

For immune system monitoring in large-scale studies at the single-cell resolution using CyTOF, (semi-)automated computational methods are applied for annotating live cells of mixed cell types. Here, we show that the live cell pool can be highly enriched with undefined heterogeneous cells, i.e., 'ungated' cells, and that current semi-automated approaches ignore their modeling resulting in misclassified annotations.

RESULT

We introduce 'CyAnno', a novel semi-automated approach for deconvoluting the unlabeled cytometry dataset based on a machine learning framework utilizing manually gated training data that allows the integrative modeling of 'gated' cell types and the 'ungated' cells. By applying this framework on several CyTOF datasets, we demonstrated that including the 'ungated' cells can lead to a significant increase in the precision of the 'gated' cell types prediction. CyAnno can be used to identify even a single cell type, including rare cells, with higher efficacy than current state-of-the-art semi-automated approaches.

AVAILABILITY

The CyAnno is available as a python script with a user-manual and sample dataset at https://github.com/abbioinfo/CyAnno.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Censcyt: censored covariates in differential abundance analysis in cytometry

BACKGROUND

Innovations in single cell technologies have lead to a flurry of datasets and computational tools to process and interpret them, including analyses of cell composition changes and transition in cell states. The diffcyt workflow for differential discovery in cytometry data consist of several steps, including preprocessing, cell population identification and differential testing for an association with a binary or continuous covariate. However, the commonly measured quantity of survival time in clinical studies often results in a censored covariate where classical differential testing is inapplicable.

RESULTS

To overcome this limitation, multiple methods to directly include censored covariates in differential abundance analysis were examined with the use of simulation studies and a case study. Results show that multiple imputation based methods offer on-par performance with the Cox proportional hazards model in terms of sensitivity and error control, while offering flexibility to account for covariates. The tested methods are implemented in the R package censcyt as an extension of diffcyt and are available at https://bioconductor.org/packages/censcyt .

CONCLUSION

Methods for the direct inclusion of a censored variable as a predictor in GLMMs are a valid alternative to classical survival analysis methods, such as the Cox proportional hazard model, while allowing for more flexibility in the differential analysis.

Data-Driven Modeling of Pregnancy-Related Complications

A healthy pregnancy depends on complex interrelated biological adaptations involving placentation, maternal immune responses, and hormonal homeostasis. Recent advances in high-throughput technologies have provided access to multiomics biological data that, combined with clinical and social data, can provide a deeper understanding of normal and abnormal pregnancies. Integration of these heterogeneous datasets using state-of-the-art machine-learning methods can enable the prediction of short- and long-term health trajectories for a mother and offspring and the development of treatments to prevent or minimize complications. We review advanced machine-learning methods that could: provide deeper biological insights into a pregnancy not yet unveiled by current methodologies; clarify the etiologies and heterogeneity of pathologies that affect a pregnancy; and suggest the best approaches to address disparities in outcomes affecting vulnerable populations.

Single-cell peripheral immunoprofiling of Alzheimer's and Parkinson's diseases

Peripheral blood mononuclear cells (PBMCs) may provide insight into the pathogenesis of Alzheimer's disease (AD) or Parkinson's disease (PD). We investigated PBMC samples from 132 well-characterized research participants using seven canonical immune stimulants, mass cytometric identification of 35 PBMC subsets, and single-cell quantification of 15 intracellular signaling markers, followed by machine learning model development to increase predictive power. From these, three main intracellular signaling pathways were identified specifically in PBMC subsets from people with AD versus controls: reduced activation of PLCγ2 across many cell types and stimulations and selectively variable activation of STAT1 and STAT5, depending on stimulant and cell type. Our findings functionally buttress the now multiply-validated observation that a rare coding variant in PLCG2 is associated with a decreased risk of AD. Together, these data suggest enhanced PLCγ2 activity as a potential new therapeutic target for AD with a readily accessible pharmacodynamic biomarker.

Integration of mechanistic immunological knowledge into a machine learning pipeline improves predictions

The dense network of interconnected cellular signalling responses that are quantifiable in peripheral immune cells provides a wealth of actionable immunological insights. Although high-throughput single-cell profiling techniques, including polychromatic flow and mass cytometry, have matured to a point that enables detailed immune profiling of patients in numerous clinical settings, the limited cohort size and high dimensionality of data increase the possibility of false-positive discoveries and model overfitting. We introduce a generalizable machine learning platform, the immunological Elastic-Net (iEN), which incorporates immunological knowledge directly into the predictive models. Importantly, the algorithm maintains the exploratory nature of the high-dimensional dataset, allowing for the inclusion of immune features with strong predictive capabilities even if not consistent with prior knowledge. In three independent studies our method demonstrates improved predictions for clinically relevant outcomes from mass cytometry data generated from whole blood, as well as a large simulated dataset. The iEN is available under an open-source licence.

Preferential inhibition of adaptive immune system dynamics by glucocorticoids in patients after acute surgical trauma

Glucocorticoids (GC) are a controversial yet commonly used intervention in the clinical management of acute inflammatory conditions, including sepsis or traumatic injury. In the context of major trauma such as surgery, concerns have been raised regarding adverse effects from GC, thereby necessitating a better understanding of how GCs modulate the immune response. Here we report the results of a randomized controlled trial (NCT02542592) in which we employ a high-dimensional mass cytometry approach to characterize innate and adaptive cell signaling dynamics after a major surgery (primary outcome) in patients treated with placebo or methylprednisolone (MP). A robust, unsupervised bootstrap clustering of immune cell subsets coupled with random forest analysis shows profound (AUC = 0.92, p-value = 3.16E-8) MP-induced alterations of immune cell signaling trajectories, particularly in the adaptive compartments. By contrast, key innate signaling responses previously associated with pain and functional recovery after surgery, including STAT3 and CREB phosphorylation, are not affected by MP. These results imply cell-specific and pathway-specific effects of GCs, and also prompt future studies to examine GCs' effects on clinical outcomes likely dependent on functional adaptive immune responses.