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.

Multimodal charting of molecular and functional cell states via in situ electro-sequencing.

Paired mapping of single-cell gene expression and electrophysiology is essential to understand gene-to-function relationships in electrogenic tissues. Here, we developed in situ electro-sequencing (electro-seq) that combines flexible bioelectronics with in situ RNA sequencing to stably map millisecond-timescale electrical activity and profile single-cell gene expression from the same cells across intact biological networks, including cardiac and neural patches. When applied to human-induced pluripotent stem-cell-derived cardiomyocyte patches, in situ electro-seq enabled multimodal in situ analysis of cardiomyocyte electrophysiology and gene expression at the cellular level, jointly defining cell states and developmental trajectories. Using machine-learning-based cross-modal analysis, in situ electro-seq identified gene-to-electrophysiology relationships throughout cardiomyocyte development and accurately reconstructed the evolution of gene expression profiles based on long-term stable electrical measurements. In situ electro-seq could be applicable to create spatiotemporal multimodal maps in electrogenic tissues, potentiating the discovery of cell types and gene programs responsible for electrophysiological function and dysfunction.

Systematic identification and characterization of genes in the regulation and biogenesis of photosynthetic machinery.

Photosynthesis is central to food production and the Earth's biogeochemistry, yet the molecular basis for its regulation remains poorly understood. Here, using high-throughput genetics in the model eukaryotic alga Chlamydomonas reinhardtii, we identify with high confidence (false discovery rate [FDR] < 0.11) 70 poorly characterized genes required for photosynthesis. We then enable the functional characterization of these genes by providing a resource of proteomes of mutant strains, each lacking one of these genes. The data allow assignment of 34 genes to the biogenesis or regulation of one or more specific photosynthetic complexes. Further analysis uncovers biogenesis/regulatory roles for at least seven proteins, including five photosystem I mRNA maturation factors, the chloroplast translation factor MTF1, and the master regulator PMR1, which regulates chloroplast genes via nuclear-expressed factors. Our work provides a rich resource identifying regulatory and functional genes and placing them into pathways, thereby opening the door to a system-level understanding of photosynthesis.

Early Alzheimer's disease pathology in human cortex involves transient cell states.

Cellular perturbations underlying Alzheimer's disease (AD) are primarily studied in human postmortem samples and model organisms. Here, we generated a single-nucleus atlas from a rare cohort of cortical biopsies from living individuals with varying degrees of AD pathology. We next performed a systematic cross-disease and cross-species integrative analysis to identify a set of cell states that are specific to early AD pathology. These changes-which we refer to as the early cortical amyloid response-were prominent in neurons, wherein we identified a transitional hyperactive state preceding the loss of excitatory neurons, which we confirmed by acute slice physiology on independent biopsy specimens. Microglia overexpressing neuroinflammatory-related processes also expanded as AD pathology increased. Finally, both oligodendrocytes and pyramidal neurons upregulated genes associated with ß-amyloid production and processing during this early hyperactive phase. Our integrative analysis provides an organizing framework for targeting circuit dysfunction, neuroinflammation, and amyloid production early in AD pathogenesis.

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.

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.

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.

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.

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 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.

Human neutralizing antibodies against SARS-CoV-2 require intact Fc effector functions for optimal therapeutic protection.

SARS-CoV-2 has caused the global COVID-19 pandemic. Although passively delivered neutralizing antibodies against SARS-CoV-2 show promise in clinical trials, their mechanism of action in vivo is incompletely understood. Here, we define correlates of protection of neutralizing human monoclonal antibodies (mAbs) in SARS-CoV-2-infected animals. Whereas Fc effector functions are dispensable when representative neutralizing mAbs are administered as prophylaxis, they are required for optimal protection as therapy. When given after infection, intact mAbs reduce SARS-CoV-2 burden and lung disease in mice and hamsters better than loss-of-function Fc variant mAbs. Fc engagement of neutralizing antibodies mitigates inflammation and improves respiratory mechanics, and transcriptional profiling suggests these phenotypes are associated with diminished innate immune signaling and preserved tissue repair. Immune cell depletions establish that neutralizing mAbs require monocytes and CD8+ T cells for optimal clinical and virological benefit. Thus, potently neutralizing mAbs utilize Fc effector functions during therapy to mitigate lung infection and disease.

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.