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.

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.

Changes in adenoviral chromatin organization precede early gene activation upon infection.

Within the virion, adenovirus DNA associates with the virus-encoded, protamine-like structural protein pVII. Whether this association is organized, and how genome packaging changes during infection and subsequent transcriptional activation is currently unclear. Here, we combined RNA-seq, MNase-seq, ChIP-seq, and single genome imaging during early adenovirus infection to unveil the structure- and time-resolved dynamics of viral chromatin changes as well as their correlation with gene transcription. Our MNase mapping data indicates that the adenoviral genome is arranged in precisely positioned nucleoprotein particles with nucleosome-like characteristics, that we term adenosomes. We identified 238 adenosomes that are positioned by a DNA sequence code and protect about 60-70 bp of DNA. The incoming adenoviral genome is more accessible at early gene loci that undergo additional chromatin de-condensation upon infection. Histone H3.3 containing nucleosomes specifically replaces pVII at distinct genomic sites and at the transcription start sites of early genes. Acetylation of H3.3 is predominant at the transcription start sites and precedes transcriptional activation. Based on our results, we propose a central role for the viral pVII nucleoprotein architecture, which is required for the dynamic structural changes during early infection, including the regulation of nucleosome assembly prior to transcription initiation. Our study thus may aid the rational development of recombinant adenoviral vectors exhibiting sustained expression in gene therapy.

Engineering an inducible leukemia-associated fusion protein enables large-scale ex vivo production of functional human phagocytes.

Ex vivo expansion of human CD34+ hematopoietic stem and progenitor cells remains a challenge due to rapid differentiation after detachment from the bone marrow niche. In this study, we assessed the capacity of an inducible fusion protein to enable sustained ex vivo proliferation of hematopoietic precursors and their capacity to differentiate into functional phagocytes. We fused the coding sequences of an FK506-Binding Protein 12 (FKBP12)-derived destabilization domain (DD) to the myeloid/lymphoid lineage leukemia/eleven nineteen leukemia (MLL-ENL) fusion gene to generate the fusion protein DD-MLL-ENL and retrovirally expressed the protein switch in human CD34+ progenitors. Using Shield1, a chemical inhibitor of DD fusion protein degradation, we established large-scale and long-term expansion of late monocytic precursors. Upon Shield1 removal, the cells lost self-renewal capacity and spontaneously differentiated, even after 2.5 y of continuous ex vivo expansion. In the absence of Shield1, stimulation with IFN-γ, LPS, and GM-CSF triggered terminal differentiation. Gene expression analysis of the obtained phagocytes revealed marked similarity with naïve monocytes. In functional assays, the novel phagocytes migrated toward CCL2, attached to VCAM-1 under shear stress, produced reactive oxygen species, and engulfed bacterial particles, cellular particles, and apoptotic cells. Finally, we demonstrated Fcγ receptor recognition and phagocytosis of opsonized lymphoma cells in an antibody-dependent manner. Overall, we have established an engineered protein that, as a single factor, is useful for large-scale ex vivo production of human phagocytes. Such adjustable proteins have the potential to be applied as molecular tools to produce functional immune cells for experimental cell-based approaches.

Biosensors for inflammation as a strategy to engineer regulatory T cells for cell therapy.

Engineered regulatory T cell (Treg cell) therapy is a promising strategy to treat patients suffering from inflammatory diseases, autoimmunity, and transplant rejection. However, in many cases, disease-related antigens that can be targeted by Treg cells are not available. In this study, we introduce a class of synthetic biosensors, named artificial immune receptors (AIRs), for murine and human Treg cells. AIRs consist of three domains: (a) extracellular binding domain of a tumor necrosis factor (TNF)-receptor superfamily member, (b) intracellular costimulatory signaling domain of CD28, and (c) T cell receptor signaling domain of CD3-ζ chain. These AIR receptors equip Treg cells with an inflammation-sensing machinery and translate this environmental information into a CD3-ζ chain-dependent TCR-activation program. Different AIRs were generated, recognizing the inflammatory ligands of the TNF-receptor superfamily, including LIGHT, TNFα, and TNF-like ligand 1A (TL1A), leading to activation, differentiation, and proliferation of AIR-Treg cells. In a graft-versus-host disease model, Treg cells expressing lymphotoxin ß receptor-AIR, which can be activated by the ligand LIGHT, protect significantly better than control Treg cells. Expression and signaling of the corresponding human AIR in human Treg cells prove that this concept can be translated. Engineering Treg cells that target inflammatory ligands leading to TCR signaling and activation might be used as a Treg cell-based therapy approach for a broad range of inflammation-driven diseases.

D-2-hydroxyglutarate supports a tolerogenic phenotype with lowered major histocompatibility class II expression in non-malignant dendritic cells and acute myeloid leukemia cells.

