Dendritic cells (DCs), monocytes, and macrophages are leukocytes with critical roles in immunity and tolerance. The DC network is evolutionarily conserved; the homologs of human tissue CD141(hi)XCR1⁺ CLEC9A⁺ DCs and CD1c⁺ DCs are murine CD103⁺ DCs and CD64⁻ CD11b⁺ DCs. In addition, human tissues also contain CD14⁺ cells, currently designated as DCs, with an as-yet unknown murine counterpart. Here we have demonstrated that human dermal CD14⁺ cells are a tissue-resident population of monocyte-derived macrophages with a short half-life of <6 days. The decline and reconstitution kinetics of human blood CD14⁺ monocytes and dermal CD14⁺ cells in vivo supported their precursor-progeny relationship. The murine homologs of human dermal CD14⁺ cells are CD11b⁺ CD64⁺ monocyte-derived macrophages. Human and mouse monocytes and macrophages were defined by highly conserved gene transcripts, which were distinct from DCs. The demonstration of monocyte-derived macrophages in the steady state in human tissue supports a conserved organization of human and mouse mononuclear phagocyte system.

The homozygous K108E mutation of interferon regulatory factor 8 (IRF8) is reported to cause dendritic cell (DC) and monocyte deficiency. However, more widespread immune dysfunction is predicted from the multiple roles ascribed to IRF8 in immune cell development and function.

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.

Rationale: Obesity affects 40% of U.S. adults, is associated with a proinflammatory state, and presents a significant risk factor for the development of severe coronavirus disease (COVID-19). To date, there is limited information on how obesity might affect immune cell responses in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Objectives: To determine the impact of obesity on respiratory tract immunity in COVID-19 across the human lifespan. Methods: We analyzed single-cell transcriptomes from BAL in three ventilated adult cohorts with (n = 24) or without (n = 9) COVID-19 from nasal immune cells in children with (n = 14) or without (n = 19) COVID-19, and from peripheral blood mononuclear cells in an independent adult COVID-19 cohort (n = 42), comparing obese and nonobese subjects. Measurements and Main Results: Surprisingly, we found that obese adult subjects had attenuated lung immune or inflammatory responses in SARS-CoV-2 infection, with decreased expression of IFN-α, IFN-γ, and TNF-α (tumor necrosis factor α) response gene signatures in almost all lung epithelial and immune cell subsets, and lower expression of IFNG and TNF in specific lung immune cells. Peripheral blood immune cells in an independent adult cohort showed a similar but less marked reduction in type-I IFN and IFNγ response genes, as well as decreased serum IFNα, in obese patients with SARS-CoV-2. Nasal immune cells from obese children with COVID-19 also showed reduced enrichment of IFN-α and IFN-γ response genes. Conclusions: These findings show blunted tissue immune responses in obese patients with COVID-19, with implications for treatment stratification, supporting the specific application of inhaled recombinant type-I IFNs in this vulnerable subset.

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.

Congenital or acquired cellular deficiencies in humans have the potential to reveal much about normal hematopoiesis and immune function. We show that a recently described syndrome of monocytopenia, B and NK lymphoid deficiency additionally includes the near absence of dendritic cells. Four subjects showed severe depletion of the peripheral blood HLA-DR(+) lineage(-) compartment, with virtually no CD123(+) or CD11c(+) dendritic cells (DCs) and very few CD14(+) or CD16(+) monocytes. The only remaining HLA-DR(+) lineage(-) cells were circulating CD34(+) progenitor cells. Dermal CD14(+) and CD1a(+) DC were also absent, consistent with their dependence on blood-derived precursors. In contrast, epidermal Langerhans cells and tissue macrophages were largely preserved. Combined loss of peripheral DCs, monocytes, and B and NK lymphocytes was mirrored in the bone marrow by complete absence of multilymphoid progenitors and depletion of granulocyte-macrophage progenitors. Depletion of the HLA-DR(+) peripheral blood compartment was associated with elevated serum fms-like tyrosine kinase ligand and reduced circulating CD4(+)CD25(hi)FoxP3(+) T cells, supporting a role for DC in T reg cell homeostasis.

