Single-cell chromatin accessibility and transposable element landscapes reveal shared features of tissue-residing immune cells.

Tissue adaptation is required for regulatory T (Treg) cell function within organs. Whether this program shares aspects with other tissue-localized immune populations is unclear. Here, we analyzed single-cell chromatin accessibility data, including the transposable element (TE) landscape of CD45+ immune cells from colon, skin, adipose tissue, and spleen. We identified features of organ-specific tissue adaptation across different immune cells. Focusing on tissue Treg cells, we found conservation of the Treg tissue adaptation program in other tissue-localized immune cells, such as amphiregulin-producing T helper (Th)17 cells. Accessible TEs can act as regulatory elements, but their contribution to tissue adaptation is not understood. TE landscape analysis revealed an enrichment of specific transcription factor binding motifs in TE regions within accessible chromatin peaks. TEs, specifically from the LTR family, were located in enhancer regions and associated with tissue adaptation. These findings broaden our understanding of immune tissue residency and provide an important step toward organ-specific immune interventions.

Integrated single-cell profiling dissects cell-state-specific enhancer landscapes of human tumor-infiltrating CD8+ T cells.

Despite extensive studies on the chromatin landscape of exhausted T cells, the transcriptional wiring underlying the heterogeneous functional and dysfunctional states of human tumor-infiltrating lymphocytes (TILs) is incompletely understood. Here, we identify gene-regulatory landscapes in a wide breadth of functional and dysfunctional CD8+ TIL states covering four cancer entities using single-cell chromatin profiling. We map enhancer-promoter interactions in human TILs by integrating single-cell chromatin accessibility with single-cell RNA-seq data from tumor-entity-matching samples and prioritize cell-state-specific genes by super-enhancer analysis. Besides revealing entity-specific chromatin remodeling in exhausted TILs, our analyses identify a common chromatin trajectory to TIL dysfunction and determine key enhancers, transcriptional regulators, and deregulated genes involved in this process. Finally, we validate enhancer regulation at immunotherapeutically relevant loci by targeting non-coding regulatory elements with potent CRISPR activators and repressors. In summary, our study provides a framework for understanding and manipulating cell-state-specific gene-regulatory cues from human tumor-infiltrating lymphocytes.

Single-cell chromatin accessibility landscape identifies tissue repair program in human regulatory T cells.

Murine regulatory T (Treg) cells in tissues promote tissue homeostasis and regeneration. We sought to identify features that characterize human Treg cells with these functions in healthy tissues. Single-cell chromatin accessibility profiles of murine and human tissue Treg cells defined a conserved, microbiota-independent tissue-repair Treg signature with a prevailing footprint of the transcription factor BATF. This signature, combined with gene expression profiling and TCR fate mapping, identified a population of tissue-like Treg cells in human peripheral blood that expressed BATF, chemokine receptor CCR8 and HLA-DR. Human BATF+CCR8+ Treg cells from normal skin and adipose tissue shared features with nonlymphoid T follicular helper-like (Tfh-like) cells, and induction of a Tfh-like differentiation program in naive human Treg cells partially recapitulated tissue Treg regenerative characteristics, including wound healing potential. Human BATF+CCR8+ Treg cells from healthy tissue share features with tumor-resident Treg cells, highlighting the importance of understanding the context-specific functions of these cells.

Precursors for Nonlymphoid-Tissue Treg Cells Reside in Secondary Lymphoid Organs and Are Programmed by the Transcription Factor BATF.

Specialized regulatory T (Treg) cells accumulate and perform homeostatic and regenerative functions in nonlymphoid tissues. Whether common precursors for nonlymphoid-tissue Treg cells exist and how they differentiate remain elusive. Using transcription factor nuclear factor, interleukin 3 regulated (Nfil3) reporter mice and single-cell RNA-sequencing (scRNA-seq), we identified two precursor stages of interleukin 33 (IL-33) receptor ST2-expressing nonlymphoid tissue Treg cells, which resided in the spleen and lymph nodes. Global chromatin profiling of nonlymphoid tissue Treg cells and the two precursor stages revealed a stepwise acquisition of chromatin accessibility and reprogramming toward the nonlymphoid-tissue Treg cell phenotype. Mechanistically, we identified and validated the transcription factor Batf as the driver of the molecular tissue program in the precursors. Understanding this tissue development program will help to harness regenerative properties of tissue Treg cells for therapy.

