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.

Spatial multiomics map of trophoblast development in early pregnancy.

The relationship between the human placenta-the extraembryonic organ made by the fetus, and the decidua-the mucosal layer of the uterus, is essential to nurture and protect the fetus during pregnancy. Extravillous trophoblast cells (EVTs) derived from placental villi infiltrate the decidua, transforming the maternal arteries into high-conductance vessels1. Defects in trophoblast invasion and arterial transformation established during early pregnancy underlie common pregnancy disorders such as pre-eclampsia2. Here we have generated a spatially resolved multiomics single-cell atlas of the entire human maternal-fetal interface including the myometrium, which enables us to resolve the full trajectory of trophoblast differentiation. We have used this cellular map to infer the possible transcription factors mediating EVT invasion and show that they are preserved in in vitro models of EVT differentiation from primary trophoblast organoids3,4 and trophoblast stem cells5. We define the transcriptomes of the final cell states of trophoblast invasion: placental bed giant cells (fused multinucleated EVTs) and endovascular EVTs (which form plugs inside the maternal arteries). We predict the cell-cell communication events contributing to trophoblast invasion and placental bed giant cell formation, and model the dual role of interstitial EVTs and endovascular EVTs in mediating arterial transformation during early pregnancy. Together, our data provide a comprehensive analysis of postimplantation trophoblast differentiation that can be used to inform the design of experimental models of the human placenta in early pregnancy.

Single-cell roadmap of human gonadal development.

Gonadal development is a complex process that involves sex determination followed by divergent maturation into either testes or ovaries1. Historically, limited tissue accessibility, a lack of reliable in vitro models and critical differences between humans and mice have hampered our knowledge of human gonadogenesis, despite its importance in gonadal conditions and infertility. Here, we generated a comprehensive map of first- and second-trimester human gonads using a combination of single-cell and spatial transcriptomics, chromatin accessibility assays and fluorescent microscopy. We extracted human-specific regulatory programmes that control the development of germline and somatic cell lineages by profiling equivalent developmental stages in mice. In both species, we define the somatic cell states present at the time of sex specification, including the bipotent early supporting population that, in males, upregulates the testis-determining factor SRY and sPAX8s, a gonadal lineage located at the gonadal-mesonephric interface. In females, we resolve the cellular and molecular events that give rise to the first and second waves of granulosa cells that compartmentalize the developing ovary to modulate germ cell differentiation. In males, we identify human SIGLEC15+ and TREM2+ fetal testicular macrophages, which signal to somatic cells outside and inside the developing testis cords, respectively. This study provides a comprehensive spatiotemporal map of human and mouse gonadal differentiation, which can guide in vitro gonadogenesis.

Spatial genomics maps the structure, nature and evolution of cancer clones.

Genome sequencing of cancers often reveals mosaics of different subclones present in the same tumour1-3. Although these are believed to arise according to the principles of somatic evolution, the exact spatial growth patterns and underlying mechanisms remain elusive4,5. Here, to address this need, we developed a workflow that generates detailed quantitative maps of genetic subclone composition across whole-tumour sections. These provide the basis for studying clonal growth patterns, and the histological characteristics, microanatomy and microenvironmental composition of each clone. The approach rests on whole-genome sequencing, followed by highly multiplexed base-specific in situ sequencing, single-cell resolved transcriptomics and dedicated algorithms to link these layers. Applying the base-specific in situ sequencing workflow to eight tissue sections from two multifocal primary breast cancers revealed intricate subclonal growth patterns that were validated by microdissection. In a case of ductal carcinoma in situ, polyclonal neoplastic expansions occurred at the macroscopic scale but segregated within microanatomical structures. Across the stages of ductal carcinoma in situ, invasive cancer and lymph node metastasis, subclone territories are shown to exhibit distinct transcriptional and histological features and cellular microenvironments. These results provide examples of the benefits afforded by spatial genomics for deciphering the mechanisms underlying cancer evolution and microenvironmental ecology.

Cells of the human intestinal tract mapped across space and time.