D-2-hydroxyglutarate (D-2-HG) accumulates in patients with acute myeloid leukemia (AML) with mutated isocitrate dehydrogenase (IDH) and in other malignancies. D-2-HG suppresses antitumor T-cell immunity but little is known about potential effects on non-malignant myeloid cells. Here we show that D-2-HG impairs human but not murine dendritic cell differentiation, resulting in a tolerogenic phenotype with low major histocompatibility class II expression. In line with this, IDH-mutated AML blasts exhibited lower expression of HLA-DP and were less susceptible to lysis by HLA-DP-specific T cells. Interestingly, besides its expected impact on DNA demethylation, D-2-HG reprogrammed metabolism towards increased lactate production in dendritic cells and AML. Vitamin C accelerated DNA demethylation, but only the combination of vitamin C and glycolytic inhibition lowered lactate levels and supported major histocompatibility complex class II expression. Our results indicate an unexpected link between the immunosuppressive metabolites 2-HG and lactic acid and suggest a potentially novel therapeutic strategy with combinations of anti-glycolytic drugs and epigenetic modulators (hypomethylating agents) or other therapeutics for the treatment of AML.

Allele-specific expression of GATA2 due to epigenetic dysregulation in CEBPA double-mutant AML.

Transcriptional deregulation is a central event in the development of acute myeloid leukemia (AML). To identify potential disturbances in gene regulation, we conducted an unbiased screen of allele-specific expression (ASE) in 209 AML cases. The gene encoding GATA binding protein 2 (GATA2) displayed ASE more often than any other myeloid- or cancer-related gene. GATA2 ASE was strongly associated with CEBPA double mutations (DMs), with 95% of cases presenting GATA2 ASE. In CEBPA DM AML with GATA2 mutations, the mutated allele was preferentially expressed. We found that GATA2 ASE was a somatic event lost in complete remission, supporting the notion that it plays a role in CEBPA DM AML. Acquisition of GATA2 ASE involved silencing of 1 allele via promoter methylation and concurrent overactivation of the other allele, thereby preserving expression levels. Notably, promoter methylation was also lost in remission along with GATA2 ASE. In summary, we propose that GATA2 ASE is acquired by epigenetic mechanisms and is a prerequisite for the development of AML with CEBPA DMs. This finding constitutes a novel example of an epigenetic hit cooperating with a genetic hit in the pathogenesis of AML.

Aberrant DNA methylation patterns in microsatellite stable human colorectal cancers define a new marker panel for the CpG island methylator phenotype.

A distinct group of colorectal carcinomas (CRCs) referred to as the "CpG island methylator phenotype" (CIMP) shows an extremely high incidence of de novo DNA methylation and may share common pathological, clinical or molecular features. However, there is limited consensus about which CpG islands (CGIs) define a CIMP, particularly in microsatellite stable (MSS) carcinomas. To study this phenotype in a systematic manner, we analyzed genome-wide CGI DNA methylation profiles of 19 MSS CRC using methyl-CpG immunoprecipitation (MCIp) and hybridization on 244K CGI oligonucleotide microarrays, determined KRAS and BRAF mutation status and compared disease-related DNA methylation changes to chromosomal instability as detected by microarray-based comparative genomic hybridization. Results were validated using mass spectrometry analysis of bisulfite-converted DNA at a subset of 76 individual CGIs in 120 CRC and 43 matched normal tissue samples. Both genome-wide profiling and CpG methylation fine mapping segregated a group of CRC showing pronounced and frequent de novo DNA methylation of a distinct group of CGIs that only partially overlapped with previously established classifiers. The CIMP group defined in our study revealed significant association with colon localization, either KRAS or BRAF mutation, and mostly minor chromosomal losses but no association with known histopathological features. Our data provide a basis for defining novel marker panels that may enable a more reliable classification of CIMP in all CRCs, independently of the MS status.

Physiological levels of 25-hydroxyvitamin D3 induce a suppressive CD4+ T cell phenotype not reflected in the epigenetic landscape.

1,25-dihydroxyvitamin D3 (1,25(OH)2 D3 ), the active metabolite of vitamin D3 has a strong impact on the differentiation and function of immune cells. Here we analysed the influence of its precursor 25-hydroxyvitamin D3 (25(OH)D3 ) on the differentiation of human CD4+ T cells applying physiological concentrations in vitro. Our data show that 25(OH)D3 is converted to its active form 1,25(OH)2 D3 by T cells, which in turn supports FOXP3, CD25 and CTLA-4 expression and inhibits IFN-γ production. These changes were not reflected in the demethylation of the respective promoters. Furthermore, we investigated the impact of vitamin D3 metabolites under induced Treg (iTreg) polarization conditions using TGF-ß. Surprisingly, no additive effect but a decreased percentage of FOXP3 expressing cells was observed. However, the combination of 25(OH)D3 or 1,25(OH)2 D3 together with TGF-ß further upregulated CD25 and CTLA-4 and significantly increased soluble CTLA-4 and IL-10 secretion whereas IFN-γ expression of iTreg was decreased. Our data suggest that physiological levels of 25(OH)D3 act as potent modulator of human CD4+ T cells and autocrine or paracrine production of 1,25(OH)2 D3 by T cells might be crucial for the local regulation of an adaptive immune response. However, since no epigenetic changes are detected by 25(OH)D3 a rather transient phenotype is induced.