Structural cells are key regulators of organ-specific immune responses.

The mammalian immune system implements a remarkably effective set of mechanisms for fighting pathogens1. Its main components are haematopoietic immune cells, including myeloid cells that control innate immunity, and lymphoid cells that constitute adaptive immunity2. However, immune functions are not unique to haematopoietic cells, and many other cell types display basic mechanisms of pathogen defence3-5. To advance our understanding of immunology outside the haematopoietic system, here we systematically investigate the regulation of immune genes in the three major types of structural cells: epithelium, endothelium and fibroblasts. We characterize these cell types across twelve organs in mice, using cellular phenotyping, transcriptome sequencing, chromatin accessibility profiling and epigenome mapping. This comprehensive dataset revealed complex immune gene activity and regulation in structural cells. The observed patterns were highly organ-specific and seem to modulate the extensive interactions between structural cells and haematopoietic immune cells. Moreover, we identified an epigenetically encoded immune potential in structural cells under tissue homeostasis, which was triggered in response to systemic viral infection. This study highlights the prevalence and organ-specific complexity of immune gene activity in non-haematopoietic structural cells, and it provides a high-resolution, multi-omics atlas of the epigenetic and transcriptional networks that regulate structural cells in the mouse.

Genome-wide cooperation of EMT transcription factor ZEB1 with YAP and AP-1 in breast cancer.

Invasion, metastasis and therapy resistance are the major cause of cancer-associated deaths, and the EMT-inducing transcription factor ZEB1 is a crucial stimulator of these processes. While work on ZEB1 has mainly focused on its role as a transcriptional repressor, it can also act as a transcriptional activator. To further understand these two modes of action, we performed a genome-wide ZEB1 binding study in triple-negative breast cancer cells. We identified ZEB1 as a novel interactor of the AP-1 factors FOSL1 and JUN and show that, together with the Hippo pathway effector YAP, they form a transactivation complex, predominantly activating tumour-promoting genes, thereby synergising with its function as a repressor of epithelial genes. High expression of ZEB1, YAP, FOSL1 and JUN marks the aggressive claudin-low subtype of breast cancer, indicating the translational relevance of our findings. Thus, our results link critical tumour-promoting transcription factors: ZEB1, AP-1 and Hippo pathway factors. Disturbing their molecular interaction may provide a promising treatment option for aggressive cancer types.

Targeting clinical epigenetic reprogramming for chemoprevention of metabolic and viral hepatocellular carcinoma.

OBJECTIVE: Hepatocellular carcinoma (HCC) is the fastest-growing cause of cancer-related mortality with chronic viral hepatitis and non-alcoholic steatohepatitis (NASH) as major aetiologies. Treatment options for HCC are unsatisfactory and chemopreventive approaches are absent. Chronic hepatitis C (CHC) results in epigenetic alterations driving HCC risk and persisting following cure. Here, we aimed to investigate epigenetic modifications as targets for liver cancer chemoprevention. DESIGN: Liver tissues from patients with NASH and CHC were analysed by ChIP-Seq (H3K27ac) and RNA-Seq. The liver disease-specific epigenetic and transcriptional reprogramming in patients was modelled in a liver cell culture system. Perturbation studies combined with a targeted small molecule screen followed by in vivo and ex vivo validation were used to identify chromatin modifiers and readers for HCC chemoprevention. RESULTS: In patients, CHC and NASH share similar epigenetic and transcriptomic modifications driving cancer risk. Using a cell-based system modelling epigenetic modifications in patients, we identified chromatin readers as targets to revert liver gene transcription driving clinical HCC risk. Proof-of-concept studies in a NASH-HCC mouse model showed that the pharmacological inhibition of chromatin reader bromodomain 4 inhibited liver disease progression and hepatocarcinogenesis by restoring transcriptional reprogramming of the genes that were epigenetically altered in patients. CONCLUSION: Our results unravel the functional relevance of metabolic and virus-induced epigenetic alterations for pathogenesis of HCC development and identify chromatin readers as targets for chemoprevention in patients with chronic liver diseases.

Combined chemosensitivity and chromatin profiling prioritizes drug combinations in CLL.