As the first line of defence against pathogens, cells mount an innate immune response, which varies widely from cell to cell. The response must be potent but carefully controlled to avoid self-damage. How these constraints have shaped the evolution of innate immunity remains poorly understood. Here we characterize the innate immune response's transcriptional divergence between species and variability in expression among cells. Using bulk and single-cell transcriptomics in fibroblasts and mononuclear phagocytes from different species, challenged with immune stimuli, we map the architecture of the innate immune response. Transcriptionally diverging genes, including those that encode cytokines and chemokines, vary across cells and have distinct promoter structures. Conversely, genes that are involved in the regulation of this response, such as those that encode transcription factors and kinases, are conserved between species and display low cell-to-cell variability in expression. We suggest that this expression pattern, which is observed across species and conditions, has evolved as a mechanism for fine-tuned regulation to achieve an effective but balanced response.

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.

T cell pathology in the skin leads to monocyte influx, but we have little understanding of the fate of recruited cells within the diseased niche, or the long-term impact on cutaneous immune homeostasis. By combining a murine model of acute graft-versus-host disease (aGVHD) with analysis of patient samples, we demonstrate that pathology initiates dermis-specific macrophage differentiation and show that aGVHD-primed macrophages continue to dominate the dermal compartment at the relative expense of quiescent MHCIIint cells. Exposure of the altered dermal niche to topical haptens after disease resolution results in hyper-activation of regulatory T cells (Treg), but local breakdown in tolerance. Disease-imprinted macrophages express increased IL-1β and are predicted to elicit altered TNF superfamily interactions with cutaneous Treg, and we demonstrate the direct loss of T cell regulation within the resolved skin. Thus, T cell pathology leaves an immunological scar in the skin marked by failure to re-set immune homeostasis.

Conjunctival epithelial cells, which express viral-entry receptors angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine type 2 (TMPRSS2), constitute the largest exposed epithelium of the ocular surface tissue and may represent a relevant viral-entry route. To address this question, we generated an organotypic air-liquid-interface model of conjunctival epithelium, composed of basal, suprabasal, and superficial epithelial cells, and fibroblasts, which could be maintained successfully up to day 75 of differentiation. Using single-cell RNA sequencing (RNA-seq), with complementary imaging and virological assays, we observed that while all conjunctival cell types were permissive to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome expression, a productive infection did not ensue. The early innate immune response to SARS-CoV-2 infection in conjunctival cells was characterised by a robust autocrine and paracrine NF-κB activity, without activation of antiviral interferon signalling. Collectively, these data enrich our understanding of SARS-CoV-2 infection at the human ocular surface, with potential implications for the design of preventive strategies and conjunctival transplantation.

The extraembryonic yolk sac (YS) ensures delivery of nutritional support and oxygen to the developing embryo but remains ill-defined in humans. We therefore assembled a comprehensive multiomic reference of the human YS from 3 to 8 postconception weeks by integrating single-cell protein and gene expression data. Beyond its recognized role as a site of hematopoiesis, we highlight roles in metabolism, coagulation, vascular development, and hematopoietic regulation. We reconstructed the emergence and decline of YS hematopoietic stem and progenitor cells from hemogenic endothelium and revealed a YS-specific accelerated route to macrophage production that seeds developing organs. The multiorgan functions of the YS are superseded as intraembryonic organs develop, effecting a multifaceted relay of vital functions as pregnancy proceeds.

Non-lymphoid tissues (NLTs) harbor a pool of adaptive immune cells with largely unexplored phenotype and development. We used single-cell RNA-seq to characterize 35,000 CD4+ regulatory (Treg) and memory (Tmem) T cells in mouse skin and colon, their respective draining lymph nodes (LNs) and spleen. In these tissues, we identified Treg cell subpopulations with distinct degrees of NLT phenotype. Subpopulation pseudotime ordering and gene kinetics were consistent in recruitment to skin and colon, yet the initial NLT-priming in LNs and the final stages of NLT functional adaptation reflected tissue-specific differences. Predicted kinetics were recapitulated using an in vivo melanoma-induction model, validating key regulators and receptors. Finally, we profiled human blood and NLT Treg and Tmem cells, and identified cross-mammalian conserved tissue signatures. In summary, we describe the relationship between Treg cell heterogeneity and recruitment to NLTs through the combined use of computational prediction and in vivo validation.

