How stem cells give rise to epidermis is unclear despite the crucial role the epidermis plays in barrier and appendage formation. Here we use single cell-RNA sequencing to interrogate basal stem cell heterogeneity of human interfollicular epidermis and find four spatially distinct stem cell populations at the top and bottom of rete ridges and transitional positions between the basal and suprabasal epidermal layers. Cell-cell communication modeling suggests that basal cell populations serve as crucial signaling hubs to maintain epidermal communication. Combining pseudotime, RNA velocity, and cellular entropy analyses point to a hierarchical differentiation lineage supporting multi-stem cell interfollicular epidermal homeostasis models and suggest that transitional basal stem cells are stable states essential for proper stratification. Finally, alterations in differentially expressed transitional basal stem cell genes result in severe thinning of human skin equivalents, validating their essential role in epidermal homeostasis and reinforcing the critical nature of basal stem cell heterogeneity.

Supplementary data are available at Bioinformatics online.

The use of single-cell transcriptomics has become a major approach to delineate cell subpopulations and the transitions between them. While various computational tools using different mathematical methods have been developed to infer clusters, marker genes, and cell lineage, none yet integrate these within a mathematical framework to perform multiple tasks coherently. Such coherence is critical for the inference of cell–cell communication, a major remaining challenge. Here, we present similarity matrix-based optimization for single-cell data analysis (SoptSC), in which unsupervised clustering, pseudotemporal ordering, lineage inference, and marker gene identification are inferred via a structured cell-to-cell similarity matrix. SoptSC then predicts cell–cell communication networks, enabling reconstruction of complex cell lineages that include feedback or feedforward interactions. Application of SoptSC to early embryonic development, epidermal regeneration, and hematopoiesis demonstrates robust identification of subpopulations, lineage relationships, and pseudotime, and prediction of pathway-specific cell communication patterns regulating processes of development and differentiation.

Highlights d scRNA-seq identifies four epidermal basal cell states in homeostatic adult skin d Computational analysis supports a ''hierarchical'' model of epidermal homeostasis d Basal cell states are metabolically distinct and spatially partitioned in wounded skin d Epidermal basal cells show enhanced cell fate and state plasticity during wound healing

As our catalog of cell states expands, appropriate characterization of these states and the transitions between them is crucial. Here we discuss the roles of intermediate cell states (ICSs) in this growing collection. We begin with definitions and discuss evidence for the existence of ICSs and their relevance in various tissues. We then provide a list of possible functions for ICSs with examples. Finally, we describe means by which ICSs and their functional roles can be identified from single-cell data or predicted from models.

Highlights d scRNA-seq identifies four epidermal basal cell states in homeostatic adult skin d Computational analysis supports a ''hierarchical'' model of epidermal homeostasis d Basal cell states are metabolically distinct and spatially partitioned in wounded skin d Epidermal basal cells show enhanced cell fate and state plasticity during wound healing

Chronic inflammation is a serious risk factor for cancer; however, the routes from inflammation to cancer are poorly understood. On the basis of the processes implicated by frequently mutated genes associated with inflammation and cancer in three organs (stomach, colon, and liver) extracted from the Gene Expression Omnibus, The Cancer Genome Atlas, and Gene Ontology databases, we present a multiscale model of the long-term evolutionary dynamics leading from inflammation to tumorigenesis. The model incorporates cross-talk among interactions on several scales, including responses to DNA damage, gene mutation, cell-cycle behavior, population dynamics, inflammation, and metabolism-immune balance. Model simulations revealed two stages of inflammation-induced tumorigenesis: a precancerous state and tumorigenesis. The precancerous state was mainly caused by mutations in the cell proliferation pathway; the transition from the precancerous to tumorigenic states was induced by mutations in pathways associated with apoptosis, differentiation, and metabolism-immune balance. We identified opposing effects of inflammation on tumorigenesis. Mild inflammation removed cells with DNA damage through DNA damage-induced cell death, whereas severe inflammation accelerated accumulation of mutations and hence promoted tumorigenesis. These results provide insight into the evolutionary dynamics of inflammation-induced tumorigenesis and highlight the combinatorial effects of inflammation and metabolism-immune balance. This approach establishes methods for quantifying cancer risk, for the discovery of driver pathways in inflammation-induced tumorigenesis, and has direct relevance for early detection and prevention and development of new treatment regimes. .

Chronic myeloid leukemia (CML) is a blood disease that disrupts normal function of the hematopoietic system. Despite the great progress made in terms of molecular therapies for CML, there remain large gaps in our understanding. By comparing mathematical models that describe CML progression and etiology we sought to identify those models that provide the best description of disease dynamics and their underlying mechanisms. Data for two clinical outcomesdisease remission or relapse-are considered, and we investigate these using Bayesian inference techniques throughout. We find that it is not possible to choose between the models based on fits to the data alone; however, by studying model predictions we can discard models that fail to take niche effects into account. More detailed analysis of the remaining models reveals mechanistic differences: for one model, leukemia stem cell dynamics determine the disease outcome; and for the other model disease progression is determined at the stage of progenitor cells, in particular by differences in progenitor death rates. This analysis also reveals distinct transient dynamics that will be experimentally accessible, but are currently at the limits of what is possible to measure. To resolve these differences we need to be able to probe the hematopoietic stem cell niche directly. Our analysis highlights the importance of further mapping of the bone marrow hematopoietic niche microenvironment as the "ecological" interactions between cells in this niche appear to be intricately linked to disease outcome.ematopoiesis is the intricate process that provides the vast diversity of cells found in the blood and immune system of vertebrate organisms. To generate great numbers of very different cellsranging from erythrocytes to highly specialized lymphocytes-a delicate balance of activity of transcriptional programs and hematopoietic growth factors is required (1). Within this system, hematopoietic stem cells (HSCs) give rise to (lymphoid and myeloid) progenitors, which in turn through further differentiation steps yield the panoply of cells found in the blood (2). Given the complexity of the hematopoietic system it is remarkable how robustly this process works in general.The cellular and molecular processes involved in healthy hematopoiesis must be disrupted in disease (3). There is however solid evidence that disease processes operate in parallel to the normal hematopoietic system. Thus, a branch of the hematopoietic system will propagate deficient cells and generate the ensuing pathology, but other cells of the system will continue to operate alongside (4-7). This is, of course, akin to problems encountered in ecology. Based on this analogy we here attempt to understand the interplay between the healthy and chronic myeloid leukemia (CML) systems. This is a particularly helpful perspective if interactions between healthy and diseased cell types can affect the disease progression, as is now widely believed.Interactions will include shared reliance on nutrients, cofactors, and physiologi...