Search

The liver has been studied extensively due to the broad number of diseases affecting its vital functions. However, therapeutic advances have been hampered by the lack of knowledge concerning human hepatic development. Here, we addressed this limitation by describing the developmental trajectories of different cell types that make up the human liver at single-cell resolution. These transcriptomic analyses revealed that sequential cell-to-cell interactions direct functional maturation of hepatocytes, with non-parenchymal cells playing essential roles during organogenesis. We utilized this information to derive bipotential hepatoblast organoids and then exploited this model system to validate the importance of signalling pathways in hepatocyte and cholangiocyte specification. Further insights into hepatic maturation also enabled the identification of stage-specific transcription factors to improve the functionality of hepatocyte-like cells generated from human pluripotent stem cells. Thus, our study establishes a platform to investigate the basic mechanisms directing human liver development and to produce cell types for clinical applications., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell-cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions., Competing Interests: Conflict of interest: F.J. Theis reports personal fees and nonfinancial support from Cellarity and Dermagnostix, during the conduct of the study. J. Schniering was supported by an ERS/EU RESPIRE4 Marie Skłlodowska-Curie Postdoctoral Fellowship. J. Lundeberg reports personal fees from 10x Genomics, outside the submitted work. P. Powell and J. Denning are employees of the European Lung Foundation. W. Timens reports personal fees from Merck Sharp Dohme and Bristol-Myers-Squibb, outside the submitted work. M. Nilsson reports personal fees from 10x Genomics, outside the submitted work. G.H. Koppelman reports grants from the Netherlands Lung Foundation, GlaxoSmithKline, Vertex, TEVA, UBBO EMMIUS Foundation and European Union (H2020 grant), outside the submitted work, and has participated in advisory board meetings to GlaxoSmithKline and PURE-IMS, outside the submitted work. M. van den Berge reports research grants paid to UMCG from GlaxoSmithKline, Genentech, Roche and Novartis, outside the submitted work. M.C. Nawijn reports grants from the European Commission, Chan Zuckerberg Initiative and Netherlands Lung Foundation, during the conduct of the study; grants from GlaxoSmithKline, outside the submitted work. All other authors do not have a conflict of interest., (Copyright ©The authors 2022.)

Pathogenic variants in genes that cause dilated cardiomyopathy (DCM) and arrhythmogenic cardiomyopathy (ACM) convey high risks for the development of heart failure through unknown mechanisms. Using single-nucleus RNA sequencing, we characterized the transcriptome of 880,000 nuclei from 18 control and 61 failing, nonischemic human hearts with pathogenic variants in DCM and ACM genes or idiopathic disease. We performed genotype-stratified analyses of the ventricular cell lineages and transcriptional states. The resultant DCM and ACM ventricular cell atlas demonstrated distinct right and left ventricular responses, highlighting genotype-associated pathways, intercellular interactions, and differential gene expression at single-cell resolution. Together, these data illuminate both shared and distinct cellular and molecular architectures of human heart failure and suggest candidate therapeutic targets.

Human cerebral cancers are known to contain cell types resembling the varying stages of neural development. However, the basis of this association remains unclear. Here, we map the development of mouse cerebrum across the developmental time-course, from embryonic day 12.5 to postnatal day 365, performing single-cell transcriptomics on >100,000 cells. By comparing this reference atlas to single-cell data from >100 glial tumours of the adult and paediatric human cerebrum, we find that tumour cells have an expression signature that overlaps with temporally restricted, embryonic radial glial precursors (RGPs) and their immediate sublineages. Further, we demonstrate that prenatal transformation of RGPs in a genetic mouse model gives rise to adult cerebral tumours that show an embryonic/juvenile RGP identity. Together, these findings implicate the acquisition of embryonic-like states in the genesis of adult glioma, providing insight into the origins of human glioma, and identifying specific developmental cell types for therapeutic targeting., (© 2022. The Author(s).)

B cells play a significant role in established Rheumatoid Arthritis (RA). However, it is unclear to what extent differentiated B cells are present in joint tissue already at the onset of disease. Here, we studied synovial biopsies (n = 8) captured from untreated patients at time of diagnosis. 3414 index-sorted B cells underwent RNA sequencing and paired tissue pieces were subjected to spatial transcriptomics (n = 4). We performed extensive bioinformatics analyses to dissect the local B cell composition. Select plasma cell immunoglobulin sequences were expressed as monoclonal antibodies and tested by ELISA. Memory and plasma cells were found irrespective of autoantibody status of the patients. Double negative memory B cells were prominent, but did not display a distinct transcriptional profile. The tissue architecture implicate both local B cell maturation via T cell help and plasma cell survival niches with a strong CXCL12-CXCR4 axis. The immunoglobulin sequence analyses revealed clonality between the memory B and plasma cell pools further supporting local maturation. One of the plasma cell-derived antibodies displayed citrulline autoreactivity, demonstrating local autoreactive plasma cell differentiation in joint biopsies captured from untreated early RA. Hence, plasma cell niches are not a consequence of chronic inflammation, but are already present at the time of diagnosis., (© 2022. The Author(s).)

