Adam L. MacLean

state to early-response state, proliferate at the juncture of these two states, or become growth arrested before differentiating into spinous cells. Wound healing induces plasticity manifested by dynamic basal-spinous interconversions at multiple basal transcriptional states. Our study provides a systematic view of epidermal cellular dynamics, supporting a revised "hierarchical-lineage" model of homeostasis.

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.

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.

Motivation: Single-cell RNA-sequencing (scRNA-seq) offers unprecedented resolution for studying cellular decision-making processes. Robust inference of cell state transition paths and probabilities is an important yet challenging step in the analysis of these data. Results: Here we present scEpath, an algorithm that calculates energy landscapes and probabilistic directed graphs in order to reconstruct developmental trajectories. We quantify the energy landscape using 'single-cell energy' and distance-based measures, and find that the combination of these enables robust inference of the transition probabilities and lineage relationships between cell states. We also identify marker genes and gene expression patterns associated with cell state transitions. Our approach produces pseudotemporal orderings that are-in combination-more robust and accurate than current methods, and offers higher resolution dynamics of the cell state transitions, leading to new insight into key transition events during differentiation and development. Moreover, scEpath is robust to variation in the size of the input gene set, and is broadly unsupervised, requiring few parameters to be set by the user. Applications of scEpath led to the identification of a cell-cell communication network implicated in early human embryo development, and novel transcription factors important for myoblast differentiation. scEpath allows us to identify common and specific temporal dynamics and transcriptional factor programs along branched lineages, as well as the transition probabilities that control cell fates. Availability and implementation: A MATLAB package of scEpath is available at https://github.com/sqjin/scEpath. Supplementary information: Supplementary data are available at Bioinformatics online.

Abstract 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. Cancer Res; 77(22); 6429–41. ©2017 AACR.