Systematic discovery and perturbation of regulatory genes in human T cells reveals the architecture of immune networks.

Gene regulatory networks ensure that important genes are expressed at precise levels. When gene expression is sufficiently perturbed, it can lead to disease. To understand how gene expression disruptions percolate through a network, we must first map connections between regulatory genes and their downstream targets. However, we lack comprehensive knowledge of the upstream regulators of most genes. Here, we developed an approach for systematic discovery of upstream regulators of critical immune factors-IL2RA, IL-2 and CTLA4-in primary human T cells. Then, we mapped the network of the target genes of these regulators and putative cis-regulatory elements using CRISPR perturbations, RNA-seq and ATAC-seq. These regulators form densely interconnected networks with extensive feedback loops. Furthermore, this network is enriched for immune-associated disease variants and genes. These results provide insight into how immune-associated disease genes are regulated in T cells and broader principles about the structure of human gene regulatory networks.

NEAT-seq: simultaneous profiling of intra-nuclear proteins, chromatin accessibility and gene expression in single cells.

In this work, we describe NEAT-seq (sequencing of nuclear protein epitope abundance, chromatin accessibility and the transcriptome in single cells), enabling interrogation of regulatory mechanisms spanning the central dogma. We apply this technique to profile CD4 memory T cells using a panel of master transcription factors (TFs) that drive T cell subsets and identify examples of TFs with regulatory activity gated by transcription, translation and regulation of chromatin binding. We also link a noncoding genome-wide association study single-nucleotide polymorphism (SNP) within a GATA motif to a putative target gene, using NEAT-seq data to internally validate SNP impact on GATA3 regulation.

The dynamic, combinatorial cis-regulatory lexicon of epidermal differentiation.

Transcription factors bind DNA sequence motif vocabularies in cis-regulatory elements (CREs) to modulate chromatin state and gene expression during cell state transitions. A quantitative understanding of how motif lexicons influence dynamic regulatory activity has been elusive due to the combinatorial nature of the cis-regulatory code. To address this, we undertook multiomic data profiling of chromatin and expression dynamics across epidermal differentiation to identify 40,103 dynamic CREs associated with 3,609 dynamically expressed genes, then applied an interpretable deep-learning framework to model the cis-regulatory logic of chromatin accessibility. This analysis framework identified cooperative DNA sequence rules in dynamic CREs regulating synchronous gene modules with diverse roles in skin differentiation. Massively parallel reporter assay analysis validated temporal dynamics and cooperative cis-regulatory logic. Variants linked to human polygenic skin disease were enriched in these time-dependent combinatorial motif rules. This integrative approach shows the combinatorial cis-regulatory lexicon of epidermal differentiation and represents a general framework for deciphering the organizational principles of the cis-regulatory code of dynamic gene regulation.

Integrated single-cell transcriptomics and epigenomics reveals strong germinal center-associated etiology of autoimmune risk loci.

The germinal center (GC) response is critical for both effective adaptive immunity and establishing peripheral tolerance by limiting autoreactive B cells. Dysfunction in these processes can lead to defective immune responses to infection or contribute to autoimmune disease. To understand the gene regulatory principles underlying the GC response, we generated a single-cell transcriptomic and epigenomic atlas of the human tonsil, a widely studied and representative lymphoid tissue. We characterize diverse immune cell subsets and build a trajectory of dynamic gene expression and transcription factor activity during B cell activation, GC formation, and plasma cell differentiation. We subsequently leverage cell type­specific transcriptomic and epigenomic maps to interpret potential regulatory impact of genetic variants implicated in autoimmunity, revealing that many exhibit their greatest regulatory potential in GC-associated cellular populations. These included gene loci linked with known roles in GC biology (IL21, IL21R, IL4R, and BCL6) and transcription factors regulating B cell differentiation (POU2AF1 and HHEX). Together, these analyses provide a powerful new cell type­resolved resource for the interpretation of cellular and genetic causes underpinning autoimmune disease.

Chromatin and gene-regulatory dynamics of the developing human cerebral cortex at single-cell resolution.

Genetic perturbations of cortical development can lead to neurodevelopmental disease, including autism spectrum disorder (ASD). To identify genomic regions crucial to corticogenesis, we mapped the activity of gene-regulatory elements generating a single-cell atlas of gene expression and chromatin accessibility both independently and jointly. This revealed waves of gene regulation by key transcription factors (TFs) across a nearly continuous differentiation trajectory, distinguished the expression programs of glial lineages, and identified lineage-determining TFs that exhibited strong correlation between linked gene-regulatory elements and expression levels. These highly connected genes adopted an active chromatin state in early differentiating cells, consistent with lineage commitment. Base-pair-resolution neural network models identified strong cell-type-specific enrichment of noncoding mutations predicted to be disruptive in a cohort of ASD individuals and identified frequently disrupted TF binding sites. This approach illustrates how cell-type-specific mapping can provide insights into the programs governing human development and disease.

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.

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.

