EmptyDropsMultiome discriminates real cells from background in single-cell multiomics assays.

Multiomic droplet-based technologies allow different molecular modalities, such as chromatin accessibility and gene expression (scATAC-seq and scRNA-seq), to be probed in the same nucleus. We develop EmptyDropsMultiome, an approach that distinguishes true nuclei-containing droplets from background. Using simulations, we show that EmptyDropsMultiome has higher statistical power and accuracy than existing approaches, including CellRanger-arc and EmptyDrops. On real datasets, we observe that CellRanger-arc misses more than half of the nuclei identified by EmptyDropsMultiome and, moreover, is biased against certain cell types, some of which have a retrieval rate lower than 20%.

The Virtual Child.

SUMMARY: We are building the world's first Virtual Child-a computer model of normal and cancerous human development at the level of each individual cell. The Virtual Child will "develop cancer" that we will subject to unlimited virtual clinical trials that pinpoint, predict, and prioritize potential new treatments, bringing forward the day when no child dies of cancer, giving each one the opportunity to lead a full and healthy life.

A single-cell atlas enables mapping of homeostatic cellular shifts in the adult human breast.

Here we use single-cell RNA sequencing to compile a human breast cell atlas assembled from 55 donors that had undergone reduction mammoplasties or risk reduction mastectomies. From more than 800,000 cells we identified 41 cell subclusters across the epithelial, immune and stromal compartments. The contribution of these different clusters varied according to the natural history of the tissue. Age, parity and germline mutations, known to modulate the risk of developing breast cancer, affected the homeostatic cellular state of the breast in different ways. We found that immune cells from BRCA1 or BRCA2 carriers had a distinct gene expression signature indicative of potential immune exhaustion, which was validated by immunohistochemistry. This suggests that immune-escape mechanisms could manifest in non-cancerous tissues very early during tumor initiation. This atlas is a rich resource that can be used to inform novel approaches for early detection and prevention of breast cancer.

The dynamic genetic determinants of increased transcriptional divergence in spermatids.

Cis-genetic effects are key determinants of transcriptional divergence in discrete tissues and cell types. However, how cis- and trans-effects act across continuous trajectories of cellular differentiation in vivo is poorly understood. Here, we quantify allele-specific expression during spermatogenic differentiation at single-cell resolution in an F1 hybrid mouse system, allowing for the comprehensive characterisation of cis- and trans-genetic effects, including their dynamics across cellular differentiation. Collectively, almost half of the genes subject to genetic regulation show evidence for dynamic cis-effects that vary during differentiation. Our system also allows us to robustly identify dynamic trans-effects, which are less pervasive than cis-effects. In aggregate, genetic effects were strongest in round spermatids, which parallels their increased transcriptional divergence we identified between species. Our approach provides a comprehensive quantification of the variability of genetic effects in vivo, and demonstrates a widely applicable strategy to dissect the impact of regulatory variants on gene regulation in dynamic systems.

scSNV-seq: high-throughput phenotyping of single nucleotide variants by coupled single-cell genotyping and transcriptomics.

CRISPR screens with single-cell transcriptomic readouts are a valuable tool to understand the effect of genetic perturbations including single nucleotide variants (SNVs) associated with diseases. Interpretation of these data is currently limited as genotypes cannot be accurately inferred from guide RNA identity alone. scSNV-seq overcomes this limitation by coupling single-cell genotyping and transcriptomics of the same cells enabling accurate and high-throughput screening of SNVs. Analysis of variants across the JAK1 gene with scSNV-seq demonstrates the importance of determining the precise genetic perturbation and accurately classifies clinically observed missense variants into three functional categories: benign, loss of function, and separation of function.

Stabilized mosaic single-cell data integration using unshared features.

