Discovery of sparse, reliable omic biomarkers with Stabl.

Adoption of high-content omic technologies in clinical studies, coupled with computational methods, has yielded an abundance of candidate biomarkers. However, translating such findings into bona fide clinical biomarkers remains challenging. To facilitate this process, we introduce Stabl, a general machine learning method that identifies a sparse, reliable set of biomarkers by integrating noise injection and a data-driven signal-to-noise threshold into multivariable predictive modeling. Evaluation of Stabl on synthetic datasets and five independent clinical studies demonstrates improved biomarker sparsity and reliability compared to commonly used sparsity-promoting regularization methods while maintaining predictive performance; it distills datasets containing 1,400-35,000 features down to 4-34 candidate biomarkers. Stabl extends to multi-omic integration tasks, enabling biological interpretation of complex predictive models, as it hones in on a shortlist of proteomic, metabolomic and cytometric events predicting labor onset, microbial biomarkers of pre-term birth and a pre-operative immune signature of post-surgical infections. Stabl is available at https://github.com/gregbellan/Stabl .

Spatial subsetting enables integrative modeling of oral squamous cell carcinoma multiplex imaging data.

Oral squamous cell carcinoma (OSCC), a prevalent and aggressive neoplasm, poses a significant challenge due to poor prognosis and limited prognostic biomarkers. Leveraging highly multiplexed imaging mass cytometry, we investigated the tumor immune microenvironment (TIME) in OSCC biopsies, characterizing immune cell distribution and signaling activity at the tumor-invasive front. Our spatial subsetting approach standardized cellular populations by tissue zone, improving feature reproducibility and revealing TIME patterns accompanying loss-of-differentiation. Employing a machine-learning pipeline combining reliable feature selection with multivariable modeling, we achieved accurate histological grade classification (AUC = 0.88). Three model features correlated with clinical outcomes in an independent cohort: granulocyte MAPKAPK2 signaling at the tumor front, stromal CD4+ memory T cell size, and the distance of fibroblasts from the tumor border. This study establishes a robust modeling framework for distilling complex imaging data, uncovering sentinel characteristics of the OSCC TIME to facilitate prognostic biomarkers discovery for recurrence risk stratification and immunomodulatory therapy development.

Longitudinal clinical phenotyping of post COVID condition in Mexican adults recovering from severe COVID-19: a prospective cohort study.

Introduction: Few studies have evaluated the presence of Post COVID-19 conditions (PCC) in people from Latin America, a region that has been heavily afflicted by the COVID-19 pandemic. In this study, we describe the frequency, co-occurrence, predictors, and duration of 23 symptoms in a cohort of Mexican patients with PCC. Methods: We prospectively enrolled and followed adult patients hospitalized for severe COVID-19 at a tertiary care centre in Mexico City. The incidence of PCC symptoms was determined using questionnaires. Unsupervised clustering of PCC symptom co-occurrence and Kaplan-Meier analyses of symptom persistence were performed. The effect of baseline clinical characteristics was evaluated using Cox regression models and reported with hazard ratios (HR). Results: We found that amongst 192 patients with PCC, respiratory problems were the most prevalent and commonly co-occurred with functional activity impairment. 56% had ≥5 persistent symptoms. Symptom persistence probability at 360 days 0.78. Prior SARS-CoV-2 vaccination and infection during the Delta variant wave were associated with a shorter duration of PCC. Male sex was associated with a shorter duration of functional activity impairment and respiratory symptoms. Hypertension and diabetes were associated with a longer duration of functional impairment. Previous vaccination accelerated PCC recovery. Discussion: In our cohort, PCC symptoms were frequent (particularly respiratory and neurocognitive ones) and persistent. Importantly, prior SARS-CoV-2 vaccination resulted in a shorter duration of PCC.

Expanded vacuum-stable gels for multiplexed high-resolution spatial histopathology.

