Most axonal mitochondria in cortical pyramidal neurons lack mitochondrial DNA and consume ATP.

In neurons of the mammalian central nervous system (CNS), axonal mitochondria are thought to be indispensable for supplying ATP during energy-consuming processes such as neurotransmitter release. Here, we demonstrate using multiple, independent, in vitro and in vivo approaches that the majority (~80-90%) of axonal mitochondria in cortical pyramidal neurons (CPNs), lack mitochondrial DNA (mtDNA). Using dynamic, optical imaging analysis of genetically encoded sensors for mitochondrial matrix ATP and pH, we demonstrate that in axons of CPNs, but not in their dendrites, mitochondrial complex V (ATP synthase) functions in a reverse way, consuming ATP and protruding H+ out of the matrix to maintain mitochondrial membrane potential. Our results demonstrate that in mammalian CPNs, axonal mitochondria do not play a major role in ATP supply, despite playing other functions critical to regulating neurotransmission such as Ca2+ buffering.

A Roadmap for the Human Gut Cell Atlas.

The number of studies investigating the human gastrointestinal tract using various single-cell profiling methods has increased substantially in the past few years. Although this increase provides a unique opportunity for the generation of the first comprehensive Human Gut Cell Atlas (HGCA), there remains a range of major challenges ahead. Above all, the ultimate success will largely depend on a structured and coordinated approach that aligns global efforts undertaken by a large number of research groups. In this Roadmap, we discuss a comprehensive forward-thinking direction for the generation of the HGCA on behalf of the Gut Biological Network of the Human Cell Atlas. Based on the consensus opinion of experts from across the globe, we outline the main requirements for the first complete HGCA by summarizing existing data sets and highlighting anatomical regions and/or tissues with limited coverage. We provide recommendations for future studies and discuss key methodologies and the importance of integrating the healthy gut atlas with related diseases and gut organoids. Importantly, we critically overview the computational tools available and provide recommendations to overcome key challenges.

CRISPRmap: Sequencing-free optical pooled screens mapping multi-omic phenotypes in cells and tissue.

Pooled genetic screens are powerful tools to study gene function in a high-throughput manner. Typically, sequencing-based screens require cell lysis, which limits the examination of critical phenotypes such as cell morphology, protein subcellular localization, and cell-cell/tissue interactions. In contrast, emerging optical pooled screening methods enable the investigation of these spatial phenotypes in response to targeted CRISPR perturbations. In this study, we report a multi-omic optical pooled CRISPR screening method, which we have named CRISPRmap. Our method combines a novel in situ CRISPR guide identifying barcode readout approach with concurrent multiplexed immunofluorescence and in situ RNA detection. CRISPRmap barcodes are detected and read out through combinatorial hybridization of DNA oligos, enhancing barcode detection efficiency, while reducing both dependency on third party proprietary sequencing reagents and assay cost. Notably, we conducted a multi-omic base-editing screen in a breast cancer cell line on core DNA damage repair genes involved in the homologous recombination and Fanconi anemia pathways investigating how nucleotide variants in those genes influence DNA damage signaling and cell cycle regulation following treatment with ionizing radiation or DNA damaging agents commonly used for cancer therapy. Approximately a million cells were profiled with our multi-omic approach, providing a comprehensive phenotypic assessment of the functional consequences of the studied variants. CRISPRmap enabled us to pinpoint likely-pathogenic patient-derived mutations that were previously classified as variants of unknown clinical significance. Furthermore, our approach effectively distinguished barcodes of a pooled library in tumor tissue, and we coupled it with cell-type and molecular phenotyping by cyclic immunofluorescence. Multi-omic spatial analysis of how CRISPR-perturbed cells respond to various environmental cues in the tissue context offers the potential to significantly expand our understanding of tissue biology in both health and disease.

Promotion of cholangiocarcinoma growth by diverse cancer-associated fibroblast subpopulations.

Cancer-associated fibroblasts (CAF) are a poorly characterized cell population in the context of liver cancer. Our study investigates CAF functions in intrahepatic cholangiocarcinoma (ICC), a highly desmoplastic liver tumor. Genetic tracing, single-cell RNA sequencing, and ligand-receptor analyses uncovered hepatic stellate cells (HSC) as the main source of CAF and HSC-derived CAF as the dominant population interacting with tumor cells. In mice, CAF promotes ICC progression, as revealed by HSC-selective CAF depletion. In patients, a high panCAF signature is associated with decreased survival and increased recurrence. Single-cell RNA sequencing segregates CAF into inflammatory and growth factor-enriched (iCAF) and myofibroblastic (myCAF) subpopulations, displaying distinct ligand-receptor interactions. myCAF-expressed hyaluronan synthase 2, but not type I collagen, promotes ICC. iCAF-expressed hepatocyte growth factor enhances ICC growth via tumor-expressed MET, thus directly linking CAF to tumor cells. In summary, our data demonstrate promotion of desmoplastic ICC growth by therapeutically targetable CAF subtype-specific mediators, but not by type I collagen.