The development of single-cell and spatial transcriptomics methods was instrumental in the conception of the Human Cell Atlas initiative, which aims to generate an integrated map of all cells across the human body. These technology advances are bringing increasing depth and resolution to maps of human organs and tissues, as well as our understanding of individual human cell types. Commonalities as well as tissue-specific features of primary and supportive cell types across human organs are beginning to emerge from these human tissue maps. In this Review, we highlight key biological insights obtained from cross-tissue studies into epithelial, fibroblast, vascular and immune cells based on single-cell gene expression data in humans and contrast it with mechanisms reported in mice., (© 2022. Springer Nature Limited.)

Single-cell genomics studies have decoded the immune cell composition of several human prenatal organs but were limited in describing the developing immune system as a distributed network across tissues. We profiled nine prenatal tissues combining single-cell RNA sequencing, antigen-receptor sequencing, and spatial transcriptomics to reconstruct the developing human immune system. This revealed the late acquisition of immune-effector functions by myeloid and lymphoid cell subsets and the maturation of monocytes and T cells before peripheral tissue seeding. Moreover, we uncovered system-wide blood and immune cell development beyond primary hematopoietic organs, characterized human prenatal B1 cells, and shed light on the origin of unconventional T cells. Our atlas provides both valuable data resources and biological insights that will facilitate cell engineering, regenerative medicine, and disease understanding.

γδ T cells are generally considered innate-like lymphocytes, however, an "adaptive-like" γδ compartment has now emerged. To understand transcriptional regulation of adaptive γδ T cell immunobiology, we combined single-cell transcriptomics, T cell receptor (TCR)-clonotype assignment, ATAC-seq, and immunophenotyping. We show that adult Vδ1 + T cells segregate into TCF7 + LEF1 + Granzyme B neg (T naive ) or T-bet + Eomes + BLIMP-1 + Granzyme B + (T effector ) transcriptional subtypes, with clonotypically expanded TCRs detected exclusively in T effector cells. Transcriptional reprogramming mirrors changes within CD8 + αβ T cells following antigen-specific maturation and involves chromatin remodeling, enhancing cytokine production and cytotoxicity. Consistent with this, in vitro TCR engagement induces comparable BLIMP-1, Eomes, and T-bet expression in naive Vδ1 + and CD8 + T cells. Finally, both human cytomegalovirus and Plasmodium falciparum infection in vivo drive adaptive Vδ1 T cell differentiation from T naive to T effector transcriptional status, alongside clonotypic expansion. Contrastingly, semi-invariant Vγ9 + Vδ2 + T cells exhibit a distinct "innate-effector" transcriptional program established by early childhood. In summary, adaptive-like γδ subsets undergo a pathogen-driven differentiation process analogous to conventional CD8 + T cells., Competing Interests: Declaration of interests M.J.T.S. has been employed by 10× Genomics since April 2018; this employment had no bearing on this work. The other authors declare no competing financial interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Despite their crucial role in health and disease, our knowledge of immune cells within human tissues remains limited. We surveyed the immune compartment of 16 tissues from 12 adult donors by single-cell RNA sequencing and VDJ sequencing generating a dataset of ~360,000 cells. To systematically resolve immune cell heterogeneity across tissues, we developed CellTypist, a machine learning tool for rapid and precise cell type annotation. Using this approach, combined with detailed curation, we determined the tissue distribution of finely phenotyped immune cell types, revealing hitherto unappreciated tissue-specific features and clonal architecture of T and B cells. Our multitissue approach lays the foundation for identifying highly resolved immune cell types by leveraging a common reference dataset, tissue-integrated expression analysis, and antigen receptor sequencing.

KMT2A-rearranged infant ALL is an aggressive childhood leukemia with poor prognosis. Here, we investigated the developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia (B-ALL) using bulk messenger RNA (mRNA) meta-analysis and examination of single lymphoblast transcriptomes against a developing bone marrow reference. KMT2A-rearranged infant B-ALL was uniquely dominated by an early lymphocyte precursor (ELP) state, whereas less adverse NUTM1-rearranged infant ALL demonstrated signals of later developing B cells, in line with most other childhood B-ALLs. We compared infant lymphoblasts with ELP cells and revealed that the cancer harbored hybrid myeloid-lymphoid features, including nonphysiological antigen combinations potentially targetable to achieve cancer specificity. We validated surface coexpression of exemplar combinations by flow cytometry. Through analysis of shared mutations in separate leukemias from a child with infant KMT2A-rearranged B-ALL relapsing as AML, we established that KMT2A rearrangement occurred in very early development, before hematopoietic specification, emphasizing that cell of origin cannot be inferred from the transcriptional state., (© 2022. The Author(s).)

