Bacterial droplet-based single-cell RNA-seq reveals antibiotic-associated heterogeneous cellular states.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis.

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

Discovering dominant tumor immune archetypes in a pan-cancer census.

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

Vision-dependent specification of cell types and function in the developing cortex.

The role of postnatal experience in sculpting cortical circuitry, while long appreciated, is poorly understood at the level of cell types. We explore this in the mouse primary visual cortex (V1) using single-nucleus RNA sequencing, visual deprivation, genetics, and functional imaging. We find that vision selectively drives the specification of glutamatergic cell types in upper layers (L) (L2/3/4), while deeper-layer glutamatergic, GABAergic, and non-neuronal cell types are established prior to eye opening. L2/3 cell types form an experience-dependent spatial continuum defined by the graded expression of ∼200 genes, including regulators of cell adhesion and synapse formation. One of these genes, Igsf9b, a vision-dependent gene encoding an inhibitory synaptic cell adhesion molecule, is required for the normal development of binocular responses in L2/3. In summary, vision preferentially regulates the development of upper-layer glutamatergic cell types through the regulation of cell-type-specific gene expression programs.

Cellular architecture of human brain metastases.

Brain metastasis (BrM) is the most common form of brain cancer, characterized by neurologic disability and an abysmal prognosis. Unfortunately, our understanding of the biology underlying human BrMs remains rudimentary. Here, we present an integrative analysis of >100,000 malignant and non-malignant cells from 15 human parenchymal BrMs, generated by single-cell transcriptomics, mass cytometry, and complemented with mouse model- and in silico approaches. We interrogated the composition of BrM niches, molecularly defined the blood-tumor interface, and revealed stromal immunosuppressive states enriched with infiltrated T cells and macrophages. Specific single-cell interrogation of metastatic tumor cells provides a framework of 8 functional cell programs that coexist or anticorrelate. Collectively, these programs delineate two functional BrM archetypes, one proliferative and the other inflammatory, that are evidently shaped through tumor-immune interactions. Our resource provides a foundation to understand the molecular basis of BrM in patients with tumor cell-intrinsic and host environmental traits.

Chimeric contribution of human extended pluripotent stem cells to monkey embryos ex vivo.

Interspecies chimera formation with human pluripotent stem cells (hPSCs) represents a necessary alternative to evaluate hPSC pluripotency in vivo and might constitute a promising strategy for various regenerative medicine applications, including the generation of organs and tissues for transplantation. Studies using mouse and pig embryos suggest that hPSCs do not robustly contribute to chimera formation in species evolutionarily distant to humans. We studied the chimeric competency of human extended pluripotent stem cells (hEPSCs) in cynomolgus monkey (Macaca fascicularis) embryos cultured ex vivo. We demonstrate that hEPSCs survived, proliferated, and generated several peri- and early post-implantation cell lineages inside monkey embryos. We also uncovered signaling events underlying interspecific crosstalk that may help shape the unique developmental trajectories of human and monkey cells within chimeric embryos. These results may help to better understand early human development and primate evolution and develop strategies to improve human chimerism in evolutionarily distant species.

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.

Molecular topography of an entire nervous system.

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

EASI-FISH for thick tissue defines lateral hypothalamus spatio-molecular organization.

Determining the spatial organization and morphological characteristics of molecularly defined cell types is a major bottleneck for characterizing the architecture underpinning brain function. We developed Expansion-Assisted Iterative Fluorescence In Situ Hybridization (EASI-FISH) to survey gene expression in brain tissue, as well as a turnkey computational pipeline to rapidly process large EASI-FISH image datasets. EASI-FISH was optimized for thick brain sections (300 µm) to facilitate reconstruction of spatio-molecular domains that generalize across brains. Using the EASI-FISH pipeline, we investigated the spatial distribution of dozens of molecularly defined cell types in the lateral hypothalamic area (LHA), a brain region with poorly defined anatomical organization. Mapping cell types in the LHA revealed nine spatially and molecularly defined subregions. EASI-FISH also facilitates iterative reanalysis of scRNA-seq datasets to determine marker-genes that further dissociated spatial and morphological heterogeneity. The EASI-FISH pipeline democratizes mapping molecularly defined cell types, enabling discoveries about brain organization.

Antigen dominance hierarchies shape TCF1+ progenitor CD8 T cell phenotypes in tumors.

