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.

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.

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.

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.

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.

An immune clock of human pregnancy.

The maintenance of pregnancy relies on finely tuned immune adaptations. We demonstrate that these adaptations are precisely timed, reflecting an immune clock of pregnancy in women delivering at term. Using mass cytometry, the abundance and functional responses of all major immune cell subsets were quantified in serial blood samples collected throughout pregnancy. Cell signaling-based Elastic Net, a regularized regression method adapted from the elastic net algorithm, was developed to infer and prospectively validate a predictive model of interrelated immune events that accurately captures the chronology of pregnancy. Model components highlighted existing knowledge and revealed previously unreported biology, including a critical role for the interleukin-2-dependent STAT5ab signaling pathway in modulating T cell function during pregnancy. These findings unravel the precise timing of immunological events occurring during a term pregnancy and provide the analytical framework to identify immunological deviations implicated in pregnancy-related pathologies.

Deep Immune Profiling of an Arginine-Enriched Nutritional Intervention in Patients Undergoing Surgery.

Application of high-content immune profiling technologies has enormous potential to advance medicine. Whether these technologies reveal pertinent biology when implemented in interventional clinical trials is an important question. The beneficial effects of preoperative arginine-enriched dietary supplements (AES) are highly context specific, as they reduce infection rates in elective surgery, but possibly increase morbidity in critically ill patients. This study combined single-cell mass cytometry with the multiplex analysis of relevant plasma cytokines to comprehensively profile the immune-modifying effects of this much-debated intervention in patients undergoing surgery. An elastic net algorithm applied to the high-dimensional mass cytometry dataset identified a cross-validated model consisting of 20 interrelated immune features that separated patients assigned to AES from controls. The model revealed wide-ranging effects of AES on innate and adaptive immune compartments. Notably, AES increased STAT1 and STAT3 signaling responses in lymphoid cell subsets after surgery, consistent with enhanced adaptive mechanisms that may protect against postsurgical infection. Unexpectedly, AES also increased ERK and P38 MAPK signaling responses in monocytic myeloid-derived suppressor cells, which was paired with their pronounced expansion. These results provide novel mechanistic arguments as to why AES may exert context-specific beneficial or adverse effects in patients with critical illness. This study lays out an analytical framework to distill high-dimensional datasets gathered in an interventional clinical trial into a fairly simple model that converges with known biology and provides insight into novel and clinically relevant cellular mechanisms.

Multicenter Systems Analysis of Human Blood Reveals Immature Neutrophils in Males and During Pregnancy.

Despite clear differences in immune system responses and in the prevalence of autoimmune diseases between males and females, there is little understanding of the processes involved. In this study, we identified a gene signature of immature-like neutrophils, characterized by the overexpression of genes encoding for several granule-containing proteins, which was found at higher levels (up to 3-fold) in young (20-30 y old) but not older (60 to >89 y old) males compared with females. Functional and phenotypic characterization of peripheral blood neutrophils revealed more mature and responsive neutrophils in young females, which also exhibited an elevated capacity in neutrophil extracellular trap formation at baseline and upon microbial or sterile autoimmune stimuli. The expression levels of the immature-like neutrophil signature increased linearly with pregnancy, an immune state of increased susceptibility to certain infections. Using mass cytometry, we also find increased frequencies of immature forms of neutrophils in the blood of women during late pregnancy. Thus, our findings show novel sex differences in innate immunity and identify a common neutrophil signature in males and in pregnant women.

Expression of specific inflammasome gene modules stratifies older individuals into two extreme clinical and immunological states.

Low-grade, chronic inflammation has been associated with many diseases of aging, but the mechanisms responsible for producing this inflammation remain unclear. Inflammasomes can drive chronic inflammation in the context of an infectious disease or cellular stress, and they trigger the maturation of interleukin-1ß (IL-1ß). Here we find that the expression of specific inflammasome gene modules stratifies older individuals into two extremes: those with constitutive expression of IL-1ß, nucleotide metabolism dysfunction, elevated oxidative stress, high rates of hypertension and arterial stiffness; and those without constitutive expression of IL-1ß, who lack these characteristics. Adenine and N4-acetylcytidine, nucleotide-derived metabolites that are detectable in the blood of the former group, prime and activate the NLRC4 inflammasome, induce the production of IL-1ß, activate platelets and neutrophils and elevate blood pressure in mice. In individuals over 85 years of age, the elevated expression of inflammasome gene modules was associated with all-cause mortality. Thus, targeting inflammasome components may ameliorate chronic inflammation and various other age-associated conditions.

Mapping the Fetomaternal Peripheral Immune System at Term Pregnancy.

Preterm labor and infections are the leading causes of neonatal deaths worldwide. During pregnancy, immunological cross talk between the mother and her fetus is critical for the maintenance of pregnancy and the delivery of an immunocompetent neonate. A precise understanding of healthy fetomaternal immunity is the important first step to identifying dysregulated immune mechanisms driving adverse maternal or neonatal outcomes. This study combined single-cell mass cytometry of paired peripheral and umbilical cord blood samples from mothers and their neonates with a graphical approach developed for the visualization of high-dimensional data to provide a high-resolution reference map of the cellular composition and functional organization of the healthy fetal and maternal immune systems at birth. The approach enabled mapping of known phenotypical and functional characteristics of fetal immunity (including the functional hyperresponsiveness of CD4+ and CD8+ T cells and the global blunting of innate immune responses). It also allowed discovery of new properties that distinguish the fetal and maternal immune systems. For example, examination of paired samples revealed differences in endogenous signaling tone that are unique to a mother and her offspring, including increased ERK1/2, MAPK-activated protein kinase 2, rpS6, and CREB phosphorylation in fetal Tbet+CD4+ T cells, CD8+ T cells, B cells, and CD56loCD16+ NK cells and decreased ERK1/2, MAPK-activated protein kinase 2, and STAT1 phosphorylation in fetal intermediate and nonclassical monocytes. This highly interactive functional map of healthy fetomaternal immunity builds the core reference for a growing data repository that will allow inferring deviations from normal associated with adverse maternal and neonatal outcomes.