The intestinal immune system represents the largest collection of immune cells in the body and is continually exposed to antigens from food and the microbiota. Here we discuss the contribution of single-cell transcriptomics in shaping our understanding of this complex system. We consider the impact on resolving early intestine development, engagement with the neighbouring microbiota, diversity of intestinal immune cells, compartmentalisation within the intestines and interactions with non-immune cells. Finally, we offer a perspective on open questions about gut immunity that evolving single-cell technologies are well placed to address., (© 2021. The Author(s).)

Effective clinical management of acute dengue virus (DENV) infection relies on the timing of suitable treatments during the disease progression. We analyzed single-cell transcriptomic profiles of the peripheral blood mononuclear cell samples from two DENV patients, collected daily during acute phase and also at convalescence. Key immune cell types demonstrated different dynamic responses over the course of the infection. On the day before defervescence (Day -1), we observed the peak expression of several prominent genes in the adaptive immunological pathways. We also characterized unique effector T cell clusters that expressed skin-homing signature genes at Day -1, whereas upregulation of skin and gut homing genes was also observed in plasma cells and plasmablasts during the febrile period. This work provides an overview of unique molecular dynamics that signify the entry of the critical phase, and the findings could improve the patient management of DENV infection., Competing Interests: In the past 3 years, Sarah A. Teichmann has consulted for Genentech and Roche and sits on scientific advisory boards for Qiagen, Foresite Labs, Biogen, and GlaxoSmithKline and is a co-founder and equity holder of Transition Bio. The other authors declare no competing interests., (© 2022 The Authors.)

Cutaneous T cell lymphoma (CTCL) represents a heterogeneous group of non-Hodgkin lymphoma distinguished by the presence of clonal malignant T cells. The heterogeneity of malignant T cells and the complex tumor microenvironment remain poorly characterized. With single-cell RNA analysis and bulk whole-exome sequencing on 19 skin lesions from 15 CTCL patients, we decipher the intra-tumor and inter-lesion diversity of CTCL patients and propose a multi-step tumor evolution model. We further establish a subtyping scheme based on the molecular features of malignant T cells and their pro-tumorigenic microenvironments: the T CyEM group, demonstrating a cytotoxic effector memory T cell phenotype, shows more M2 macrophages infiltration, while the T CM group, featured by a central memory T cell phenotype and adverse patient outcome, is infiltrated by highly exhausted CD8 + reactive T cells, B cells and Tregs with suppressive activities. Our results establish a solid basis for understanding the nature of CTCL and pave the way for future precision medicine for CTCL patients., (© 2022. The Author(s).)

The SARS-CoV-2 Omicron BA.1 variant emerged in 2021 1 and has multiple mutations in its spike protein 2 . Here we show that the spike protein of Omicron has a higher affinity for ACE2 compared with Delta, and a marked change in its antigenicity increases Omicron's evasion of therapeutic monoclonal and vaccine-elicited polyclonal neutralizing antibodies after two doses. mRNA vaccination as a third vaccine dose rescues and broadens neutralization. Importantly, the antiviral drugs remdesivir and molnupiravir retain efficacy against Omicron BA.1. Replication was similar for Omicron and Delta virus isolates in human nasal epithelial cultures. However, in lung cells and gut cells, Omicron demonstrated lower replication. Omicron spike protein was less efficiently cleaved compared with Delta. The differences in replication were mapped to the entry efficiency of the virus on the basis of spike-pseudotyped virus assays. The defect in entry of Omicron pseudotyped virus to specific cell types effectively correlated with higher cellular RNA expression of TMPRSS2, and deletion of TMPRSS2 affected Delta entry to a greater extent than Omicron. Furthermore, drug inhibitors targeting specific entry pathways 3 demonstrated that the Omicron spike inefficiently uses the cellular protease TMPRSS2, which promotes cell entry through plasma membrane fusion, with greater dependency on cell entry through the endocytic pathway. Consistent with suboptimal S1/S2 cleavage and inability to use TMPRSS2, syncytium formation by the Omicron spike was substantially impaired compared with the Delta spike. The less efficient spike cleavage of Omicron at S1/S2 is associated with a shift in cellular tropism away from TMPRSS2-expressing cells, with implications for altered pathogenesis., (© 2022. The Author(s).)

Current computational workflows for comparative analyses of single-cell datasets typically use discrete clusters as input when testing for differential abundance among experimental conditions. However, clusters do not always provide the appropriate resolution and cannot capture continuous trajectories. Here we present Milo, a scalable statistical framework that performs differential abundance testing by assigning cells to partially overlapping neighborhoods on a k-nearest neighbor graph. Using simulations and single-cell RNA sequencing (scRNA-seq) data, we show that Milo can identify perturbations that are obscured by discretizing cells into clusters, that it maintains false discovery rate control across batch effects and that it outperforms alternative differential abundance testing strategies. Milo identifies the decline of a fate-biased epithelial precursor in the aging mouse thymus and identifies perturbations to multiple lineages in human cirrhotic liver. As Milo is based on a cell-cell similarity structure, it might also be applicable to single-cell data other than scRNA-seq. Milo is provided as an open-source R software package at https://github.com/MarioniLab/miloR ., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