CD8 T cell responses against different tumor neoantigens occur simultaneously, yet little is known about the interplay between responses and its impact on T cell function and tumor control. In mouse lung adenocarcinoma, we found that immunodominance is established in tumors, wherein CD8 T cell expansion is predominantly driven by the antigen that most stably binds MHC. T cells responding to subdominant antigens were enriched for a TCF1+ progenitor phenotype correlated with response to immune checkpoint blockade (ICB) therapy. However, the subdominant T cell response did not preferentially benefit from ICB due to a dysfunctional subset of TCF1+ cells marked by CCR6 and Tc17 differentiation. Analysis of human samples and sequencing datasets revealed that CCR6+ TCF1+ cells exist across human cancers and are not correlated with ICB response. Vaccination eliminated CCR6+ TCF1+ cells and dramatically improved the subdominant response, highlighting a strategy to optimally engage concurrent neoantigen responses against tumors.

Stem-like intestinal Th17 cells give rise to pathogenic effector T cells during autoimmunity.

While intestinal Th17 cells are critical for maintaining tissue homeostasis, recent studies have implicated their roles in the development of extra-intestinal autoimmune diseases including multiple sclerosis. However, the mechanisms by which tissue Th17 cells mediate these dichotomous functions remain unknown. Here, we characterized the heterogeneity, plasticity, and migratory phenotypes of tissue Th17 cells in vivo by combined fate mapping with profiling of the transcriptomes and TCR clonotypes of over 84,000 Th17 cells at homeostasis and during CNS autoimmune inflammation. Inter- and intra-organ single-cell analyses revealed a homeostatic, stem-like TCF1+ IL-17+ SLAMF6+ population that traffics to the intestine where it is maintained by the microbiota, providing a ready reservoir for the IL-23-driven generation of encephalitogenic GM-CSF+ IFN-γ+ CXCR6+ T cells. Our study defines a direct in vivo relationship between IL-17+ non-pathogenic and GM-CSF+ and IFN-γ+ pathogenic Th17 populations and provides a mechanism by which homeostatic intestinal Th17 cells direct extra-intestinal autoimmune disease.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Single-Cell Transcriptome Atlas of Murine Endothelial Cells.

The heterogeneity of endothelial cells (ECs) across tissues remains incompletely inventoried. We constructed an atlas of >32,000 single-EC transcriptomes from 11 mouse tissues and identified 78 EC subclusters, including Aqp7+ intestinal capillaries and angiogenic ECs in healthy tissues. ECs from brain/testis, liver/spleen, small intestine/colon, and skeletal muscle/heart pairwise expressed partially overlapping marker genes. Arterial, venous, and lymphatic ECs shared more markers in more tissues than did heterogeneous capillary ECs. ECs from different vascular beds (arteries, capillaries, veins, lymphatics) exhibited transcriptome similarity across tissues, but the tissue (rather than the vessel) type contributed to the EC heterogeneity. Metabolic transcriptome analysis revealed a similar tissue-grouping phenomenon of ECs and heterogeneous metabolic gene signatures in ECs between tissues and between vascular beds within a single tissue in a tissue-type-dependent pattern. The EC atlas taxonomy enabled identification of EC subclusters in public scRNA-seq datasets and provides a powerful discovery tool and resource value.

Large-Scale Topological Changes Restrain Malignant Progression in Colorectal Cancer.

Widespread changes to DNA methylation and chromatin are well documented in cancer, but the fate of higher-order chromosomal structure remains obscure. Here we integrated topological maps for colon tumors and normal colons with epigenetic, transcriptional, and imaging data to characterize alterations to chromatin loops, topologically associated domains, and large-scale compartments. We found that spatial partitioning of the open and closed genome compartments is profoundly compromised in tumors. This reorganization is accompanied by compartment-specific hypomethylation and chromatin changes. Additionally, we identify a compartment at the interface between the canonical A and B compartments that is reorganized in tumors. Remarkably, similar shifts were evident in non-malignant cells that have accumulated excess divisions. Our analyses suggest that these topological changes repress stemness and invasion programs while inducing anti-tumor immunity genes and may therefore restrain malignant progression. Our findings call into question the conventional view that tumor-associated epigenomic alterations are primarily oncogenic.

Multi-Omics Resolves a Sharp Disease-State Shift between Mild and Moderate COVID-19.

