Dynamic epigenetic enhancer signatures reveal key transcription factors associated with monocytic differentiation states

Down-regulation of PLCγ2-β-catenin pathway promotes activation and expansion of myeloid-derived suppressor cells in cancer

Myeloid-derived suppressor cells (MDSCs) favor tumor promotion, mainly by suppressing antitumor T cell responses in many cancers. Although the mechanism of T cell inhibition is established, the pathways leading to MDSC accumulation in bone marrow and secondary lymphoid organs of tumor-bearing hosts remain unclear. We demonstrate that down-regulation of PLCγ2 signaling in MDSCs is responsible for their aberrant expansion during tumor progression. PLCγ2(-/-) MDSCs show stronger immune-suppressive activity against CD8(+) T cells than WT MDSCs and potently promote tumor growth when adoptively transferred into WT mice. Mechanistically, PLCγ2(-/-) MDSCs display reduced β-catenin levels, and restoration of β-catenin expression decreases their expansion and tumor growth. Consistent with a negative role for β-catenin in MDSCs, its deletion in the myeloid population leads to MDSC accumulation and supports tumor progression, whereas expression of β-catenin constitutively active reduces MDSC numbers and protects from tumor growth. Further emphasizing the clinical relevance of these findings, MDSCs isolated from pancreatic cancer patients show reduced p-PLCγ2 and β-catenin levels compared with healthy controls, similar to tumor-bearing mice. Thus, for the first time, we demonstrate that down-regulation of PLCγ2-β-catenin pathway occurs in mice and humans and leads to MDSC-mediated tumor expansion, raising concerns about the efficacy of systemic β-catenin blockade as anti-cancer therapy.

Lysine-specific demethylase 1 represses THP-1 monocyte-to-macrophage differentiation

OBJECTIVE

To investigate the role of lysine-specific demethylase 1 (LSD1) in the process of THP-1 monocyte-to-macrophage differentiation.

METHODS

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and Western blotting were performed to analyze the expression of LSD1 and interleukin-6 (IL-6) in THP-1 monocytes and THP-1-derived macrophages. Chromatin immunoprecipitation (ChIP) assay was applied to detect the occupancy of LSD1 and H3K4 methylation at IL-6 promoter during THP-1 monocyte-to-macrophage differentiation. IL-6 mRNA level and H3K4 methylation at IL-6 promoter were analyzed using qRT-PCR and ChIP assay in LSD1-knockdown THP-1 cells treated with 12-O-tetradecanoylphorbol-13-acetate (TPA) for 0, 4, 8, 12, and 24 hours. Fluorescence activated flow cytometry was performed to reveal the percentage of macrophages differentiated from THP-1 monocytes.

RESULTS

The expression of LSD1 reduced during THP-1 monocyte-to-macrophage differentiation (P<0.01). LSD1 occupancy decreased and H3K4 methylation increased at IL-6 promoter during the differentiation. With knockdown of LSD1, H3K4 methylation at IL-6 promoter was found increased after TPA treatment at different times points (all P<0.05, except 24 hours). The percentage of macrophages increased significantly in the THP-1 cells with LSD1 knockdown (P<0.05).

CONCLUSIONS

LSD1 is repressed during the monocyte-to-macrophage differentiation of THP-1 cells. Suppression of LSD1-mediated H3K4 demethylation may be required for THP-1 monocyte-to-macrophage differentiation.

Synergistic activation of inflammatory cytokine genes by interferon-γ-induced chromatin remodeling and toll-like receptor signaling

Synergistic activation of inflammatory cytokine genes by interferon-γ (IFN-γ) and Toll-like receptor (TLR) signaling is important for innate immunity and inflammatory disease pathogenesis. Enhancement of TLR signaling, a previously proposed mechanism, is insufficient to explain strong synergistic activation of cytokine production in human macrophages. Rather, we found that IFN-γ induced sustained occupancy of transcription factors STAT1, IRF-1, and associated histone acetylation at promoters and enhancers at the TNF, IL6, and IL12B loci. This priming of chromatin did not activate transcription but greatly increased and prolonged recruitment of TLR4-induced transcription factors and RNA polymerase II to gene promoters and enhancers. Priming sensitized cytokine transcription to suppression by Jak inhibitors. Genome-wide analysis revealed pervasive priming of regulatory elements by IFN-γ and linked coordinate priming of promoters and enhancers with synergistic induction of transcription. Our results provide a synergy mechanism whereby IFN-γ creates a primed chromatin environment to augment TLR-induced gene transcription.

Identification and characterization of regulatory elements in the promoter of ACVR1, the gene mutated in Fibrodysplasia Ossificans Progressiva

Background