It is not fully understood why COVID-19 is typically milder in children 1-3 . Here, to examine the differences between children and adults in their response to SARS-CoV-2 infection, we analysed paediatric and adult patients with COVID-19 as well as healthy control individuals (total n = 93) using single-cell multi-omic profiling of matched nasal, tracheal, bronchial and blood samples. In the airways of healthy paediatric individuals, we observed cells that were already in an interferon-activated state, which after SARS-CoV-2 infection was further induced especially in airway immune cells. We postulate that higher paediatric innate interferon responses restrict viral replication and disease progression. The systemic response in children was characterized by increases in naive lymphocytes and a depletion of natural killer cells, whereas, in adults, cytotoxic T cells and interferon-stimulated subpopulations were significantly increased. We provide evidence that dendritic cells initiate interferon signalling in early infection, and identify epithelial cell states associated with COVID-19 and age. Our matching nasal and blood data show a strong interferon response in the airways with the induction of systemic interferon-stimulated populations, which were substantially reduced in paediatric patients. Together, we provide several mechanisms that explain the milder clinical syndrome observed in children., (© 2021. The Author(s).)

Multimodal data is rapidly growing in many fields of science and engineering, including single-cell biology. We introduce MultiMAP, a novel algorithm for dimensionality reduction and integration. MultiMAP can integrate any number of datasets, leverages features not present in all datasets, is not restricted to a linear mapping, allows the user to specify the influence of each dataset, and is extremely scalable to large datasets. We apply MultiMAP to single-cell transcriptomics, chromatin accessibility, methylation, and spatial data and show that it outperforms current approaches. On a new thymus dataset, we use MultiMAP to integrate cells along a temporal trajectory. This enables quantitative comparison of transcription factor expression and binding site accessibility over the course of T cell differentiation, revealing patterns of expression versus binding site opening kinetics., (© 2021. The Author(s).)

The endometrium, the mucosal lining of the uterus, undergoes dynamic changes throughout the menstrual cycle in response to ovarian hormones. We have generated dense single-cell and spatial reference maps of the human uterus and three-dimensional endometrial organoid cultures. We dissect the signaling pathways that determine cell fate of the epithelial lineages in the lumenal and glandular microenvironments. Our benchmark of the endometrial organoids reveals the pathways and cell states regulating differentiation of the secretory and ciliated lineages both in vivo and in vitro. In vitro downregulation of WNT or NOTCH pathways increases the differentiation efficiency along the secretory and ciliated lineages, respectively. We utilize our cellular maps to deconvolute bulk data from endometrial cancers and endometriotic lesions, illuminating the cell types dominating in each of these disorders. These mechanistic insights provide a platform for future development of treatments for common conditions including endometriosis and endometrial carcinoma., (© 2021. The Author(s).)

Massive single-cell profiling efforts have accelerated our discovery of the cellular composition of the human body while at the same time raising the need to formalize this new knowledge. Here, we discuss current efforts to harmonize and integrate different sources of annotations of cell types and states into a reference cell ontology. We illustrate with examples how a unified ontology can consolidate and advance our understanding of cell types across scientific communities and biological domains., (© 2021. Springer Nature Limited.)

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function 1-3 . Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi 4-6 . However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation 7 , are invisible with such approaches 8 . Here we developed immunoGAM, an extension of genome architecture mapping (GAM) 2,9 , to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation 2,10 . We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions., (© 2021. The Author(s).)

The Human Reference Atlas (HRA) aims to map all of the cells of the human body to advance biomedical research and clinical practice. This Perspective presents collaborative work by members of 16 international consortia on two essential and interlinked parts of the HRA: (1) three-dimensional representations of anatomy that are linked to (2) tables that name and interlink major anatomical structures, cell types, plus biomarkers (ASCT+B). We discuss four examples that demonstrate the practical utility of the HRA., (© 2021. Springer Nature Limited.)

Haematopoiesis in the bone marrow (BM) maintains blood and immune cell production throughout postnatal life. Haematopoiesis first emerges in human BM at 11-12 weeks after conception 1,2 , yet almost nothing is known about how fetal BM (FBM) evolves to meet the highly specialized needs of the fetus and newborn. Here we detail the development of FBM, including stroma, using multi-omic assessment of mRNA and multiplexed protein epitope expression. We find that the full blood and immune cell repertoire is established in FBM in a short time window of 6-7 weeks early in the second trimester. FBM promotes rapid and extensive diversification of myeloid cells, with granulocytes, eosinophils and dendritic cell subsets emerging for the first time. The substantial expansion of B lymphocytes in FBM contrasts with fetal liver at the same gestational age. Haematopoietic progenitors from fetal liver, FBM and cord blood exhibit transcriptional and functional differences that contribute to tissue-specific identity and cellular diversification. Endothelial cell types form distinct vascular structures that we show are regionally compartmentalized within FBM. Finally, we reveal selective disruption of B lymphocyte, erythroid and myeloid development owing to a cell-intrinsic differentiation bias as well as extrinsic regulation through an altered microenvironment in Down syndrome (trisomy 21)., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

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., (© 2021. The Author(s).)

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development., (© 2021. Springer Nature Limited.)

