Pan-cancer profiling of tumor-infiltrating natural killer cells through transcriptional reference mapping.

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions.

Patient-derived organoid biobank identifies epigenetic dysregulation of intestinal epithelial MHC-I as a novel mechanism in severe Crohn's Disease.

OBJECTIVE: Epigenetic mechanisms, including DNA methylation (DNAm), have been proposed to play a key role in Crohn's disease (CD) pathogenesis. However, the specific cell types and pathways affected as well as their potential impact on disease phenotype and outcome remain unknown. We set out to investigate the role of intestinal epithelial DNAm in CD pathogenesis. DESIGN: We generated 312 intestinal epithelial organoids (IEOs) from mucosal biopsies of 168 patients with CD (n=72), UC (n=23) and healthy controls (n=73). We performed genome-wide molecular profiling including DNAm, bulk as well as single-cell RNA sequencing. Organoids were subjected to gene editing and the functional consequences of DNAm changes evaluated using an organoid-lymphocyte coculture and a nucleotide-binding oligomerisation domain, leucine-rich repeat and CARD domain containing 5 (NLRC5) dextran sulphate sodium (DSS) colitis knock-out mouse model. RESULTS: We identified highly stable, CD-associated loss of DNAm at major histocompatibility complex (MHC) class 1 loci including NLRC5 and cognate gene upregulation. Single-cell RNA sequencing of primary mucosal tissue and IEOs confirmed the role of NLRC5 as transcriptional transactivator in the intestinal epithelium. Increased mucosal MHC-I and NLRC5 expression in adult and paediatric patients with CD was validated in additional cohorts and the functional role of MHC-I highlighted by demonstrating a relative protection from DSS-mediated mucosal inflammation in NLRC5-deficient mice. MHC-I DNAm in IEOs showed a significant correlation with CD disease phenotype and outcomes. Application of machine learning approaches enabled the development of a disease prognostic epigenetic molecular signature. CONCLUSIONS: Our study has identified epigenetically regulated intestinal epithelial MHC-I as a novel mechanism in CD pathogenesis.

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.

Enhanced CD95 and interleukin 18 signalling accompany T cell receptor Vß21.3+ activation in multi-inflammatory syndrome in children.

Multisystem inflammatory syndrome in children is a post-infectious presentation SARS-CoV-2 associated with expansion of the T cell receptor Vß21.3+ T-cell subgroup. Here we apply muti-single cell omics to compare the inflammatory process in children with acute respiratory COVID-19 and those presenting with non SARS-CoV-2 infections in children. Here we show that in Multi-Inflammatory Syndrome in Children (MIS-C), the natural killer cell and monocyte population demonstrate heightened CD95 (Fas) and Interleuking 18 receptor expression. Additionally, TCR Vß21.3+ CD4+ T-cells exhibit skewed differentiation towards T helper 1, 17 and regulatory T cells, with increased expression of the co-stimulation receptors ICOS, CD28 and interleukin 18 receptor. We observe no functional evidence for NLRP3 inflammasome pathway overactivation, though MIS-C monocytes show elevated active caspase 8. This, coupled with raised IL18 mRNA expression in CD16- NK cells on single cell RNA sequencing analysis, suggests interleukin 18 and CD95 signalling may trigger activation of TCR Vß21.3+ T-cells in MIS-C, driven by increased IL-18 production from activated monocytes and CD16- Natural Killer cells.

Spatial genomics: mapping human steatotic liver disease.

Metabolic dysfunction-associated steatotic liver disease (MASLD, formerly known as non-alcoholic fatty liver disease) is a leading cause of chronic liver disease worldwide. MASLD can progress to metabolic dysfunction-associated steatohepatitis (MASH, formerly known as non-alcoholic steatohepatitis) with subsequent liver cirrhosis and hepatocellular carcinoma formation. The advent of current technologies such as single-cell and single-nuclei RNA sequencing have transformed our understanding of the liver in homeostasis and disease. The next frontier is contextualizing this single-cell information in its native spatial orientation. This understanding will markedly accelerate discovery science in hepatology, resulting in a further step-change in our knowledge of liver biology and pathobiology. In this Review, we discuss up-to-date knowledge of MASLD development and progression and how the burgeoning field of spatial genomics is driving exciting new developments in our understanding of human liver disease pathogenesis and therapeutic target identification.

Immune microniches shape intestinal Treg function.

The intestinal immune system is highly adapted to maintaining tolerance to the commensal microbiota and self-antigens while defending against invading pathogens1,2. Recognizing how the diverse network of local cells establish homeostasis and maintains it in the complex immune environment of the gut is critical to understanding how tolerance can be re-established following dysfunction, such as in inflammatory disorders. Although cell and molecular interactions that control T regulatory (Treg) cell development and function have been identified3,4, less is known about the cellular neighbourhoods and spatial compartmentalization that shapes microorganism-reactive Treg cell function. Here we used in vivo live imaging, photo-activation-guided single-cell RNA sequencing5-7 and spatial transcriptomics to follow the natural history of T cells that are reactive towards Helicobacter hepaticus through space and time in the settings of tolerance and inflammation. Although antigen stimulation can occur anywhere in the tissue, the lamina propria-but not embedded lymphoid aggregates-is the key microniche that supports effector Treg (eTreg) cell function. eTreg cells are stable once their niche is established; however, unleashing inflammation breaks down compartmentalization, leading to dominance of CD103+SIRPα+ dendritic cells in the lamina propria. We identify and validate the putative tolerogenic interaction between CD206+ macrophages and eTreg cells in the lamina propria and identify receptor-ligand pairs that are likely to govern the interaction. Our results reveal a spatial mechanism of tolerance in the lamina propria and demonstrate how knowledge of local interactions may contribute to the next generation of tolerance-inducing therapies.

