Single-cell RNA-sequencing technologies suffer from many sources of technical noise, including under-sampling of mRNA molecules, often termed ‘dropout’, which can severely obscure important gene-gene relationships. To address this, we developed MAGIC (Markov Affinity-based Graph Imputation of Cells), a method that shares information across similar cells, via data diffusion, to denoise the cell count matrix and fill in missing transcripts. We validate MAGIC on several biological systems and find it effective at recovering gene-gene relationships and additional structures. MAGIC reveals a phenotypic continuum, with the majority of cells residing in intermediate states that display stem-like signatures and uncovers known and previously uncharacterized regulatory interactions, demonstrating that our approach can successfully uncover regulatory relations without perturbations.

Summary To delineate the ontogeny of the mammalian endoderm, we generated 112,217 single-cell transcriptomes representing all endoderm populations within the mouse embryo until midgestation. By employing graph-based approaches, we modelled differentiating cells for spatio-temporal characterization of developmental trajectories and defined the transcriptional architecture that accompanies the emergence of the first (primitive or extra-embryonic) endodermal population and its sister pluripotent (embryonic) epiblast lineage. We uncovered a relationship between descendants of these two lineages, whereby epiblast cells differentiate into endoderm at two distinct time-points, before and during gastrulation. Trajectories of endoderm cells were mapped as they acquired embryonic versus extra-embryonic fates, and as they spatially converged within the nascent gut endoderm; revealing them to be globally similar but retaining aspects of their lineage history. We observed the regionalized identity of cells along the anterior-posterior axis of the emergent gut tube, reflecting their embryonic or extra-embryonic origin, and their coordinate patterning into organ-specific territories.

No abstract

Immune checkpoint blockade therapy targets T cell-negative costimulatory molecules such as cytotoxic T lymphocyte antigen-4 (CTLA-4) and programmed cell death-1 (PD-1). Combination anti–CTLA-4 and anti–PD-1 blockade therapy has enhanced efficacy, but it remains unclear through what mechanisms such effects are mediated. A critical question is whether combination therapy targets and modulates the same T cell populations as monotherapies. Using a mass cytometry-based systems approach, we comprehensively profiled the response of T cell populations to monotherapy and combination anti–CTLA-4 plus anti–PD-1 therapy in syngeneic murine tumors and clinical samples. Most effects of monotherapies were additive in the context of combination therapy; however, multiple combination therapy-specific effects were observed. Highly phenotypically exhausted cluster of differentiation 8 (CD8) T cells expand in frequency following anti–PD-1 monotherapy but not combination therapy, while activated terminally differentiated effector CD8 T cells expand only following combination therapy. Combination therapy also led to further increased frequency of T helper type 1 (Th1)-like CD4 effector T cells even though anti–PD-1 monotherapy is not sufficient to do so. Mass cytometry analyses of peripheral blood from melanoma patients treated with immune checkpoint blockade therapies similarly revealed mostly additive effects on the frequencies of T cell subsets along with unique modulation of terminally differentiated effector CD8 T cells by combination ipilimumab plus nivolumab therapy. Together, these findings indicate that dual blockade of CTLA-4 and PD-1 therapy is sufficient to induce unique cellular responses compared with either monotherapy.

Metastasis-initiating cells with stem-like properties drive cancer lethality, yet their origins and relationship to primary-tumorinitiating stem cells are not known. We show that L1CAM + cells in human colorectal cancer (CRC) have metastasis-initiating capacity, and we define their relationship to tissue regeneration. L1CAM is not expressed in the homeostatic intestinal epithelium, but is induced and required for epithelial regeneration following colitis and in CRC organoid growth. By using human tissues and mouse models, we show that L1CAM is dispensable for adenoma initiation but required for orthotopic carcinoma propagation, liver metastatic colonization and chemoresistance. L1CAM high cells partially overlap with LGR5 high stem-like cells in human CRC organoids. Disruption of intercellular epithelial contacts causes E-cadherin-REST transcriptional derepression of L1CAM, switching chemoresistant CRC progenitors from an L1CAM low to an L1CAM high state. Thus, L1CAM dependency emerges in regenerative intestinal cells when epithelial integrity is lost, a phenotype of wound healing deployed in metastasisinitiating cells.

Purpose: Human papillomavirus (HPV)-negative head and neck squamous cell carcinomas (HNSCC) commonly bear disruptive mutations in TP53, resulting in treatment resistance. In these patients, direct targeting of p53 has not been successful, but synthetic lethal approaches have promise. Although Aurora A kinase (AURKA) is overexpressed and an oncogenic driver, its inhibition has only modest clinical effects in HPV-negative HNSCC. We explored a novel combination of AURKA and WEE1 inhibition to overcome intrinsic resistance to AURKA inhibition. Experimental Design: AURKA protein expression was determined by fluorescence-based automated quantitative analysis of patient specimens and correlated with survival. We evaluated treatment with the AURKA inhibitor alisertib (MLN8237) and the WEE1 inhibitor adavosertib (AZD1775), alone or in combination, using in vitro and in vivo HNSCC models. Results: Elevated nuclear AURKA correlated with worse survival among p16(−) HNSCC patients. Alisertib caused spindle defects, G2/M arrest and inhibitory CDK1 phosphorylation, and cytostasis in TP53 mutant HNSCC FaDu and UNC7 cells. Addition of adavosertib to alisertib instead triggered mitotic entry and mitotic catastrophe. Moreover, in FaDu and Detroit 562 xenografts, this combination demonstrated synergistic effects on tumor growth and extended overall survival compared to either vehicle or single agent treatment. Conclusions: Combinatorial treatment with adavosertib and alisertib leads to synergistic antitumor effects in in vitro and in vivo HNSCC models. These findings suggest a novel rational combination, providing a promising therapeutic avenue for TP53-mutated cancers.

Metacells are cell groupings derived from single-cell sequencing data that represent highly granular, distinct cell states. Here we present single-cell aggregation of cell states (SEACells), an algorithm for identifying metacells that overcome the sparsity of single-cell data while retaining heterogeneity obscured by traditional cell clustering. SEACells outperforms existing algorithms in identifying comprehensive, compact and well-separated metacells in both RNA and assay for transposase-accessible chromatin (ATAC) modalities across datasets with discrete cell types and continuous trajectories. We demonstrate the use of SEACells to improve gene–peak associations, compute ATAC gene scores and infer the activities of critical regulators during differentiation. Metacell-level analysis scales to large datasets and is particularly well suited for patient cohorts, where per-patient aggregation provides more robust units for data integration. We use our metacells to reveal expression dynamics and gradual reconfiguration of the chromatin landscape during hematopoietic differentiation and to uniquely identify CD4 T cell differentiation and activation states associated with disease onset and severity in a Coronavirus Disease 2019 (COVID-19) patient cohort.

Highlights d Genetic variation in evolutionarily distant mice regulates gene expression in cis d Distinct TF families mediate stable and transient chromatinlevel responses d T-box and Runx family TFs cooperate in stable chromatin remodeling in CD8 T cells d Stability is associated with congruent remodeling of neighboring chromatin sites