Landscape of stimulation-responsive chromatin across diverse human immune cells.

A hallmark of the immune system is the interplay among specialized cell types transitioning between resting and stimulated states. The gene regulatory landscape of this dynamic system has not been fully characterized in human cells. Here we collected assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing data under resting and stimulated conditions for up to 32 immune cell populations. Stimulation caused widespread chromatin remodeling, including response elements shared between stimulated B and T cells. Furthermore, several autoimmune traits showed significant heritability in stimulation-responsive elements from distinct cell types, highlighting the importance of these cell states in autoimmunity. Allele-specific read mapping identified variants that alter chromatin accessibility in particular conditions, allowing us to observe evidence of function for a candidate causal variant that is undetected by existing large-scale studies in resting cells. Our results provide a resource of chromatin dynamics and highlight the need to characterize the effects of genetic variation in stimulated cells.

BLIMP1 Induces Transient Metastatic Heterogeneity in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most metastatic and deadly cancers. Despite the clinical significance of metastatic spread, our understanding of molecular mechanisms that drive PDAC metastatic ability remains limited. By generating a genetically engineered mouse model of human PDAC, we uncover a transient subpopulation of cancer cells with exceptionally high metastatic ability. Global gene expression profiling and functional analyses uncovered the transcription factor BLIMP1 as a driver of PDAC metastasis. The highly metastatic PDAC subpopulation is enriched for hypoxia-induced genes, and hypoxia-mediated induction of BLIMP1 contributes to the regulation of a subset of hypoxia-associated gene expression programs. These findings support a model in which upregulation of BLIMP1 links microenvironmental cues to a metastatic stem cell character.Significance: PDAC is an almost uniformly lethal cancer, largely due to its tendency for metastasis. We define a highly metastatic subpopulation of cancer cells, uncover a key transcriptional regulator of metastatic ability, and define hypoxia as an important factor within the tumor microenvironment that increases metastatic proclivity. Cancer Discov; 7(10); 1184-99. ©2017 AACR.See related commentary by Vakoc and Tuveson, p. 1067This article is highlighted in the In This Issue feature, p. 1047.

An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues.

We present Omni-ATAC, an improved ATAC-seq protocol for chromatin accessibility profiling that works across multiple applications with substantial improvement of signal-to-background ratio and information content. The Omni-ATAC protocol generates chromatin accessibility profiles from archival frozen tissue samples and 50-µm sections, revealing the activities of disease-associated DNA elements in distinct human brain structures. The Omni-ATAC protocol enables the interrogation of personal regulomes in tissue context and translational studies.

Single-cell epigenomic variability reveals functional cancer heterogeneity.

BACKGROUND: Cell-to-cell heterogeneity is a major driver of cancer evolution, progression, and emergence of drug resistance. Epigenomic variation at the single-cell level can rapidly create cancer heterogeneity but is difficult to detect and assess functionally. RESULTS: We develop a strategy to bridge the gap between measurement and function in single-cell epigenomics. Using single-cell chromatin accessibility and RNA-seq data in K562 leukemic cells, we identify the cell surface marker CD24 as co-varying with chromatin accessibility changes linked to GATA transcription factors in single cells. Fluorescence-activated cell sorting of CD24 high versus low cells prospectively isolated GATA1 and GATA2 high versus low cells. GATA high versus low cells express differential gene regulatory networks, differential sensitivity to the drug imatinib mesylate, and differential self-renewal capacity. Lineage tracing experiments show that GATA/CD24hi cells have the capability to rapidly reconstitute the heterogeneity within the entire starting population, suggesting that GATA expression levels drive a phenotypically relevant source of epigenomic plasticity. CONCLUSION: Single-cell chromatin accessibility can guide prospective characterization of cancer heterogeneity. Epigenomic subpopulations in cancer impact drug sensitivity and the clonal dynamics of cancer evolution.

Nrf2 is not required for epithelial prohibitin-dependent attenuation of experimental colitis.

