DURIAN: an integrative deconvolution and imputation method for robust signaling analysis of single-cell transcriptomics data.

Single-cell RNA sequencing trades read-depth for dimensionality, often leading to loss of critical signaling gene information that is typically present in bulk data sets. We introduce DURIAN (Deconvolution and mUltitask-Regression-based ImputAtioN), an integrative method for recovery of gene expression in single-cell data. Through systematic benchmarking, we demonstrate the accuracy, robustness and empirical convergence of DURIAN using both synthetic and published data sets. We show that use of DURIAN improves single-cell clustering, low-dimensional embedding, and recovery of intercellular signaling networks. Our study resolves several inconsistent results of cell-cell communication analysis using single-cell or bulk data independently. The method has broad application in biomarker discovery and cell signaling analysis using single-cell transcriptomics data sets.

Noncatalytic PTEN missense mutation predisposes to organ-selective cancer development in vivo.

Inactivation of phosphatase and tensin homology deleted on chromosome 10 (PTEN) is linked to increased PI3K-AKT signaling, enhanced organismal growth, and cancer development. Here we generated and analyzed Pten knock-in mice harboring a C2 domain missense mutation at phenylalanine 341 (Pten(FV)), found in human cancer. Despite having reduced levels of PTEN protein, homozygous Pten(FV/FV) embryos have intact AKT signaling, develop normally, and are carried to term. Heterozygous Pten(FV/+) mice develop carcinoma in the thymus, stomach, adrenal medulla, and mammary gland but not in other organs typically sensitive to Pten deficiency, including the thyroid, prostate, and uterus. Progression to carcinoma in sensitive organs ensues in the absence of overt AKT activation. Carcinoma in the uterus, a cancer-resistant organ, requires a second clonal event associated with the spontaneous activation of AKT and downstream signaling. In summary, this PTEN noncatalytic missense mutation exposes a core tumor suppressor function distinct from inhibition of canonical AKT signaling that predisposes to organ-selective cancer development in vivo.

Cell lineage and communication network inference via optimization for single-cell transcriptomics.

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.

Sequential emergence and contraction of epithelial subtypes in the prenatal human choroid plexus revealed by a stem cell model.

Despite the major roles of choroid plexus epithelial cells (CPECs) in brain homeostasis and repair, their developmental lineage and diversity remain undefined. In simplified differentiations from human pluripotent stem cells, derived CPECs (dCPECs) displayed canonical properties and dynamic multiciliated phenotypes that interacted with Aß uptake. Single dCPEC transcriptomes over time correlated well with human organoid and fetal CPECs, while pseudotemporal and cell cycle analyses highlighted the direct CPEC origin from neuroepithelial cells. In addition, time series analyses defined metabolic (type 1) and ciliogenic dCPECs (type 2) at early timepoints, followed by type 1 diversification into anabolic-secretory (type 1a) and catabolic-absorptive subtypes (type 1b) as type 2 cells contracted. These temporal patterns were then confirmed in independent derivations and mapped to prenatal stages using human tissues. In addition to defining the prenatal lineage of human CPECs, these findings suggest new dynamic models of ChP support for the developing human brain.

Single-cell analysis of human basal cell carcinoma reveals novel regulators of tumor growth and the tumor microenvironment.

How basal cell carcinoma (BCC) interacts with its tumor microenvironment to promote growth is unclear. We use singe-cell RNA sequencing to define the human BCC ecosystem and discriminate between normal and malignant epithelial cells. We identify spatial biomarkers of tumors and their surrounding stroma that reinforce the heterogeneity of each tissue type. Combining pseudotime, RNA velocity-PAGA, cellular entropy, and regulon analysis in stromal cells reveals a cancer-specific rewiring of fibroblasts, where STAT1, TGF-ß, and inflammatory signals induce a noncanonical WNT5A program that maintains the stromal inflammatory state. Cell-cell communication modeling suggests that tumors respond to the sudden burst of fibroblast-specific inflammatory signaling pathways by producing heat shock proteins, whose expression we validated in situ. Last, dose-dependent treatment with an HSP70 inhibitor suppresses in vitro vismodegib-resistant BCC cell growth, Hedgehog signaling, and in vivo tumor growth in a BCC mouse model, validating HSP70's essential role in tumor growth and reinforcing the critical nature of tumor microenvironment cross-talk in BCC progression.

Modeling the effects of EMT-immune dynamics on carcinoma disease progression.

During progression from carcinoma in situ to an invasive tumor, the immune system is engaged in complex sets of interactions with various tumor cells. Tumor cell plasticity alters disease trajectories via epithelial-to-mesenchymal transition (EMT). Several of the same pathways that regulate EMT are involved in tumor-immune interactions, yet little is known about the mechanisms and consequences of crosstalk between these regulatory processes. Here we introduce a multiscale evolutionary model to describe tumor-immune-EMT interactions and their impact on epithelial cancer progression from in situ to invasive disease. Through simulation of patient cohorts in silico, the model predicts that a controllable region maximizes invasion-free survival. This controllable region depends on properties of the mesenchymal tumor cell phenotype: its growth rate and its immune-evasiveness. In light of the model predictions, we analyze EMT-inflammation-associated data from The Cancer Genome Atlas, and find that association with EMT worsens invasion-free survival probabilities. This result supports the predictions of the model, and leads to the identification of genes that influence outcomes in bladder and uterine cancer, including FGF pathway members. These results suggest new means to delay disease progression, and demonstrate the importance of studying cancer-immune interactions in light of EMT.

Single-cell transcriptomic analysis of zebrafish cranial neural crest reveals spatiotemporal regulation of lineage decisions during development.

Neural crest (NC) cells migrate throughout vertebrate embryos to give rise to a huge variety of cell types, but when and where lineages emerge and their regulation remain unclear. We have performed single-cell RNA sequencing (RNA-seq) of cranial NC cells from the first pharyngeal arch in zebrafish over several stages during migration. Computational analysis combining pseudotime and real-time data reveals that these NC cells first adopt a transitional state, becoming specified mid-migration, with the first lineage decisions being skeletal and pigment, followed by neural and glial progenitors. In addition, by computationally integrating these data with RNA-seq data from a transgenic Wnt reporter line, we identify gene cohorts with similar temporal responses to Wnts during migration and show that one, Atp6ap2, is required for melanocyte differentiation. Together, our results show that cranial NC cell lineages arise progressively and uncover a series of spatially restricted cell interactions likely to regulate such cell-fate decisions.