PU.1 and BCL11B sequentially cooperate with RUNX1 to anchor mSWI/SNF to poise the T cell effector landscape.

Adaptive immunity relies on specialized effector functions elicited by lymphocytes, yet how antigen recognition activates appropriate effector responses through nonspecific signaling intermediates is unclear. Here we examined the role of chromatin priming in specifying the functional outputs of effector T cells and found that most of the cis-regulatory landscape active in effector T cells was poised early in development before the expression of the T cell antigen receptor. We identified two principal mechanisms underpinning this poised landscape: the recruitment of the nucleosome remodeler mammalian SWItch/Sucrose Non-Fermentable (mSWI/SNF) by the transcription factors RUNX1 and PU.1 to establish chromatin accessibility at T effector loci; and a 'relay' whereby the transcription factor BCL11B succeeded PU.1 to maintain occupancy of the chromatin remodeling complex mSWI/SNF together with RUNX1, after PU.1 silencing during lineage commitment. These mechanisms define modes by which T cells acquire the potential to elicit specialized effector functions early in their ontogeny and underscore the importance of integrating extrinsic cues to the developmentally specified intrinsic program.

Single-cell expression analyses during cellular reprogramming reveal an early stochastic and a late hierarchic phase.

During cellular reprogramming, only a small fraction of cells become induced pluripotent stem cells (iPSCs). Previous analyses of gene expression during reprogramming were based on populations of cells, impeding single-cell level identification of reprogramming events. We utilized two gene expression technologies to profile 48 genes in single cells at various stages during the reprogramming process. Analysis of early stages revealed considerable variation in gene expression between cells in contrast to late stages. Expression of Esrrb, Utf1, Lin28, and Dppa2 is a better predictor for cells to progress into iPSCs than expression of the previously suggested reprogramming markers Fbxo15, Fgf4, and Oct4. Stochastic gene expression early in reprogramming is followed by a late hierarchical phase with Sox2 being the upstream factor in a gene expression hierarchy. Finally, downstream factors derived from the late phase, which do not include Oct4, Sox2, Klf4, c-Myc, and Nanog, can activate the pluripotency circuitry.

Chromatin accessibility and the regulatory epigenome.

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.

Alleviating Cell Lysate-Induced Inhibition to Enable RT-PCR from Single Cells in Picoliter-Volume Double Emulsion Droplets.

Microfluidic droplet assays enable single-cell polymerase chain reaction (PCR) and sequencing analyses at unprecedented scales, with most methods encapsulating cells within nanoliter-sized single emulsion droplets (water-in-oil). Encapsulating cells within picoliter double emulsion (DE) (water-in-oil-in-water) allows sorting droplets with commercially available fluorescence-activated cell sorter (FACS) machines, making it possible to isolate single cells based on phenotypes of interest for downstream analyses. However, sorting DE droplets with standard cytometers requires small droplets that can pass FACS nozzles. This poses challenges for molecular biology, as prior reports suggest that reverse transcription (RT) and PCR amplification cannot proceed efficiently at volumes below 1 nL due to cell lysate-induced inhibition. To overcome this limitation, we used a plate-based RT-PCR assay designed to mimic reactions in picoliter droplets to systematically quantify and ameliorate the inhibition. We find that RT-PCR is blocked by lysate-induced cleavage of nucleic acid probes and primers, which can be efficiently alleviated through heat lysis. We further show that the magnitude of inhibition depends on the cell type, but that RT-PCR can proceed in low-picoscale reaction volumes for most mouse and human cell lines tested. Finally, we demonstrate one-step RT-PCR from single cells in 20 pL DE droplets with fluorescence quantifiable via FACS. These results open up new avenues for improving picoscale droplet RT-PCR reactions and expanding microfluidic droplet-based single-cell analysis technologies.

Transcriptional profiling of cells sorted by RNA abundance.

