The Myeloid Cell Compartment-Cell by Cell.

Myeloid cells are a major cellular compartment of the immune system comprising monocytes, dendritic cells, tissue macrophages, and granulocytes. Models of cellular ontogeny, activation, differentiation, and tissue-specific functions of myeloid cells have been revisited during the last years with surprising results; for example, most tissue macrophages are yolk sac derived, monocytes and macrophages follow a multidimensional model of activation, and tissue signals have a significant impact on the functionality of all these cells. While these exciting results have brought these cells back to center stage, their enormous plasticity and heterogeneity, during both homeostasis and disease, are far from understood. At the same time, the ongoing revolution in single-cell genomics, with single-cell RNA sequencing (scRNA-seq) leading the way, promises to change this. Prevailing models of hematopoiesis with distinct intermediates are challenged by scRNA-seq data suggesting more continuous developmental trajectories in the myeloid cell compartment. Cell subset structures previously defined by protein marker expression need to be revised based on unbiased analyses of scRNA-seq data. Particularly in inflammatory conditions, myeloid cells exhibit substantially vaster heterogeneity than previously anticipated, and work performed within large international projects, such as the Human Cell Atlas, has already revealed novel tissue macrophage subsets. Based on these exciting developments, we propose the next steps to a full understanding of the myeloid cell compartment in health and diseases.

Systems Immunology: Learning the Rules of the Immune System.

Given the many cell types and molecular components of the human immune system, along with vast variations across individuals, how should we go about developing causal and predictive explanations of immunity? A central strategy in human studies is to leverage natural variation to find relationships among variables, including DNA variants, epigenetic states, immune phenotypes, clinical descriptors, and others. Here, we focus on how natural variation is used to find patterns, infer principles, and develop predictive models for two areas: (a) immune cell activation-how single-cell profiling boosts our ability to discover immune cell types and states-and (b) antigen presentation and recognition-how models can be generated to predict presentation of antigens on MHC molecules and their detection by T cell receptors. These are two examples of a shift in how we find the drivers and targets of immunity, especially in the human system in the context of health and disease.

Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes.

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

Neutrophil profiling illuminates anti-tumor antigen-presenting potency.

Neutrophils, the most abundant and efficient defenders against pathogens, exert opposing functions across cancer types. However, given their short half-life, it remains challenging to explore how neutrophils adopt specific fates in cancer. Here, we generated and integrated single-cell neutrophil transcriptomes from 17 cancer types (225 samples from 143 patients). Neutrophils exhibited extraordinary complexity, with 10 distinct states including inflammation, angiogenesis, and antigen presentation. Notably, the antigen-presenting program was associated with favorable survival in most cancers and could be evoked by leucine metabolism and subsequent histone H3K27ac modification. These neutrophils could further invoke both (neo)antigen-specific and antigen-independent T cell responses. Neutrophil delivery or a leucine diet fine-tuned the immune balance to enhance anti-PD-1 therapy in various murine cancer models. In summary, these data not only indicate the neutrophil divergence across cancers but also suggest therapeutic opportunities such as antigen-presenting neutrophil delivery.

A multimodal zebrafish developmental atlas reveals the state-transition dynamics of late-vertebrate pluripotent axial progenitors.

Elucidating organismal developmental processes requires a comprehensive understanding of cellular lineages in the spatial, temporal, and molecular domains. In this study, we introduce Zebrahub, a dynamic atlas of zebrafish embryonic development that integrates single-cell sequencing time course data with lineage reconstructions facilitated by light-sheet microscopy. This atlas offers high-resolution and in-depth molecular insights into zebrafish development, achieved through the sequencing of individual embryos across ten developmental stages, complemented by reconstructions of cellular trajectories. Zebrahub also incorporates an interactive tool to navigate the complex cellular flows and lineages derived from light-sheet microscopy data, enabling in silico fate-mapping experiments. To demonstrate the versatility of our multimodal resource, we utilize Zebrahub to provide fresh insights into the pluripotency of neuro-mesodermal progenitors (NMPs) and the origins of a joint kidney-hemangioblast progenitor population.

Molecular and spatial signatures of mouse brain aging at single-cell resolution.

The diversity and complex organization of cells in the brain have hindered systematic characterization of age-related changes in its cellular and molecular architecture, limiting our ability to understand the mechanisms underlying its functional decline during aging. Here, we generated a high-resolution cell atlas of brain aging within the frontal cortex and striatum using spatially resolved single-cell transcriptomics and quantified changes in gene expression and spatial organization of major cell types in these regions over the mouse lifespan. We observed substantially more pronounced changes in cell state, gene expression, and spatial organization of non-neuronal cells over neurons. Our data revealed molecular and spatial signatures of glial and immune cell activation during aging, particularly enriched in the subcortical white matter, and identified both similarities and notable differences in cell-activation patterns induced by aging and systemic inflammatory challenge. These results provide critical insights into age-related decline and inflammation in the brain.