Currently available single-cell omics technologies capture many unique features with different biological information content. Data integration aims to place cells, captured with different technologies, onto a common embedding to facilitate downstream analytical tasks. Current horizontal data integration techniques use a set of common features, thereby ignoring non-overlapping features and losing information. Here we introduce StabMap, a mosaic data integration technique that stabilizes mapping of single-cell data by exploiting the non-overlapping features. StabMap first infers a mosaic data topology based on shared features, then projects all cells onto supervised or unsupervised reference coordinates by traversing shortest paths along the topology. We show that StabMap performs well in various simulation contexts, facilitates 'multi-hop' mosaic data integration where some datasets do not share any features and enables the use of spatial gene expression features for mapping dissociated single-cell data onto a spatial transcriptomic reference.

Tracking early mammalian organogenesis - prediction and validation of differentiation trajectories at whole organism scale.

Early organogenesis represents a key step in animal development, during which pluripotent cells diversify to initiate organ formation. Here, we sampled 300,000 single-cell transcriptomes from mouse embryos between E8.5 and E9.5 in 6-h intervals and combined this new dataset with our previous atlas (E6.5-E8.5) to produce a densely sampled timecourse of >400,000 cells from early gastrulation to organogenesis. Computational lineage reconstruction identified complex waves of blood and endothelial development, including a new programme for somite-derived endothelium. We also dissected the E7.5 primitive streak into four adjacent regions, performed scRNA-seq and predicted cell fates computationally. Finally, we defined developmental state/fate relationships by combining orthotopic grafting, microscopic analysis and scRNA-seq to transcriptionally determine cell fates of grafted primitive streak regions after 24 h of in vitro embryo culture. Experimentally determined fate outcomes were in good agreement with computationally predicted fates, demonstrating how classical grafting experiments can be revisited to establish high-resolution cell state/fate relationships. Such interdisciplinary approaches will benefit future studies in developmental biology and guide the in vitro production of cells for organ regeneration and repair.

A human embryonic limb cell atlas resolved in space and time.

Human limbs emerge during the fourth post-conception week as mesenchymal buds, which develop into fully formed limbs over the subsequent months1. This process is orchestrated by numerous temporally and spatially restricted gene expression programmes, making congenital alterations in phenotype common2. Decades of work with model organisms have defined the fundamental mechanisms underlying vertebrate limb development, but an in-depth characterization of this process in humans has yet to be performed. Here we detail human embryonic limb development across space and time using single-cell and spatial transcriptomics. We demonstrate extensive diversification of cells from a few multipotent progenitors to myriad differentiated cell states, including several novel cell populations. We uncover two waves of human muscle development, each characterized by different cell states regulated by separate gene expression programmes, and identify musculin (MSC) as a key transcriptional repressor maintaining muscle stem cell identity. Through assembly of multiple anatomically continuous spatial transcriptomic samples using VisiumStitcher, we map cells across a sagittal section of a whole fetal hindlimb. We reveal a clear anatomical segregation between genes linked to brachydactyly and polysyndactyly, and uncover transcriptionally and spatially distinct populations of the mesenchyme in the autopod. Finally, we perform single-cell RNA sequencing on mouse embryonic limbs to facilitate cross-species developmental comparison, finding substantial homology between the two species.

Early human lung immune cell development and its role in epithelial cell fate.

Studies of human lung development have focused on epithelial and mesenchymal cell types and function, but much less is known about the developing lung immune cells, even though the airways are a major site of mucosal immunity after birth. An unanswered question is whether tissue-resident immune cells play a role in shaping the tissue as it develops in utero. Here, we profiled human embryonic and fetal lung immune cells using scRNA-seq, smFISH, and immunohistochemistry. At the embryonic stage, we observed an early wave of innate immune cells, including innate lymphoid cells, natural killer cells, myeloid cells, and lineage progenitors. By the canalicular stage, we detected naive T lymphocytes expressing high levels of cytotoxicity genes and the presence of mature B lymphocytes, including B-1 cells. Our analysis suggests that fetal lungs provide a niche for full B cell maturation. Given the presence and diversity of immune cells during development, we also investigated their possible effect on epithelial maturation. We found that IL-1ß drives epithelial progenitor exit from self-renewal and differentiation to basal cells in vitro. In vivo, IL-1ß-producing myeloid cells were found throughout the lung and adjacent to epithelial tips, suggesting that immune cells may direct human lung epithelial development.