The origin of human metaplastic states and their propensity for cancer is poorly understood. Barrett's esophagus is a common metaplastic condition that increases the risk for esophageal adenocarcinoma, and its cellular origin is enigmatic. To address this, we harvested tissues spanning the gastroesophageal junction from healthy and diseased donors, including isolation of esophageal submucosal glands. A combination of single-cell transcriptomic profiling, in silico lineage tracing from methylation, open chromatin and somatic mutation analyses, and functional studies in organoid models showed that Barrett's esophagus originates from gastric cardia through c-MYC and HNF4A-driven transcriptional programs. Furthermore, our data indicate that esophageal adenocarcinoma likely arises from undifferentiated Barrett's esophagus cell types even in the absence of a pathologically identifiable metaplastic precursor, illuminating early detection strategies., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Background: Cellular diversity of the lung endothelium has not been systematically characterized in humans. We provide a reference atlas of human lung endothelial cells (ECs) to facilitate a better understanding of the phenotypic diversity and composition of cells comprising the lung endothelium., Methods: We reprocessed human control single-cell RNA sequencing (scRNAseq) data from 6 datasets. EC populations were characterized through iterative clustering with subsequent differential expression analysis. Marker genes were validated by fluorescent microscopy and in situ hybridization. scRNAseq of primary lung ECs cultured in vitro was performed. The signaling network between different lung cell types was studied. For cross-species analysis or disease relevance, we applied the same methods to scRNAseq data obtained from mouse lungs or from human lungs with pulmonary hypertension., Results: Six lung scRNAseq datasets were reanalyzed and annotated to identify >15 000 vascular EC cells from 73 individuals. Differential expression analysis of EC revealed signatures corresponding to endothelial lineage, including panendothelial, panvascular, and subpopulation-specific marker gene sets. Beyond the broad cellular categories of lymphatic, capillary, arterial, and venous ECs, we found previously indistinguishable subpopulations; among venous EC, we identified 2 previously indistinguishable populations: pulmonary-venous ECs (COL15A1 neg ) localized to the lung parenchyma and systemic-venous ECs (COL15A1 pos ) localized to the airways and the visceral pleura; among capillary ECs, we confirmed their subclassification into recently discovered aerocytes characterized by EDNRB , SOSTDC1 , and TBX2 and general capillary EC. We confirmed that all 6 endothelial cell types, including the systemic-venous ECs and aerocytes, are present in mice and identified endothelial marker genes conserved in humans and mice. Ligand-receptor connectome analysis revealed important homeostatic crosstalk of EC with other lung resident cell types. scRNAseq of commercially available primary lung ECs demonstrated a loss of their native lung phenotype in culture. scRNAseq revealed that endothelial diversity is maintained in pulmonary hypertension. Our article is accompanied by an online data mining tool (www.LungEndothelialCellAtlas.com)., Conclusions: Our integrated analysis provides a comprehensive and well-crafted reference atlas of ECs in the normal lung and confirms and describes in detail previously unrecognized endothelial populations across a large number of humans and mice.

Comprehensively characterizing the cellular composition and organization of tissues has been a long-term scientific challenge that has limited our ability to study fundamental and clinical aspects of human physiology. The Human Cell Atlas (HCA) is a global collaborative effort to create a reference map of all human cells as a basis for both understanding human health and diagnosing, monitoring, and treating disease. Many aspects of the HCA are analogous to the Human Genome Project (HGP), whose completion presents a major milestone in modern biology. To commemorate the HGP's 20-year anniversary of completion, we discuss the launch of the HCA in light of the HGP, and highlight recent progress by the HCA consortium., Competing Interests: Declaration of Interests In the last 3 years, S.A.T. has consulted for Genentech and Roche, and is a member of SABs of Biogen, GlaxoSmithKline, and Foresite Labs. A.R. is a cofounder and equity holder of Celsius Therapeutics, an equity holder of Immunitas and was an SAB member of Neogene Therapeutics, Thermo Fisher Scientific, Asimov, and Syros Pharmaceuticals until July 31, 2020. Since August 1, 2020, A.R. is an employee of Genentech, a member of the Roche group. A.R. is an inventor on multiple patents to the Broad Institute in the area of single cell genomics. R.G.H.L. has no interests to declare., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Tumor cells may share some patterns of gene expression with their cell of origin, providing clues into the differentiation state and origin of cancer. Here, we study the differentiation state and cellular origin of 1300 childhood and adult kidney tumors. Using single cell mRNA reference maps of normal tissues, we quantify reference "cellular signals" in each tumor. Quantifying global differentiation, we find that childhood tumors exhibit fetal cellular signals, replacing the presumption of "fetalness" with a quantitative measure of immaturity. By contrast, in adult cancers our assessment refutes the suggestion of dedifferentiation towards a fetal state in most cases. We find an intimate connection between developmental mesenchymal populations and childhood renal tumors. We demonstrate the diagnostic potential of our approach with a case study of a cryptic renal tumor. Our findings provide a cellular definition of human renal tumors through an approach that is broadly applicable to human cancer.