Tissue-resident memory and circulating T cells are early responders to pre-surgical cancer immunotherapy.

Neoadjuvant immune checkpoint blockade has shown promising clinical activity. Here, we characterized early kinetics in tumor-infiltrating and circulating immune cells in oral cancer patients treated with neoadjuvant anti-PD-1 or anti-PD-1/CTLA-4 in a clinical trial (NCT02919683). Tumor-infiltrating CD8 T cells that clonally expanded during immunotherapy expressed elevated tissue-resident memory and cytotoxicity programs, which were already active prior to therapy, supporting the capacity for rapid response. Systematic target discovery revealed that treatment-expanded tumor T cell clones in responding patients recognized several self-antigens, including the cancer-specific antigen MAGEA1. Treatment also induced a systemic immune response characterized by expansion of activated T cells enriched for tumor-infiltrating T cell clonotypes, including both pre-existing and emergent clonotypes undetectable prior to therapy. The frequency of activated blood CD8 T cells, notably pre-treatment PD-1-positive KLRG1-negative T cells, was strongly associated with intra-tumoral pathological response. These results demonstrate how neoadjuvant checkpoint blockade induces local and systemic tumor immunity.

Single-cell approaches in human microbiome research.

Microbial culturing and meta-omic profiling technologies have significantly advanced our understanding of the taxonomic and functional variation of the human microbiome and its impact on host processes. The next increase in resolution will come by understanding the role of low-abundant and less-prevalent bacteria and the study of individual cell behaviors that underlie the complexity of microbial ecosystems. To this aim, single-cell techniques are being rapidly developed to isolate, culture, and characterize the genomes and transcriptomes of individual microbes in complex communities. Here, we discuss how these single-cell technologies are providing unique insights into the biology and behavior of human microbiomes.

Tissue-resident FOLR2+ macrophages associate with CD8+ T cell infiltration in human breast cancer.

Macrophage infiltration is a hallmark of solid cancers, and overall macrophage infiltration correlates with lower patient survival and resistance to therapy. Tumor-associated macrophages, however, are phenotypically and functionally heterogeneous. Specific subsets of tumor-associated macrophage might be endowed with distinct roles on cancer progression and antitumor immunity. Here, we identify a discrete population of FOLR2+ tissue-resident macrophages in healthy mammary gland and breast cancer primary tumors. FOLR2+ macrophages localize in perivascular areas in the tumor stroma, where they interact with CD8+ T cells. FOLR2+ macrophages efficiently prime effector CD8+ T cells ex vivo. The density of FOLR2+ macrophages in tumors positively correlates with better patient survival. This study highlights specific roles for tumor-associated macrophage subsets and paves the way for subset-targeted therapeutic interventions in macrophages-based cancer therapies.

Quiescent cancer cells resist T cell attack by forming an immunosuppressive niche.

Immunotherapy is a promising treatment for triple-negative breast cancer (TNBC), but patients relapse, highlighting the need to understand the mechanisms of resistance. We discovered that in primary breast cancer, tumor cells that resist T cell attack are quiescent. Quiescent cancer cells (QCCs) form clusters with reduced immune infiltration. They also display superior tumorigenic capacity and higher expression of chemotherapy resistance and stemness genes. We adapted single-cell RNA-sequencing with precise spatial resolution to profile infiltrating cells inside and outside the QCC niche. This transcriptomic analysis revealed hypoxia-induced programs and identified more exhausted T cells, tumor-protective fibroblasts, and dysfunctional dendritic cells inside clusters of QCCs. This uncovered differential phenotypes in infiltrating cells based on their intra-tumor location. Thus, QCCs constitute immunotherapy-resistant reservoirs by orchestrating a local hypoxic immune-suppressive milieu that blocks T cell function. Eliminating QCCs holds the promise to counteract immunotherapy resistance and prevent disease recurrence in TNBC.

Mapping information-rich genotype-phenotype landscapes with genome-scale Perturb-seq.

A central goal of genetics is to define the relationships between genotypes and phenotypes. High-content phenotypic screens such as Perturb-seq (CRISPR-based screens with single-cell RNA-sequencing readouts) enable massively parallel functional genomic mapping but, to date, have been used at limited scales. Here, we perform genome-scale Perturb-seq targeting all expressed genes with CRISPR interference (CRISPRi) across >2.5 million human cells. We use transcriptional phenotypes to predict the function of poorly characterized genes, uncovering new regulators of ribosome biogenesis (including CCDC86, ZNF236, and SPATA5L1), transcription (C7orf26), and mitochondrial respiration (TMEM242). In addition to assigning gene function, single-cell transcriptional phenotypes allow for in-depth dissection of complex cellular phenomena-from RNA processing to differentiation. We leverage this ability to systematically identify genetic drivers and consequences of aneuploidy and to discover an unanticipated layer of stress-specific regulation of the mitochondrial genome. Our information-rich genotype-phenotype map reveals a multidimensional portrait of gene and cellular function.