Graft-versus-host disease (GVHD) is a life-threatening complication of allogeneic stem cell transplantation induced by the influx of donor-derived effector T cells (TE) into peripheral tissues. Current treatment strategies rely on targeting systemic T cells; however, the precise location and nature of instructions that program TE to become pathogenic and trigger injury are unknown. We therefore used weighted gene coexpression network analysis to construct an unbiased spatial map of TE differentiation during the evolution of GVHD and identified wide variation in effector programs in mice and humans according to location. Idiosyncrasy of effector programming in affected organs did not result from variation in T cell receptor repertoire or the selection of optimally activated TE. Instead, TE were reprogrammed by tissue-autonomous mechanisms in target organs for site-specific proinflammatory functions that were highly divergent from those primed in lymph nodes. In the skin, we combined the correlation-based network with a module-based differential expression analysis and showed that Langerhans cells provided in situ instructions for a Notch-dependent T cell gene cluster critical for triggering local injury. Thus, the principal determinant of TE pathogenicity in GVHD is the final destination, highlighting the need for target organ-specific approaches to block immunopathology while avoiding global immune suppression.

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.

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.

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.

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells and multipotent progenitors (HSC/MPPs) but remains poorly defined in humans. Here, using single-cell transcriptome profiling of approximately 140,000 liver and 74,000 skin, kidney and yolk sac cells, we identify the repertoire of human blood and immune cells during development. We infer differentiation trajectories from HSC/MPPs and evaluate the influence of the tissue microenvironment on blood and immune cell development. We reveal physiological erythropoiesis in fetal skin and the presence of mast cells, natural killer and innate lymphoid cell precursors in the yolk sac. We demonstrate a shift in the haemopoietic composition of fetal liver during gestation away from being predominantly erythroid, accompanied by a parallel change in differentiation potential of HSC/MPPs, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a reference for harnessing the therapeutic potential of HSC/MPPs.

The development of the human immune system during embryonic and fetal life has historically been difficult to research due to limited access to human tissue. Experimental animal models have been widely used to study development but cellular and molecular programmes may not be conserved across species. The advent of multiomic single-cell technologies and an increase in human developmental tissue biobank resources have facilitated single-cell multiomic studies focused on human immune development. A critical question in the near future is "How do we best reconcile scientific findings across multiple omic modalities, developmental time, and organismic space?" In this review, we discuss the application of single-cell multiomic technologies to unravel the major cellular lineages in the prenatal human immune system. We also identify key areas where the combined power of multiomics technologies can be leveraged to address specific immunological gaps in our current knowledge and explore new research horizons in human development.

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.

Rare hematopoietic stem and progenitor cell (HSPC) pools outside the bone marrow (BM) contribute to blood production in stress and disease but remain ill-defined. Although nonmobilized peripheral blood (PB) is routinely sampled for clinical management, the diagnosis and monitoring potential of PB HSPCs remain untapped, as no healthy PB HSPC baseline has been reported. Here we comprehensively delineate human extramedullary HSPC compartments comparing spleen, PB, and mobilized PB to BM using single-cell RNA-sequencing and/or functional assays. We uncovered HSPC features shared by extramedullary tissues and others unique to PB. First, in contrast to actively dividing BM HSPCs, we found no evidence of substantial ongoing hematopoiesis in extramedullary tissues at steady state but report increased splenic HSPC proliferative output during stress erythropoiesis. Second, extramedullary hematopoietic stem cells/multipotent progenitors (HSCs/MPPs) from spleen, PB, and mobilized PB share a common transcriptional signature and increased abundance of lineage-primed subsets compared with BM. Third, healthy PB HSPCs display a unique bias toward erythroid-megakaryocytic differentiation. At the HSC/MPP level, this is functionally imparted by a subset of phenotypic CD71+ HSCs/MPPs, exclusively producing erythrocytes and megakaryocytes, highly abundant in PB but rare in other adult tissues. Finally, the unique erythroid-megakaryocytic-skewing of PB is perturbed with age in essential thrombocythemia and β-thalassemia. Collectively, we identify extramedullary lineage-primed HSPC reservoirs that are nonproliferative in situ and report involvement of splenic HSPCs during demand-adapted hematopoiesis. Our data also establish aberrant composition and function of circulating HSPCs as potential clinical indicators of BM dysfunction.