Author Correction: Nuclei multiplexing with barcoded antibodies for single-nucleus genomics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Nuclei multiplexing with barcoded antibodies for single-nucleus genomics.

Single-nucleus RNA-seq (snRNA-seq) enables the interrogation of cellular states in complex tissues that are challenging to dissociate or are frozen, and opens the way to human genetics studies, clinical trials, and precise cell atlases of large organs. However, such applications are currently limited by batch effects, processing, and costs. Here, we present an approach for multiplexing snRNA-seq, using sample-barcoded antibodies to uniquely label nuclei from distinct samples. Comparing human brain cortex samples profiled with or without hashing antibodies, we demonstrate that nucleus hashing does not significantly alter recovered profiles. We develop DemuxEM, a computational tool that detects inter-sample multiplets and assigns singlets to their sample of origin, and validate its accuracy using sex-specific gene expression, species-mixing and natural genetic variation. Our approach will facilitate tissue atlases of isogenic model organisms or from multiple biopsies or longitudinal samples of one donor, and large-scale perturbation screens.

Geometry-dependent functional changes in iPSC-derived cardiomyocytes probed by functional imaging and RNA sequencing.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) are a promising platform for cardiac studies in vitro, and possibly for tissue repair in humans. However, hiPSC-CM cells tend to retain morphology, metabolism, patterns of gene expression, and electrophysiology similar to that of embryonic cardiomyocytes. We grew hiPSC-CM in patterned islands of different sizes and shapes, and measured the effect of island geometry on action potential waveform and calcium dynamics using optical recordings of voltage and calcium from 970 islands of different sizes. hiPSC-CM in larger islands showed electrical and calcium dynamics indicative of greater functional maturity. We then compared transcriptional signatures of the small and large islands against a developmental time course of cardiac differentiation. Although island size had little effect on expression of most genes whose levels differed between hiPSC-CM and adult primary CM, we identified a subset of genes for which island size drove the majority (58%) of the changes associated with functional maturation. Finally, we patterned hiPSC-CM on islands with a variety of shapes to probe the relative contributions of soluble factors, electrical coupling, and direct cell-cell contacts to the functional maturation. Collectively, our data show that optical electrophysiology is a powerful tool for assaying hiPSC-CM maturation, and that island size powerfully drives activation of a subset of genes involved in cardiac maturation.

Protein C receptor (PROCR) is a negative regulator of Th17 pathogenicity.

Th17 cells are key players in defense against pathogens and maintaining tissue homeostasis, but also act as critical drivers of autoimmune diseases. Based on single-cell RNA-seq profiling of pathogenic versus nonpathogenic Th17 cells, we identified protein C receptor (PROCR) as a cell surface molecule expressed in covariance with the regulatory module of Th17 cells. Although PROCR expression in T cells was controlled by the cooperative action of the Th17 lineage-specific transcription factors RORγt, IRF4, and STAT3, PROCR negatively regulated Th17 differentiation. CD4+ T cells from PROCR low expressor mutant mice readily differentiated into Th17 cells, whereas addition of the PROCR ligand, activated protein C, inhibited Th17 differentiation in vitro. In addition, PROCR acted as a negative regulator of Th17 pathogenicity in that it down-regulated expression of several pathogenic signature genes, including IL-1 and IL-23 receptors. Furthermore, T cell-specific deficiency of PROCR resulted in the exacerbation of experimental autoimmune encephalomyelitis (EAE) and higher frequencies of Th17 cell in vivo, indicating that PROCR also inhibits pathogenicity of Th17 cells in vivo. PROCR thus does not globally inhibit Th17 responses but could be targeted to selectively inhibit proinflammatory Th17 cells.

Identification and characterization of latency-associated peptide-expressing γδ T cells.