A cellular and spatial map of the choroid plexus across brain ventricles and ages.

The choroid plexus (ChP) in each brain ventricle produces cerebrospinal fluid (CSF) and forms the blood-CSF barrier. Here, we construct a single-cell and spatial atlas of each ChP in the developing, adult, and aged mouse brain. We delineate diverse cell types, subtypes, cell states, and expression programs in epithelial and mesenchymal cells across ages and ventricles. In the developing ChP, we predict a common progenitor pool for epithelial and neuronal cells, validated by lineage tracing. Epithelial and fibroblast cells show regionalized expression by ventricle, starting at embryonic stages and persisting with age, with a dramatic transcriptional shift with maturation, and a smaller shift in each aged cell type. With aging, epithelial cells upregulate host-defense programs, and resident macrophages upregulate interleukin-1ß (IL-1ß) signaling genes. Our atlas reveals cellular diversity, architecture and signaling across ventricles during development, maturation, and aging of the ChP-brain barrier.

A taxonomy of transcriptomic cell types across the isocortex and hippocampal formation.

The isocortex and hippocampal formation (HPF) in the mammalian brain play critical roles in perception, cognition, emotion, and learning. We profiled ∼1.3 million cells covering the entire adult mouse isocortex and HPF and derived a transcriptomic cell-type taxonomy revealing a comprehensive repertoire of glutamatergic and GABAergic neuron types. Contrary to the traditional view of HPF as having a simpler cellular organization, we discover a complete set of glutamatergic types in HPF homologous to all major subclasses found in the six-layered isocortex, suggesting that HPF and the isocortex share a common circuit organization. We also identify large-scale continuous and graded variations of cell types along isocortical depth, across the isocortical sheet, and in multiple dimensions in hippocampus and subiculum. Overall, our study establishes a molecular architecture of the mammalian isocortex and hippocampal formation and begins to shed light on its underlying relationship with the development, evolution, connectivity, and function of these two brain structures.

Microenvironment drives cell state, plasticity, and drug response in pancreatic cancer.

Prognostically relevant RNA expression states exist in pancreatic ductal adenocarcinoma (PDAC), but our understanding of their drivers, stability, and relationship to therapeutic response is limited. To examine these attributes systematically, we profiled metastatic biopsies and matched organoid models at single-cell resolution. In vivo, we identify a new intermediate PDAC transcriptional cell state and uncover distinct site- and state-specific tumor microenvironments (TMEs). Benchmarking models against this reference map, we reveal strong culture-specific biases in cancer cell transcriptional state representation driven by altered TME signals. We restore expression state heterogeneity by adding back in vivo-relevant factors and show plasticity in culture models. Further, we prove that non-genetic modulation of cell state can strongly influence drug responses, uncovering state-specific vulnerabilities. This work provides a broadly applicable framework for aligning cell states across in vivo and ex vivo settings, identifying drivers of transcriptional plasticity and manipulating cell state to target associated vulnerabilities.

A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells.

Tumor-infiltrating myeloid cells (TIMs) are key regulators in tumor progression, but the similarity and distinction of their fundamental properties across different tumors remain elusive. Here, by performing a pan-cancer analysis of single myeloid cells from 210 patients across 15 human cancer types, we identified distinct features of TIMs across cancer types. Mast cells in nasopharyngeal cancer were found to be associated with better prognosis and exhibited an anti-tumor phenotype with a high ratio of TNF+/VEGFA+ cells. Systematic comparison between cDC1- and cDC2-derived LAMP3+ cDCs revealed their differences in transcription factors and external stimulus. Additionally, pro-angiogenic tumor-associated macrophages (TAMs) were characterized with diverse markers across different cancer types, and the composition of TIMs appeared to be associated with certain features of somatic mutations and gene expressions. Our results provide a systematic view of the highly heterogeneous TIMs and suggest future avenues for rational, targeted immunotherapies.

Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has high relapse and low 5-year survival rates. Single-cell profiling in relapsed HCC may aid in the design of effective anticancer therapies, including immunotherapies. We profiled the transcriptomes of ∼17,000 cells from 18 primary or early-relapse HCC cases. Early-relapse tumors have reduced levels of regulatory T cells, increased dendritic cells (DCs), and increased infiltrated CD8+ T cells, compared with primary tumors, in two independent cohorts. Remarkably, CD8+ T cells in recurrent tumors overexpressed KLRB1 (CD161) and displayed an innate-like low cytotoxic state, with low clonal expansion, unlike the classical exhausted state observed in primary HCC. The enrichment of these cells was associated with a worse prognosis. Differential gene expression and interaction analyses revealed potential immune evasion mechanisms in recurrent tumor cells that dampen DC antigen presentation and recruit innate-like CD8+ T cells. Our comprehensive picture of the HCC ecosystem provides deeper insights into immune evasion mechanisms associated with tumor relapse.