The Bruton tyrosine kinase (BTK) inhibitor ibrutinib has substantially improved therapeutic options for chronic lymphocytic leukemia (CLL). Although ibrutinib is not curative, it has a profound effect on CLL cells and may create new pharmacologically exploitable vulnerabilities. To identify such vulnerabilities, we developed a systematic approach that combines epigenome profiling (charting the gene-regulatory basis of cell state) with single-cell chemosensitivity profiling (quantifying cell-type-specific drug response) and bioinformatic data integration. By applying our method to a cohort of matched patient samples collected before and during ibrutinib therapy, we identified characteristic ibrutinib-induced changes that provide a starting point for the rational design of ibrutinib combination therapies. Specifically, we observed and validated preferential sensitivity to proteasome, PLK1, and mTOR inhibitors during ibrutinib treatment. More generally, our study establishes a broadly applicable method for investigating treatment-specific vulnerabilities by integrating the complementary perspectives of epigenetic cell states and phenotypic drug responses in primary patient samples.

HCV-Induced Epigenetic Changes Associated With Liver Cancer Risk Persist After Sustained Virologic Response.

BACKGROUND & AIMS: Chronic hepatitis C virus (HCV) infection is an important risk factor for hepatocellular carcinoma (HCC). Despite effective antiviral therapies, the risk for HCC is decreased but not eliminated after a sustained virologic response (SVR) to direct-acting antiviral (DAA) agents, and the risk is higher in patients with advanced fibrosis. We investigated HCV-induced epigenetic alterations that might affect risk for HCC after DAA treatment in patients and mice with humanized livers. METHODS: We performed genome-wide ChIPmentation-based ChIP-Seq and RNA-seq analyses of liver tissues from 6 patients without HCV infection (controls), 18 patients with chronic HCV infection, 8 patients with chronic HCV infection cured by DAA treatment, 13 patients with chronic HCV infection cured by interferon therapy, 4 patients with chronic hepatitis B virus infection, and 7 patients with nonalcoholic steatohepatitis in Europe and Japan. HCV-induced epigenetic modifications were mapped by comparative analyses with modifications associated with other liver disease etiologies. uPA/SCID mice were engrafted with human hepatocytes to create mice with humanized livers and given injections of HCV-infected serum samples from patients; mice were given DAAs to eradicate the virus. Pathways associated with HCC risk were identified by integrative pathway analyses and validated in analyses of paired HCC tissues from 8 patients with an SVR to DAA treatment of HCV infection. RESULTS: We found chronic HCV infection to induce specific genome-wide changes in H3K27ac, which correlated with changes in expression of mRNAs and proteins. These changes persisted after an SVR to DAAs or interferon-based therapies. Integrative pathway analyses of liver tissues from patients and mice with humanized livers demonstrated that HCV-induced epigenetic alterations were associated with liver cancer risk. Computational analyses associated increased expression of SPHK1 with HCC risk. We validated these findings in an independent cohort of patients with HCV-related cirrhosis (n = 216), a subset of which (n = 21) achieved viral clearance. CONCLUSIONS: In an analysis of liver tissues from patients with and without an SVR to DAA therapy, we identified epigenetic and gene expression alterations associated with risk for HCC. These alterations might be targeted to prevent liver cancer in patients treated for HCV infection.

Pooled CRISPR screening with single-cell transcriptome readout.

CRISPR-based genetic screens are accelerating biological discovery, but current methods have inherent limitations. Widely used pooled screens are restricted to simple readouts including cell proliferation and sortable marker proteins. Arrayed screens allow for comprehensive molecular readouts such as transcriptome profiling, but at much lower throughput. Here we combine pooled CRISPR screening with single-cell RNA sequencing into a broadly applicable workflow, directly linking guide RNA expression to transcriptome responses in thousands of individual cells. Our method for CRISPR droplet sequencing (CROP-seq) enables pooled CRISPR screens with single-cell transcriptome resolution, which will facilitate high-throughput functional dissection of complex regulatory mechanisms and heterogeneous cell populations.

A promoter-level mammalian expression atlas.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

An atlas of active enhancers across human cell types and tissues.

Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.

High-throughput ChIPmentation: freely scalable, single day ChIPseq data generation from very low cell-numbers.