Methyl donor micronutrients, CD40-ligand methylation and disease activity in systemic lupus erythematosus: A cross-sectional association study.

OBJECTIVE: Hypomethylation of CD40-ligand (CD40L) in T-cells is associated with increased disease activity in systemic lupus erythematosus (SLE). We therefore investigated possible associations of dietary methyl donors and products with CD40L methylation status in SLE. METHODS: Food frequency questionnaires were employed to calculate methyl donor micronutrients in 61 female SLE patients (age 45.7 ± 12.0 years, disease duration 16.2 ± 8.4 years) and compared to methylation levels of previously identified key DNA methylation sites (CpG17 and CpG22) within CD40L promotor of T-cells using quantitative DNA methylation analysis on the EpiTYPER mass spectrometry platform. Disease activity was assessed by SLE Disease Activity Index (SLEDAI). Linear regression modelling was used. P values were adjusted according to Benjamini & Hochberg. RESULTS: Amongst the micronutrients assessed (g per day), methionine and cysteine were associated with methylation of CpG17 (ß = 5.0 (95%CI: 0.6-9.4), p = 0.04; and ß = 2.4 (0.6-4.1), p = 0.02, respectively). Methionine, choline, and cysteine were additionally associated with the mean methylation of the entire CD40L (ß = 9.5 (1.0-18.0), p = 0.04; ß = 1.6 (0.4-3.0), p = 0.04; and ß = 4.3 (0.9-7.7), p = 0.02, respectively). Associations of the SLEDAI with hypomethylation were confirmed for CpG17 (ß=-32.6 (-60.6 to -4.6), p = 0.04) and CpG22 (ß=-38.3 (-61.2 to -15.4), p = 0.004), but not the mean methylation of CD40L. Dietary products with the highest impact on methylation included meat, ice cream, white bread, and cooked potatoes. CONCLUSIONS: Dietary methyl donors may influence DNA methylation levels and thereby disease activity in SLE.

Associations of site-specific CD4+-T-cell hypomethylation within CD40-ligand promotor and enhancer regions with disease activity of women with systemic lupus erythematosus.

OBJECTIVE: To comprehensively assess associations of site-specific CD4+-T-cell hypomethylation of the CD40-Ligand gene (CD40L) with disease activity of women with systemic lupus erythematosus (SLE). METHODS: CpG-sites within the DNA of the promotor and two enhancer regions (n = 22) of CD40L were identified and numbered consecutively. The rate of methylated DNA in isolated CD4+-T-cells of women with SLE were quantified for each methylation site by MALDI-TOF. Disease activity was assessed by SLE Disease Activity Index (SLEDAI). Associations of site-specific methylation rates with the SLEDAI scores were assessed by linear regression modelling. P values were adjusted according to Bonferroni-Holm as indicated. RESULTS: 60 female SLE patients participated in the study (age 45.7 ± 11.1 years, disease duration 17.0 ± 8.3 years). Significant associations to the SLEDAI were noted for CpG22 hypomethylation of the promotor (ß = -40.1, p = 0.017, adjusted p = 0.027), trends were noted for CpG17 hypomethylation of the promotor (ß = -30.5, p = 0.032, adjusted p = 0.6), and for CpG11 hypermethylation of the second enhancer (ß = 15.0, p = 0.046, adjusted p = 0.8). CONCLUSION: Site-specific hypomethylation of the CD40L promotor in CD4+-T-cells show associations with disease activity in female SLE patients.

Temporal autoregulation during human PU.1 locus SubTAD formation.

Epigenetic control of gene expression occurs within discrete spatial chromosomal units called topologically associating domains (TADs), but the exact spatial requirements of most genes are unknown; this is of particular interest for genes involved in cancer. We therefore applied high-resolution chromosomal conformation capture sequencing to map the three-dimensional (3D) organization of the human locus encoding the key myeloid transcription factor PU.1 in healthy monocytes and acute myeloid leukemia (AML) cells. We identified a dynamic ∼75-kb unit (SubTAD) as the genomic region in which spatial interactions between PU.1 gene regulatory elements occur during myeloid differentiation and are interrupted in AML. Within this SubTAD, proper initiation of the spatial chromosomal interactions requires PU.1 autoregulation and recruitment of the chromatin-adaptor protein LDB1 (LIM domain-binding protein 1). However, once these spatial interactions have occurred, LDB1 stabilizes them independently of PU.1 autoregulation. Thus, our data support that PU.1 autoregulates its expression in a "hit-and-run" manner by initiating stable chromosomal loops that result in a transcriptionally active chromatin architecture.

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.