Spatiotemporal analysis of human intestinal development at single-cell resolution.

Development of the human intestine is not well understood. Here, we link single-cell RNA sequencing and spatial transcriptomics to characterize intestinal morphogenesis through time. We identify 101 cell states including epithelial and mesenchymal progenitor populations and programs linked to key morphogenetic milestones. We describe principles of crypt-villus axis formation; neural, vascular, mesenchymal morphogenesis, and immune population of the developing gut. We identify the differentiation hierarchies of developing fibroblast and myofibroblast subtypes and describe diverse functions for these including as vascular niche cells. We pinpoint the origins of Peyer's patches and gut-associated lymphoid tissue (GALT) and describe location-specific immune programs. We use our resource to present an unbiased analysis of morphogen gradients that direct sequential waves of cellular differentiation and define cells and locations linked to rare developmental intestinal disorders. We compile a publicly available online resource, spatio-temporal analysis resource of fetal intestinal development (STAR-FINDer), to facilitate further work.

A transcriptional rheostat couples past activity to future sensory responses.

Animals traversing different environments encounter both stable background stimuli and novel cues, which are thought to be detected by primary sensory neurons and then distinguished by downstream brain circuits. Here, we show that each of the ∼1,000 olfactory sensory neuron (OSN) subtypes in the mouse harbors a distinct transcriptome whose content is precisely determined by interactions between its odorant receptor and the environment. This transcriptional variation is systematically organized to support sensory adaptation: expression levels of more than 70 genes relevant to transforming odors into spikes continuously vary across OSN subtypes, dynamically adjust to new environments over hours, and accurately predict acute OSN-specific odor responses. The sensory periphery therefore separates salient signals from predictable background via a transcriptional rheostat whose moment-to-moment state reflects the past and constrains the future; these findings suggest a general model in which structured transcriptional variation within a cell type reflects individual experience.

Single-cell protein activity analysis identifies recurrence-associated renal tumor macrophages.

Clear cell renal carcinoma (ccRCC) is a heterogeneous disease with a variable post-surgical course. To assemble a comprehensive ccRCC tumor microenvironment (TME) atlas, we performed single-cell RNA sequencing (scRNA-seq) of hematopoietic and non-hematopoietic subpopulations from tumor and tumor-adjacent tissue of treatment-naive ccRCC resections. We leveraged the VIPER algorithm to quantitate single-cell protein activity and validated this approach by comparison to flow cytometry. The analysis identified key TME subpopulations, as well as their master regulators and candidate cell-cell interactions, revealing clinically relevant populations, undetectable by gene-expression analysis. Specifically, we uncovered a tumor-specific macrophage subpopulation characterized by upregulation of TREM2/APOE/C1Q, validated by spatially resolved, quantitative multispectral immunofluorescence. In a large clinical validation cohort, these markers were significantly enriched in tumors from patients who recurred following surgery. The study thus identifies TREM2/APOE/C1Q-positive macrophage infiltration as a potential prognostic biomarker for ccRCC recurrence, as well as a candidate therapeutic target.

Identification of a Master Regulator of Differentiation in Toxoplasma.

Toxoplasma gondii chronically infects a quarter of the world's population, and its recrudescence can cause life-threatening disease in immunocompromised individuals and recurrent ocular lesions in the immunocompetent. Acute-stage tachyzoites differentiate into chronic-stage bradyzoites, which form intracellular cysts resistant to immune clearance and existing therapies. The molecular basis of this differentiation is unknown, despite being efficiently triggered by stresses in culture. Through Cas9-mediated screening and single-cell profiling, we identify a Myb-like transcription factor (BFD1) necessary for differentiation in cell culture and in mice. BFD1 accumulates during stress and its synthetic expression is sufficient to drive differentiation. Consistent with its function as a transcription factor, BFD1 binds the promoters of many stage-specific genes and represents a counterpoint to the ApiAP2 factors that dominate our current view of parasite gene regulation. BFD1 provides a genetic switch to study and control Toxoplasma differentiation and will inform prevention and treatment of chronic infections.