Cellular organization and functions encompass multiple scales in vivo. Emerging high-plex imaging technologies are limited in resolving subcellular biomolecular features. Expansion Microscopy (ExM) and related techniques physically expand samples for enhanced spatial resolution, but are challenging to be combined with high-plex imaging technologies to enable integrative multiscaled tissue biology insights. Here, we introduce Expand and comPRESS hydrOgels (ExPRESSO), an ExM framework that allows high-plex protein staining, physical expansion, and removal of water, while retaining the lateral tissue expansion. We demonstrate ExPRESSO imaging of archival clinical tissue samples on Multiplexed Ion Beam Imaging and Imaging Mass Cytometry platforms, with detection capabilities of > 40 markers. Application of ExPRESSO on archival human lymphoid and brain tissues resolved tissue architecture at the subcellular level, particularly that of the blood-brain barrier. ExPRESSO hence provides a platform for extending the analysis compatibility of hydrogel-expanded biospecimens to mass spectrometry, with minimal modifications to protocols and instrumentation.

Stabl: sparse and reliable biomarker discovery in predictive modeling of high-dimensional omic data.

High-content omic technologies coupled with sparsity-promoting regularization methods (SRM) have transformed the biomarker discovery process. However, the translation of computational results into a clinical use-case scenario remains challenging. A rate-limiting step is the rigorous selection of reliable biomarker candidates among a host of biological features included in multivariate models. We propose Stabl, a machine learning framework that unifies the biomarker discovery process with multivariate predictive modeling of clinical outcomes by selecting a sparse and reliable set of biomarkers. Evaluation of Stabl on synthetic datasets and four independent clinical studies demonstrates improved biomarker sparsity and reliability compared to commonly used SRMs at similar predictive performance. Stabl readily extends to double- and triple-omics integration tasks and identifies a sparser and more reliable set of biomarkers than those selected by state-of-the-art early- and late-fusion SRMs, thereby facilitating the biological interpretation and clinical translation of complex multi-omic predictive models. The complete package for Stabl is available online at https://github.com/gregbellan/Stabl.

Differential dynamics of bone graft transplantation and mesenchymal stem cell therapy during bone defect healing in a murine critical size defect.

Background: A critical size bone defect is a clinical scenario in which bone is lost or excised due to trauma, infection, tumor, or other causes, and cannot completely heal spontaneously. The most common treatment for this condition is autologous bone grafting to the defect site. However, autologous bone graft is often insufficient in quantity or quality for transplantation to these large defects. Recently, tissue engineering methods using mesenchymal stem cells (MSCs) have been proposed as an alternative treatment. However, the underlying biological principles and optimal techniques for tissue regeneration of bone using stem cell therapy have not been completely elucidated. Methods: In this study, we compare the early cellular dynamics of healing between bone graft transplantation and MSC therapy in a murine chronic femoral critical-size bone defect. We employ high-dimensional mass cytometry to provide a comprehensive view of the differences in cell composition, stem cell functionality, and immunomodulatory activity between these two treatment methods one week after transplantation. Results: We reveal distinct cell compositions among tissues from bone defect sites compared with original bone graft, show active recruitment of MSCs to the bone defect sites, and demonstrate the phenotypic diversity of macrophages and T cells in each group that may affect the clinical outcome. Conclusion: Our results provide critical data and future directions on the use of MSCs for treating critical size defects to regenerate bone.Translational Potential of this article: This study showed systematic comparisons of the cellular and immunomodulatory profiles among different interventions to improve the healing of the critical-size bone defect. The results provided potential strategies for designing robust therapeutic interventions for the unmet clinical need of treating critical-size bone defects.

Integrated plasma proteomic and single-cell immune signaling network signatures demarcate mild, moderate, and severe COVID-19.

The biological determinants underlying the range of coronavirus 2019 (COVID-19) clinical manifestations are not fully understood. Here, over 1,400 plasma proteins and 2,600 single-cell immune features comprising cell phenotype, endogenous signaling activity, and signaling responses to inflammatory ligands are cross-sectionally assessed in peripheral blood from 97 patients with mild, moderate, and severe COVID-19 and 40 uninfected patients. Using an integrated computational approach to analyze the combined plasma and single-cell proteomic data, we identify and independently validate a multi-variate model classifying COVID-19 severity (multi-class area under the curve [AUC]training = 0.799, p = 4.2e-6; multi-class AUCvalidation = 0.773, p = 7.7e-6). Examination of informative model features reveals biological signatures of COVID-19 severity, including the dysregulation of JAK/STAT, MAPK/mTOR, and nuclear factor κB (NF-κB) immune signaling networks in addition to recapitulating known hallmarks of COVID-19. These results provide a set of early determinants of COVID-19 severity that may point to therapeutic targets for prevention and/or treatment of COVID-19 progression.

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.