The cellular landscape of the human intestinal tract is dynamic throughout life, developing in utero and changing in response to functional requirements and environmental exposures. Here, to comprehensively map cell lineages, we use single-cell RNA sequencing and antigen receptor analysis of almost half a million cells from up to 5 anatomical regions in the developing and up to 11 distinct anatomical regions in the healthy paediatric and adult human gut. This reveals the existence of transcriptionally distinct BEST4 epithelial cells throughout the human intestinal tract. Furthermore, we implicate IgG sensing as a function of intestinal tuft cells. We describe neural cell populations in the developing enteric nervous system, and predict cell-type-specific expression of genes associated with Hirschsprung's disease. Finally, using a systems approach, we identify key cell players that drive the formation of secondary lymphoid tissue in early human development. We show that these programs are adopted in inflammatory bowel disease to recruit and retain immune cells at the site of inflammation. This catalogue of intestinal cells will provide new insights into cellular programs in development, homeostasis and disease.

Cells of the adult human heart.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

Combinatorial temporal patterning in progenitors expands neural diversity.

Human outer subventricular zone (OSVZ) neural progenitors and Drosophila type II neuroblasts both generate intermediate neural progenitors (INPs) that populate the adult cerebral cortex or central complex, respectively. It is unknown whether INPs simply expand or also diversify neural cell types. Here we show that Drosophila INPs sequentially generate distinct neural subtypes, that INPs sequentially express Dichaete, Grainy head and Eyeless transcription factors, and that these transcription factors are required for the production of distinct neural subtypes. Moreover, parental type II neuroblasts also sequentially express transcription factors and generate different neuronal/glial progeny over time, providing a second temporal identity axis. We conclude that neuroblast and INP temporal patterning axes act together to generate increased neural diversity within the adult central complex; OSVZ progenitors may use similar mechanisms to increase neural diversity in the human brain.

Generation of Mammalian Astrocyte Functional Heterogeneity.

Mammalian astrocytes have regional roles within the brain parenchyma. Indeed, the notion that astrocytes are molecularly heterogeneous could help explain how the central nervous system (CNS) retains embryonic positional information through development into specialized regions into adulthood. A growing body of evidence supports the concept of morphological and molecular differences between astrocytes in different brain regions, which might relate to their derivation from regionally patterned radial glia and/or local neuron inductive cues. Here, we review evidence for regionally encoded functions of astrocytes to provide an integrated concept on lineage origins and heterogeneity to understand regional brain organization, as well as emerging technologies to identify and further investigate novel roles for astrocytes.

Distinct lung cell signatures define the temporal evolution of diffuse alveolar damage in fatal COVID-19.

BACKGROUND: Lung damage in severe COVID-19 is highly heterogeneous however studies with dedicated spatial distinction of discrete temporal phases of diffuse alveolar damage (DAD) and alternate lung injury patterns are lacking. Existing studies have also not accounted for progressive airspace obliteration in cellularity estimates. We used an imaging mass cytometry (IMC) analysis with an airspace correction step to more accurately identify the cellular immune response that underpins the heterogeneity of severe COVID-19 lung disease. METHODS: Lung tissue was obtained at post-mortem from severe COVID-19 deaths. Pathologist-selected regions of interest (ROIs) were chosen by light microscopy representing the patho-evolutionary spectrum of DAD and alternate disease phenotypes were selected for comparison. Architecturally normal SARS-CoV-2-positive lung tissue and tissue from SARS-CoV-2-negative donors served as controls. ROIs were stained for 40 cellular protein markers and ablated using IMC before segmented cells were classified. Cell populations corrected by ROI airspace and their spatial relationships were compared across lung injury patterns. FINDINGS: Forty patients (32M:8F, age: 22-98), 345 ROIs and >900k single cells were analysed. DAD progression was marked by airspace obliteration and significant increases in mononuclear phagocytes (MnPs), T and B lymphocytes and significant decreases in alveolar epithelial and endothelial cells. Neutrophil populations proved stable overall although several interferon-responding subsets demonstrated expansion. Spatial analysis revealed immune cell interactions occur prior to microscopically appreciable tissue injury. INTERPRETATION: The immunopathogenesis of severe DAD in COVID-19 lung disease is characterised by sustained increases in MnPs and lymphocytes with key interactions occurring even prior to lung injury is established. FUNDING: UK Research and Innovation/Medical Research Council through the UK Coronavirus Immunology Consortium, Barbour Foundation, General Sir John Monash Foundation, Newcastle University, JGW Patterson Foundation, Wellcome Trust.