Analysis of the human monocyte-derived macrophage transcriptome and response to lipopolysaccharide provides new insights into genetic aetiology of inflammatory bowel disease.

The FANTOM5 consortium utilised cap analysis of gene expression (CAGE) to provide an unprecedented insight into transcriptional regulation in human cells and tissues. In the current study, we have used CAGE-based transcriptional profiling on an extended dense time course of the response of human monocyte-derived macrophages grown in macrophage colony-stimulating factor (CSF1) to bacterial lipopolysaccharide (LPS). We propose that this system provides a model for the differentiation and adaptation of monocytes entering the intestinal lamina propria. The response to LPS is shown to be a cascade of successive waves of transient gene expression extending over at least 48 hours, with hundreds of positive and negative regulatory loops. Promoter analysis using motif activity response analysis (MARA) identified some of the transcription factors likely to be responsible for the temporal profile of transcriptional activation. Each LPS-inducible locus was associated with multiple inducible enhancers, and in each case, transient eRNA transcription at multiple sites detected by CAGE preceded the appearance of promoter-associated transcripts. LPS-inducible long non-coding RNAs were commonly associated with clusters of inducible enhancers. We used these data to re-examine the hundreds of loci associated with susceptibility to inflammatory bowel disease (IBD) in genome-wide association studies. Loci associated with IBD were strongly and specifically (relative to rheumatoid arthritis and unrelated traits) enriched for promoters that were regulated in monocyte differentiation or activation. Amongst previously-identified IBD susceptibility loci, the vast majority contained at least one promoter that was regulated in CSF1-dependent monocyte-macrophage transitions and/or in response to LPS. On this basis, we concluded that IBD loci are strongly-enriched for monocyte-specific genes, and identified at least 134 additional candidate genes associated with IBD susceptibility from reanalysis of published GWA studies. We propose that dysregulation of monocyte adaptation to the environment of the gastrointestinal mucosa is the key process leading to inflammatory bowel disease.

Shared activity patterns arising at genetic susceptibility loci reveal underlying genomic and cellular architecture of human disease.

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

Controlled re-activation of epigenetically silenced Tet promoter-driven transgene expression by targeted demethylation.

Faithful expression of transgenes in cell cultures and mice is often challenged by locus dependent epigenetic silencing. We investigated silencing of Tet-controlled expression cassettes within the mouse ROSA26 locus. We observed pronounced DNA methylation of the Tet promoter concomitant with loss of expression in mES cells as well as in differentiated cells and transgenic animals. Strikingly, the ROSA26 promoter remains active and methylation free indicating that this silencing mechanism specifically affects the transgene, but does not spread to the host's chromosomal neighborhood. To reactivate Tet cassettes a synthetic fusion protein was constructed and expressed in silenced cells. This protein includes the enzymatic domains of ten eleven translocation methylcytosine dioxygenase 1 (TET-1) as well as the Tet repressor DNA binding domain. Expression of the synthetic fusion protein and Doxycycline treatment allowed targeted demethylation of the Tet promoter in the ROSA26 locus and in another genomic site, rescuing transgene expression in cells and transgenic mice. Thus, inducible, reversible and site-specific epigenetic modulation is a promising strategy for reactivation of silenced transgene expression, independent of the integration site.

An autonomous CEBPA enhancer specific for myeloid-lineage priming and neutrophilic differentiation.

Neutrophilic differentiation is dependent on CCAAT enhancer-binding protein α (C/EBPα), a transcription factor expressed in multiple organs including the bone marrow. Using functional genomic technologies in combination with clustered regularly-interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 genome editing and in vivo mouse modeling, we show that CEBPA is located in a 170-kb topological-associated domain that contains 14 potential enhancers. Of these, 1 enhancer located +42 kb from CEBPA is active and engages with the CEBPA promoter in myeloid cells only. Germ line deletion of the homologous enhancer in mice in vivo reduces Cebpa levels exclusively in hematopoietic stem cells (HSCs) and myeloid-primed progenitor cells leading to severe defects in the granulocytic lineage, without affecting any other Cebpa-expressing organ studied. The enhancer-deleted progenitor cells lose their myeloid transcription program and are blocked in differentiation. Deletion of the enhancer also causes loss of HSC maintenance. We conclude that a single +42-kb enhancer is essential for CEBPA expression in myeloid cells only.

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.