A Peripheral Immune Signature of Labor Induction.

Approximately 1 in 4 pregnant women in the United States undergo labor induction. The onset and establishment of labor, particularly induced labor, is a complex and dynamic process influenced by multiple endocrine, inflammatory, and mechanical factors as well as obstetric and pharmacological interventions. The duration from labor induction to the onset of active labor remains unpredictable. Moreover, prolonged labor is associated with severe complications for the mother and her offspring, most importantly chorioamnionitis, uterine atony, and postpartum hemorrhage. While maternal immune system adaptations that are critical for the maintenance of a healthy pregnancy have been previously characterized, the role of the immune system during the establishment of labor is poorly understood. Understanding maternal immune adaptations during labor initiation can have important ramifications for predicting successful labor induction and labor complications in both induced and spontaneous types of labor. The aim of this study was to characterize labor-associated maternal immune system dynamics from labor induction to the start of active labor. Serial blood samples from fifteen participants were collected immediately prior to labor induction (baseline) and during the latent phase until the start of active labor. Using high-dimensional mass cytometry, a total of 1,059 single-cell immune features were extracted from each sample. A multivariate machine-learning method was employed to characterize the dynamic changes of the maternal immune system after labor induction until the establishment of active labor. A cross-validated linear sparse regression model (least absolute shrinkage and selection operator, LASSO) predicted the minutes since induction of labor with high accuracy (R = 0.86, p = 6.7e-15, RMSE = 277 min). Immune features most informative for the model included STAT5 signaling in central memory CD8+ T cells and pro-inflammatory STAT3 signaling responses across multiple adaptive and innate immune cell subsets. Our study reports a peripheral immune signature of labor induction, and provides important insights into biological mechanisms that may ultimately predict labor induction success as well as complications, thereby facilitating clinical decision-making to improve maternal and fetal well-being.

Single-Cell Analysis of the Neonatal Immune System Across the Gestational Age Continuum.

Although most causes of death and morbidity in premature infants are related to immune maladaptation, the premature immune system remains poorly understood. We provide a comprehensive single-cell depiction of the neonatal immune system at birth across the spectrum of viable gestational age (GA), ranging from 25 weeks to term. A mass cytometry immunoassay interrogated all major immune cell subsets, including signaling activity and responsiveness to stimulation. An elastic net model described the relationship between GA and immunome (R=0.85, p=8.75e-14), and unsupervised clustering highlighted previously unrecognized GA-dependent immune dynamics, including decreasing basal MAP-kinase/NFκB signaling in antigen presenting cells; increasing responsiveness of cytotoxic lymphocytes to interferon-α; and decreasing frequency of regulatory and invariant T cells, including NKT-like cells and CD8+CD161+ T cells. Knowledge gained from the analysis of the neonatal immune landscape across GA provides a mechanistic framework to understand the unique susceptibility of preterm infants to both hyper-inflammatory diseases and infections.

Integrated trajectories of the maternal metabolome, proteome, and immunome predict labor onset.

