Increasing Fluid Viscosity Ensures Consistent Single-Cell Encapsulation.

High-throughput single-cell analysis typically relies on the isolation of cells of interest in separate compartments for subsequent phenotypic or genotypic characterization. Using microfluidics, this is achieved by isolating individual cells in microdroplets or microwells. However, due to cell-to-cell variability in size, shape, and density, the cell capture efficiencies may vary significantly. This variability can negatively impact the measurements and introduce undesirable artifacts when trying to isolate and characterize heterogeneous cell populations. In this study, we show that single-cell isolation biases in microfluidics can be circumvented by increasing the viscosity of fluids in which cells are dispersed. At a viscosity of 40-50 cP (cP), the cell sedimentation is effectively reduced, resulting in a steady cell flow inside the microfluidics chip and consistent encapsulation in water-in-oil droplets over extended periods of time. This approach allows nearly all cells in a sample to be isolated with the same efficiency, irrespective of their type. Our results show that increased fluid viscosity, rather than cell-adjusted density, provides a more reliable approach to mitigate single-cell isolation biases.

Activation of NLRP3 Inflammasome by Virus-Like Particles of Human Polyomaviruses in Macrophages.

Viral antigens can activate phagocytes, inducing inflammation, but the mechanisms are barely explored. The aim of this study is to investigate how viral oligomeric proteins of different structures induce inflammatory response in macrophages. Human THP-1 cell line was used to prepare macrophages that were treated with filamentous nucleocapsid-like particles (NLPs) of paramyxoviruses and spherical virus-like particles (VLPs) of human polyomaviruses. The effects of viral proteins on cell viability, pro-inflammatory cytokines' production, and NLRP3 inflammasome activation were investigated. Filamentous NLPs did not induce inflammation while spherical VLPs mediated inflammatory response followed by NLRP3 inflammasome activation. Inhibitors of cathepsins and K+ efflux decreased IL-1ß release and cell death, indicating a complex inflammasome activation process. A similar activation pattern was observed in primary human macrophages. Single-cell RNAseq analysis of THP-1 cells revealed several cell activation states different in inflammation-related genes. This study provides new insights into the interaction of viral proteins with immune cells and suggests that structural properties of oligomeric proteins may define cell activation pathways.

Macrophage-Targeted Therapy Unlocks Antitumoral Cross-talk between IFNγ-Secreting Lymphocytes and IL12-Producing Dendritic Cells.

Macrophages often abound within tumors, express colony-stimulating factor 1 receptor (CSF1R), and are linked to adverse patient survival. Drugs blocking CSF1R signaling have been used to suppress tumor-promoting macrophage responses; however, their mechanisms of action remain incompletely understood. Here, we assessed the lung tumor immune microenvironment in mice treated with BLZ945, a prototypical small-molecule CSF1R inhibitor, using single-cell RNA sequencing and mechanistic validation approaches. We showed that tumor control was not caused by CSF1R+ cell depletion; instead, CSF1R targeting reshaped the CSF1R+ cell landscape, which unlocked cross-talk between antitumoral CSF1R- cells. These cells included IFNγ-producing natural killer and T cells, and an IL12-producing dendritic cell subset, denoted as DC3, which were all necessary for CSF1R inhibitor-mediated lung tumor control. These data indicate that CSF1R targeting can activate a cardinal cross-talk between cells that are not macrophages and that are essential to mediate the effects of T cell-targeted immunotherapies and promote antitumor immunity.See related Spotlight by Burrello and de Visser, p. 4.

Resident Kupffer cells and neutrophils drive liver toxicity in cancer immunotherapy.

Immunotherapy is revolutionizing cancer treatment but is often restricted by toxicities. What distinguishes adverse events from concomitant antitumor reactions is poorly understood. Here, using anti-CD40 treatment in mice as a model of TH1-promoting immunotherapy, we showed that liver macrophages promoted local immune-related adverse events. Mechanistically, tissue-resident Kupffer cells mediated liver toxicity by sensing lymphocyte-derived IFN-γ and subsequently producing IL-12. Conversely, dendritic cells were dispensable for toxicity but drove tumor control. IL-12 and IFN-γ were not toxic themselves but prompted a neutrophil response that determined the severity of tissue damage. We observed activation of similar inflammatory pathways after anti-PD-1 and anti-CTLA-4 immunotherapies in mice and humans. These findings implicated macrophages and neutrophils as mediators and effectors of aberrant inflammation in TH1-promoting immunotherapy, suggesting distinct mechanisms of toxicity and antitumor immunity.

