The CD4 Versus CD8 T Cell Fate Decision: A Multiomics-Informed Perspective.

The choice of developing thymocytes to become CD8+ cytotoxic or CD4+ helper T cells has been intensely studied, but many of the underlying mechanisms remain to be elucidated. Recent multiomics approaches have provided much higher resolution analysis of gene expression in developing thymocytes than was previously achievable, thereby offering a fresh perspective on this question. Focusing on our recent studies using CITE-seq (cellular indexing of transcriptomes and epitopes) analyses of mouse thymocytes, we present a detailed timeline of RNA and protein expression changes during CD8 versus CD4 T cell differentiation. We also revisit our current understanding of the links between T cell receptor signaling and expression of the lineage-defining transcription factors ThPOK and RUNX3. Finally, we propose a sequential selection model to explain the tight linkage between MHC-I versus MHC-II recognition and T cell lineage choice. This model incorporates key aspects of previously proposed kinetic signaling, instructive, and stochastic/selection models.

Insights Gained from Single-Cell Analysis of Immune Cells in the Tumor Microenvironment.

Understanding tumor immune microenvironments is critical for identifying immune modifiers of cancer progression and developing cancer immunotherapies. Recent applications of single-cell RNA sequencing (scRNA-seq) in dissecting tumor microenvironments have brought important insights into the biology of tumor-infiltrating immune cells, including their heterogeneity, dynamics, and potential roles in both disease progression and response to immune checkpoint inhibitors and other immunotherapies. This review focuses on the advances in knowledge of tumor immune microenvironments acquired from scRNA-seq studies across multiple types of human tumors, with a particular emphasis on the study of phenotypic plasticity and lineage dynamics of immune cells in the tumor environment. We also discuss several imminent questions emerging from scRNA-seq observations and their potential solutions on the horizon.

Microglia and Central Nervous System-Associated Macrophages-From Origin to Disease Modulation.

The immune system of the central nervous system (CNS) consists primarily of innate immune cells. These are highly specialized macrophages found either in the parenchyma, called microglia, or at the CNS interfaces, such as leptomeningeal, perivascular, and choroid plexus macrophages. While they were primarily thought of as phagocytes, their function extends well beyond simple removal of cell debris during development and diseases. Brain-resident innate immune cells were found to be plastic, long-lived, and host to an outstanding number of risk genes for multiple pathologies. As a result, they are now considered the most suitable targets for modulating CNS diseases. Additionally, recent single-cell technologies enhanced our molecular understanding of their origins, fates, interactomes, and functional cell statesduring health and perturbation. Here, we review the current state of our understanding and challenges of the myeloid cell biology in the CNS and treatment options for related diseases.

Immunology in the Era of Single-Cell Technologies.

Immune cells are characterized by diversity, specificity, plasticity, and adaptability-properties that enable them to contribute to homeostasis and respond specifically and dynamically to the many threats encountered by the body. Single-cell technologies, including the assessment of transcriptomics, genomics, and proteomics at the level of individual cells, are ideally suited to studying these properties of immune cells. In this review we discuss the benefits of adopting single-cell approaches in studying underappreciated qualities of immune cells and highlight examples where these technologies have been critical to advancing our understanding of the immune system in health and disease.

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.

Private information leakage from single-cell count matrices.

The increase in publicly available human single-cell datasets, encompassing millions of cells from many donors, has significantly enhanced our understanding of complex biological processes. However, the accessibility of these datasets raises significant privacy concerns. Due to the inherent noise in single-cell measurements and the scarcity of population-scale single-cell datasets, recent private information quantification studies have focused on bulk gene expression data sharing. To address this gap, we demonstrate that individuals in single-cell gene expression datasets are vulnerable to linking attacks, where attackers can infer their sensitive phenotypic information using publicly available tissue or cell-type-specific expression quantitative trait loci (eQTLs) information. We further develop a method for genotype prediction and genotype-phenotype linking that remains effective without relying on eQTL information. We show that variants from one study can be exploited to uncover private information about individuals in another study.

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.

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.

Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging.

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

The future of rapid and automated single-cell data analysis using reference mapping.

As the number of single-cell datasets continues to grow rapidly, workflows that map new data to well-curated reference atlases offer enormous promise for the biological community. In this perspective, we discuss key computational challenges and opportunities for single-cell reference-mapping algorithms. We discuss how mapping algorithms will enable the integration of diverse datasets across disease states, molecular modalities, genetic perturbations, and diverse species and will eventually replace manual and laborious unsupervised clustering pipelines.

A programmable reaction-diffusion system for spatiotemporal cell signaling circuit design.

Cells self-organize molecules in space and time to generate complex behaviors, but we lack synthetic strategies for engineering spatiotemporal signaling. We present a programmable reaction-diffusion platform for designing protein oscillations, patterns, and circuits in mammalian cells using two bacterial proteins, MinD and MinE (MinDE). MinDE circuits act like "single-cell radios," emitting frequency-barcoded fluorescence signals that can be spectrally isolated and analyzed using digital signal processing tools. We define how to genetically program these signals and connect their spatiotemporal dynamics to cell biology using engineerable protein-protein interactions. This enabled us to construct sensitive reporter circuits that broadcast endogenous cell signaling dynamics on a frequency-barcoded imaging channel and to build control signal circuits that synthetically pattern activities in the cell, such as protein condensate assembly and actin filamentation. Our work establishes a paradigm for visualizing, probing, and engineering cellular activities at length and timescales critical for biological function.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders.

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

The life-saving benefit of dexamethasone in severe COVID-19 is linked to a reversal of monocyte dysregulation.