Human skeletal muscle aging atlas.

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

Site-specific binding of a PPR protein defines and stabilizes 5' and 3' mRNA termini in chloroplasts.

Chloroplast mRNA populations are characterized by overlapping transcripts derived by processing from polycistronic precursors. The mechanisms and functional significance of these processing events are poorly understood. We describe a pentatricopeptide repeat (PPR) protein, PPR10, whose binding defines mRNA segments derived from two transcription units in maize chloroplasts. PPR10 interacts in vivo and in vitro with two intergenic RNA regions of similar sequence. The processed 5' and 3' RNA termini in these regions overlap by approximately 25 nucleotides. The PPR10-binding sites map precisely to these overlapping sequences, and PPR10 is required specifically for the accumulation of RNAs with these termini. These findings show that PPR10 serves as a barrier to RNA decay from either the 5' or 3' direction and that a bound protein provides an alternative to an RNA hairpin as a barrier to 3' exonucleases. The results imply that protein 'caps' at both 5' and 3' ends can define the termini of chloroplast mRNA segments. These results, together with recent insights into bacterial RNA decay, suggest a unifying model for the biogenesis of chloroplast transcript populations and for the determinants of chloroplast mRNA stability.

Identification of a protein expression signature distinguishing early from organising diffuse alveolar damage in COVID-19 patients.

Diffuse alveolar damage (DAD) is the histological expression of acute respiratory distress syndrome and characterises lung pathology due to infection with SARS-CoV-2, and other respiratory pathogens of clinical significance. DAD reflects a time-dependent immunopathological process, progressing from an early/exudative stage through to an organising/fibrotic stage, yet within an individual these different stages of DAD may coexist. Understanding the progression of DAD is central to the development of new therapeutics to limit progressive lung damage. Here, we applied highly multiplexed spatial protein profiling to autopsy lung tissues derived from 27 patients who died from COVID-19 and identified a protein signature (ARG1, CD127, GZMB, IDO1, Ki67, phospho-PRAS40 (T246) and VISTA) that distinguishes early DAD from late DAD with good predictive accuracy. These proteins warrant further investigation as potential regulators of DAD progression.

A spatially resolved atlas of the human lung characterizes a gland-associated immune niche.

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

Cell2location maps fine-grained cell types in spatial transcriptomics.

Spatial transcriptomic technologies promise to resolve cellular wiring diagrams of tissues in health and disease, but comprehensive mapping of cell types in situ remains a challenge. Here we present сell2location, a Bayesian model that can resolve fine-grained cell types in spatial transcriptomic data and create comprehensive cellular maps of diverse tissues. Cell2location accounts for technical sources of variation and borrows statistical strength across locations, thereby enabling the integration of single-cell and spatial transcriptomics with higher sensitivity and resolution than existing tools. We assessed cell2location in three different tissues and show improved mapping of fine-grained cell types. In the mouse brain, we discovered fine regional astrocyte subtypes across the thalamus and hypothalamus. In the human lymph node, we spatially mapped a rare pre-germinal center B cell population. In the human gut, we resolved fine immune cell populations in lymphoid follicles. Collectively, our results present сell2location as a versatile analysis tool for mapping tissue architectures in a comprehensive manner.

Two CRM protein subfamilies cooperate in the splicing of group IIB introns in chloroplasts.

Chloroplast genomes in angiosperms encode approximately 20 group II introns, approximately half of which are classified as subgroup IIB. The splicing of all but one of the subgroup IIB introns requires a heterodimer containing the peptidyl-tRNA hydrolase homolog CRS2 and one of two closely related proteins, CAF1 or CAF2, that harbor a recently recognized RNA binding domain called the CRM domain. Two CRS2/CAF-dependent introns require, in addition, a CRM domain protein called CFM2 that is only distantly related to CAF1 and CAF2. Here, we show that CFM3, a close relative of CFM2, associates in vivo with those CRS2/CAF-dependent introns that are not CFM2 ligands. Mutant phenotypes in rice and Arabidopsis support a role for CFM3 in the splicing of most of the introns with which it associates. These results show that either CAF1 or CAF2 and either CFM2 or CFM3 simultaneously bind most chloroplast subgroup IIB introns in vivo, and that the CAF and CFM subunits play nonredundant roles in splicing. These results suggest that the expansion of the CRM protein family in plants resulted in two subfamilies that play different roles in group II intron splicing, with further diversification within a subfamily to accommodate multiple intron ligands.