Identification of a targetable KRAS-mutant epithelial population in non-small cell lung cancer.

Lung cancer is the leading cause of cancer deaths. Tumor heterogeneity, which hampers development of targeted therapies, was herein deconvoluted via single cell RNA sequencing in aggressive human adenocarcinomas (carrying Kras-mutations) and comparable murine model. We identified a tumor-specific, mutant-KRAS-associated subpopulation which is conserved in both human and murine lung cancer. We previously reported a key role for the oncogene BMI-1 in adenocarcinomas. We therefore investigated the effects of in vivo PTC596 treatment, which affects BMI-1 activity, in our murine model. Post-treatment, MRI analysis showed decreased tumor size, while single cell transcriptomics concomitantly detected near complete ablation of the mutant-KRAS-associated subpopulation, signifying the presence of a pharmacologically targetable, tumor-associated subpopulation. Our findings therefore hold promise for the development of a targeted therapy for KRAS-mutant adenocarcinomas.

Tumor-Promoting Ly-6G+ SiglecFhigh Cells Are Mature and Long-Lived Neutrophils.

Myeloid cells co-expressing the markers CD11b, Ly-6G, and SiglecF can be found in large numbers in murine lung adenocarcinomas and accelerate cancer growth by fostering tumor cell invasion, angiogenesis, and immunosuppression; however, some of these cells' fundamental features remain unexplored. Here, we show that tumor-infiltrating CD11b+ Ly-6G+ SiglecFhigh cells are bona fide mature neutrophils and therefore differ from other myeloid cells, including SiglecFhigh eosinophils, SiglecFhigh macrophages, and CD11b+ Ly-6G+ myeloid-derived suppressor cells. We further show that SiglecFhigh neutrophils gradually accumulate in growing tumors, where they can live for several days; this lifespan is in marked contrast to that of their SiglecFlow counterparts and neutrophils in general, which live for several hours only. Together, these findings reveal distinct attributes for tumor-promoting SiglecFhigh neutrophils and help explain their deleterious accumulation in the tumor bed.

A large pool of actively cycling progenitors orchestrates self-renewal and injury repair of an ectodermal appendage.

The classical model of tissue renewal posits that small numbers of quiescent stem cells (SCs) give rise to proliferating transit-amplifying cells before terminal differentiation. However, many organs house pools of SCs with proliferative and differentiation potentials that diverge from this template. Resolving SC identity and organization is therefore central to understanding tissue renewal. Here, using a combination of single-cell RNA sequencing (scRNA-seq), mouse genetics and tissue injury approaches, we uncover cellular hierarchies and mechanisms that underlie the maintenance and repair of the continuously growing mouse incisor. Our results reveal that, during homeostasis, a group of actively cycling epithelial progenitors generates enamel-producing ameloblasts and adjacent layers of non-ameloblast cells. After injury, tissue repair was achieved through transient increases in progenitor-cell proliferation and through direct conversion of Notch1-expressing cells to ameloblasts. We elucidate epithelial SC identity, position and function, providing a mechanistic basis for the homeostasis and repair of a fast-turnover ectodermal appendage.

Single-Cell Transcriptomics of Human and Mouse Lung Cancers Reveals Conserved Myeloid Populations across Individuals and Species.

Tumor-infiltrating myeloid cells (TIMs) comprise monocytes, macrophages, dendritic cells, and neutrophils, and have emerged as key regulators of cancer growth. These cells can diversify into a spectrum of states, which might promote or limit tumor outgrowth but remain poorly understood. Here, we used single-cell RNA sequencing (scRNA-seq) to map TIMs in non-small-cell lung cancer patients. We uncovered 25 TIM states, most of which were reproducibly found across patients. To facilitate translational research of these populations, we also profiled TIMs in mice. In comparing TIMs across species, we identified a near-complete congruence of population structures among dendritic cells and monocytes; conserved neutrophil subsets; and species differences among macrophages. By contrast, myeloid cell population structures in patients' blood showed limited overlap with those of TIMs. This study determines the lung TIM landscape and sets the stage for future investigations into the potential of TIMs as immunotherapy targets.