Here, we describe protocols for the preparation and dissociation of human fetal and pediatric intestinal tissue to a high-viability epithelial single-cell suspension. This epithelium-enriched single-cell suspension can then be used to generate single-cell RNA sequencing data as well as to create human intestinal organoids from both the fetal and pediatric intestine. Finally, this protocol details the dissociation of the intestinal organoids for use in single-cell analysis or passaging of organoids. For complete details on the use and execution of this protocol, please refer to Elmentaite et al. (2020)., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

Existing protocols for full-length single-cell RNA sequencing produce libraries of high complexity (thousands of distinct genes) with outstanding sensitivity and specificity of transcript quantification. These full-length libraries have the advantage of allowing probing of transcript isoforms, are informative regarding single-nucleotide polymorphisms and allow assembly of the VDJ region of the T- and B-cell-receptor sequences. Since full-length protocols are mostly plate-based at present, they are also suited to profiling cell types where cell numbers are limiting, such as rare cell types during development. A disadvantage of these methods has been the scalability and cost of the experiments, which has limited their popularity as compared with droplet-based and nanowell approaches. Here, we describe an automated protocol for full-length single-cell RNA sequencing, including both an in-house automated Smart-seq2 protocol and a commercial kit-based workflow. The protocols take 3-5 d to complete, depending on the number of plates processed in a batch. We discuss these two protocols in terms of ease of use, equipment requirements, running time, cost per sample and sequencing quality. By benchmarking the lysis buffers, reverse transcription enzymes and their combinations, we have optimized the in-house automated protocol to dramatically reduce its cost. An automated setup can be adopted easily by a competent researcher with basic laboratory skills and no prior automation experience. These pipelines have been employed successfully for several research projects allied with the Human Cell Atlas initiative ( www.humancellatlas.org ).

Germinal centers (GCs) are anatomic structures where B cells undergo affinity maturation, leading to production of high-affinity antibodies. The balance between T follicular helper (T FH ) and regulatory (T FR ) cells is critical for adequate control of GC responses. The study of human T FH and T FR cell development has been hampered because of the lack of in vitro assays reproducing in vivo biology, along with difficult access to healthy human lymphoid tissues. We used a single-cell transcriptomics approach to study the maturation of T FH and T FR cells isolated from human blood, iliac lymph nodes (LNs), and tonsils. As independent tissues have distinct proportions of follicular T cells in different maturation states, we leveraged the heterogeneity to reconstruct the maturation trajectory for human T FH and T FR cells. We found that the dominant maturation of T FR cells follows a bifurcated trajectory from precursor T reg cells, with one arm of the bifurcation leading to blood T FR cells and the other leading to the most mature GC T FR cells. Overall, our data provide a comprehensive resource for the transcriptomics of different follicular T cell populations and their dynamic relationship across different tissues., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The transcription factor Rora has been shown to be important for the development of ILC2 and the regulation of ILC3, macrophages and Treg cells. Here we investigate the role of Rora across CD4+ T cells in general, but with an emphasis on Th2 cells, both in vitro as well as in the context of several in vivo type 2 infection models. We dissect the function of Rora using overexpression and a CD4-conditional Rora-knockout mouse, as well as a RORA-reporter mouse. We establish the importance of Rora in CD4+ T cells for controlling lung inflammation induced by Nippostrongylus brasiliensis infection, and have measured the effect on downstream genes using RNA-seq. Using a systematic stimulation screen of CD4+ T cells, coupled with RNA-seq, we identify upstream regulators of Rora, most importantly IL-33 and CCL7. Our data suggest that Rora is a negative regulator of the immune system, possibly through several downstream pathways, and is under control of the local microenvironment., Competing Interests: The authors have read the journal’s policy and have the following competing interests: JF was an employee of Microsoft Research Cambridge at the time the study was conducted, but is no longer an employee of the company. EN is an employee of Aleph Labs, but was not employed by the company at the time the study was conducted. GK is an employee of AstraZeneca, but was not employed by the company at the time the study was conducted. This does not alter our adherence to PLOS ONE policies on sharing data and materials. There are no patents, products in development or marketed products associated with this research to declare.