BACKGROUND: Chromatin immunoprecipitation coupled to sequencing (ChIP-seq) is widely used to map histone modifications and transcription factor binding on a genome-wide level. RESULTS: We present high-throughput ChIPmentation (HT-ChIPmentation) that eliminates the need for DNA purification prior to library amplification and reduces reverse-crosslinking time from hours to minutes. CONCLUSIONS: The resulting workflow is easily established, extremely rapid, and compatible with requirements for very low numbers of FACS sorted cells, high-throughput applications and single day data generation.

Epigenetic mechanisms regulating T-cell responses.

During the last decade, advances in sequencing technologies allowed production of a wealth of information on epigenetic modifications in T cells. Epigenome maps, in combination with mechanistic studies, have demonstrated that T cells undergo extensive epigenome remodeling in response to signals, which has a strong effect on phenotypic stability and function of lymphocytes. In this review we focus on DNA methylation, histone modifications, and chromatin structure as important epigenetic mechanisms involved in controlling T-cell responses. In particular, we discuss epigenetic processes in light of the development, activation, and differentiation of CD4+ T helper (TH), regulatory T, and CD8+ T cells. As central aspects of the adaptive immune system, we review mechanisms that ensure molecular memory, stability, plasticity, and exhaustion of T cells. We further discuss the effect of the tissue environment on imprinting T-cell epigenomes with potential implications for immunotherapy.

ChIPmentation: fast, robust, low-input ChIP-seq for histones and transcription factors.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is widely used to map histone marks and transcription factor binding throughout the genome. Here we present ChIPmentation, a method that combines chromatin immunoprecipitation with sequencing library preparation by Tn5 transposase ('tagmentation'). ChIPmentation introduces sequencing-compatible adaptors in a single-step reaction directly on bead-bound chromatin, which reduces time, cost and input requirements, thus providing a convenient and broadly useful alternative to existing ChIP-seq protocols.

Genome-wide identification of hypoxia-inducible factor-1 and -2 binding sites in hypoxic human macrophages alternatively activated by IL-10.

Macrophages (MΦ) often accumulate in hypoxic areas, where they significantly influence disease progression. Anti-inflammatory cytokines, such as IL-10, generate alternatively activated macrophages that support tumor growth. To understand how alternative activation affects the transcriptional profile of hypoxic macrophages, we globally mapped binding sites of hypoxia-inducible factor (HIF)-1α and HIF-2α in primary human monocyte-derived macrophages prestimulated with IL-10. 713 HIF-1 and 795 HIF-2 binding sites were identified under hypoxia. Pretreatment with IL-10 altered the binding pattern, with 120 new HIF-1 and 188 new HIF-2 binding sites emerging. HIF-1 binding was most prominent in promoters, while HIF-2 binding was more abundant in enhancer regions. Comparison of ChIP-seq data obtained in other cells revealed a highly cell type specific binding of HIF. In MΦ HIF binding occurred preferentially in already active enhancers or promoters. To assess the roles of HIF on gene expression, primary human macrophages were treated with siRNA against HIF-1α or HIF-2α, followed by genome-wide gene expression analysis. Comparing mRNA expression to the HIF binding profile revealed a significant enrichment of hypoxia-inducible genes previously identified by ChIP-seq. Analysis of gene expression under hypoxia alone and hypoxia/IL-10 showed the enhanced induction of a set of genes including PLOD2 and SLC2A3, while another group including KDM3A and ADM remained unaffected or was reduced by IL-10. Taken together IL-10 influences the DNA binding pattern of HIF and the level of gene induction.

The enhancer and promoter landscape of human regulatory and conventional T-cell subpopulations.

CD4(+)CD25(+)FOXP3(+) human regulatory T cells (Tregs) are essential for self-tolerance and immune homeostasis. Here, we describe the promoterome of CD4(+)CD25(high)CD45RA(+) naïve and CD4(+)CD25(high)CD45RA(-) memory Tregs and their CD25(-) conventional T-cell (Tconv) counterparts both before and after in vitro expansion by cap analysis of gene expression (CAGE) adapted to single-molecule sequencing (HeliScopeCAGE). We performed comprehensive comparative digital gene expression analyses and revealed novel transcription start sites, of which several were validated as alternative promoters of known genes. For all in vitro expanded subsets, we additionally generated global maps of poised and active enhancer elements marked by histone H3 lysine 4 monomethylation and histone H3 lysine 27 acetylation, describe their cell type-specific motif signatures, and evaluate the role of candidate transcription factors STAT5, FOXP3, RUNX1, and ETS1 in both Treg- and Tconv-specific enhancer architectures. Network analyses of gene expression data revealed additional candidate transcription factors contributing to cell type specificity and a transcription factor network in Tregs that is dominated by FOXP3 interaction partners and targets. In summary, we provide a comprehensive and easily accessible resource of gene expression and gene regulation in human Treg and Tconv subpopulations.