γδ T cells are a subset of lymphocytes specialized in protecting the host against pathogens and tumours. Here we describe a subset of regulatory γδ T cells that express the latency-associated peptide (LAP), a membrane-bound TGF-ß1. Thymic CD27+IFN-γ+CCR9+α4ß7+TCRγδ+ cells migrate to the periphery, particularly to Peyer's patches and small intestine lamina propria, where they upregulate LAP, downregulate IFN-γ via ATF-3 expression and acquire a regulatory phenotype. TCRγδ+LAP+ cells express antigen presentation molecules and function as antigen presenting cells that induce CD4+Foxp3+ regulatory T cells, although TCRγδ+LAP+ cells do not themselves express Foxp3. Identification of TCRγδ+LAP+ regulatory cells provides an avenue for understanding immune regulation and biologic processes linked to intestinal function and disease.

CD5L/AIM Regulates Lipid Biosynthesis and Restrains Th17 Cell Pathogenicity.

Th17 cells play a critical role in host defense against extracellular pathogens and tissue homeostasis but can induce autoimmunity. The mechanisms implicated in balancing "pathogenic" and "non-pathogenic" Th17 cell states remain largely unknown. We used single-cell RNA-seq to identify CD5L/AIM as a regulator expressed in non-pathogenic, but not in pathogenic Th17 cells. Although CD5L does not affect Th17 differentiation, it is a functional switch that regulates the pathogenicity of Th17 cells. Loss of CD5L converts non-pathogenic Th17 cells into pathogenic cells that induce autoimmunity. CD5L mediates this effect by modulating the intracellular lipidome, altering fatty acid composition and restricting cholesterol biosynthesis and, thus, ligand availability for Rorγt, the master transcription factor of Th17 cells. Our study identifies CD5L as a critical regulator of the Th17 cell functional state and highlights the importance of lipid metabolism in balancing immune protection and disease induced by T cells.

Single-Cell Genomics Unveils Critical Regulators of Th17 Cell Pathogenicity.

Extensive cellular heterogeneity exists within specific immune-cell subtypes classified as a single lineage, but its molecular underpinnings are rarely characterized at a genomic scale. Here, we use single-cell RNA-seq to investigate the molecular mechanisms governing heterogeneity and pathogenicity of Th17 cells isolated from the central nervous system (CNS) and lymph nodes (LN) at the peak of autoimmune encephalomyelitis (EAE) or differentiated in vitro under either pathogenic or non-pathogenic polarization conditions. Computational analysis relates a spectrum of cellular states in vivo to in-vitro-differentiated Th17 cells and unveils genes governing pathogenicity and disease susceptibility. Using knockout mice, we validate four new genes: Gpr65, Plzp, Toso, and Cd5l (in a companion paper). Cellular heterogeneity thus informs Th17 function in autoimmunity and can identify targets for selective suppression of pathogenic Th17 cells while potentially sparing non-pathogenic tissue-protective ones.

Single-cell RNA-seq reveals dynamic paracrine control of cellular variation.

High-throughput single-cell transcriptomics offers an unbiased approach for understanding the extent, basis and function of gene expression variation between seemingly identical cells. Here we sequence single-cell RNA-seq libraries prepared from over 1,700 primary mouse bone-marrow-derived dendritic cells spanning several experimental conditions. We find substantial variation between identically stimulated dendritic cells, in both the fraction of cells detectably expressing a given messenger RNA and the transcript's level within expressing cells. Distinct gene modules are characterized by different temporal heterogeneity profiles. In particular, a 'core' module of antiviral genes is expressed very early by a few 'precocious' cells in response to uniform stimulation with a pathogenic component, but is later activated in all cells. By stimulating cells individually in sealed microfluidic chambers, analysing dendritic cells from knockout mice, and modulating secretion and extracellular signalling, we show that this response is coordinated by interferon-mediated paracrine signalling from these precocious cells. Notably, preventing cell-to-cell communication also substantially reduces variability between cells in the expression of an early-induced 'peaked' inflammatory module, suggesting that paracrine signalling additionally represses part of the inflammatory program. Our study highlights the importance of cell-to-cell communication in controlling cellular heterogeneity and reveals general strategies that multicellular populations can use to establish complex dynamic responses.

Somatic mutation as a mechanism of Wnt/ß-catenin pathway activation in CLL.