Background: Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs., Objectives: To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis., Results: We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state., Conclusions: Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC., (© 2021 The Authors. Journal of Thrombosis and Haemostasis published by Wiley Periodicals LLC on behalf of International Society on Thrombosis and Haemostasis.)

Analysis of human blood immune cells provides insights into the coordinated response to viral infections such as severe acute respiratory syndrome coronavirus 2, which causes coronavirus disease 2019 (COVID-19). We performed single-cell transcriptome, surface proteome and T and B lymphocyte antigen receptor analyses of over 780,000 peripheral blood mononuclear cells from a cross-sectional cohort of 130 patients with varying severities of COVID-19. We identified expansion of nonclassical monocytes expressing complement transcripts (CD16 + C1QA/B/C + ) that sequester platelets and were predicted to replenish the alveolar macrophage pool in COVID-19. Early, uncommitted CD34 + hematopoietic stem/progenitor cells were primed toward megakaryopoiesis, accompanied by expanded megakaryocyte-committed progenitors and increased platelet activation. Clonally expanded CD8 + T cells and an increased ratio of CD8 + effector T cells to effector memory T cells characterized severe disease, while circulating follicular helper T cells accompanied mild disease. We observed a relative loss of IgA2 in symptomatic disease despite an overall expansion of plasmablasts and plasma cells. Our study highlights the coordinated immune response that contributes to COVID-19 pathogenesis and reveals discrete cellular components that can be targeted for therapy.

Despite signs of infection-including taste loss, dry mouth and mucosal lesions such as ulcerations, enanthema and macules-the involvement of the oral cavity in coronavirus disease 2019 (COVID-19) is poorly understood. To address this, we generated and analyzed two single-cell RNA sequencing datasets of the human minor salivary glands and gingiva (9 samples, 13,824 cells), identifying 50 cell clusters. Using integrated cell normalization and annotation, we classified 34 unique cell subpopulations between glands and gingiva. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral entry factors such as ACE2 and TMPRSS members were broadly enriched in epithelial cells of the glands and oral mucosae. Using orthogonal RNA and protein expression assessments, we confirmed SARS-CoV-2 infection in the glands and mucosae. Saliva from SARS-CoV-2-infected individuals harbored epithelial cells exhibiting ACE2 and TMPRSS expression and sustained SARS-CoV-2 infection. Acellular and cellular salivary fractions from asymptomatic individuals were found to transmit SARS-CoV-2 ex vivo. Matched nasopharyngeal and saliva samples displayed distinct viral shedding dynamics, and salivary viral burden correlated with COVID-19 symptoms, including taste loss. Upon recovery, this asymptomatic cohort exhibited sustained salivary IgG antibodies against SARS-CoV-2. Collectively, these data show that the oral cavity is an important site for SARS-CoV-2 infection and implicate saliva as a potential route of SARS-CoV-2 transmission.

Pregnenolone (P5) promotes prostate cancer cell growth, and de novo synthesis of intratumoural P5 is a potential cause of development of castration resistance. Immune cells can also synthesize P5 de novo . Despite its biological importance, little is known about P5's mode of actions, which appears to be context dependent and pleiotropic. A comprehensive proteome-wide spectrum of P5-binding proteins that are involved in its trafficking and functionality remains unknown. Here, we describe an approach that integrates chemical biology for probe synthesis with chemoproteomics to map P5-protein interactions in live prostate cancer cells and murine CD8 + T cells. We subsequently identified P5-binding proteins potentially involved in P5-trafficking and in P5's non-genomic action that may drive the promotion of castrate-resistance prostate cancer and regulate CD8 + T cell function. We envisage that this methodology could be employed for other steroids to map their interactomes directly in a broad range of living cells, tissues, and organisms., Competing Interests: The authors declare no competing financial interests., (© 2021 The Authors.)

Organoid technology holds great promise for regenerative medicine but has not yet been applied to humans. We address this challenge using cholangiocyte organoids in the context of cholangiopathies, which represent a key reason for liver transplantation. Using single-cell RNA sequencing, we show that primary human cholangiocytes display transcriptional diversity that is lost in organoid culture. However, cholangiocyte organoids remain plastic and resume their in vivo signatures when transplanted back in the biliary tree. We then utilize a model of cell engraftment in human livers undergoing ex vivo normothermic perfusion to demonstrate that this property allows extrahepatic organoids to repair human intrahepatic ducts after transplantation. Our results provide proof of principle that cholangiocyte organoids can be used to repair human biliary epithelium., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