Quantitative Clonal Analysis and Single-Cell Transcriptomics Reveal Division Kinetics, Hierarchy, and Fate of Oral Epithelial Progenitor Cells.

The oral mucosa is one of the most rapidly dividing tissues in the body and serves as a barrier to physical and chemical insults from mastication, food, and microorganisms. Breakdown of this barrier can lead to significant morbidity and potentially life-threatening infections for patients. Determining the identity and organization of oral epithelial progenitor cells (OEPCs) is therefore paramount to understanding their roles in homeostasis and disease. Using lineage tracing and label retention experiments, we show that rapidly dividing OEPCs are located broadly within the basal layer of the mucosa throughout the oral cavity. Quantitative clonal analysis demonstrated that OEPCs undergo population-asymmetrical divisions following neutral drift dynamics and that they respond to chemotherapy-induced damage by altering daughter cell fates. Finally, using single-cell RNA-seq, we establish the basal layer population structure and propose a model that defines the organization of cells within the basal layer.

Molecular diversification of regulatory T cells in nonlymphoid tissues.

Foxp3+CD4+ regulatory T cells (Tregs) accumulate in certain nonlymphoid tissues, where they control diverse aspects of organ homeostasis. Populations of tissue Tregs, as they have been termed, have transcriptomes distinct from those of their counterparts in lymphoid organs and other nonlymphoid tissues. We examined the diversification of Tregs in visceral adipose tissue, skeletal muscle, and the colon vis-à-vis lymphoid organs from the same individuals. The unique transcriptomes of the various tissue Treg populations resulted from layering of tissue-restricted open chromatin regions over regions already open in the spleen, the latter tagged by super-enhancers and particular histone marks. The binding motifs for a small number of transcription factor (TF) families were repeatedly enriched within the accessible chromatin stretches of Tregs in the three nonlymphoid tissues. However, a bioinformatically and experimentally validated transcriptional network, constructed by integrating chromatin accessibility and single-cell transcriptomic data, predicted reliance on different TF family members in the different tissues. The network analysis also revealed that tissue-restricted and broadly acting TFs were integrated into feed-forward loops to enforce tissue-specific gene expression in nonlymphoid-tissue Tregs. Overall, this study provides a framework for understanding the epigenetic dynamics of T cells operating in nonlymphoid tissues, which should inform strategies for specifically targeting them.

A single-cell atlas of the airway epithelium reveals the CFTR-rich pulmonary ionocyte.

The functions of epithelial tissues are dictated by the types, abundance and distribution of the differentiated cells they contain. Attempts to restore tissue function after damage require knowledge of how physiological tasks are distributed among cell types, and how cell states vary between homeostasis, injury-repair and disease. In the conducting airway, a heterogeneous basal cell population gives rise to specialized luminal cells that perform mucociliary clearance1. Here we perform single-cell profiling of human bronchial epithelial cells and mouse tracheal epithelial cells to obtain a comprehensive census of cell types in the conducting airway and their behaviour in homeostasis and regeneration. Our analysis reveals cell states that represent known and novel cell populations, delineates their heterogeneity and identifies distinct differentiation trajectories during homeostasis and tissue repair. Finally, we identified a novel, rare cell type that we call the 'pulmonary ionocyte', which co-expresses FOXI1, multiple subunits of the vacuolar-type H+-ATPase (V-ATPase) and CFTR, the gene that is mutated in cystic fibrosis. Using immunofluorescence, modulation of signalling pathways and electrophysiology, we show that Notch signalling is necessary and FOXI1 expression is sufficient to drive the production of the pulmonary ionocyte, and that the pulmonary ionocyte is a major source of CFTR activity in the conducting airway epithelium.

Publisher Correction: Single-cell gene expression reveals a landscape of regulatory T cell phenotypes shaped by the TCR.

In the version of this article initially published, the Supplementary Note was missing. The Supplementary Note has now been provided online and is cited in the Methods section of the article. The error has been corrected in the HTML and PDF version of the article.

Publisher Correction: Single-cell analysis of experience-dependent transcriptomic states in the mouse visual cortex.