Dexamethasone is a life-saving treatment for severe COVID-19, yet its mechanism of action is unknown, and many patients deteriorate or die despite timely treatment initiation. Here, we identify dexamethasone treatment-induced cellular and molecular changes associated with improved survival in COVID-19 patients. We observed a reversal of transcriptional hallmark signatures in monocytes associated with severe COVID-19 and the induction of a monocyte substate characterized by the expression of glucocorticoid-response genes. These molecular responses to dexamethasone were detected in circulating and pulmonary monocytes, and they were directly linked to survival. Monocyte single-cell RNA sequencing (scRNA-seq)-derived signatures were enriched in whole blood transcriptomes of patients with fatal outcome in two independent cohorts, highlighting the potential for identifying non-responders refractory to dexamethasone. Our findings link the effects of dexamethasone to specific immunomodulation and reversal of monocyte dysregulation, and they highlight the potential of single-cell omics for monitoring in vivo target engagement of immunomodulatory drugs and for patient stratification for precision medicine approaches.

Neoadjuvant PARPi or chemotherapy in ovarian cancer informs targeting effector Treg cells for homologous-recombination-deficient tumors.

Homologous recombination deficiency (HRD) is prevalent in cancer, sensitizing tumor cells to poly (ADP-ribose) polymerase (PARP) inhibition. However, the impact of HRD and related therapies on the tumor microenvironment (TME) remains elusive. Our study generates single-cell gene expression and T cell receptor profiles, along with validatory multimodal datasets from >100 high-grade serous ovarian cancer (HGSOC) samples, primarily from a phase II clinical trial (NCT04507841). Neoadjuvant monotherapy with the PARP inhibitor (PARPi) niraparib achieves impressive 62.5% and 73.6% response rates per RECIST v.1.1 and GCIG CA125, respectively. We identify effector regulatory T cells (eTregs) as key responders to HRD and neoadjuvant therapies, co-occurring with other tumor-reactive T cells, particularly terminally exhausted CD8+ T cells (Tex). TME-wide interferon signaling correlates with cancer cells upregulating MHC class II and co-inhibitory ligands, potentially driving Treg and Tex fates. Depleting eTregs in HRD mouse models, with or without PARP inhibition, significantly suppresses tumor growth without observable toxicities, underscoring the potential of eTreg-focused therapeutics for HGSOC and other HRD-related tumors.

Time-resolved single-cell transcriptomics defines immune trajectories in glioblastoma.

Deciphering the cell-state transitions underlying immune adaptation across time is fundamental for advancing biology. Empirical in vivo genomic technologies that capture cellular dynamics are currently lacking. We present Zman-seq, a single-cell technology recording transcriptomic dynamics across time by introducing time stamps into circulating immune cells, tracking them in tissues for days. Applying Zman-seq resolved cell-state and molecular trajectories of the dysfunctional immune microenvironment in glioblastoma. Within 24 hours of tumor infiltration, cytotoxic natural killer cells transitioned to a dysfunctional program regulated by TGFB1 signaling. Infiltrating monocytes differentiated into immunosuppressive macrophages, characterized by the upregulation of suppressive myeloid checkpoints Trem2, Il18bp, and Arg1, over 36 to 48 hours. Treatment with an antagonistic anti-TREM2 antibody reshaped the tumor microenvironment by redirecting the monocyte trajectory toward pro-inflammatory macrophages. Zman-seq is a broadly applicable technology, enabling empirical measurements of differentiation trajectories, which can enhance the development of more efficacious immunotherapies.

A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Cellular architecture of evolving neuroinflammatory lesions and multiple sclerosis pathology.

Multiple sclerosis (MS) is a neurological disease characterized by multifocal lesions and smoldering pathology. Although single-cell analyses provided insights into cytopathology, evolving cellular processes underlying MS remain poorly understood. We investigated the cellular dynamics of MS by modeling temporal and regional rates of disease progression in mouse experimental autoimmune encephalomyelitis (EAE). By performing single-cell spatial expression profiling using in situ sequencing (ISS), we annotated disease neighborhoods and found centrifugal evolution of active lesions. We demonstrated that disease-associated (DA)-glia arise independently of lesions and are dynamically induced and resolved over the disease course. Single-cell spatial mapping of human archival MS spinal cords confirmed the differential distribution of homeostatic and DA-glia, enabled deconvolution of active and inactive lesions into sub-compartments, and identified new lesion areas. By establishing a spatial resource of mouse and human MS neuropathology at a single-cell resolution, our study unveils the intricate cellular dynamics underlying MS.

Method of moments framework for differential expression analysis of single-cell RNA sequencing data.

Differential expression analysis of single-cell RNA sequencing (scRNA-seq) data is central for characterizing how experimental factors affect the distribution of gene expression. However, distinguishing between biological and technical sources of cell-cell variability and assessing the statistical significance of quantitative comparisons between cell groups remain challenging. We introduce Memento, a tool for robust and efficient differential analysis of mean expression, variability, and gene correlation from scRNA-seq data, scalable to millions of cells and thousands of samples. We applied Memento to 70,000 tracheal epithelial cells to identify interferon-responsive genes, 160,000 CRISPR-Cas9 perturbed T cells to reconstruct gene-regulatory networks, 1.2 million peripheral blood mononuclear cells (PBMCs) to map cell-type-specific quantitative trait loci (QTLs), and the 50-million-cell CELLxGENE Discover corpus to compare arbitrary cell groups. In all cases, Memento identified more significant and reproducible differences in mean expression compared with existing methods. It also identified differences in variability and gene correlation that suggest distinct transcriptional regulation mechanisms imparted by perturbations.

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.