A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Defining the Design Principles of Skin Epidermis Postnatal Growth.

During embryonic and postnatal development, organs and tissues grow steadily to achieve their final size at the end of puberty. However, little is known about the cellular dynamics that mediate postnatal growth. By combining in vivo clonal lineage tracing, proliferation kinetics, single-cell transcriptomics, and in vitro micro-pattern experiments, we resolved the cellular dynamics taking place during postnatal skin epidermis expansion. Our data revealed that harmonious growth is engineered by a single population of developmental progenitors presenting a fixed fate imbalance of self-renewing divisions with an ever-decreasing proliferation rate. Single-cell RNA sequencing revealed that epidermal developmental progenitors form a more uniform population compared with adult stem and progenitor cells. Finally, we found that the spatial pattern of cell division orientation is dictated locally by the underlying collagen fiber orientation. Our results uncover a simple design principle of organ growth where progenitors and differentiated cells expand in harmony with their surrounding tissues.

CD4+ T cell activation distinguishes response to anti-PD-L1+anti-CTLA4 therapy from anti-PD-L1 monotherapy.

Cancer patients often receive a combination of antibodies targeting programmed death-ligand 1 (PD-L1) and cytotoxic T lymphocyte antigen-4 (CTLA4). We conducted a window-of-opportunity study in head and neck squamous cell carcinoma (HNSCC) to examine the contribution of anti-CTLA4 to anti-PD-L1 therapy. Single-cell profiling of on- versus pre-treatment biopsies identified T cell expansion as an early response marker. In tumors, anti-PD-L1 triggered the expansion of mostly CD8+ T cells, whereas combination therapy expanded both CD4+ and CD8+ T cells. Such CD4+ T cells exhibited an activated T helper 1 (Th1) phenotype. CD4+ and CD8+ T cells co-localized with and were surrounded by dendritic cells expressing T cell homing factors or antibody-producing plasma cells. T cell receptor tracing suggests that anti-CTLA4, but not anti-PD-L1, triggers the trafficking of CD4+ naive/central-memory T cells from tumor-draining lymph nodes (tdLNs), via blood, to the tumor wherein T cells acquire a Th1 phenotype. Thus, CD4+ T cell activation and recruitment from tdLNs are hallmarks of early response to anti-PD-L1 plus anti-CTLA4 in HNSCC.

Toward Minimal Residual Disease-Directed Therapy in Melanoma.

Many patients with advanced cancers achieve dramatic responses to a panoply of therapeutics yet retain minimal residual disease (MRD), which ultimately results in relapse. To gain insights into the biology of MRD, we applied single-cell RNA sequencing to malignant cells isolated from BRAF mutant patient-derived xenograft melanoma cohorts exposed to concurrent RAF/MEK-inhibition. We identified distinct drug-tolerant transcriptional states, varying combinations of which co-occurred within MRDs from PDXs and biopsies of patients on treatment. One of these exhibited a neural crest stem cell (NCSC) transcriptional program largely driven by the nuclear receptor RXRG. An RXR antagonist mitigated accumulation of NCSCs in MRD and delayed the development of resistance. These data identify NCSCs as key drivers of resistance and illustrate the therapeutic potential of MRD-directed therapy. They also highlight how gene regulatory network architecture reprogramming may be therapeutically exploited to limit cellular heterogeneity, a key driver of disease progression and therapy resistance.

A Single-Cell Transcriptome Atlas of the Aging Drosophila Brain.