Estimating the time of delivery is of high clinical importance because pre- and postterm deviations are associated with complications for the mother and her offspring. However, current estimations are inaccurate. As pregnancy progresses toward labor, major transitions occur in fetomaternal immune, metabolic, and endocrine systems that culminate in birth. The comprehensive characterization of maternal biology that precedes labor is key to understanding these physiological transitions and identifying predictive biomarkers of delivery. Here, a longitudinal study was conducted in 63 women who went into labor spontaneously. More than 7000 plasma analytes and peripheral immune cell responses were analyzed using untargeted mass spectrometry, aptamer-based proteomic technology, and single-cell mass cytometry in serial blood samples collected during the last 100 days of pregnancy. The high-dimensional dataset was integrated into a multiomic model that predicted the time to spontaneous labor [R = 0.85, 95% confidence interval (CI) [0.79 to 0.89], P = 1.2 × 10-40, N = 53, training set; R = 0.81, 95% CI [0.61 to 0.91], P = 3.9 × 10-7, N = 10, independent test set]. Coordinated alterations in maternal metabolome, proteome, and immunome marked a molecular shift from pregnancy maintenance to prelabor biology 2 to 4 weeks before delivery. A surge in steroid hormone metabolites and interleukin-1 receptor type 4 that preceded labor coincided with a switch from immune activation to regulation of inflammatory responses. Our study lays the groundwork for developing blood-based methods for predicting the day of labor, anchored in mechanisms shared in preterm and term pregnancies.

Integrated plasma proteomic and single-cell immune signaling network signatures demarcate mild, moderate, and severe COVID-19.

The biological determinants of the wide spectrum of COVID-19 clinical manifestations are not fully understood. Here, over 1400 plasma proteins and 2600 single-cell immune features comprising cell phenotype, basal signaling activity, and signaling responses to inflammatory ligands were assessed in peripheral blood from patients with mild, moderate, and severe COVID-19, at the time of diagnosis. Using an integrated computational approach to analyze the combined plasma and single-cell proteomic data, we identified and independently validated a multivariate model classifying COVID-19 severity (multi-class AUCtraining = 0.799, p-value = 4.2e-6; multi-class AUCvalidation = 0.773, p-value = 7.7e-6). Features of this high-dimensional model recapitulated recent COVID-19 related observations of immune perturbations, and revealed novel biological signatures of severity, including the mobilization of elements of the renin-angiotensin system and primary hemostasis, as well as dysregulation of JAK/STAT, MAPK/mTOR, and NF-κB immune signaling networks. These results provide a set of early determinants of COVID-19 severity that may point to therapeutic targets for the prevention of COVID-19 progression.

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.

Author Correction: Preferential inhibition of adaptive immune system dynamics by glucocorticoids in patients after acute surgical trauma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Systematic Immunophenotyping Reveals Sex-Specific Responses After Painful Injury in Mice.

Many diseases display unequal prevalence between sexes. The sex-specific immune response to both injury and persistent pain remains underexplored and would inform treatment paradigms. We utilized high-dimensional mass cytometry to perform a comprehensive analysis of phenotypic and functional immune system differences between male and female mice after orthopedic injury. Multivariate modeling of innate and adaptive immune cell responses after injury using an elastic net algorithm, a regularized regression method, revealed sex-specific divergence at 12 h and 7 days after injury with a stronger immune response to injury in females. At 12 h, females upregulated STAT3 signaling in neutrophils but downregulated STAT1 and STAT6 signals in T regulatory cells, suggesting a lack of engagement of immune suppression pathways by females. Furthermore, at 7 days females upregulated MAPK pathways (p38, ERK, NFkB) in CD4T memory cells, setting up a possible heightened immune memory of painful injury. Taken together, our findings provide the first comprehensive and functional analysis of sex-differences in the immune response to painful injury.

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.

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