A spatial human thymus cell atlas mapped to a continuous tissue axis.

T cells develop from circulating precursors, which enter the thymus and migrate throughout specialised sub-compartments to support maturation and selection. This process starts already in early fetal development and is highly active until the involution of the thymus in adolescence. To map the micro-anatomical underpinnings of this process in pre- vs. post-natal states, we undertook a spatially resolved analysis and established a new quantitative morphological framework for the thymus, the Cortico-Medullary Axis. Using this axis in conjunction with the curation of a multimodal single-cell, spatial transcriptomics and high-resolution multiplex imaging atlas, we show that canonical thymocyte trajectories and thymic epithelial cells are highly organised and fully established by post-conception week 12, pinpoint TEC progenitor states, find that TEC subsets and peripheral tissue genes are associated with Hassall's Corpuscles and uncover divergence in the pace and drivers of medullary entry between CD4 vs. CD8 T cell lineages. These findings are complemented with a holistic toolkit for spatial analysis and annotation, providing a basis for a detailed understanding of T lymphocyte development.

Cytotoxic T lymphocytes require transcription for infiltration but not target cell lysis.

Effector cytotoxic T lymphocytes (CTLs) are critical for ridding the body of infected or cancerous cells. CTL T cell receptor (TCR) ligation not only drives the delivery and release of cytolytic granules but also initiates a new wave of transcription. In order to address whether TCR-induced transcriptomic changes impact the ability of CTLs to kill, we asked which genes are expressed immediately after CTLs encounter targets and how CTL responses change when inhibiting transcription. Our data demonstrate that while expression of cytokines/chemokines and transcriptional machinery depend on transcription, cytotoxic protein expression and cytolytic activity are relatively robust to transcription blockade, with CTLs lysing nearby target cells for several hours after actinomycin D treatment. Monitoring CTL movement among target cells after inhibiting transcription demonstrates an infiltration defect that is not rectified by provision of exogenous cytokine/chemokine gradients, indicating a cell-intrinsic transcriptional requirement for infiltration. Together, our results reveal differential molecular control of CTL functions, separating recruitment and infiltration from cytolysis.

Precise identification of cell states altered in disease using healthy single-cell references.

Joint analysis of single-cell genomics data from diseased tissues and a healthy reference can reveal altered cell states. We investigate whether integrated collections of data from healthy individuals (cell atlases) are suitable references for disease-state identification and whether matched control samples are needed to minimize false discoveries. We demonstrate that using a reference atlas for latent space learning followed by differential analysis against matched controls leads to improved identification of disease-associated cells, especially with multiple perturbed cell types. Additionally, when an atlas is available, reducing control sample numbers does not increase false discovery rates. Jointly analyzing data from a COVID-19 cohort and a blood cell atlas, we improve detection of infection-related cell states linked to distinct clinical severities. Similarly, we studied disease states in pulmonary fibrosis using a healthy lung atlas, characterizing two distinct aberrant basal states. Our analysis provides guidelines for designing disease cohort studies and optimizing cell atlas use.

Time space and single-cell resolved tissue lineage trajectories and laterality of body plan at gastrulation.