Age-specific nasal epithelial responses to SARS-CoV-2 infection.

Children infected with SARS-CoV-2 rarely progress to respiratory failure. However, the risk of mortality in infected people over 85 years of age remains high. Here we investigate differences in the cellular landscape and function of paediatric (<12 years), adult (30-50 years) and older adult (>70 years) ex vivo cultured nasal epithelial cells in response to infection with SARS-CoV-2. We show that cell tropism of SARS-CoV-2, and expression of ACE2 and TMPRSS2 in nasal epithelial cell subtypes, differ between age groups. While ciliated cells are viral replication centres across all age groups, a distinct goblet inflammatory subtype emerges in infected paediatric cultures and shows high expression of interferon-stimulated genes and incomplete viral replication. In contrast, older adult cultures infected with SARS-CoV-2 show a proportional increase in basaloid-like cells, which facilitate viral spread and are associated with altered epithelial repair pathways. We confirm age-specific induction of these cell types by integrating data from in vivo COVID-19 studies and validate that our in vitro model recapitulates early epithelial responses to SARS-CoV-2 infection.

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.

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.

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 spatially resolved atlas of the human lung characterizes a gland-associated immune niche.

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

Obesity Is Associated with Attenuated Tissue Immunity in COVID-19.

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.

Impact of the Human Cell Atlas on medicine.

Single-cell atlases promise to provide a 'missing link' between genes, diseases and therapies. By identifying the specific cell types, states, programs and contexts where disease-implicated genes act, we will understand the mechanisms of disease at the cellular and tissue levels and can use this understanding to develop powerful disease diagnostics; identify promising new drug targets; predict their efficacy, toxicity and resistance mechanisms; and empower new kinds of therapies, from cancer therapies to regenerative medicine. Here, we lay out a vision for the potential of cell atlases to impact the future of medicine, and describe how advances over the past decade have begun to realize this potential in common complex diseases, infectious diseases (including COVID-19), rare diseases and cancer.

A human fetal lung cell atlas uncovers proximal-distal gradients of differentiation and key regulators of epithelial fates.

We present a multiomic cell atlas of human lung development that combines single-cell RNA and ATAC sequencing, high-throughput spatial transcriptomics, and single-cell imaging. Coupling single-cell methods with spatial analysis has allowed a comprehensive cellular survey of the epithelial, mesenchymal, endothelial, and erythrocyte/leukocyte compartments from 5-22 post-conception weeks. We identify previously uncharacterized cell states in all compartments. These include developmental-specific secretory progenitors and a subtype of neuroendocrine cell related to human small cell lung cancer. Our datasets are available through our web interface (https://lungcellatlas.org). To illustrate its general utility, we use our cell atlas to generate predictions about cell-cell signaling and transcription factor hierarchies which we rigorously test using organoid models.

Mapping single-cell transcriptomes in the intra-tumoral and associated territories of kidney cancer.

Tumor behavior is intricately dependent on the oncogenic properties of cancer cells and their multi-cellular interactions. To understand these dependencies within the wider microenvironment, we studied over 270,000 single-cell transcriptomes and 100 microdissected whole exomes from 12 patients with kidney tumors, prior to validation using spatial transcriptomics. Tissues were sampled from multiple regions of the tumor core, the tumor-normal interface, normal surrounding tissues, and peripheral blood. We find that the tissue-type location of CD8+ T cell clonotypes largely defines their exhaustion state with intra-tumoral spatial heterogeneity that is not well explained by somatic heterogeneity. De novo mutation calling from single-cell RNA-sequencing data allows us to broadly infer the clonality of stromal cells and lineage-trace myeloid cell development. We report six conserved meta-programs that distinguish tumor cell function, and find an epithelial-mesenchymal transition meta-program highly enriched at the tumor-normal interface that co-localizes with IL1B-expressing macrophages, offering a potential therapeutic target.

SOX9 maintains human foetal lung tip progenitor state by enhancing WNT and RTK signalling.

The balance between self-renewal and differentiation in human foetal lung epithelial progenitors controls the size and function of the adult organ. Moreover, progenitor cell gene regulation networks are employed by both regenerating and malignant lung cells, where modulators of their effects could potentially be of therapeutic value. Details of the molecular networks controlling human lung progenitor self-renewal remain unknown. We performed the first CRISPRi screen in primary human lung organoids to identify transcription factors controlling progenitor self-renewal. We show that SOX9 promotes proliferation of lung progenitors and inhibits precocious airway differentiation. Moreover, by identifying direct transcriptional targets using Targeted DamID, we place SOX9 at the centre of a transcriptional network, which amplifies WNT and RTK signalling to stabilise the progenitor cell state. In addition, the proof-of-principle CRISPRi screen and Targeted DamID tools establish a new workflow for using primary human organoids to elucidate detailed functional mechanisms underlying normal development and disease.

Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies.

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.

A brain precursor atlas reveals the acquisition of developmental-like states in adult cerebral tumours.

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.

Mapping the developing human immune system across organs.

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.

Single-cell transcriptomics links malignant T cells to the tumor immune landscape in cutaneous T cell lymphoma.

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 TCyEM group, demonstrating a cytotoxic effector memory T cell phenotype, shows more M2 macrophages infiltration, while the TCM 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.