We present INSIGHT [isothermal NASBA (nucleic acid sequence-based amplification) sequencing-based high-throughput test], a two-stage coronavirus disease 2019 testing strategy, using a barcoded isothermal NASBA reaction. It combines point-of-care diagnosis with next-generation sequencing, aiming to achieve population-scale testing. Stage 1 allows a quick decentralized readout for early isolation of presymptomatic or asymptomatic patients. It gives results within 1 to 2 hours, using either fluorescence detection or a lateral flow readout, while simultaneously incorporating sample-specific barcodes. The same reaction products from potentially hundreds of thousands of samples can then be pooled and used in a highly multiplexed sequencing-based assay in stage 2. This second stage confirms the near-patient testing results and facilitates centralized data collection. The 95% limit of detection is <50 copies of viral RNA per reaction. INSIGHT is suitable for further development into a rapid home-based, point-of-care assay and is potentially scalable to the population level., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

Protective humoral memory forms in secondary lymphoid organs where B cells undergo affinity maturation and differentiation into memory or plasma cells. Here, we provide a comprehensive roadmap of human B cell maturation with single-cell transcriptomics matched with bulk and single-cell antibody repertoires to define gene expression, antibody repertoires, and clonal sharing of B cell states at single-cell resolution, including memory B cell heterogeneity that reflects diverse functional and signaling states. We reconstruct gene expression dynamics during B cell activation to reveal a pre-germinal center state primed to undergo class switch recombination and dissect how antibody class-dependent gene expression in germinal center and memory B cells is linked with a distinct transcriptional wiring with potential to influence their fate and function. Our analyses reveal the dynamic cellular states that shape human B cell-mediated immunity and highlight how antibody isotype may play a role during their antibody-based selection., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Neuroblastoma is a childhood cancer that resembles developmental stages of the neural crest. It is not established what developmental processes neuroblastoma cancer cells represent. Here, we sought to reveal the phenotype of neuroblastoma cancer cells by comparing cancer ( n = 19,723) with normal fetal adrenal single-cell transcriptomes ( n = 57,972). Our principal finding was that the neuroblastoma cancer cell resembled fetal sympathoblasts, but no other fetal adrenal cell type. The sympathoblastic state was a universal feature of neuroblastoma cells, transcending cell cluster diversity, individual patients, and clinical phenotypes. We substantiated our findings in 650 neuroblastoma bulk transcriptomes and by integrating canonical features of the neuroblastoma genome with transcriptional signals. Overall, our observations indicate that a pan-neuroblastoma cancer cell state exists, which may be attractive for novel immunotherapeutic and targeted avenues., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

The skin confers biophysical and immunological protection through a complex cellular network established early in embryonic development. We profiled the transcriptomes of more than 500,000 single cells from developing human fetal skin, healthy adult skin, and adult skin with atopic dermatitis and psoriasis. We leveraged these datasets to compare cell states across development, homeostasis, and disease. Our analysis revealed an enrichment of innate immune cells in skin during the first trimester and clonal expansion of disease-associated lymphocytes in atopic dermatitis and psoriasis. We uncovered and validated in situ a reemergence of prenatal vascular endothelial cell and macrophage cellular programs in atopic dermatitis and psoriasis lesional skin. These data illustrate the dynamism of cutaneous immunity and provide opportunities for targeting pathological developmental programs in inflammatory skin diseases., (Copyright © 2021, American Association for the Advancement of Science.)

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., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

The thymus is a primary lymphoid organ, essential for T cell maturation and selection. There has been long-standing interest in processes underpinning thymus generation and the potential to manipulate it clinically, because alterations of thymus development or function can result in severe immunodeficiency and autoimmunity. Here, we identify epithelial-mesenchymal hybrid cells, capable of long-term expansion in vitro, and able to reconstitute an anatomic phenocopy of the native thymus, when combined with thymic interstitial cells and a natural decellularised extracellular matrix (ECM) obtained by whole thymus perfusion. This anatomical human thymus reconstruction is functional, as judged by its capacity to support mature T cell development in vivo after transplantation into humanised immunodeficient mice. These findings establish a basis for dissecting the cellular and molecular crosstalk between stroma, ECM and thymocytes, and offer practical prospects for treating congenital and acquired immunological diseases.

The dynamics of CD4 + T cell memory development remain to be examined at genome scale. In malaria-endemic regions, antimalarial chemoprevention protects long after its cessation and associates with effects on CD4 + T cells. We applied single-cell RNA sequencing and computational modelling to track memory development during Plasmodium infection and treatment. In the absence of central memory precursors, two trajectories developed as T helper 1 (T H 1) and follicular helper T (T FH ) transcriptomes contracted and partially coalesced over three weeks. Progeny of single clones populated T H 1 and T FH trajectories, and fate-mapping suggested that there was minimal lineage plasticity. Relationships between T FH and central memory were revealed, with antimalarials modulating these responses and boosting T H 1 recall. Finally, single-cell epigenomics confirmed that heterogeneity among effectors was partially reset in memory. Thus, the effector-to-memory transition in CD4 + T cells is gradual during malaria and is modulated by antiparasitic drugs. Graphical user interfaces are presented for examining gene-expression dynamics and gene-gene correlations ( http://haquelab.mdhs.unimelb.edu.au/cd4&#95;memory/ ).

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.