High-throughput single-cell analysis technologies produce an abundance of data that is critical for profiling the heterogeneity of cellular systems. We introduce VoPo (https://github.com/stanleyn/VoPo), a machine learning algorithm for predictive modeling and comprehensive visualization of the heterogeneity captured in large single-cell datasets. In three mass cytometry datasets, with the largest measuring hundreds of millions of cells over hundreds of samples, VoPo defines phenotypically and functionally homogeneous cell populations. VoPo further outperforms state-of-the-art machine learning algorithms in classification tasks, and identified immune-correlates of clinically-relevant parameters.

Differential Dynamics of the Maternal Immune System in Healthy Pregnancy and Preeclampsia.

Preeclampsia is one of the most severe pregnancy complications and a leading cause of maternal death. However, early diagnosis of preeclampsia remains a clinical challenge. Alterations in the normal immune adaptations necessary for the maintenance of a healthy pregnancy are central features of preeclampsia. However, prior analyses primarily focused on the static assessment of select immune cell subsets have provided limited information for the prediction of preeclampsia. Here, we used a high-dimensional mass cytometry immunoassay to characterize the dynamic changes of over 370 immune cell features (including cell distribution and functional responses) in maternal blood during healthy and preeclamptic pregnancies. We found a set of eight cell-specific immune features that accurately identified patients well before the clinical diagnosis of preeclampsia (median area under the curve (AUC) 0.91, interquartile range [0.82-0.92]). Several features recapitulated previously known immune dysfunctions in preeclampsia, such as elevated pro-inflammatory innate immune responses early in pregnancy and impaired regulatory T (Treg) cell signaling. The analysis revealed additional novel immune responses that were strongly associated with, and preceded the onset of preeclampsia, notably abnormal STAT5ab signaling dynamics in CD4+T cell subsets (AUC 0.92, p = 8.0E-5). These results provide a global readout of the dynamics of the maternal immune system early in pregnancy and lay the groundwork for identifying clinically-relevant immune dysfunctions for the prediction and prevention of preeclampsia.

A year-long immune profile of the systemic response in acute stroke survivors.

Stroke is a leading cause of cognitive impairment and dementia, but the mechanisms that underlie post-stroke cognitive decline are not well understood. Stroke produces profound local and systemic immune responses that engage all major innate and adaptive immune compartments. However, whether the systemic immune response to stroke contributes to long-term disability remains ill-defined. We used a single-cell mass cytometry approach to comprehensively and functionally characterize the systemic immune response to stroke in longitudinal blood samples from 24 patients over the course of 1 year and correlated the immune response with changes in cognitive functioning between 90 and 365 days post-stroke. Using elastic net regularized regression modelling, we identified key elements of a robust and prolonged systemic immune response to ischaemic stroke that occurs in three phases: an acute phase (Day 2) characterized by increased signal transducer and activator of transcription 3 (STAT3) signalling responses in innate immune cell types, an intermediate phase (Day 5) characterized by increased cAMP response element-binding protein (CREB) signalling responses in adaptive immune cell types, and a late phase (Day 90) by persistent elevation of neutrophils, and immunoglobulin M+ (IgM+) B cells. By Day 365 there was no detectable difference between these samples and those from an age- and gender-matched patient cohort without stroke. When regressed against the change in the Montreal Cognitive Assessment scores between Days 90 and 365 after stroke, the acute inflammatory phase Elastic Net model correlated with post-stroke cognitive trajectories (r = -0.692, Bonferroni-corrected P = 0.039). The results demonstrate the utility of a deep immune profiling approach with mass cytometry for the identification of clinically relevant immune correlates of long-term cognitive trajectories.

Multiomics modeling of the immunome, transcriptome, microbiome, proteome and metabolome adaptations during human pregnancy.

Motivation: Multiple biological clocks govern a healthy pregnancy. These biological mechanisms produce immunologic, metabolomic, proteomic, genomic and microbiomic adaptations during the course of pregnancy. Modeling the chronology of these adaptations during full-term pregnancy provides the frameworks for future studies examining deviations implicated in pregnancy-related pathologies including preterm birth and preeclampsia. Results: We performed a multiomics analysis of 51 samples from 17 pregnant women, delivering at term. The datasets included measurements from the immunome, transcriptome, microbiome, proteome and metabolome of samples obtained simultaneously from the same patients. Multivariate predictive modeling using the Elastic Net (EN) algorithm was used to measure the ability of each dataset to predict gestational age. Using stacked generalization, these datasets were combined into a single model. This model not only significantly increased predictive power by combining all datasets, but also revealed novel interactions between different biological modalities. Future work includes expansion of the cohort to preterm-enriched populations and in vivo analysis of immune-modulating interventions based on the mechanisms identified. Availability and implementation: Datasets and scripts for reproduction of results are available through: https://nalab.stanford.edu/multiomics-pregnancy/. Supplementary information: Supplementary data are available at Bioinformatics online.