The diversity of cell types and regulatory states in the brain, and how these change during aging, remains largely unknown. We present a single-cell transcriptome atlas of the entire adult Drosophila melanogaster brain sampled across its lifespan. Cell clustering identified 87 initial cell clusters that are further subclustered and validated by targeted cell-sorting. Our data show high granularity and identify a wide range of cell types. Gene network analyses using SCENIC revealed regulatory heterogeneity linked to energy consumption. During aging, RNA content declines exponentially without affecting neuronal identity in old brains. This single-cell brain atlas covers nearly all cells in the normal brain and provides the tools to study cellular diversity alongside other Drosophila and mammalian single-cell datasets in our unique single-cell analysis platform: SCope (http://scope.aertslab.org). These results, together with SCope, allow comprehensive exploration of all transcriptional states of an entire aging brain.

Mutational Processes Shaping the Genome in Early Human Embryos.

In-vitro-fertilized human embryos often acquire large structural and numerical chromosomal abnormalities. Liu et al. now show that multiple smaller copy number variations may also arise in in-vivo-conceived embryos. Analysis of these variations provides insight into the DNA mutational processes occurring in early embryos and the mechanisms underlying them.

Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos.

The major and essential objective of pre-implantation development is to establish embryonic and extra-embryonic cell fates. To address when and how this fundamental process is initiated in mammals, we characterize transcriptomes of all individual cells throughout mouse pre-implantation development. This identifies targets of master pluripotency regulators Oct4 and Sox2 as being highly heterogeneously expressed between blastomeres of the 4-cell embryo, with Sox21 showing one of the most heterogeneous expression profiles. Live-cell tracking demonstrates that cells with decreased Sox21 yield more extra-embryonic than pluripotent progeny. Consistently, decreasing Sox21 results in premature upregulation of the differentiation regulator Cdx2, suggesting that Sox21 helps safeguard pluripotency. Furthermore, Sox21 is elevated following increased expression of the histone H3R26-methylase CARM1 and is lowered following CARM1 inhibition, indicating the importance of epigenetic regulation. Therefore, our results indicate that heterogeneous gene expression, as early as the 4-cell stage, initiates cell-fate decisions by modulating the balance of pluripotency and differentiation.

Methods and applications for single-cell and spatial multi-omics.

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

Dedicated macrophages organize and maintain the enteric nervous system.

Correct development and maturation of the enteric nervous system (ENS) is critical for survival1. At birth, the ENS is immature and requires considerable refinement to exert its functions in adulthood2. Here we demonstrate that resident macrophages of the muscularis externa (MMϕ) refine the ENS early in life by pruning synapses and phagocytosing enteric neurons. Depletion of MMϕ before weaning disrupts this process and results in abnormal intestinal transit. After weaning, MMϕ continue to interact closely with the ENS and acquire a neurosupportive phenotype. The latter is instructed by transforming growth factor-ß produced by the ENS; depletion of the ENS and disruption of transforming growth factor-ß signalling result in a decrease in neuron-associated MMϕ associated with loss of enteric neurons and altered intestinal transit. These findings introduce a new reciprocal cell-cell communication responsible for maintenance of the ENS and indicate that the ENS, similarly to the brain, is shaped and maintained by a dedicated population of resident macrophages that adapts its phenotype and transcriptome to the timely needs of the ENS niche.

Pharmacological targeting of netrin-1 inhibits EMT in cancer.

Epithelial-to-mesenchymal transition (EMT) regulates tumour initiation, progression, metastasis and resistance to anti-cancer therapy1-7. Although great progress has been made in understanding the role of EMT and its regulatory mechanisms in cancer, no therapeutic strategy to pharmacologically target EMT has been identified. Here we found that netrin-1 is upregulated in a primary mouse model of skin squamous cell carcinoma (SCC) exhibiting spontaneous EMT. Pharmacological inhibition of netrin-1 by administration of NP137, a netrin-1-blocking monoclonal antibody currently used in clinical trials in human cancer (ClinicalTrials.gov identifier NCT02977195 ), decreased the proportion of EMT tumour cells in skin SCC, decreased the number of metastases and increased the sensitivity of tumour cells to chemotherapy. Single-cell RNA sequencing revealed the presence of different EMT states, including epithelial, early and late hybrid EMT, and full EMT states, in control SCC. By contrast, administration of NP137 prevented the progression of cancer cells towards a late EMT state and sustained tumour epithelial states. Short hairpin RNA knockdown of netrin-1 and its receptor UNC5B in EPCAM+ tumour cells inhibited EMT in vitro in the absence of stromal cells and regulated a common gene signature that promotes tumour epithelial state and restricts EMT. To assess the relevance of these findings to human cancers, we treated mice transplanted with the A549 human cancer cell line-which undergoes EMT following TGFß1 administration8,9-with NP137. Netrin-1 inhibition decreased EMT in these transplanted A549 cells. Together, our results identify a pharmacological strategy for targeting EMT in cancer, opening up novel therapeutic interventions for anti-cancer therapy.

Mechanisms of stretch-mediated skin expansion at single-cell resolution.

The ability of the skin to grow in response to stretching has been exploited in reconstructive surgery1. Although the response of epidermal cells to stretching has been studied in vitro2,3, it remains unclear how mechanical forces affect their behaviour in vivo. Here we develop a mouse model in which the consequences of stretching on skin epidermis can be studied at single-cell resolution. Using a multidisciplinary approach that combines clonal analysis with quantitative modelling and single-cell RNA sequencing, we show that stretching induces skin expansion by creating a transient bias in the renewal activity of epidermal stem cells, while a second subpopulation of basal progenitors remains committed to differentiation. Transcriptional and chromatin profiling identifies how cell states and gene-regulatory networks are modulated by stretching. Using pharmacological inhibitors and mouse mutants, we define the step-by-step mechanisms that control stretch-mediated tissue expansion at single-cell resolution in vivo.

Heterotypic cell-cell communication regulates glandular stem cell multipotency.

Glandular epithelia, including the mammary and prostate glands, are composed of basal cells (BCs) and luminal cells (LCs)1,2. Many glandular epithelia develop from multipotent basal stem cells (BSCs) that are replaced in adult life by distinct pools of unipotent stem cells1,3-8. However, adult unipotent BSCs can reactivate multipotency under regenerative conditions and upon oncogene expression3,9-13. This suggests that an active mechanism restricts BSC multipotency under normal physiological conditions, although the nature of this mechanism is unknown. Here we show that the ablation of LCs reactivates the multipotency of BSCs from multiple epithelia both in vivo in mice and in vitro in organoids. Bulk and single-cell RNA sequencing revealed that, after LC ablation, BSCs activate a hybrid basal and luminal cell differentiation program before giving rise to LCs-reminiscent of the genetic program that regulates multipotency during embryonic development7. By predicting ligand-receptor pairs from single-cell data14, we find that TNF-which is secreted by LCs-restricts BC multipotency under normal physiological conditions. By contrast, the Notch, Wnt and EGFR pathways were activated in BSCs and their progeny after LC ablation; blocking these pathways, or stimulating the TNF pathway, inhibited regeneration-induced BC multipotency. Our study demonstrates that heterotypic communication between LCs and BCs is essential to maintain lineage fidelity in glandular epithelial stem cells.

LifeTime and improving European healthcare through cell-based interceptive medicine.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

Adipose tissue macrophage dysfunction is associated with a breach of vascular integrity in NASH.

BACKGROUND & AIMS: In non-alcoholic fatty liver disease (NAFLD), monocytes infiltrate visceral adipose tissue promoting local and hepatic inflammation. However, it remains unclear what drives inflammation and how the immune landscape in adipose tissue differs across the NAFLD severity spectrum. We aimed to assess adipose tissue macrophage (ATM) heterogeneity in a NAFLD cohort. METHODS: Visceral adipose tissue macrophages from lean and obese patients, stratified by NAFLD phenotypes, underwent single-cell RNA sequencing. Adipose tissue vascular integrity and breaching was assessed on a protein level via immunohistochemistry and immunofluorescence to determine targets of interest. RESULTS: We discovered multiple ATM populations, including resident vasculature-associated macrophages (ResVAMs) and distinct metabolically active macrophages (MMacs). Using trajectory analysis, we show that ResVAMs and MMacs are replenished by a common transitional macrophage (TransMac) subtype and that, during NASH, MMacs are not effectively replenished by TransMac precursors. We postulate an accessory role for MMacs and ResVAMs in protecting the adipose tissue vascular barrier, since they both interact with endothelial cells and localize around the vasculature. However, across the NAFLD severity spectrum, alterations occur in these subsets that parallel an adipose tissue vasculature breach characterized by albumin extravasation into the perivascular tissue. CONCLUSIONS: NAFLD-related macrophage dysfunction coincides with a loss of adipose tissue vascular integrity, providing a plausible mechanism by which tissue inflammation is perpetuated in adipose tissue and downstream in the liver. IMPACT AND IMPLICATIONS: Our study describes for the first time the myeloid cell landscape in human visceral adipose tissue at single-cell level within a cohort of well-characterized patients with non-alcoholic fatty liver disease. We report unique non-alcoholic steatohepatitis-specific transcriptional changes within metabolically active macrophages (MMacs) and resident vasculature-associated macrophages (ResVAMs) and we demonstrate their spatial location surrounding the vasculature. These dysfunctional transcriptional macrophage states coincided with the loss of adipose tissue vascular integrity, providing a plausible mechanism by which tissue inflammation is perpetuated in adipose tissue and downstream in the liver. Our study provides a theoretical basis for new therapeutic strategies to be directed towards reinstating the endogenous metabolic, homeostatic and cytoprotective functions of ResVAMs and MMacs, including their role in protecting vascular integrity.

Publisher Correction: A TRP channel trio mediates acute noxious heat sensing.

In this Letter, the trace is missing in Fig. 1e. This error has been corrected online.

A TRP channel trio mediates acute noxious heat sensing.

Acute pain represents a crucial alarm signal to protect us from injury. Whereas the nociceptive neurons that convey pain signals were described more than a century ago, the molecular sensors that detect noxious thermal or mechanical insults have yet to be fully identified. Here we show that acute noxious heat sensing in mice depends on a triad of transient receptor potential (TRP) ion channels: TRPM3, TRPV1, and TRPA1. We found that robust somatosensory heat responsiveness at the cellular and behavioural levels is observed only if at least one of these TRP channels is functional. However, combined genetic or pharmacological elimination of all three channels largely and selectively prevents heat responses in both isolated sensory neurons and rapidly firing C and Aδ sensory nerve fibres that innervate the skin. Strikingly, Trpv1-/-Trpm3-/-Trpa1-/- triple knockout (TKO) mice lack the acute withdrawal response to noxious heat that is necessary to avoid burn injury, while showing normal nociceptive responses to cold or mechanical stimuli and a preserved preference for moderate temperatures. These findings indicate that the initiation of the acute heat-evoked pain response in sensory nerve endings relies on three functionally redundant TRP channels, representing a fault-tolerant mechanism to avoid burn injury.

Identification of the tumour transition states occurring during EMT.

In cancer, the epithelial-to-mesenchymal transition (EMT) is associated with tumour stemness, metastasis and resistance to therapy. It has recently been proposed that, rather than being a binary process, EMT occurs through distinct intermediate states. However, there is no direct in vivo evidence for this idea. Here we screen a large panel of cell surface markers in skin and mammary primary tumours, and identify the existence of multiple tumour subpopulations associated with different EMT stages: from epithelial to completely mesenchymal states, passing through intermediate hybrid states. Although all EMT subpopulations presented similar tumour-propagating cell capacity, they displayed differences in cellular plasticity, invasiveness and metastatic potential. Their transcriptional and epigenetic landscapes identify the underlying gene regulatory networks, transcription factors and signalling pathways that control these different EMT transition states. Finally, these tumour subpopulations are localized in different niches that differentially regulate EMT transition states.

Erratum: Pluripotent state transitions coordinate morphogenesis in mouse and human embryos.

This corrects the article DOI: 10.1038/nature24675.

Expansion microscopy allows high resolution single cell analysis of epigenetic readers.

Interactions between epigenetic readers and histone modifications play a pivotal role in gene expression regulation and aberrations can enact etiopathogenic roles in both developmental and acquired disorders like cancer. Typically, epigenetic interactions are studied by mass spectrometry or chromatin immunoprecipitation sequencing. However, in these methods, spatial information is completely lost. Here, we devise an expansion microscopy based method, termed Expansion Microscopy for Epigenetics or ExEpi, to preserve spatial information and improve resolution. We calculated relative co-localization ratios for two epigenetic readers, lens epithelium derived growth factor (LEDGF) and bromodomain containing protein 4 (BRD4), with marks for heterochromatin (H3K9me3 and H3K27me3) and euchromatin (H3K36me2, H3K36me3 and H3K9/14ac). ExEpi confirmed their preferred epigenetic interactions, showing co-localization for LEDGF with H3K36me3/me2 and for BRD4 with H3K9/14ac. Moreover addition of JQ1, a known BET-inhibitor, abolished BRD4 interaction with H3K9/14ac with an IC50 of 137 nM, indicating ExEpi could serve as a platform for epigenetic drug discovery. Since ExEpi retains spatial information, the nuclear localization of marks and readers was determined, which is one of the main advantages of ExEpi. The heterochromatin mark, H3K9me3, is located in the nuclear rim whereas LEDGF co-localization with H3K36me3 and BRD4 co-localization with H3K9/14ac occur further inside the nucleus.

Single-cell genome-wide concurrent haplotyping and copy-number profiling through genotyping-by-sequencing.

Single-cell whole-genome haplotyping allows simultaneous detection of haplotypes associated with monogenic diseases, chromosome copy-numbering and subsequently, has revealed mosaicism in embryos and embryonic stem cells. Methods, such as karyomapping and haplarithmisis, were deployed as a generic and genome-wide approach for preimplantation genetic testing (PGT) and are replacing traditional PGT methods. While current methods primarily rely on single-nucleotide polymorphism (SNP) array, we envision sequencing-based methods to become more accessible and cost-efficient. Here, we developed a novel sequencing-based methodology to haplotype and copy-number profile single cells. Following DNA amplification, genomic size and complexity is reduced through restriction enzyme digestion and DNA is genotyped through sequencing. This single-cell genotyping-by-sequencing (scGBS) is the input for haplarithmisis, an algorithm we previously developed for SNP array-based single-cell haplotyping. We established technical parameters and developed an analysis pipeline enabling accurate concurrent haplotyping and copy-number profiling of single cells. We demonstrate its value in human blastomere and trophectoderm samples as application for PGT for monogenic disorders. Furthermore, we demonstrate the method to work in other species through analyzing blastomeres of bovine embryos. Our scGBS method opens up the path for single-cell haplotyping of any species with diploid genomes and could make its way into the clinic as a PGT application.