The ACVR1 gene encodes a type I receptor for bone morphogenetic proteins (BMPs). Mutations in the ACVR1 gene are associated with Fibrodysplasia Ossificans Progressiva (FOP), a rare and extremely disabling disorder characterized by congenital malformation of the great toes and progressive heterotopic endochondral ossification in muscles and other non-skeletal tissues. Several aspects of FOP pathophysiology are still poorly understood, including mechanisms regulating ACVR1 expression. This work aimed to identify regulatory elements that control ACVR1 gene transcription.

Methods and results

We first characterized the structure and composition of human ACVR1 gene transcripts by identifying the transcription start site, and then characterized a 2.9 kb upstream region. This region showed strong activating activity when tested by reporter gene assays in transfected cells. We identified specific elements within the 2.9 kb region that are important for transcription factor binding using deletion constructs, co-transfection experiments with plasmids expressing selected transcription factors, site-directed mutagenesis of consensus binding-site sequences, and by protein/DNA binding assays. We also characterized a GC-rich minimal promoter region containing binding sites for the Sp1 transcription factor.

Conclusions

Our results showed that several transcription factors such as Egr-1, Egr-2, ZBTB7A/LRF, and Hey1, regulate the ACVR1 promoter by binding to the -762/-308 region, which is essential to confer maximal transcriptional activity. The Sp1 transcription factor acts at the most proximal promoter segment upstream of the transcription start site. We observed significant differences in different cell types suggesting tissue specificity of transcriptional regulation. These findings provide novel insights into the molecular mechanisms that regulate expression of the ACVR1 gene and that could be targets of new strategies for future therapeutic treatments.

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.

Enhancer transcripts mark active estrogen receptor binding sites

We have integrated and analyzed a large number of data sets from a variety of genomic assays using a novel computational pipeline to provide a global view of estrogen receptor 1 (ESR1; a.k.a. ERα) enhancers in MCF-7 human breast cancer cells. Using this approach, we have defined a class of primary transcripts (eRNAs) that are transcribed uni- or bidirectionally from estrogen receptor binding sites (ERBSs) with an average transcription unit length of ∼3–5 kb. The majority are up-regulated by short treatments with estradiol (i.e., 10, 25, or 40 min) with kinetics that precede or match the induction of the target genes. The production of eRNAs at ERBSs is strongly correlated with the enrichment of a number of genomic features that are associated with enhancers (e.g., H3K4me1, H3K27ac, EP300/CREBBP, RNA polymerase II, open chromatin architecture), as well as enhancer looping to target gene promoters. In the absence of eRNA production, strong enrichment of these features is not observed, even though ESR1 binding is evident. We find that flavopiridol, a CDK9 inhibitor that blocks transcription elongation, inhibits eRNA production but does not affect other molecular indicators of enhancer activity, suggesting that eRNA production occurs after the assembly of active enhancers. Finally, we show that an enhancer transcription “signature” based on GRO-seq data can be used for de novo enhancer prediction across cell types. Together, our studies shed new light on the activity of ESR1 at its enhancer sites and provide new insights about enhancer function.

Predicting cell types and genetic variations contributing to disease by combining GWAS and epigenetic data

Genome-wide association studies (GWASs) identify single nucleotide polymorphisms (SNPs) that are enriched in individuals suffering from a given disease. Most disease-associated SNPs fall into non-coding regions, so that it is not straightforward to infer phenotype or function; moreover, many SNPs are in tight genetic linkage, so that a SNP identified as associated with a particular disease may not itself be causal, but rather signify the presence of a linked SNP that is functionally relevant to disease pathogenesis. Here, we present an analysis method that takes advantage of the recent rapid accumulation of epigenomics data to address these problems for some SNPs. Using asthma as a prototypic example; we show that non-coding disease-associated SNPs are enriched in genomic regions that function as regulators of transcription, such as enhancers and promoters. Identifying enhancers based on the presence of the histone modification marks such as H3K4me1 in different cell types, we show that the location of enhancers is highly cell-type specific. We use these findings to predict which SNPs are likely to be directly contributing to disease based on their presence in regulatory regions, and in which cell types their effect is expected to be detectable. Moreover, we can also predict which cell types contribute to a disease based on overlap of the disease-associated SNPs with the locations of enhancers present in a given cell type. Finally, we suggest that it will be possible to re-analyze GWAS studies with much higher power by limiting the SNPs considered to those in coding or regulatory regions of cell types relevant to a given disease.

Cell-Type Specific Determinants of NRAMP1 Expression in Professional Phagocytes