Transcription and enhancer profiling in human monocyte subsets.

Human blood monocytes comprise at least 3 subpopulations that differ in phenotype and function. Here, we present the first in-depth regulome analysis of human classical (CD14(++)CD16(-)), intermediate (CD14(+)CD16(+)), and nonclassical (CD14(dim)CD16(+)) monocytes. Cap analysis of gene expression adapted to Helicos single-molecule sequencing was used to map transcription start sites throughout the genome in all 3 subsets. In addition, global maps of H3K4me1 and H3K27ac deposition were generated for classical and nonclassical monocytes defining enhanceosomes of the 2 major subsets. We identified differential regulatory elements (including promoters and putative enhancers) that were associated with subset-specific motif signatures corresponding to different transcription factor activities and exemplarily validated novel downstream enhancer elements at the CD14 locus. In addition to known subset-specific features, pathway analysis revealed marked differences in metabolic gene signatures. Whereas classical monocytes expressed higher levels of genes involved in carbohydrate metabolism, priming them for anaerobic energy production, nonclassical monocytes expressed higher levels of oxidative pathway components and showed a higher mitochondrial routine activity. Our findings describe promoter/enhancer landscapes and provide novel insights into the specific biology of human monocyte subsets.

Comprehensive genome and epigenome characterization of CHO cells in response to evolutionary pressures and over time.

The most striking characteristic of CHO cells is their adaptability, which enables efficient production of proteins as well as growth under a variety of culture conditions, but also results in genomic and phenotypic instability. To investigate the relative contribution of genomic and epigenetic modifications towards phenotype evolution, comprehensive genome and epigenome data are presented for six related CHO cell lines, both in response to perturbations (different culture conditions and media as well as selection of a specific phenotype with increased transient productivity) and in steady state (prolonged time in culture under constant conditions). Clear transitions were observed in DNA-methylation patterns upon each perturbation, while few changes occurred over time under constant conditions. Only minor DNA-methylation changes were observed between exponential and stationary growth phase; however, throughout a batch culture the histone modification pattern underwent continuous adaptation. Variation in genome sequence between the six cell lines on the level of SNPs, InDels, and structural variants is high, both upon perturbation and under constant conditions over time. The here presented comprehensive resource may open the door to improved control and manipulation of gene expression during industrial bioprocesses based on epigenetic mechanisms. Biotechnol. Bioeng. 2016;113: 2241-2253. © 2016 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.

Mechanisms of in vivo binding site selection of the hematopoietic master transcription factor PU.1.

The transcription factor PU.1 is crucial for the development of many hematopoietic lineages and its binding patterns significantly change during differentiation processes. However, the 'rules' for binding or not-binding of potential binding sites are only partially understood. To unveil basic characteristics of PU.1 binding site selection in different cell types, we studied the binding properties of PU.1 during human macrophage differentiation. Using in vivo and in vitro binding assays, as well as computational prediction, we show that PU.1 selects its binding sites primarily based on sequence affinity, which results in the frequent autonomous binding of high affinity sites in DNase I inaccessible regions (25-45% of all occupied sites). Increasing PU.1 concentrations and the availability of cooperative transcription factor interactions during lineage differentiation both decrease affinity thresholds for in vivo binding and fine-tune cell type-specific PU.1 binding, which seems to be largely independent of DNA methylation. Occupied sites were predominantly detected in active chromatin domains, which are characterized by higher densities of PU.1 recognition sites and neighboring motifs for cooperative transcription factors. Our study supports a model of PU.1 binding control that involves motif-binding affinity, PU.1 concentration, cooperativeness with neighboring transcription factor sites and chromatin domain accessibility, which likely applies to all PU.1 expressing cells.