We present an integrated analysis of the clinical measurements, immune cells, and plasma multi-omics of 139 COVID-19 patients representing all levels of disease severity, from serial blood draws collected during the first week of infection following diagnosis. We identify a major shift between mild and moderate disease, at which point elevated inflammatory signaling is accompanied by the loss of specific classes of metabolites and metabolic processes. Within this stressed plasma environment at moderate disease, multiple unusual immune cell phenotypes emerge and amplify with increasing disease severity. We condensed over 120,000 immune features into a single axis to capture how different immune cell classes coordinate in response to SARS-CoV-2. This immune-response axis independently aligns with the major plasma composition changes, with clinical metrics of blood clotting, and with the sharp transition between mild and moderate disease. This study suggests that moderate disease may provide the most effective setting for therapeutic intervention.

Origins and Proliferative States of Human Oligodendrocyte Precursor Cells.

Human cerebral cortex size and complexity has increased greatly during evolution. While increased progenitor diversity and enhanced proliferative potential play important roles in human neurogenesis and gray matter expansion, the mechanisms of human oligodendrogenesis and white matter expansion remain largely unknown. Here, we identify EGFR-expressing "Pre-OPCs" that originate from outer radial glial cells (oRGs) and undergo mitotic somal translocation (MST) during division. oRG-derived Pre-OPCs provide an additional source of human cortical oligodendrocyte precursor cells (OPCs) and define a lineage trajectory. We further show that human OPCs undergo consecutive symmetric divisions to exponentially increase the progenitor pool size. Additionally, we find that the OPC-enriched gene, PCDH15, mediates daughter cell repulsion and facilitates proliferation. These findings indicate properties of OPC derivation, proliferation, and dispersion important for human white matter expansion and myelination.

Type V Collagen in Scar Tissue Regulates the Size of Scar after Heart Injury.

Scar tissue size following myocardial infarction is an independent predictor of cardiovascular outcomes, yet little is known about factors regulating scar size. We demonstrate that collagen V, a minor constituent of heart scars, regulates the size of heart scars after ischemic injury. Depletion of collagen V led to a paradoxical increase in post-infarction scar size with worsening of heart function. A systems genetics approach across 100 in-bred strains of mice demonstrated that collagen V is a critical driver of postinjury heart function. We show that collagen V deficiency alters the mechanical properties of scar tissue, and altered reciprocal feedback between matrix and cells induces expression of mechanosensitive integrins that drive fibroblast activation and increase scar size. Cilengitide, an inhibitor of specific integrins, rescues the phenotype of increased post-injury scarring in collagen-V-deficient mice. These observations demonstrate that collagen V regulates scar size in an integrin-dependent manner.

Niche-Selective Inhibition of Pathogenic Th17 Cells by Targeting Metabolic Redundancy.

Targeting glycolysis has been considered therapeutically intractable owing to its essential housekeeping role. However, the context-dependent requirement for individual glycolytic steps has not been fully explored. We show that CRISPR-mediated targeting of glycolysis in T cells in mice results in global loss of Th17 cells, whereas deficiency of the glycolytic enzyme glucose phosphate isomerase (Gpi1) selectively eliminates inflammatory encephalitogenic and colitogenic Th17 cells, without substantially affecting homeostatic microbiota-specific Th17 cells. In homeostatic Th17 cells, partial blockade of glycolysis upon Gpi1 inactivation was compensated by pentose phosphate pathway flux and increased mitochondrial respiration. In contrast, inflammatory Th17 cells experience a hypoxic microenvironment known to limit mitochondrial respiration, which is incompatible with loss of Gpi1. Our study suggests that inhibiting glycolysis by targeting Gpi1 could be an effective therapeutic strategy with minimum toxicity for Th17-mediated autoimmune diseases, and, more generally, that metabolic redundancies can be exploited for selective targeting of disease processes.

CD81 Controls Beige Fat Progenitor Cell Growth and Energy Balance via FAK Signaling.

Adipose tissues dynamically remodel their cellular composition in response to external cues by stimulating beige adipocyte biogenesis; however, the developmental origin and pathways regulating this process remain insufficiently understood owing to adipose tissue heterogeneity. Here, we employed single-cell RNA-seq and identified a unique subset of adipocyte progenitor cells (APCs) that possessed the cell-intrinsic plasticity to give rise to beige fat. This beige APC population is proliferative and marked by cell-surface proteins, including PDGFRα, Sca1, and CD81. Notably, CD81 is not only a beige APC marker but also required for de novo beige fat biogenesis following cold exposure. CD81 forms a complex with αV/ß1 and αV/ß5 integrins and mediates the activation of integrin-FAK signaling in response to irisin. Importantly, CD81 loss causes diet-induced obesity, insulin resistance, and adipose tissue inflammation. These results suggest that CD81 functions as a key sensor of external inputs and controls beige APC proliferation and whole-body energy homeostasis.