The Natural resistance-associated macrophage protein 1 (Nramp1 or Solute carrier 11 member 1, Slc11a1) transports divalent metals across the membrane of late endosomes and lysosomes in professional phagocytes. Nramp1 represents an ancient eukaryotic cell-autonomous defense whereas the gene duplication that yielded Nramp1 and Nramp2 predated the origin of Sarcopterygians (lobe-finned fishes and tetrapods). SLC11A1 genetic polymorphisms associated with human resistance to tuberculosis consist of potential regulatory variants. Herein, current knowledge of the regulation of SLC11A1 gene expression is reviewed and comprehensive analysis of ENCODE data available for hematopoietic cell-types suggests a hypothesis for the regulation of SLC11A1 expression during myeloid development and phagocyte functional polarization. SLC11A1 is part of a 34.6 kb CTCF-insulated locus scattered with predicted regulatory elements: a 3' enhancer, a large 5' enhancer domain and four elements spread around the transcription start site (TSS), including several C/EBP and PU.1 sites. SLC11A1 locus ends appear mobilized by ETS-related factors early during myelopoiesis; activation of both 5' and 3' enhancers in myelo-monocytic cells correlate with transcription factor binding at the TSS. Characterizing the corresponding cis/trans determinants functionally will establish the mechanisms involved and possibly reveal genetic variation that impacts susceptibility to infectious or immune diseases.

PU.1 target genes undergo Tet2-coupled demethylation and DNMT3b-mediated methylation in monocyte-to-osteoclast differentiation

BACKGROUND

DNA methylation is a key epigenetic mechanism for driving and stabilizing cell-fate decisions. Local deposition and removal of DNA methylation are tightly coupled with transcription factor binding, although the relationship varies with the specific differentiation process. Conversion of monocytes to osteoclasts is a unique terminal differentiation process within the hematopoietic system. This differentiation model is relevant to autoimmune disease and cancer, and there is abundant knowledge on the sets of transcription factors involved.

RESULTS

Here we focused on DNA methylation changes during osteoclastogenesis. Hypermethylation and hypomethylation changes took place in several thousand genes, including all relevant osteoclast differentiation and function categories. Hypomethylation occurred in association with changes in 5-hydroxymethylcytosine, a proposed intermediate toward demethylation. Transcription factor binding motif analysis revealed an over-representation of PU.1, NF-κB, and AP-1 (Jun/Fos) binding motifs in genes undergoing DNA methylation changes. Among these, only PU.1 motifs were significantly enriched in both hypermethylated and hypomethylated genes; ChIP-seq data analysis confirmed its association to both gene sets. Moreover, PU.1 interacts with both DNMT3b and TET2, suggesting its participation in driving hypermethylation and hydroxymethylation-mediated hypomethylation. Consistent with this, siRNA-mediated PU.1 knockdown in primary monocytes impaired the acquisition of DNA methylation and expression changes, and reduced the association of TET2 and DNMT3b at PU.1 targets during osteoclast differentiation.

CONCLUSIONS

The work described here identifies key changes in DNA methylation during monocyte-to-osteoclast differentiation and reveals novel roles for PU.1 in this process.

Epigenetic regulation of macrophage polarization and function

Macrophage polarization refers to development of a specific phenotype important for tissue homeostasis or host defense in response to environmental cues. Environmental factors that induce macrophage polarization include cytokines and microbial factors produced by pathogens or commensal microbiota. Signaling pathways utilized by these polarizing factors have been well characterized, but it is less clear how signals are converted into complex and sustained patterns of gene expression, and how macrophages are reprogrammed during polarization to alter their responses to subsequent environmental challenges. Emerging evidence, reviewed here, suggests an important role for epigenetic mechanisms in modulating and transmitting signals during macrophage polarization and reprogramming. Deeper understanding of epigenetic regulation of macrophage phenotype will enable development of gene-specific therapeutic approaches to enhance host defense while preserving tissue integrity and preventing chronic inflammatory diseases.

Transcription factor binding at enhancers: shaping a genomic regulatory landscape in flux

The mammalian genome is packed tightly in the nucleus of the cell. This packing is primarily facilitated by histone proteins and results in an ordered organization of the genome in chromosome territories that can be roughly divided in heterochromatic and euchromatic domains. On top of this organization several distinct gene regulatory elements on the same chromosome or other chromosomes are thought to dynamically communicate via chromatin looping. Advances in genome-wide technologies have revealed the existence of a plethora of these regulatory elements in various eukaryotic genomes. These regulatory elements are defined by particular in vitro assays as promoters, enhancers, insulators, and boundary elements. However, recent studies indicate that the in vivo distinction between these elements is often less strict. Regulatory elements are bound by a mixture of common and lineage-specific transcription factors which mediate the long-range interactions between these elements. Inappropriate modulation of the binding of these transcription factors can alter the interactions between regulatory elements, which in turn leads to aberrant gene expression with disease as an ultimate consequence. Here we discuss the bi-modal behavior of regulatory elements that act in cis (with a focus on enhancers), how their activity is modulated by transcription factor binding and the effect this has on gene regulation.