We have developed a quantitative technique for sorting cells on the basis of endogenous RNA abundance, with a molecular resolution of 10-20 transcripts. We demonstrate efficient and unbiased RNA extraction from transcriptionally sorted cells and report a high-fidelity transcriptome measurement of mouse induced pluripotent stem cells (iPSCs) isolated from a heterogeneous reprogramming culture. This method is broadly applicable to profiling transcriptionally distinct cellular states without requiring antibodies or transgenic fluorescent proteins.

A versatile genome-scale PCR-based pipeline for high-definition DNA FISH.

We developed a cost-effective genome-scale PCR-based method for high-definition DNA FISH (HD-FISH). We visualized gene loci with diffraction-limited resolution, chromosomes as spot clusters and single genes together with transcripts by combining HD-FISH with single-molecule RNA FISH. We provide a database of over 4.3 million primer pairs targeting the human and mouse genomes that is readily usable for rapid and flexible generation of probes.

Malaria-driven expansion of adaptive-like functional CD56-negative NK cells correlates with clinical immunity to malaria.

Natural killer (NK) cells likely play an important role in immunity to malaria, but the effect of repeated malaria on NK cell responses remains unclear. Here, we comprehensively profiled the NK cell response in a cohort of 264 Ugandan children. Repeated malaria exposure was associated with expansion of an atypical, CD56neg population of NK cells that differed transcriptionally, epigenetically, and phenotypically from CD56dim NK cells, including decreased expression of PLZF and the Fc receptor γ-chain, increased histone methylation, and increased protein expression of LAG-3, KIR, and LILRB1. CD56neg NK cells were highly functional and displayed greater antibody-dependent cellular cytotoxicity than CD56dim NK cells. Higher frequencies of CD56neg NK cells were associated with protection against symptomatic malaria and high parasite densities. After marked reductions in malaria transmission, frequencies of these cells rapidly declined, suggesting that continuous exposure to Plasmodium falciparum is required to maintain this modified, adaptive-like NK cell subset.

Single-cell transcriptomic analysis of the adult mouse spinal cord reveals molecular diversity of autonomic and skeletal motor neurons.