Automation of Multiplexed RNAscope Single-Molecule Fluorescent In Situ Hybridization and Immunohistochemistry for Spatial Tissue Mapping.

In situ transcriptomic methods hold immense promise for spatially resolved mapping of cell types across human tissues. Here, we describe a protocol for automated single-molecule fluorescent in situ hybridization (smFISH) on standard histology sections at high throughput. We focus on the RNAscope smFISH assay that combines branched DNA amplification with tyramide signal amplification (TSA) to obtain high signal-to-noise ratio for tissue imaging. We describe the use of the robotic Leica BOND RX system for automation of liquid handling and the combination of the RNAscope assay with TSA-based immunohistochemistry without the need for specialized demultiplexed imaging.

Single-Cell Sequencing of Developing Human Gut Reveals Transcriptional Links to Childhood Crohn's Disease.

Human gut development requires the orchestrated interaction of differentiating cell types. Here, we generate an in-depth single-cell map of the developing human intestine at 6-10 weeks post-conception. Our analysis reveals the transcriptional profile of cycling epithelial precursor cells; distinct from LGR5-expressing cells. We propose that these cells may contribute to differentiated cell subsets via the generation of LGR5-expressing stem cells and receive signals from surrounding mesenchymal cells. Furthermore, we draw parallels between the transcriptomes of ex vivo tissues and in vitro fetal organoids, revealing the maturation of organoid cultures in a dish. Lastly, we compare scRNA-seq profiles from pediatric Crohn's disease epithelium alongside matched healthy controls to reveal disease-associated changes in the epithelial composition. Contrasting these with the fetal profiles reveals the re-activation of fetal transcription factors in Crohn's disease. Our study provides a resource available at www.gutcellatlas.org, and underscores the importance of unraveling fetal development in understanding disease.

Astrocyte layers in the mammalian cerebral cortex revealed by a single-cell in situ transcriptomic map.

Although the cerebral cortex is organized into six excitatory neuronal layers, it is unclear whether glial cells show distinct layering. In the present study, we developed a high-content pipeline, the large-area spatial transcriptomic (LaST) map, which can quantify single-cell gene expression in situ. Screening 46 candidate genes for astrocyte diversity across the mouse cortex, we identified superficial, mid and deep astrocyte identities in gradient layer patterns that were distinct from those of neurons. Astrocyte layer features, established in the early postnatal cortex, mostly persisted in adult mouse and human cortex. Single-cell RNA sequencing and spatial reconstruction analysis further confirmed the presence of astrocyte layers in the adult cortex. Satb2 and Reeler mutations that shifted neuronal post-mitotic development were sufficient to alter glial layering, indicating an instructive role for neuronal cues. Finally, astrocyte layer patterns diverged between mouse cortical regions. These findings indicate that excitatory neurons and astrocytes are organized into distinct lineage-associated laminae.

Aging Neural Progenitors Lose Competence to Respond to Mitogenic Notch Signaling.

Drosophila neural stem cells (neuroblasts) are a powerful model system for investigating stem cell self-renewal, specification of temporal identity, and progressive restriction in competence. Notch signaling is a conserved cue that is an important determinant of cell fate in many contexts across animal development; for example, mammalian T cell differentiation in the thymus and neuroblast specification in Drosophila are both regulated by Notch signaling. However, Notch also functions as a mitogen, and constitutive Notch signaling potentiates T cell leukemia as well as Drosophila neuroblast tumors. While the role of Notch signaling has been studied in these and other cell types, it remains unclear how stem cells and progenitors change competence to respond to Notch over time. Notch is required in type II neuroblasts for normal development of their transit amplifying progeny, intermediate neural progenitors (INPs). Here, we find that aging INPs lose competence to respond to constitutively active Notch signaling. Moreover, we show that reducing the levels of the old INP temporal transcription factor Eyeless/Pax6 allows Notch signaling to promote the de-differentiation of INP progeny into ectopic INPs, thereby creating a proliferative mass of ectopic progenitors in the brain. These findings provide a new system for studying progenitor competence and identify a novel role for the conserved transcription factor Eyeless/Pax6 in blocking Notch signaling during development.