Inflammatory bowel disease is associated with increased reactive oxygen species (ROS) and decreased antioxidant response in the intestinal mucosa. Expression of the mitochondrial protein prohibitin (PHB) is also decreased during intestinal inflammation. Our previous study showed that genetic restoration of colonic epithelial PHB expression [villin-PHB transgenic (PHB Tg) mice] attenuated dextran sodium sulfate (DSS)-induced colitis/oxidative stress and sustained expression of colonic nuclear factor erythroid 2-related factor 2 (Nrf2), a cytoprotective transcription factor. This study investigated the role of Nrf2 in mediating PHB-induced protection against colitis and expression of the antioxidant response element (ARE)-regulated antioxidant genes heme oxygenase-1 (HO-1) and NAD(P)H quinone oxidoreductase-1 (NQO-1). PHB-transfected Caco-2-BBE human intestinal epithelial cells maintained increased ARE activation and decreased intracellular ROS levels compared with control vector-transfected cells during Nrf2 knockdown by small interfering RNA. Treatment with the ERK inhibitor PD-98059 decreased PHB-induced ARE activation, suggesting that ERK constitutes a significant portion of PHB-mediated ARE activation in Caco-2-BBE cells. PHB Tg, Nrf2(-/-), and PHB Tg/Nrf2(-/-) mice were treated with DSS or 2,4,6-trinitrobenzene sulfonic acid (TNBS), and inflammation and expression of HO-1 and NQO-1 were assessed. PHB Tg/Nrf2(-/-) mice mimicked PHB Tg mice, with attenuated DSS- or TNBS-induced colitis and induction of colonic HO-1 and NQO-1 expression, despite deletion of Nrf2. PHB Tg/Nrf2(-/-) mice exhibited increased activation of ERK during colitis. Our results suggest that maintaining expression of intestinal epithelial cell PHB, which is decreased during colitis, reduces the severity of inflammation and increases colonic levels of the antioxidants HO-1 and NQO-1 via a mechanism independent of Nrf2.

Prohibitin attenuates colitis-associated tumorigenesis in mice by modulating p53 and STAT3 apoptotic responses.

Although inflammatory bowel disease is associated with higher risk of colorectal cancer, the precise pathogenic mechanisms underlying this association are not completely understood. Prohibitin 1 (PHB), a protein implicated in the regulation of proliferation, apoptosis, and transcription, is decreased in intestinal inflammation. In this study, we have established a key function for PHB in mediating colitis-associated cancer. Wild-type and transgenic (Tg) mice specifically overexpressing PHB in intestinal epithelial cells were subjected to a classical two-stage protocol of colitis-associated carcinogenesis. In addition, wild-type and p53 null human cell models were used to assess PHB interaction with STAT3 and p53. Wild-type mice exhibited decreased mucosal PHB protein expression during colitis-associated carcinogenesis. Tg mice exhibited decreased susceptibility in a manner associated with increased apoptosis, p53, Bax, and Bad expression plus decreased Bcl-xL and Bcl-2 expression. PHB overexpression in wild-type but not p53 null human cells increased expression of Bax, Bad, and caspase-3 cleavage. In wild-type p53 cells, PHB overexpression decreased basal and interleukin-6-induced STAT3 activation and expression of the STAT3 responsive genes Bcl-xL and Bcl-2. PHB coimmunoprecipitated with phospho-STAT3 in addition to p53 in cultured cell lysates and colon mucosa. This is the first study to show interaction between PHB and STAT3 in vivo. In summary, our findings suggest that PHB protects against colitis-associated cancer by modulating p53- and STAT3-mediated apoptosis. Modulation of PHB expression in intestinal epithelial cells may offer a potential therapeutic approach to prevent colitis-associated carcinogenesis.

Prohibitin 1 modulates mitochondrial stress-related autophagy in human colonic epithelial cells.

INTRODUCTION: Autophagy is an adaptive response to extracellular and intracellular stress by which cytoplasmic components and organelles, including damaged mitochondria, are degraded to promote cell survival and restore cell homeostasis. Certain genes involved in autophagy confer susceptibility to Crohn's disease. Reactive oxygen species and pro-inflammatory cytokines such as tumor necrosis factor α (TNFα), both of which are increased during active inflammatory bowel disease, promote cellular injury and autophagy via mitochondrial damage. Prohibitin (PHB), which plays a role in maintaining normal mitochondrial respiratory function, is decreased during active inflammatory bowel disease. Restoration of colonic epithelial PHB expression protects mice from experimental colitis and combats oxidative stress. In this study, we investigated the potential role of PHB in modulating mitochondrial stress-related autophagy in intestinal epithelial cells. METHODS: We measured autophagy activation in response to knockdown of PHB expression by RNA interference in Caco2-BBE and HCT116 WT and p53 null cells. The effect of exogenous PHB expression on TNFα- and IFNγ-induced autophagy was assessed. Autophagy was inhibited using Bafilomycin A(1) or siATG16L1 during PHB knockdown and the affect on intracellular oxidative stress, mitochondrial membrane potential, and cell viability were determined. The requirement of intracellular ROS in siPHB-induced autophagy was assessed using the ROS scavenger N-acetyl-L-cysteine. RESULTS: TNFα and IFNγ-induced autophagy inversely correlated with PHB protein expression. Exogenous PHB expression reduced basal autophagy and TNFα-induced autophagy. Gene silencing of PHB in epithelial cells induces mitochondrial autophagy via increased intracellular ROS. Inhibition of autophagy during PHB knockdown exacerbates mitochondrial depolarization and reduces cell viability. CONCLUSIONS: Decreased PHB levels coupled with dysfunctional autophagy renders intestinal epithelial cells susceptible to mitochondrial damage and cytotoxicity. Repletion of PHB may represent a therapeutic approach to combat oxidant and cytokine-induced mitochondrial damage in diseases such as inflammatory bowel disease.