The spinal cord is a fascinating structure that is responsible for coordinating movement in vertebrates. Spinal motor neurons control muscle activity by transmitting signals from the spinal cord to diverse peripheral targets. In this study, we profiled 43,890 single-nucleus transcriptomes from the adult mouse spinal cord using fluorescence-activated nuclei sorting to enrich for motor neuron nuclei. We identified 16 sympathetic motor neuron clusters, which are distinguishable by spatial localization and expression of neuromodulatory signaling genes. We found surprising skeletal motor neuron heterogeneity in the adult spinal cord, including transcriptional differences that correlate with electrophysiologically and spatially distinct motor pools. We also provide evidence for a novel transcriptional subpopulation of skeletal motor neuron (γ*). Collectively, these data provide a single-cell transcriptional atlas ( http://spinalcordatlas.org ) for investigating the organizing molecular logic of adult motor neuron diversity, as well as the cellular and molecular basis of motor neuron function in health and disease.

Double emulsion flow cytometry with high-throughput single droplet isolation and nucleic acid recovery.

Droplet microfluidics has made large impacts in diverse areas such as enzyme evolution, chemical product screening, polymer engineering, and single-cell analysis. However, while droplet reactions have become increasingly sophisticated, phenotyping droplets by a fluorescent signal and sorting them to isolate individual variants-of-interest at high-throughput remains challenging. Here, we present sdDE-FACS (s[combining low line]ingle d[combining low line]roplet D[combining low line]ouble E[combining low line]mulsion-FACS), a new method that uses a standard flow cytometer to phenotype, select, and isolate individual double emulsion droplets of interest. Using a 130 µm nozzle at high sort frequency (12-14 kHz), we demonstrate detection of droplet fluorescence signals with a dynamic range spanning 5 orders of magnitude and robust post-sort recovery of intact double emulsion (DE) droplets using 2 commercially-available FACS instruments. We report the first demonstration of single double emulsion droplet isolation with post-sort recovery efficiencies >70%, equivalent to the capabilities of single-cell FACS. Finally, we establish complete downstream recovery of nucleic acids from single, sorted double emulsion droplets via qPCR with little to no cross-contamination. sdDE-FACS marries the full power of droplet microfluidics with flow cytometry to enable a variety of new droplet assays, including rare variant isolation and multiparameter single-cell analysis.

Uncovering a macrophage transcriptional program by integrating evidence from motif scanning and expression dynamics.

Macrophages are versatile immune cells that can detect a variety of pathogen-associated molecular patterns through their Toll-like receptors (TLRs). In response to microbial challenge, the TLR-stimulated macrophage undergoes an activation program controlled by a dynamically inducible transcriptional regulatory network. Mapping a complex mammalian transcriptional network poses significant challenges and requires the integration of multiple experimental data types. In this work, we inferred a transcriptional network underlying TLR-stimulated murine macrophage activation. Microarray-based expression profiling and transcription factor binding site motif scanning were used to infer a network of associations between transcription factor genes and clusters of co-expressed target genes. The time-lagged correlation was used to analyze temporal expression data in order to identify potential causal influences in the network. A novel statistical test was developed to assess the significance of the time-lagged correlation. Several associations in the resulting inferred network were validated using targeted ChIP-on-chip experiments. The network incorporates known regulators and gives insight into the transcriptional control of macrophage activation. Our analysis identified a novel regulator (TGIF1) that may have a role in macrophage activation.

Single-cell multiomic analysis identifies regulatory programs in mixed-phenotype acute leukemia.

Identifying the causes of human diseases requires deconvolution of abnormal molecular phenotypes spanning DNA accessibility, gene expression and protein abundance1-3. We present a single-cell framework that integrates highly multiplexed protein quantification, transcriptome profiling and analysis of chromatin accessibility. Using this approach, we establish a normal epigenetic baseline for healthy blood development, which we then use to deconvolve aberrant molecular features within blood from patients with mixed-phenotype acute leukemia4,5. Despite widespread epigenetic heterogeneity within the patient cohort, we observe common malignant signatures across patients as well as patient-specific regulatory features that are shared across phenotypic compartments of individual patients. Integrative analysis of transcriptomic and chromatin-accessibility maps identified 91,601 putative peak-to-gene linkages and transcription factors that regulate leukemia-specific genes, such as RUNX1-linked regulatory elements proximal to the marker gene CD69. These results demonstrate how integrative, multiomic analysis of single cells within the framework of normal development can reveal both distinct and shared molecular mechanisms of disease from patient samples.

High-throughput chromatin accessibility profiling at single-cell resolution.

Here we develop a high-throughput single-cell ATAC-seq (assay for transposition of accessible chromatin) method to measure physical access to DNA in whole cells. Our approach integrates fluorescence imaging and addressable reagent deposition across a massively parallel (5184) nano-well array, yielding a nearly 20-fold improvement in throughput (up to ~1800 cells/chip, 4-5 h on-chip processing time) and library preparation cost (~81¢ per cell) compared to prior microfluidic implementations. We apply this method to measure regulatory variation in peripheral blood mononuclear cells (PBMCs) and show robust, de novo clustering of single cells by hematopoietic cell type.

Gene expression. MicroRNA control of protein expression noise.

MicroRNAs (miRNAs) repress the expression of many genes in metazoans by accelerating messenger RNA degradation and inhibiting translation, thereby reducing the level of protein. However, miRNAs only slightly reduce the mean expression of most targeted proteins, leading to speculation about their role in the variability, or noise, of protein expression. We used mathematical modeling and single-cell reporter assays to show that miRNAs, in conjunction with increased transcription, decrease protein expression noise for lowly expressed genes but increase noise for highly expressed genes. Genes that are regulated by multiple miRNAs show more-pronounced noise reduction. We estimate that hundreds of (lowly expressed) genes in mouse embryonic stem cells have reduced noise due to substantial miRNA regulation. Our findings suggest that miRNAs confer precision to protein expression and thus offer plausible explanations for the commonly observed combinatorial targeting of endogenous genes by multiple miRNAs, as well as the preferential targeting of lowly expressed genes.