In the version of this article initially published, the x-axis labels in Fig. 3c read Vglut, Gad1/2, Aldh1l1 and Pecam1; they should have read Vglut+, Gad1/2+, Aldh1l1+ and Pecam1+. In Fig. 4, the range values were missing from the color scales; they are, from left to right, 4-15, 0-15, 4-15 and 0-15 in Fig. 4a and 4-15, 4-15 and 4-8 in Fig. 4h. In the third paragraph of the main text, the phrase reading "Previous approaches have analyzed a limited number of inhibitory cell types, thus masking the full diversity of excitatory populations" should have read "Previous approaches have analyzed a limited number of inhibitory cell types and masked the full diversity of excitatory populations." In the second paragraph of Results section "Diversity of experience-regulated ERGs," the phrase reading "thus suggesting considerable divergence within the gene expression program responding to early stimuli" should have read "thus suggesting considerable divergence within the early stimulus-responsive gene expression program." In the fourth paragraph of Results section "Excitatory neuronal LRGs," the sentence reading "The anatomical organization of these cell types into sublayers, coupled with divergent transcriptional responses to a sensory stimulus, suggested previously unappreciated functional subdivisions located within the laminae of the mouse visual cortex and resembling the cytoarchitecture in higher mammals" should have read "The anatomical organization of these cell types into sublayers, coupled with divergent transcriptional responses to a sensory stimulus, suggests previously unappreciated functional subdivisions located within the laminae of the mouse visual cortex, resembling the cytoarchitecture in higher mammals." In the last sentence of the Results, "sensory-responsive genes" should have read "sensory-stimulus-responsive genes." The errors have been corrected in the HTML and PDF versions of the article.

Single-cell gene expression reveals a landscape of regulatory T cell phenotypes shaped by the TCR.

CD4+ T regulatory cells (Treg) are central to immune homeostasis, their phenotypic heterogeneity reflecting the diverse environments and target cells that they regulate. To understand this heterogeneity, we combined single-cell RNA-seq, activation reporter and T cell receptor (TCR) analysis to profile thousands of Treg or conventional CD4+FoxP3- T cells (Tconv) from mouse lymphoid organs and human blood. Treg and Tconv pools showed areas of overlap, as resting 'furtive' Tregs with overall similarity to Tconvs or as a convergence of activated states. All Tregs expressed a small core of FoxP3-dependent transcripts, onto which additional programs were added less uniformly. Among suppressive functions, Il2ra and Ctla4 were quasiconstant, inhibitory cytokines being more sparsely distributed. TCR signal intensity did not affect resting/activated Treg proportions but molded activated Treg programs. The main lines of Treg heterogeneity in mice were strikingly conserved in human blood. These results reveal unexpected TCR-shaped states of activation, providing a framework to synthesize previous observations of Treg heterogeneity.

Population snapshots predict early haematopoietic and erythroid hierarchies.

The formation of red blood cells begins with the differentiation of multipotent haematopoietic progenitors. Reconstructing the steps of this differentiation represents a general challenge in stem-cell biology. Here we used single-cell transcriptomics, fate assays and a theory that allows the prediction of cell fates from population snapshots to demonstrate that mouse haematopoietic progenitors differentiate through a continuous, hierarchical structure into seven blood lineages. We uncovered coupling between the erythroid and the basophil or mast cell fates, a global haematopoietic response to erythroid stress and novel growth factor receptors that regulate erythropoiesis. We defined a flow cytometry sorting strategy to purify early stages of erythroid differentiation, completely isolating classically defined burst-forming and colony-forming progenitors. We also found that the cell cycle is progressively remodelled during erythroid development and during a sharp transcriptional switch that ends the colony-forming progenitor stage and activates terminal differentiation. Our work showcases the utility of linking transcriptomic data to predictive fate models, and provides insights into lineage development in vivo.

Single-cell analysis of experience-dependent transcriptomic states in the mouse visual cortex.

Activity-dependent transcriptional responses shape cortical function. However, a comprehensive understanding of the diversity of these responses across the full range of cortical cell types, and how these changes contribute to neuronal plasticity and disease, is lacking. To investigate the breadth of transcriptional changes that occur across cell types in the mouse visual cortex after exposure to light, we applied high-throughput single-cell RNA sequencing. We identified significant and divergent transcriptional responses to stimulation in each of the 30 cell types characterized, thus revealing 611 stimulus-responsive genes. Excitatory pyramidal neurons exhibited inter- and intralaminar heterogeneity in the induction of stimulus-responsive genes. Non-neuronal cells showed clear transcriptional responses that may regulate experience-dependent changes in neurovascular coupling and myelination. Together, these results reveal the dynamic landscape of the stimulus-dependent transcriptional changes occurring across cell types in the visual cortex; these changes are probably critical for cortical function and may be sites of deregulation in developmental brain disorders.