One major goal of cancer genome sequencing is to identify key genes and pathways that drive tumor pathogenesis. Although many studies have identified candidate driver genes based on recurrence of mutations in individual genes, subsets of genes with nonrecurrent mutations may also be defined as putative drivers if they affect a single biological pathway. In this fashion, we previously identified Wnt signaling as significantly mutated through large-scale massively parallel DNA sequencing of chronic lymphocytic leukemia (CLL). Here, we use a novel method of biomolecule delivery, vertical silicon nanowires, to efficiently introduce small interfering RNAs into CLL cells, and interrogate the effects of 8 of 15 mutated Wnt pathway members identified across 91 CLLs. In HEK293T cells, mutations in 2 genes did not generate functional changes, 3 led to dysregulated pathway activation, and 3 led to further activation or loss of repression of pathway activation. Silencing 4 of 8 mutated genes in CLL samples harboring the mutated alleles resulted in reduced viability compared with leukemia samples with wild-type alleles. We demonstrate that somatic mutations in CLL can generate dependence on this pathway for survival. These findings support the notion that nonrecurrent mutations at different nodes of the Wnt pathway can contribute to leukemogenesis.

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells.

Recent molecular studies have shown that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels and phenotypic output, with important functional consequences. Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs or proteins simultaneously, because genomic profiling methods could not be applied to single cells until very recently. Here we use single-cell RNA sequencing to investigate heterogeneity in the response of mouse bone-marrow-derived dendritic cells (BMDCs) to lipopolysaccharide. We find extensive, and previously unobserved, bimodal variation in messenger RNA abundance and splicing patterns, which we validate by RNA-fluorescence in situ hybridization for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, we identify a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, we show that variability in this module may be propagated through an interferon feedback circuit, involving the transcriptional regulators Stat2 and Irf7. Our study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.

Dynamic regulatory network controlling TH17 cell differentiation.

Despite their importance, the molecular circuits that control the differentiation of naive T cells remain largely unknown. Recent studies that reconstructed regulatory networks in mammalian cells have focused on short-term responses and relied on perturbation-based approaches that cannot be readily applied to primary T cells. Here we combine transcriptional profiling at high temporal resolution, novel computational algorithms, and innovative nanowire-based perturbation tools to systematically derive and experimentally validate a model of the dynamic regulatory network that controls the differentiation of mouse TH17 cells, a proinflammatory T-cell subset that has been implicated in the pathogenesis of multiple autoimmune diseases. The TH17 transcriptional network consists of two self-reinforcing, but mutually antagonistic, modules, with 12 novel regulators, the coupled action of which may be essential for maintaining the balance between TH17 and other CD4(+) T-cell subsets. Our study identifies and validates 39 regulatory factors, embeds them within a comprehensive temporal network and reveals its organizational principles; it also highlights novel drug targets for controlling TH17 cell differentiation.

Probing enzymatic activity inside living cells using a nanowire-cell "sandwich" assay.

Developing a detailed understanding of enzyme function in the context of an intracellular signal transduction pathway requires minimally invasive methods for probing enzyme activity in situ. Here, we describe a new method for monitoring enzyme activity in living cells by sandwiching live cells between two vertical silicon nanowire (NW) arrays. Specifically, we use the first NW array to immobilize the cells and then present enzymatic substrates intracellularly via the second NW array by utilizing the NWs' ability to penetrate cellular membranes without affecting cells' viability or function. This strategy, when coupled with fluorescence microscopy and mass spectrometry, enables intracellular examination of protease, phosphatase, and protein kinase activities, demonstrating the assay's potential in uncovering the physiological roles of various enzymes.

Nanowire-mediated delivery enables functional interrogation of primary immune cells: application to the analysis of chronic lymphocytic leukemia.

A circuit level understanding of immune cells and hematological cancers has been severely impeded by a lack of techniques that enable intracellular perturbation without significantly altering cell viability and function. Here, we demonstrate that vertical silicon nanowires (NWs) enable gene-specific manipulation of diverse murine and human immune cells with negligible toxicity. To illustrate the power of the technique, we then apply NW-mediated gene silencing to investigate the role of the Wnt signaling pathway in chronic lymphocytic leukemia (CLL). Remarkably, CLL-B cells from different patients exhibit tremendous heterogeneity in their response to the knockdown of a single gene, LEF1. This functional heterogeneity defines three distinct patient groups not discernible by conventional CLL cytogenetic markers and provides a prognostic indicator for patients' time to first therapy. Analyses of gene expression signatures associated with these functional patient subgroups reveal unique insights into the underlying molecular basis for disease heterogeneity. Overall, our findings suggest a functional classification that can potentially guide the selection of patient-specific therapies in CLL and highlight the opportunities for nanotechnology to drive biological inquiry.