Jellert T. Gaublomme

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.

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.

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.

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.

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 Peyers 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.

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

Kishi et al. find that protein C receptor (PROCR) is specifically expressed on the surface of Th17 cells and its loss exacerbates encephalitogenic Th17 cell responses.

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.