Understanding of the molecular drivers of lineage diversification and tissue patterning during primary germ layer development requires in-depth knowledge of the dynamic molecular trajectories of cell lineages across a series of developmental stages of gastrulation. Through computational modeling, we constructed at single-cell resolution, a spatio-temporal transcriptome of cell populations in the germ-layers of gastrula-stage mouse embryos. This molecular atlas enables the inference of molecular network activity underpinning the specification and differentiation of the germ-layer tissue lineages. Heterogeneity analysis of cellular composition at defined positions in the epiblast revealed progressive diversification of cell types. The single-cell transcriptome revealed an enhanced BMP signaling activity in the right-side mesoderm of late-gastrulation embryo. Perturbation of asymmetric BMP signaling activity at late gastrulation led to randomization of left-right molecular asymmetry in the lateral mesoderm of early-somite-stage embryo. These findings indicate the asymmetric BMP activity during gastrulation may be critical for the symmetry breaking process.

Time-resolved single-cell RNAseq profiling identifies a novel Fabp5+ subpopulation of inflammatory myeloid cells with delayed cytotoxic profile in chronic spinal cord injury.

Traumatic spinal cord injuries (SCI) are a group of highly debilitating pathologies affecting thousands annually, and adversely affecting quality of life. Currently, no fully restorative therapies exist, and SCI still results in significant personal, societal and financial burdens. Inflammation plays a major role in the evolution of SCI, with myeloid cells, including bone marrow derived macrophages (BMDMs) and microglia (MG) being primary drivers of both early secondary pathogenesis and delayed wound healing events. The precise role of myeloid cell subsets is unclear as upon crossing the blood-spinal cord barrier, infiltrating bone marrow derived macrophages (BMDMs) may take on the morphology of resident microglia, and upregulate canonical microglia markers, thus making the two populations difficult to distinguish. Here, we used time-resolved scRNAseq and transgenic fate-mapping to chart the transcriptional profiles of tissue-resident and -infiltrating myeloid cells in a mouse model of thoracic contusion SCI. Our work identifies a novel subpopulation of foam cell-like inflammatory myeloid cells with increased expression of Fatty Acid Binding Protein 5 (Fabp5) and comprise both tissue-resident and -infiltrating cells. Fabp5+ inflammatory myeloid cells display a delayed cytotoxic profile that is predominant at the lesion epicentre and extends into the chronic phase of SCI.

Unsupervised removal of systematic background noise from droplet-based single-cell experiments using CellBender.

Droplet-based single-cell assays, including single-cell RNA sequencing (scRNA-seq), single-nucleus RNA sequencing (snRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq), generate considerable background noise counts, the hallmark of which is nonzero counts in cell-free droplets and off-target gene expression in unexpected cell types. Such systematic background noise can lead to batch effects and spurious differential gene expression results. Here we develop a deep generative model based on the phenomenology of noise generation in droplet-based assays. The proposed model accurately distinguishes cell-containing droplets from cell-free droplets, learns the background noise profile and provides noise-free quantification in an end-to-end fashion. We implement this approach in the scalable and robust open-source software package CellBender. Analysis of simulated data demonstrates that CellBender operates near the theoretically optimal denoising limit. Extensive evaluations using real datasets and experimental benchmarks highlight enhanced concordance between droplet-based single-cell data and established gene expression patterns, while the learned background noise profile provides evidence of degraded or uncaptured cell types.

A transcriptional and regulatory map of mouse somite maturation.

The mammalian body plan is shaped by rhythmic segmentation of mesoderm into somites, which are transient embryonic structures that form down each side of the neural tube. We have analyzed the genome-wide transcriptional and chromatin dynamics occurring within nascent somites, from early inception of somitogenesis to the latest stages of body plan establishment. We created matched gene expression and open chromatin maps for the three leading pairs of somites at six time points during mouse embryonic development. We show that the rate of somite differentiation accelerates as development progresses. We identified a conserved maturation program followed by all somites, but somites from more developed embryos concomitantly switch on differentiation programs from derivative cell lineages soon after segmentation. Integrated analysis of the somitic transcriptional and chromatin activities identified opposing regulatory modules controlling the onset of differentiation. Our results provide a powerful, high-resolution view of the molecular genetics underlying somitic development in mammals.

Mapping interindividual dynamics of innate immune response at single-cell resolution.

