Pediatric leukemia is one of the most intensively studied cancers in modern medicine, yet the immune environment surrounding leukemic cells remains remarkably difficult to characterize in a way that is both reproducible and clinically meaningful. A new study published in BMC Bioinformatics by Seung-Hwan Kim of the Department of Biology at Fisher College takes a fresh look at this challenge by asking a deceptively simple question: does the continuous composition of immune cells in a child’s leukemia sample carry longitudinal information that coarse, cluster-based ecological labels cannot capture? The answer, according to a sophisticated Bayesian analysis of single-cell data, is a qualified but important yes.
The research, published on 16 October 2026, draws on publicly available data from the Single-cell Pediatric Cancer Atlas, or ScPCA, a resource that catalogs single-cell profiles from pediatric cancer patients. Kim reprocessed two public leukemia projects from the atlas using a common cell-type mapping strategy, ensuring that cells from different studies could be compared on the same biological vocabulary. From this harmonized foundation, a cross-project immune composition was constructed from three major compartments: T and natural killer cells, B-lineage cells, and myeloid cells. The goal was to see whether these proportions, tracked over time within individual patients, could reveal something about the trajectory of the disease that simpler categorical labels would miss.
Central to the study is the concept of ecological context. In cancer biology, the tumor microenvironment is often described in terms of discrete states, such as immune-hot versus immune-cold, or dominant cell type versus suppressed lineage. These coarse labels are convenient, but they risk discarding information. Kim’s framework treats immune composition as a continuous variable embedded in a stochastic longitudinal model, specifically an Ornstein-Uhlenbeck process, which is a mathematical model that describes a value fluctuating around a long-term mean with a tendency to revert toward that mean. In this context, the OU process models the leukemia-state trajectory over time, and the question is whether immune composition shifts the attractor, the value mu toward which the trajectory is pulled.
Among 104 compositionally eligible participants, coarse clustering with two groups was reproducible under bootstrap resampling, meaning that the split into two clusters was stable when the data were resampled repeatedly. However, this clustering showed only limited association with the project from which the data came or with the specific diagnosis. Within the 99 participants with B-cell acute lymphoblastic leukemia, the most common pediatric leukemia, the two clusters separated 73 participants with T and natural killer cell dominant profiles from 26 participants with B and myeloid shifted profiles. This separation was largely independent of which project the data originated from, suggesting a genuine biological structure rather than a batch effect.
Yet a critical complication emerged. When Kim compared cluster assignments derived from raw composition with those derived from centered log-ratio transformed composition, a standard transformation for compositional data that accounts for the constraint that proportions must sum to one, the concordance was poor. This means that the apparent ecological state a patient belongs to depends heavily on how the compositional data are represented mathematically. A patient might appear in one cluster under raw proportions and in a different cluster under log-ratio transformed proportions. This representation sensitivity undermines the reliability of coarse context labels as stable descriptors of biological reality.
Because of this instability, Kim chose to evaluate continuous ecological composition separately from coarse context labels in a corrected longitudinal analysis. The longitudinal dataset included 34 participants contributing 241 observations and 207 transitions between time points. The key finding was that higher B-lineage composition was associated with a lower effective OU attractor, the parameter mu, in both models restricted to B-lineage cells alone and models that included B plus myeloid compartments. In biological terms, this suggests that as B-lineage cells become more prevalent in the immune environment, the leukemia state trajectory is pulled toward a lower value, potentially reflecting a less aggressive or more controlled disease state.
In contrast, the coarse context effects, the discrete labels that would traditionally be used to categorize patients, were unresolved in the longitudinal models. They did not provide clear, consistent information about the direction or magnitude of the leukemia trajectory. This is a striking result because it suggests that the very labels clinicians and researchers often rely on to summarize immune context may be less informative than the underlying continuous composition they are meant to represent.
To assess whether the continuous composition genuinely adds predictive value, Kim employed leave-one-participant-out validation, a rigorous approach in which each participant is held out in turn and the model is trained on the remaining participants. In aggregate, the model including B-lineage composition modestly outperformed the OU-only model that lacked ecological covariates. However, the predictive differences were heterogeneous across participants, meaning that for some individuals the continuous composition helped substantially while for others it did not. Models that included only coarse context labels, or that combined ecology with context, provided no incremental advantage over the simpler approaches. This pattern reinforces the conclusion that the value lies in the continuous signal, not in the categorical summaries.
A secondary analysis using an ordinal-state model, which treats the leukemia state as an ordered categorical variable rather than a continuous one, revealed a directionally consistent B-associated signal. This consistency across different modeling frameworks strengthens confidence that the B-lineage association is not an artifact of a single analytical choice. Simulations using the observed irregular sampling schedules from the real data demonstrated that moderate ecological effects could be recovered even when observations were not evenly spaced in time, addressing a common concern in longitudinal single-cell studies. Posterior-predictive analyses further showed that including continuous ecological composition improved the calibration of uncertainty estimates, meaning that the model’s confidence intervals better reflected actual variability in the data.
The broader implications of this work extend beyond pediatric leukemia. The study demonstrates a general principle: stable clustering of compositional data does not necessarily imply that those clusters represent longitudinally predictive ecological states. Researchers across cancer biology, immunology, and microbial ecology often face the temptation to reduce continuous compositional information into discrete categories for ease of interpretation. Kim’s framework offers a rigorous way to test whether such reduction loses information, and it provides a template for distinguishing descriptive ecological coarse-graining from ecological covariates that genuinely add explanatory power to stochastic longitudinal models. Importantly, the author cautions that the OU mean reversion observed in the models should not be interpreted as direct evidence of stabilizing selection or of fixed leukemia ecotypes, a nuance that guards against overinterpretation of the mathematical framework. As single-cell technologies continue to generate rich longitudinal datasets, approaches like this one will be essential for extracting the full biological signal from the immune landscapes that surround cancer.
Subject of Research: Continuous immune composition as a longitudinal predictor in pediatric leukemia single-cell data
Article Title: Continuous immune composition retains longitudinal information beyond coarse ecological-context labels in pediatric leukemia
Article References: Kim, S.-H. (2026). Continuous immune composition retains longitudinal information beyond coarse ecological-context labels in pediatric leukemia. BMC Bioinformatics. https://doi.org/10.1186/s12859-026-06679-7
Image Credits: AI Generated
DOI: 10.1186/s12859-026-06679-7
Keywords: pediatric leukemia, single-cell RNA-seq, immune composition, Ornstein-Uhlenbeck process, Bayesian inference, tumor microenvironment, B-cell acute lymphoblastic leukemia, longitudinal modeling, compositional data, ScPCA, ecological modeling, BMC Bioinformatics
News Source: Nathaniel Bowman. (October 17, 2026). Continuous Immune Composition Outperforms Coarse Labels in Tracking Pediatric Leukemia Over Time. Scienmag.
