High-dimensional single-cell analysis of human natural killer cell heterogeneity.

Natural killer (NK) cells are innate lymphoid cells (ILCs) contributing to immune responses to microbes and tumors. Historically, their classification hinged on a limited array of surface protein markers. Here, we used single-cell RNA sequencing (scRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to dissect the heterogeneity of NK cells. We identified three prominent NK cell subsets in healthy human blood: NK1, NK2 and NK3, further differentiated into six distinct subgroups. Our findings delineate the molecular characteristics, key transcription factors, biological functions, metabolic traits and cytokine responses of each subgroup. These data also suggest two separate ontogenetic origins for NK cells, leading to divergent transcriptional trajectories. Furthermore, we analyzed the distribution of NK cell subsets in the lung, tonsils and intraepithelial lymphocytes isolated from healthy individuals and in 22 tumor types. This standardized terminology aims at fostering clarity and consistency in future research, thereby improving cross-study comparisons.

Multi-modal generative modeling for joint analysis of single-cell T cell receptor and gene expression data.

Recent advances in single-cell immune profiling have enabled the simultaneous measurement of transcriptome and T cell receptor (TCR) sequences, offering great potential for studying immune responses at the cellular level. However, integrating these diverse modalities across datasets is challenging due to their unique data characteristics and technical variations. Here, to address this, we develop the multimodal generative model mvTCR to fuse modality-specific information across transcriptome and TCR into a shared representation. Our analysis demonstrates the added value of multimodal over unimodal approaches to capture antigen specificity. Notably, we use mvTCR to distinguish T cell subpopulations binding to SARS-CoV-2 antigens from bystander cells. Furthermore, when combined with reference mapping approaches, mvTCR can map newly generated datasets to extensive T cell references, facilitating knowledge transfer. In summary, we envision mvTCR to enable a scalable analysis of multimodal immune profiling data and advance our understanding of immune responses.

Human SARS-CoV-2 challenge uncovers local and systemic response dynamics.

The COVID-19 pandemic is an ongoing global health threat, yet our understanding of the dynamics of early cellular responses to this disease remains limited1. Here in our SARS-CoV-2 human challenge study, we used single-cell multi-omics profiling of nasopharyngeal swabs and blood to temporally resolve abortive, transient and sustained infections in seronegative individuals challenged with pre-Alpha SARS-CoV-2. Our analyses revealed rapid changes in cell-type proportions and dozens of highly dynamic cellular response states in epithelial and immune cells associated with specific time points and infection status. We observed that the interferon response in blood preceded the nasopharyngeal response. Moreover, nasopharyngeal immune infiltration occurred early in samples from individuals with only transient infection and later in samples from individuals with sustained infection. High expression of HLA-DQA2 before inoculation was associated with preventing sustained infection. Ciliated cells showed multiple immune responses and were most permissive for viral replication, whereas nasopharyngeal T cells and macrophages were infected non-productively. We resolved 54 T cell states, including acutely activated T cells that clonally expanded while carrying convergent SARS-CoV-2 motifs. Our new computational pipeline Cell2TCR identifies activated antigen-responding T cells based on a gene expression signature and clusters these into clonotype groups and motifs. Overall, our detailed time series data can serve as a Rosetta stone for epithelial and immune cell responses and reveals early dynamic responses associated with protection against infection.

Characterising plasmacytoid and myeloid AXL+ SIGLEC-6+ dendritic cell functions and their interactions with HIV.

AXL+ Siglec-6+ dendritic cells (ASDC) are novel myeloid DCs which can be subdivided into CD11c+ and CD123+ expressing subsets. We showed for the first time that these two ASDC subsets are present in inflamed human anogenital tissues where HIV transmission occurs. Their presence in inflamed tissues was supported by single cell RNA analysis of public databases of such tissues including psoriasis diseased skin and colorectal cancer. Almost all previous studies have examined ASDCs as a combined population. Our data revealed that the two ASDC subsets differ markedly in their functions when compared with each other and to pDCs. Relative to their cell functions, both subsets of blood ASDCs but not pDCs expressed co-stimulatory and maturation markers which were more prevalent on CD11c+ ASDCs, thus inducing more T cell proliferation and activation than their CD123+ counterparts. There was also a significant polarisation of naïve T cells by both ASDC subsets toward Th2, Th9, Th22, Th17 and Treg but less toward a Th1 phenotype. Furthermore, we investigated the expression of chemokine receptors that facilitate ASDCs and pDCs migration from blood to inflamed tissues, their HIV binding receptors, and their interactions with HIV and CD4 T cells. For HIV infection, within 2 hours of HIV exposure, CD11c+ ASDCs showed a trend in more viral transfer to T cells than CD123+ ASDCs and pDCs for first phase transfer. However, for second phase transfer, CD123+ ASDCs showed a trend in transferring more HIV than CD11c+ ASDCs and there was no viral transfer from pDCs. As anogenital inflammation is a prerequisite for HIV transmission, strategies to inhibit ASDC recruitment into inflamed tissues and their ability to transmit HIV to CD4 T cells should be considered.

Inflammatory memory in psoriasis: From remission to recurrence.

The routine use of targeted systemic immunomodulatory therapies has transformed outcomes for people with severe psoriasis, with skin clearance (clinical remission) rates up to 60% at 1 year of biologic treatment. However, psoriasis may recur following drug withdrawal, and as a result, patients tend to continue receiving costly treatment indefinitely. Here, we review research into the "inflammatory memory" in resolved psoriasis skin and the potential mechanisms leading to psoriasis recurrence following drug withdrawal. Research has implicated immune cells such as tissue resident memory T cells, Langerhans cells, and dermal dendritic cells, and there is growing interest in keratinocytes and fibroblasts. A better understanding of the interactions between these cell populations, enabled by single cell technologies, will help to elucidate the events underpinning the shift from remission to recurrence. This may inform the development of personalized strategies for sustaining remission while reducing long-term drug burden.

Mapping Human Immunity and the Education of Waldeyer's Ring.

The development and deployment of single-cell genomic technologies have driven a resolution revolution in our understanding of the immune system, providing unprecedented insight into the diversity of immune cells present throughout the body and their function in health and disease. Waldeyer's ring is the collective name for the lymphoid tissue aggregations of the upper aerodigestive tract, comprising the palatine, pharyngeal (adenoids), lingual, and tubal tonsils. These tonsils are the first immune sentinels encountered by ingested and inhaled antigens and are responsible for mounting the first wave of adaptive immune response. An effective mucosal immune response is critical to neutralizing infection in the upper airway and preventing systemic spread, and dysfunctional immune responses can result in ear, nose, and throat pathologies. This review uses Waldeyer's ring to demonstrate how single-cell technologies are being applied to advance our understanding of the immune system and highlight directions for future research.

Targeted action to increase inclusion at the Wellcome Sanger Institute.

The Sanger Excellence Fellowship has been established to increase the representation of researchers with Black-heritage backgrounds at a leading research centre in the UK.

An atlas of cells in the human tonsil.

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Multiscale spatial mapping of cell populations across anatomical sites in healthy human skin and basal cell carcinoma.

Our understanding of how human skin cells differ according to anatomical site and tumour formation is limited. To address this, we have created a multiscale spatial atlas of healthy skin and basal cell carcinoma (BCC), incorporating in vivo optical coherence tomography, single-cell RNA sequencing, spatial global transcriptional profiling, and in situ sequencing. Computational spatial deconvolution and projection revealed the localisation of distinct cell populations to specific tissue contexts. Although cell populations were conserved between healthy anatomical sites and in BCC, mesenchymal cell populations including fibroblasts and pericytes retained signatures of developmental origin. Spatial profiling and in silico lineage tracing support a hair follicle origin for BCC and demonstrate that cancer-associated fibroblasts are an expansion of a POSTN+ subpopulation associated with hair follicles in healthy skin. RGS5+ pericytes are also expanded in BCC suggesting a role in vascular remodelling. We propose that the identity of mesenchymal cell populations is regulated by signals emanating from adjacent structures and that these signals are repurposed to promote the expansion of skin cancer stroma. The resource we have created is publicly available in an interactive format for the research community.

Transcriptional signals of transformation in human cancer.

BACKGROUND: As normal cells transform into cancers, their cell state changes, which may drive cancer cells into a stem-like or more primordial, foetal, or embryonic cell state. The transcriptomic profile of this final state may encode information about cancer's origin and how cancers relate to their normal cell counterparts. METHODS: Here, we used single-cell atlases to study cancer transformation in transcriptional terms. We utilised bulk transcriptomes across a wide spectrum of adult and childhood cancers, using a previously established method to interrogate their relationship to normal cell states. We extend and validate these findings using single-cell cancer transcriptomes and organ-specific atlases of colorectal and liver cancer. RESULTS: Our bulk transcriptomic data reveals that adult cancers rarely return to an embryonic state, but that a foetal state is a near-universal feature of childhood cancers. This finding was confirmed with single-cell cancer transcriptomes. CONCLUSIONS: Our findings provide a nuanced picture of transformation in human cancer, indicating cancer-specific rather than universal patterns of transformation pervade adult epithelial cancers.

Cross-Comparison of Inflammatory Skin Disease Transcriptomics Identifies PTEN as a Pathogenic Disease Classifier in Cutaneous Lupus Erythematosus.

Tissue transcriptomics is used to uncover molecular dysregulations underlying diseases. However, the majority of transcriptomics studies focus on single diseases with limited relevance for understanding the molecular relationship between diseases or for identifying disease-specific markers. In this study, we used a normalization approach to compare gene expression across nine inflammatory skin diseases. The normalized datasets were found to retain differential expression signals that allowed unsupervised disease clustering and identification of disease-specific gene signatures. Using the NS-Forest algorithm, we identified a minimal set of biomarkers and validated their use as diagnostic disease classifier. Among them, PTEN was identified as being a specific marker for cutaneous lupus erythematosus and found to be strongly expressed by lesional keratinocytes in association with pathogenic type I IFNs. In fact, PTEN facilitated the expression of IFN-ß and IFN-κ in keratinocytes by promoting activation and nuclear translocation of IRF3. Thus, cross-comparison of tissue transcriptomics is a valid strategy to establish a molecular disease classification and to identify pathogenic disease biomarkers.

Diverse macrophage populations contribute to distinct manifestations of human cutaneous graft-versus-host disease.

BACKGROUND: Graft-versus-host disease (GvHD) is a major life-threatening complication of allogeneic haematopoietic stem cell transplantation (HSCT), limiting the broad application of HSCT for haematological malignancies. Cutaneous GvHD is described as a post-transplant inflammatory reaction by skin-infiltrating donor T cells and remaining recipient tissue-resident memory T cells. Despite the major influence of lymphocytes on GvHD pathogenesis, the complex role of mononuclear phagocytes (MNPs) in tissues affected by GvHD is increasingly appreciated. OBJECTIVES: To characterize the identity, origin and functions of MNPs in patients with acute cutaneous GvHD. METHODS: Using single-cell RNA sequencing and multiplex tissue immunofluorescence, we identified an increased abundance of MNPs in skin and blood from 36 patients with acute cutaneous GvHD. In cases of sex-mismatched transplantation, we used expression of X-linked genes to detect rapid tissue adaptation of newly recruited donor MNPs resulting in similar transcriptional states of host- and donor-derived macrophages within GvHD skin lesions. RESULTS: We showed that cutaneous GvHD lesions harbour expanded CD163+ tissue-resident macrophage populations with anti-inflammatory and tissue-remodelling properties including interleukin-10 cytokine production. Cell-cell interaction analyses revealed putative signalling to strengthen regulatory T-cell responses. Notably, macrophage polarization in chronic cutaneous GvHD types was proinflammatory and drastically differed from acute GvHD, supporting the notion of distinct cellular players in different clinical GvHD subtypes. CONCLUSIONS: Overall, our data reveal a surprisingly dynamic role of MNPs after HSCT. Specific and time-resolved targeting to repolarize this cell subset may present a promising therapeutic strategy in combatting GvHD skin inflammation.

Dandelion uses the single-cell adaptive immune receptor repertoire to explore lymphocyte developmental origins.

Assessment of single-cell gene expression (single-cell RNA sequencing) and adaptive immune receptor (AIR) sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of nonproductive and partially spliced contigs. We devised a strategy to create an AIR feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double-positive T cells to mature single-positive CD4/CD8 T cells, generating predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is available at https://www.github.com/zktuong/dandelion .

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

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

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

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

Single-cell RNA sequencing of human epidermis identifies Lunatic fringe as a novel regulator of the stem cell compartment.

Single-cell RNA sequencing (scRNA-seq) of human skin provides a tool for validating observations from in vitro experimental models. By analyzing a published dataset of healthy adult epidermis, we confirm that the basal epidermal layer is heterogeneous, and three subpopulations of non-dividing cells can be distinguished. We show that Delta-like ligand 1 (DLL1) is expressed in a subset of basal cells previously identified as stem cells in cultured human keratinocytes and map the distribution of other Notch ligands and receptors to specific epidermal cell compartments. Although DLL1 is expressed at low levels, it is expressed in the same cell state as the Notch regulator, Lunatic -fringe (LFNG, O-fucosylpeptide 3-beta-N-acetylglucosaminyltransferase). Overexpression of LFNG amplifies the effects of DLL1 in cultured keratinocytes, increasing proliferation and colony-forming ability. We conclude that using scRNA-seq resources from healthy human skin not only validates previous experimental data but allows formulation of testable new hypotheses.

A Roadmap for a Consensus Human Skin Cell Atlas and Single-Cell Data Standardization.

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

Differentiation of IL-26+ TH17 intermediates into IL-17A producers via epithelial crosstalk in psoriasis.

Interleukin (IL)-26 is a TH17 cytokine with known antimicrobial and pro-inflammatory functions. However, the precise role of IL-26 in the context of pathogenic TH17 responses is unknown. Here we identify a population of blood TH17 intermediates that produce high levels of IL-26 and differentiate into IL-17A-producing TH17 cells upon TGF-ß1 exposure. By combining single cell RNA sequencing, TCR sequencing and spatial transcriptomics we show that this process occurs in psoriatic skin. In fact, IL-26+ TH17 intermediates infiltrating psoriatic skin induce TGF-ß1 expression in basal keratinocytes and thereby promote their own differentiation into IL-17A-producing cells. Thus, our study identifies IL-26-producing cells as an early differentiation stage of TH17 cells that infiltrates psoriatic skin and controls its own maturation into IL17A-producing TH17 cells, via epithelial crosstalk involving paracrine production of TGF-ß1.

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

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