Common genetic variants across individuals modulate the cellular response to pathogens and are implicated in diverse immune pathologies, yet how they dynamically alter the response upon infection is not well understood. Here, we triggered antiviral responses in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-sequencing. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), a statistical approach designed to identify nonlinear dynamic genetic effects across transcriptional trajectories of cells. This approach identified 1,275 expression quantitative trait loci (local false discovery rate 10%) that manifested during the responses, many of which were colocalized with susceptibility loci identified by genome-wide association studies of infectious and autoimmune diseases, including the OAS1 splicing quantitative trait locus in a COVID-19 susceptibility locus. In summary, our analytical approach provides a unique framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

An atlas of rabbit development as a model for single-cell comparative genomics.

Traditionally, the mouse has been the favoured vertebrate model for biomedical research, due to its experimental and genetic tractability. However, non-rodent embryological studies highlight that many aspects of early mouse development, such as its egg-cylinder gastrulation and method of implantation, diverge from other mammals, thus complicating inferences about human development. Like the human embryo, rabbits develop as a flat-bilaminar disc. Here we constructed a morphological and molecular atlas of rabbit development. We report transcriptional and chromatin accessibility profiles for over 180,000 single cells and high-resolution histology sections from embryos spanning gastrulation, implantation, amniogenesis and early organogenesis. Using a neighbourhood comparison pipeline, we compare the transcriptional landscape of rabbit and mouse at the scale of the entire organism. We characterize the gene regulatory programmes underlying trophoblast differentiation and identify signalling interactions involving the yolk sac mesothelium during haematopoiesis. We demonstrate how the combination of both rabbit and mouse atlases can be leveraged to extract new biological insights from sparse macaque and human data. The datasets and computational pipelines reported here set a framework for a broader cross-species approach to decipher early mammalian development, and are readily adaptable to deploy single-cell comparative genomics more broadly across biomedical research.

Origin and segregation of the human germline.

Human germline-soma segregation occurs during weeks 2-3 in gastrulating embryos. Although direct studies are hindered, here, we investigate the dynamics of human primordial germ cell (PGCs) specification using in vitro models with temporally resolved single-cell transcriptomics and in-depth characterisation using in vivo datasets from human and nonhuman primates, including a 3D marmoset reference atlas. We elucidate the molecular signature for the transient gain of competence for germ cell fate during peri-implantation epiblast development. Furthermore, we show that both the PGCs and amnion arise from transcriptionally similar TFAP2A-positive progenitors at the posterior end of the embryo. Notably, genetic loss of function experiments shows that TFAP2A is crucial for initiating the PGC fate without detectably affecting the amnion and is subsequently replaced by TFAP2C as an essential component of the genetic network for PGC fate. Accordingly, amniotic cells continue to emerge from the progenitors in the posterior epiblast, but importantly, this is also a source of nascent PGCs.

Organoid modeling of human fetal lung alveolar development reveals mechanisms of cell fate patterning and neonatal respiratory disease.

Variation in lung alveolar development is strongly linked to disease susceptibility. However, underlying cellular and molecular mechanisms are difficult to study in humans. We have identified an alveolar-fated epithelial progenitor in human fetal lungs, which we grow as self-organizing organoids that model key aspects of cell lineage commitment. Using this system, we have functionally validated cell-cell interactions in the developing human alveolar niche, showing that Wnt signaling from differentiating fibroblasts promotes alveolar-type-2 cell identity, whereas myofibroblasts secrete the Wnt inhibitor, NOTUM, providing spatial patterning. We identify a Wnt-NKX2.1 axis controlling alveolar differentiation. Moreover, we show that differential binding of NKX2.1 coordinates alveolar maturation, allowing us to model the effects of human genetic variation in NKX2.1 on alveolar differentiation. Our organoid system recapitulates key aspects of human fetal lung stem cell biology allowing mechanistic experiments to determine the cellular and molecular regulation of human development and disease.