Osteoblasts remotely supply lung tumors with cancer-promoting SiglecFhigh neutrophils.

Bone marrow-derived myeloid cells can accumulate within tumors and foster cancer outgrowth. Local immune-neoplastic interactions have been intensively investigated, but the contribution of the systemic host environment to tumor growth remains poorly understood. Here, we show in mice and cancer patients (n = 70) that lung adenocarcinomas increase bone stromal activity in the absence of bone metastasis. Animal studies reveal that the cancer-induced bone phenotype involves bone-resident osteocalcin-expressing (Ocn+) osteoblastic cells. These cells promote cancer by remotely supplying a distinct subset of tumor-infiltrating SiglecFhigh neutrophils, which exhibit cancer-promoting properties. Experimentally reducing Ocn+ cell numbers suppresses the neutrophil response and lung tumor outgrowth. These observations posit osteoblasts as remote regulators of lung cancer and identify SiglecFhigh neutrophils as myeloid cell effectors of the osteoblast-driven protumoral response.

Single-cell barcoding and sequencing using droplet microfluidics.

Single-cell RNA sequencing has recently emerged as a powerful tool for mapping cellular heterogeneity in diseased and healthy tissues, yet high-throughput methods are needed for capturing the unbiased diversity of cells. Droplet microfluidics is among the most promising candidates for capturing and processing thousands of individual cells for whole-transcriptome or genomic analysis in a massively parallel manner with minimal reagent use. We recently established a method called inDrops, which has the capability to index >15,000 cells in an hour. A suspension of cells is first encapsulated into nanoliter droplets with hydrogel beads (HBs) bearing barcoding DNA primers. Cells are then lysed and mRNA is barcoded (indexed) by a reverse transcription (RT) reaction. Here we provide details for (i) establishing an inDrops platform (1 d); (ii) performing hydrogel bead synthesis (4 d); (iii) encapsulating and barcoding cells (1 d); and (iv) RNA-seq library preparation (2 d). inDrops is a robust and scalable platform, and it is unique in its ability to capture and profile >75% of cells in even very small samples, on a scale of thousands or tens of thousands of cells.

End Sequence Analysis Toolkit (ESAT) expands the extractable information from single-cell RNA-seq data.

RNA-seq protocols that focus on transcript termini are well suited for applications in which template quantity is limiting. Here we show that, when applied to end-sequencing data, analytical methods designed for global RNA-seq produce computational artifacts. To remedy this, we created the End Sequence Analysis Toolkit (ESAT). As a test, we first compared end-sequencing and bulk RNA-seq using RNA from dendritic cells stimulated with lipopolysaccharide (LPS). As predicted by the telescripting model for transcriptional bursts, ESAT detected an LPS-stimulated shift to shorter 3'-isoforms that was not evident by conventional computational methods. Then, droplet-based microfluidics was used to generate 1000 cDNA libraries, each from an individual pancreatic islet cell. ESAT identified nine distinct cell types, three distinct ß-cell types, and a complex interplay between hormone secretion and vascularization. ESAT, then, offers a much-needed and generally applicable computational pipeline for either bulk or single-cell RNA end-sequencing.

DNA Nanoparticles for Improved Protein Synthesis In Vitro.

The amplification and digital quantification of single DNA molecules are important in biomedicine and diagnostics. Beyond quantifying DNA molecules in a sample, the ability to express proteins from the amplified DNA would open even broader applications in synthetic biology, directed evolution, and proteomics. Herein, a microfluidic approach is reported for the production of condensed DNA nanoparticles that can serve as efficient templates for in vitro protein synthesis. Using phi29 DNA polymerase and a multiple displacement amplification reaction, single DNA molecules were converted into DNA nanoparticles containing up to about 10(4)  clonal gene copies of the starting template. DNA nanoparticle formation was triggered by accumulation of inorganic pyrophosphate (produced during DNA synthesis) and magnesium ions from the buffer. Transcription-translation reactions performed in vitro showed that individual DNA nanoparticles can serve as efficient templates for protein synthesis in vitro.