Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 24 March 2022
Sec. Inflammation
This article is part of the Research Topic Epigenetics of the Immune Component of Inflammation View all 43 articles

BET Protein Inhibition Regulates Macrophage Chromatin Accessibility and Microbiota-Dependent Colitis

Michelle Hoffner O&#x;Connor,Michelle Hoffner O’Connor1,2Ana Berglind,Ana Berglind1,2Meaghan M. Kennedy Ng,Meaghan M. Kennedy Ng1,3Benjamin P. Keith,Benjamin P. Keith1,3Zachary J. LynchZachary J. Lynch1Matthew R. SchanerMatthew R. Schaner1Erin C. Steinbach,Erin C. Steinbach1,4Jeremy HerzogJeremy Herzog1Omar K. TradOmar K. Trad1William R. JeckWilliam R. Jeck5Janelle C. Arthur,,Janelle C. Arthur1,6,7Jeremy M. Simon,,,Jeremy M. Simon3,7,8,9R. Balfour SartorR. Balfour Sartor1Terrence S. Furey,,*&#x;Terrence S. Furey1,3,10*†Shehzad Z. Sheikh,*&#x;Shehzad Z. Sheikh1,2*†
  • 1Center for Gastrointestinal Biology and Disease, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 2Department of Genetics, Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 3Department of Genetics, Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 4Department of Medicine, Division of Rheumatology, Allergy, and Immunology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 5Department of Pathology, Duke University, Durham, NC, United States
  • 6Department of Microbiology and Immunology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 7Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 8UNC Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 9Carolina Institute for Disabilities, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
  • 10Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States

Introduction: In colitis, macrophage functionality is altered compared to normal homeostatic conditions. Loss of IL-10 signaling results in an inappropriate chronic inflammatory response to bacterial stimulation. It remains unknown if inhibition of bromodomain and extra-terminal domain (BET) proteins alters usage of DNA regulatory elements responsible for driving inflammatory gene expression. We determined if the BET inhibitor, (+)-JQ1, could suppress inflammatory activation of macrophages in Il10-/- mice.

Methods: We performed ATAC-seq and RNA-seq on Il10-/- bone marrow-derived macrophages (BMDMs) cultured in the presence and absence of lipopolysaccharide (LPS) with and without treatment with (+)-JQ1 and evaluated changes in chromatin accessibility and gene expression. Germ-free Il10-/- mice were treated with (+)-JQ1, colonized with fecal slurries and underwent histological and molecular evaluation 14-days post colonization.

Results: Treatment with (+)-JQ1 suppressed LPS-induced changes in chromatin at distal regulatory elements associated with inflammatory genes, particularly in regions that contain motifs for AP-1 and IRF transcription factors. This resulted in attenuation of inflammatory gene expression. Treatment with (+)-JQ1 in vivo resulted in a mild reduction in colitis severity as compared with vehicle-treated mice.

Conclusion: We identified the mechanism of action associated with a new class of compounds that may mitigate aberrant macrophage responses to bacteria in colitis.

Introduction

Macrophages are key innate immune cells primarily known for their ability to eliminate bacteria, foreign antigens, and debris through phagocytosis, but they also play a central role in the activation of downstream adaptive immune responses (1, 2). Macrophages have specialized roles in maintaining homeostasis in various tissues including the brain, lungs, liver, spleen, peritoneal cavity, and intestines, but these roles can vary (3). Signals from the local microenvironment determine tissue-specific macrophage functional capabilities under homeostatic conditions (4, 5). For example, macrophages found within the intestine are tolerized to the presence of the enteric microbiota, robustly performing phagocytosis when necessary without the release of potent inflammatory cytokines or recruitment of TH1 or TH17 responses (6, 7). Changes in the local environment or presence of external stimuli are sufficient to prompt rapid activation of macrophages, in which they gain additional functions required to help restore the homeostatic state (4, 5, 8).

DNA regulatory elements (DREs), including enhancers, promoters, and insulators, significantly contribute to the gene expression program and overall function in macrophages (4, 5, 8, 9). Using Chromatin Immuno-Precipitation (ChIP)-seq and Assay for Transposase-Accessible Chromatin (ATAC)-seq to examine the distribution of H3K27ac modifications and nucleosome depleted accessible chromatin, respectively, Lavin et al. determined that tissue resident macrophages have unique enhancer landscapes, resulting in the preference for a distinct set of transcription factors in a tissue-dependent manner (4). Others have demonstrated that the presence of new stimuli resulted in rapid chromatin remodeling to promote expression of relevant response genes (8, 10, 11). In particular, stimulation with lipopolysaccharide (LPS) leads to both phosphorylation of H3S28 and increased acetylation of H3K27, promoting rapid subsequent transcription of inflammatory genes, such as Il12β, in bone marrow-derived macrophages (BMDMs) (10, 11). Additionally, Ostuni et al. determined that stimulation-dependent epigenomic memory in BMDMs is established through the formation and usage of latent enhancers, identified by de novo H3K4me1 and H3K27ac marks following stimulation (8). Establishment of latent enhancers provides for a more rapid, specific response should macrophages subsequently encounter the same stimulus (8). Together these studies demonstrated that the macrophage epigenetic state, under homeostatic conditions and in response to stimuli, plays a central role in dictating macrophage function. However, epigenetic disturbances that result in dysregulated macrophage functioning are highly understudied.

The inflammatory bowel diseases (IBDs), Crohn’s disease (CD) and ulcerative colitis, are chronic inflammatory conditions of the intestine driven by exaggerated immune responses toward enteric bacteria in a genetically susceptible host (12). Macrophage loss of tolerance toward bacteria is a key event attributed to the pathogenesis of CD (1324). This is primarily driven by increased production of Il12β/IL-12p40, a subunit in the heterodimeric proteins, IL-12 and IL-23, which are responsible for downstream activation of TH1 immune responses (13, 17, 23). The Il10-/- mouse is a genotype-driven, microbiota-dependent mouse model in which loss of production of IL-10 by T regulatory cells results in inflammatory activation of macrophages leading to a classical TH1 inflammatory immune response (2528). This dysregulated inflammatory response results in major pathological changes, including increased immune infiltrate, epithelial hypoplasia, crypt abscesses, and focal ulcers throughout the caecum, colon and rectum and leaves mice vulnerable to death resulting from bacterial-induced septic shock (25). Previously, we identified altered chromatin remodeling in macrophages from both Il10-/- mice and colonic tissue isolated from CD patients (9). Using BMDMs and lamina propria macrophages (LP MΦs) isolated from wild-type (WT) and Il10-/- mice, we determined that both macrophage subsets had genome-wide chromatin accessibility changes associated with bacterial stimulation as well as genotype-specific modifications associated with the loss of IL-10. Increased accessibility was most frequently associated with LPS stimulation in BMDMs while increased accessibility in LP MΦs was most frequently associated with the loss of IL-10. Interestingly, the accessible chromatin regions associated with these two categories corresponded to a common set of genes, particularly those involved with the inflammatory response, including Nos2, Irf1, Il7r, and Rel. These data emphasize the importance of genetic predisposition as a key driver of chromatin structure that is independent of bacteria in addition to unique chromatin responses driven by bacterial stimuli. A major challenge is how to prevent the formation of an aberrant chromatin state in macrophages that promotes chronic intestinal inflammation.

The bromodomain and extra-terminal domain (BET) protein family consists of four distinct proteins known for their ability to recognize acetylated lysine residues on histone proteins, which contributes to gene regulatory activity (29). Small molecule binding of the BET proteins’ two bromodomain active sites prevents the recognition of acetylated lysines resulting in attenuated inflammatory responses (30, 31). Nicodeme et al. originally demonstrated BET protein binding increased upon LPS stimulation of BMDMs, but the addition of the small molecule, I-BET151, was sufficient to decrease BET protein binding and decrease expression of LPS responsive genes (30). Studies using other BET inhibitors, such as (+)-JQ1, confirmed that inhibition of BET proteins attenuates LPS- and/or IFN-γ-induced secretion of these cytokines in BMDMs and in vivo murine models of endotoxic shock (3134). Although BET proteins have been implicated in direct remodeling of chromatin, previous studies have not examined the full of extent of their role shaping chromatin accessibility in the LPS response (35, 36). Additionally, findings from Cheung et al. suggest that BET inhibitors are capable of limiting chronic inflammatory conditions, such as adoptive T-cell transfer-induced colitis by limiting expansion of TH1 and TH17 cell differentiation (37). However, it is unclear if the in vitro responses observed in macrophages, activated upstream of T cells in vivo, can be recapitulated and what is the mechanism by which BET inhibitors function to limit inflammatory gene expression.

In this study, we sought to determine if inhibition of BET proteins using the inhibitor (+)-JQ1 is sufficient to mitigate chromatin remodeling associated with the inflammatory state in macrophages when challenged with LPS. Our findings indicate that treatment of macrophages with (+)-JQ1 prior to LPS stimulation limited genome-wide changes in chromatin accessibility, particularly at regions distal to transcription start sites. Concomitant changes in gene expression were also observed, where combinations of treatment with (+)-JQ1 and LPS stimulation revealed 10 distinct classes of genes based on their patterns of expression. Expression of approximately 90% of LPS-induced genes identified was attenuated with (+)-JQ1 treatment. Analysis of differentially accessible chromatin regions nearby but distal to transcription start sites (TSSs) of LPS-induced genes revealed enrichment for AP-1 and IRF transcription factor binding motifs and that (+)-JQ1 treatment prevented remodeling at approximately 1,100 AP-1 and IRF motif target sequences. Finally, to evaluate efficacy of (+)-JQ1 in vivo, we treated Il10-/- with (+)-JQ1 following colonization with fecal slurries and determined (+)-JQ1 results in a mild reduction in colitis severity. Together, these data highlight that BET protein inhibition by (+)-JQ1 is sufficient to reduce LPS-induced chromatin remodeling, giving rise to attenuated inflammatory responses in Il10-/- macrophages more consistent with the function of LP MΦs.

Materials and Methods

Mice

Wild-type (WT) and Il10-/- mice on the C57BL/6J background maintained under specific pathogen free (SPF) conditions were used to generate BMDMs. Germ free Il10-/- mice (C57BL/6J background) were bred and maintained by the National Gnotobiotic Rodent Resource Center (NGRRC) and used to evaluate the clinical efficacy of (+)-JQ1.

Colonization and Treatment of Mice

8-12-week-old Il10-/- mice (C57BL/6J background; 22 male mice and 10 female mice) were removed from germ free housing and placed in isolation in a BSL2 facility for the duration of these experiments. Treatment group assignments were balanced for gender. Mice were injected with 50 mg/kg (+)-JQ1 (MedChemExpress, #HY-13030) resuspended in 20% sulfobutylether-β-cyclodextrin (SBE-β-CD) (MedChemExpress, #HY-17031) or 10% DMSO diluted in 20% SBE-β-CD. Mice were subsequently colonized with slurries made from feces of C57BL/6J mice raised under SPF conditions resuspended in pre-reduced anaerobic 1x PBS by oral and rectal swabbing. Additional injections were given on Days 3, 6, 9, and 12. Weight was measured daily through to harvest on Day 14.

Macrophage Isolation and Stimulation

BMDMs from 12-week-old male (3 Il10-/- and 4 WT) mice were harvested and cultured as previously described in biological triplicate (9). BMDMs were counted and re-plated in duplicate for RNA- and ATAC-seq experiments. Macrophages were treated with 500nM (+)-JQ1 (MedChemExpress #HY-13030), 10 ng/mL recombinant IL-10 (PeproTech, #210-10), or vehicle control for 12 hours followed by the addition of 50 ng/mL LPS (InvivoGen, #tlrl-peklps) for 4 hours. Untreated, unstimulated WT or Il10-/- BMDMs served as an additional control and remained in culture for 16 hours. Samples were collected in TRIzol or freezing media for RNA- and ATAC-seq experiments, respectively. Cell culture supernatants were collected and stored for cytokine secretion analyses using Luminex.

RNA-Seq

RNA was isolated from murine BMDMs using the Norgen Biotek Corp. Total RNA Purification Kit (Cat. 17200) according to the manufacturer’s protocols. These kits use column-based DNase treatment to eliminate DNA contamination.

RNA-seq libraries were prepared using the Illumina KAPA Stranded RNA-seq Kit with RiboErase (HMR). Paired-end (50bp) sequencing was performed on the Illumina HiSeq 4000 platform (Gene Expression Omnibus [GEO] accession no. 184563). Reads were extracted using TagDust v1.13 and aligned to the mm9 genome assemblies using STAR with default parameters (9, 38, 39). Reads were quantified using RSEM v1.2.31 with default parameters (9, 40). PCA was performed using the prcomp function in R on DESeq2 normalized variance stabilizing transformation (VST) transformed counts for the top 1,000 most variably expressed genes (41).

Pairwise differential gene expression was determined using the DESeq2 Wald test (41). Pairwise DEGs between LPS stimulated and unstimulated controls, between LPS stimulated + (+)-JQ1 treated and unstimulated controls, and LPS stimulated and LPS stimulated + (+)-JQ1 treated were determined, the latter used to indicate whether (+)-JQ1 treatment promoted increased or decreased relative gene expression. To focus on genes significantly affected by (+)-JQ1 treatment during LPS stimulation, gene expression for individual biological replicates for LPS stimulated and LPS + (+)-JQ1 conditions were normalized to the average expression across unstimulated ll10-/- BMDMs and DEGs (base mean expression ≥ 10 and FDR < 0.05) were identified using the likelihood ratio test (LRT) in DESeq2 (41). Distinct gene expression classes (I – X) were identified based on varying combinations of levels of statistical significance and log2 fold-change of individual Wald tests (Supplementary Table 1). Z-scores were calculated based on the log2 fold-change of individual biological replicates for LPS stimulated and LPS stimulated + (+)-JQ1 treated samples normalized to the average of the unstimulated and untreated Il10-/- BMDMs. These z-scores were used to generate heatmaps of the ten classes. Enriched KEGG pathways (FDR < 0.05) were identified using Enrichr (42).

ATAC-Seq

ATAC-seq was performed as previously described (43). To prepare nuclei, cells were pelleted and washed with cold PBS followed by lysis using cold lysis buffer (10mM Tris-HCl, 10mM NaCl, 3 mM MgCl2, and 0.1% NP-40). Pelleted nuclei underwent transposition using the Nextera DNA Library Prep Kit (Illumina #FC-121-1030). Samples were resuspended in the transposition reaction (12.5µL 2x TD buffer, 2 µL transposase, and 10.5 µL nuclease-free water) and incubated at 37C for 1 hour at 1000 RPM. Transposed DNA samples were purified using the Qiagen MinElute Kit (#28204) followed by amplification using 1x PCR master mix (NEB #M0541S) and 25µM Ad1_noMX and Ad2.* indexing primer for 10-14 cycles. Libraries were purified and size selected by magnetic separation using Agencourt AMPure XP magnetic beads (Beckman Coulter #A63880). Paired-end (50bp) sequencing was performed on the Illumina HiSeq 4000 platform (GEO accession no. 183564). Data was processed using PEPATAC v0.9.0 with default parameters (44). Uniquely mapped reads were aligned to the mm9 genome using bowtie2 (45). Peaks were called on individual samples using MACS2 (46, 47). ChIPSeeker was used to classify peaks as promoter or distal, where any peak that that did not fall in the -1000bp/+100bp from the TSS of mm9 reference genes was classified as distal (48). BEDTools was used to determine distal peaks +/-25kb from the TSS (50kb window) of each gene (49).

PCA was performed using the prcomp function in R on DESeq2 normalized, variance stabilizing transformation (VST) transformed counts for the top 10,000 most variably accessible peaks (41). Promoter and distal peaks were calculated and scaled with the computeMatrix tool from DeepTools v 3.5.0 and plotted using either plotHeatmap or plotProfile (50). Differentially accessible regions (DARs) across the union set of peaks was determined using DESeq2 (41). DARs with a |log2 fold-change| > 1 were identified using VennDiagram v1.6.20 (51). VST transformed counts were used for Spearman-ranked correlation analyses.

For each gene, DARs +/-25kb of the TSS were associated with the gene. Enrichments of DARs around genes in expression-defined categories (I – X) were calculated using a hypergeometric test. Known transcription factor motif analysis was performed with findMotifsGenome from HOMER using the middle 500bp window of DARs (FDR < 0.05). The annotatePeaks command was used to identify the location of enriched transcription factor motifs (FDR < 0.10) (52, 53). Motif matrices were clustered using Regulatory Sequence Analysis Tools (RSAT) using the Ncor metric for motif-to-motif similarity matrix with the lower thresholds set to 5 for width, 0.7 for Pearson correlation, and 0.4 for relative width-normalized Pearson correlation (54).

Luminex

Cytokine levels for BAFF, CRP, IL-12p70, IL-6, IL-27, LDLR, MCP-1, MCP-2, MIP-1α, and TNF-α found in cell culture medium following WT and Il10-/- BMDM stimulation were measured by Luminex (R&D) and analyzed on a Bio-plex 200 (Bio-Rad Laboratories). Cytokine levels were calculated based on the standard curve for each analyte. Samples were plated in technical duplicate and biological triplicate. Data-points are representative of the average of technical duplicates for each sample.

Histology

Slides containing sections of proximal and distal colon were prepared for H & E staining. Two independent scorers blinded to the experimental and control groups performed histological analysis using an established scoring system for evaluating goblet cell loss and submucosal infiltrate in animal models (55). Histological sub-scores were added together to generate a composite histology score.

RNA Extraction and qPCR

Total RNA was isolated from WT and Il10-/- BMDMs for all experimental conditions and caecal tissue taken from GF (Day 0) and colonized (Day 14) Il10-/- mice using the Norgen Biotek Corp. Total RNA Purification Kit (Cat. 17200). cDNA was derived from 500ng RNA by reverse transcriptase using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, 4368814) according to the manufacturer’s specifications. Quantitative real-time PCR was performed using 10ng cDNA with the PowerUp SYBR Green Master Mix (Applied Biosystems, A25741) and 10 µM of forward and reverse primers for Il12β (F: 5’ CGCAAGAAAGAAAAGATGAAGGAG 3’ R: 5’ TTGCATTGGACTTCGGTAGATG 3’), Tnf (F: 5’ ACCCTCACACTCAGATCATCTTCTC 3’ R: 5’ TGAGATCCATGCCGTTGG 3’) and β-actin (F: 5’ AGCCATGTACGTAGCCATCCAG 3’ R: 5’ TGGCGTGAGGGAGAGCATAG 3’). Reactions were performed in triplicate for each biological sample. Fold-change was calculated by determining ΔΔCt for all reactions followed by normalization to the average ΔΔCt value for unstimulated BMDMs and GF animals.

Flow Cytometry

Il10-/- BMDMs were cultured for 16 hours with 500nM (+)-JQ1. Cells were washed with cold 1x PBS followed by staining with LIVE/DEAD Fixable Blue Dead Cell Stain Kit (1:1000; Invitrogen L23105). Cell cycle changes in (+)-JQ1 treated Il10-/- BMDMs were assessed using the BD Pharmigen BrdU Flow kit according to the manufacturer (559619). Briefly, Il10-/- BMDMs were co-cultured with BrdU and treated with (+)-JQ1 for 16 hours followed by fixation and permeabilization. Cells were treated with DNase followed by staining with FITC-conjugated anti-BrdU.

Bulk LPMCs from Il10-/- mice were washed with cold 1x PBS and stained for viability followed by cell-surface marker staining for CD45 (1:100; Clone 30-F11, BioLegend), CD3ε (1:300; Clone 145-2C11, BioLegend), CD19 (1:200; Clone 6D5, BioLegend), CD11b (1:200; M1/70, BD Biosciences), CD11c (1:200; Clone N418, BioLegend), and F4/80 (1:200; Clone BM8, BioLegend) diluted in staining buffer (5% FBS/PBS). All samples were fixed with 4% PFA. Data was acquired with the FACSDIVA software using the BD LSR II and analyzed using FlowJo version 10.7.1.

Statistics

Differential analyses of ATAC-seq and RNA-seq data were performed using DESeq2 with FDR adjusted P values used to measure statistical significance (41). Spearman rank correlations of VST transformed counts were calculated using RStudio Version 1.2.5033. Known motif enrichment was determined using HOMER, which generates P values by screening a library of reliable motifs against target and background sequences. GraphPad Prism 8 for Mac (Graphpad Software Inc.) was used for all other statistical analyses. Statistical significance of change in the number of DARs with a |log2 fold-change| > 0.5 for each condition was determined using a 2x2 chi-square test. Statistical significance comparing the absolute log2 fold-change values for LPS stimulated and LPS stimulated + (+)-JQ1 treated samples were determined by Wilcoxon match-paired signed rank tests. Statistical differences in relative gene expression and cytokine secretion for the LPS and LPS + (+)-JQ1 conditions were determined using paired, parametric t-tests. Statistical differences in relative gene expression and cytokine secretion for the LPS, LPS + (+)-JQ1 and LPS + IL-10 conditions were determined using one-way ANOVA followed by post hoc analyses using Tukey’s multiple comparisons test. Paired parametric t-tests were used to determine statistical significance for (+)-JQ1 and DMSO treated BMDMs for the evaluation of the number of cells in S-phase and percentage of viable cells. All other statistics were determined using a 2-tailed unpaired, nonparametric, student’s t test. For all tests, except for motif analysis, Padj < 0.05, empirical P < 0.05, or P < 0.05 were considered statistically significant. For motif analysis, Padj < 0.1 was considered statistically significant.

Study Approval

All animal experiments were performed in accordance with protocols approved by the Institutional Animal Care and Use Committee of the University of North Carolina at Chapel Hill (19-108.0).

Results

Inhibition of BET Proteins Limits LPS-Induced Changes in Macrophage Chromatin Accessibility Genome-Wide

Genome-wide epigenetic changes, including post-translational histone modifications and chromatin accessibility, have been heavily annotated in macrophages upon LPS stimulation (810, 56). Increased BET protein binding has been associated with expression of LPS-induced inflammatory genes (3033, 37). To comprehensively determine the role of BET proteins in LPS-induced chromatin remodeling, we treated Il10-/- BMDMs with (+)-JQ1 for 12 hours prior to challenge with LPS for 4 hours and evaluated changes in chromatin accessibility by ATAC-seq (Figure 1A). Treatment with (+)-JQ1 did not affect cell viability (Supplementary Figure 1A) or cell cycle progression (Supplementary Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1 (+)-JQ1 treatment prevents LPS-induced chromatin accessibility changes in distal regions. (A) Schematic of experimental workflow. (B) Principle component analysis (PCA) of chromatin accessibility profiles for samples that remained unstimulated (triangles) or were stimulated with LPS (circles) in the absence (purple) or presence (aqua) of (+)-JQ1. Average chromatin accessibility signal profiles for (C) promoter regions (-1000/+100 bases) and (D) distal regions for unstimulated (black), LPS stimulated (purple) and LPS + (+)-JQ1 (aqua) samples. (E) Identification of regions that differentially increase (red) or decrease (blue) in accessibility during LPS stimulation with or without (+)-JQ1 as compared with unstimulated BMDMs at promoter (left) and distal (right) regions (FDR < 0.05). Red and blue hashtags indicate difference in the number of regions with increased (red) or decreased (blue) accessibility with (+)-JQ1 is statistically significant using a chi-square test. Identification of peaks shared with |log2 fold-change >1| (left) and comparison of average log2 fold-change values using Wilcoxon sum rank tests for (F) promoters increasing in accessibility, (G) promoters decreasing in accessibility, (H) distal regions increasing in accessibility, and (I) distal regions decreasing in accessibility. FDR adjusted P values determined using DESeq2. N = 3 biological replicates. N.S. stands for Not Significant.

Principle component analysis (PCA) of chromatin accessibility profiles in Il10-/- BMDMs (Figure 1B) showed the first principle component (PC1) stratified samples based on the presence or absence of LPS stimulation, demonstrating that LPS is a primary driver of chromatin accessibility changes. Samples were further segregated along the second PC based on treatment with (+)-JQ1. We note that samples treated with (+)-JQ1 prior to LPS stimulation clustered more closely to unstimulated samples along PC1 than untreated LPS stimulated samples (Figure 1B), suggesting (+)-JQ1 mitigates some LPS-induced chromatin accessibility changes in Il10-/- BMDMs.

Both LPS stimulation and the loss of IL-10 expression were previously shown to affect genome-wide chromatin accessibility in Il10-/- BMDMs (9). To determine how (+)-JQ1 modifies chromatin accessibility during LPS stimulation at both promoter (n = 36,670) and distal (n=157,614) sites genome-wide, we plotted average chromatin accessibility (Figures 1C, D) and quantified the overall number and magnitude of differentially accessible regions (DARs; FDR < 0.05) for both LPS and LPS + (+)-JQ1 conditions as compared with unstimulated Il10-/- BMDMs (Figure 1E). Average chromatin accessibility at both promoter and distal sites decreased with (+)-JQ1 treatment as compared with LPS stimulation alone (Figures 1C, D). However, the relative decrease in accessibility was greater at distal sites (Figure 1D) as compared with promoter sites (Figure 1C). The number of promoter DARS that had increased accessibility upon LPS stimulation was roughly equivalent under each condition with no statistically significant differences (Figure 1E). In contrast, (+)-JQ1 treatment prior to LPS stimulation resulted in fewer promoter DARs with decreased accessibility as compared with LPS stimulation alone (Figure 1E; p=8.54x10-8). In contrast at distal sites, both the overall number of DARs as well as the percentage of high magnitude DARs were decreased substantially with (+)-JQ1 treatment (Figure 1E; p<2.2x10-16 for both increased and decreased). Focusing on DARs with high magnitude changes at both promoter and distal sites, we found greater concordance across the two conditions when accessibility was increased by LPS stimulation than when decreased (Figures 1F–I). Even amongst these DARs, though, (+)-JQ1 treatment significantly decreased the average magnitude (|log2FC|) for all four categories (Figures 1F–I, right; Wilcoxon sum rank test, p < 0.005). Based on these data, we conclude that (+)-JQ1 treatment overall attenuates the effect of LPS on changes to genome-wide chromatin accessibility and has a greater impact on DARs at distal sites.

Exogenous IL-10 Treatment Drives Changes in Chromatin Accessibility Distinct From BET Inhibition

We previously showed that supplementation with exogenous IL-10 had minimal effects on genome-wide chromatin remodeling in Il10-/- BMDMs (9). To determine how IL-10 treatment compared to (+)-JQ1 in limiting LPS-induced chromatin remodeling, we treated Il10-/- BMDMs with IL-10 for 12 hours followed by LPS challenge for 4 hours. PCA revealed that samples treated with IL-10 clustered separately from samples treated with (+)-JQ1 and from untreated samples, regardless of LPS stimulation (Supplementary Figure 2A). IL-10 treated samples had decreased promoter accessibility compared with (+)-JQ1 treated samples, on average (Supplementary Figure 2B). In contrast, (+)-JQ1 treated samples had decreased accessibility at distal regions compared to IL-10 samples (Supplementary Figure 2C). When comparing with unstimulated and untreated Il10-/- BMDMs, IL-10 treatment prior to LPS stimulation resulted less promoter DARs as compared with LPS stimulation alone or with (+)-JQ1 treatment (Supplementary Figure 2D; p=2.003x10-15 for increased, p<2.2x10-16 for decreased). In contrast, the overall number of DARs (increased and decreased together) detected at distal sites with IL-10 treatment prior to LPS stimulation were similar to what was observed with LPS stimulation alone. When separately testing DARs with increased or decreased accessibility, differences in each were statistically different by chi-square analysis (Supplementary Figure 2D; 1.811x10-15 for increased, 2.2x10-16 for decreased). To better quantify effects of treatments on chromatin accessibility changes with LPS stimulation, we correlated normalized signals across all accessible regions between each of LPS stimulated alone, LPS + (+)-JQ1, and LPS + IL-10 with unstimulated Il10-/- BMDMs (Supplementary Figures 2E–G). We found (+)-JQ1 treated samples were the most correlated with unstimulated samples suggesting a greater impact than IL-10 on limiting chromatin remodeling. Together, these data show (+)-JQ1 better preserves a non-inflammatory chromatin state and that IL-10 signaling affects chromatin remodeling in a distinct manner.

BET Protein Inhibition Attenuates LPS-Induced Expression of Inflammatory Genes

LPS stimulation of Il10-/- BMDMs results in robust expression of key inflammatory genes, such as Il12β, that is greater than expression observed in WT BMDMs (20, 22). To better understand how BET protein inhibition affects gene expression in response to LPS stimulation, we performed RNA-seq using matched BMDM samples (Figure 1A). Similar to our initial ATAC-seq analyses, PCA on these data revealed the primary source of variation is due to LPS stimulation (PC1; Figure 2A) while (+)-JQ1 treatment orthogonally contributes to variation (PC2). Again, samples treated with (+)-JQ1 prior to LPS stimulation clustered closer to unstimulated samples along PC1 than untreated LPS stimulated samples.

FIGURE 2
www.frontiersin.org

Figure 2 (+)-JQ1 treatment attenuates LPS-induced inflammatory gene expression in Il10-/- macrophages. (A) Principle component analysis (PCA) of gene expression profiles for samples that remained unstimulated (triangles) or were stimulated with LPS (circles) in the absence (purple) or presence (aqua) of (+)-JQ1. (B) Representative heatmap of calculated z-scores based on log2 fold-change of individual replicates for the LPS and LPS + (+)-JQ1 conditions normalized to the average of unstimulated, untreated Il10-/- BMDMs. Genes were identified and divided into 2 major clusters and sub-divided into 5 classes through differential analyses outlined in Supplementary Table 1. Bar graphs represent average change in expression for LPS and LPS (+)-JQ1 conditions normalized to baseline. (C) KEGG pathway analysis for the genes identified in the 10 classes in (B). (D) Heatmaps of representative genes from Classes VI, VII, and IX that correspond to genes induced by LPS that are attenuated with (+)-JQ1. (E) Heatmaps of representative inflammatory cytokines induced by LPS that are attenuated with (+)-JQ1. Gene expression (qRT-PCR; left panels) and cytokine secretion (Luminex; right panels) analyses for (F) Il12β and IL-12p70, and (G) Tnf and TNF-α. Significance determined using paired parametric, student’s t tests. N = 3 biological replicates.

To determine which genes were affected by LPS stimulation both with and without (+)-JQ1 treatment, we performed pairwise comparisons with unstimulated samples separately for each. We identified a total of 10,656 differentially expressed genes (DEGs) present in at least one of these comparisons. We then used a likelihood ratio test (LRT) to determine more specifically which genes were significantly affected by (+)-JQ1 during LPS-stimulation, resulting in 4,257 genes.

DEGs were clustered based on the direction of effect of (+)-JQ1 treatment compared with no treatment. Cluster 1 (n = 1,536 genes) contained genes that had higher expression after LPS stimulation when treated with (+)-JQ1 compared with LPS stimulation alone, while Cluster 2 (n = 2,721 genes) contained those with lower relative expression (Figure 2B). Next, within each cluster we identified five distinct classes based on statistical significance and log2FC magnitude associated with the pair-wise differential expression analyses with unstimulated samples (Supplementary Table 1). In Cluster 1, three small classes contained genes that were significantly up-regulated by LPS with (+)-JQ1 treatment (log2FC > 1) but were either weakly induced by (Class I; n = 203 genes), down-regulated by (Class II; n = 50 genes), or showed minimal response to (Class V; n = 108 genes) LPS stimulation alone. In the other two larger classes (Class III, n = 415 genes; Class IV, n = 760 genes), (+)-JQ1 treatment also resulted in relatively higher expression compared to no treatment, but overall expression levels of these genes regardless of treatment were below baseline expression in unstimulated, untreated macrophages (Figure 2). Pathway enrichment analysis of all genes in Cluster 1 revealed a slight enrichment for Longevity Regulating, Mismatch Repair, and Systemic Lupus Erythematosus pathways, but no additional pathways (Supplementary Figure 3A).

In Cluster 2, which consists of genes whose relative expression is lower when treated with (+)-JQ1 during LPS stimulation (Figure 2B), we find nearly all (~90%) genes significantly induced during LPS stimulation alone (log2FC > 1). The most robustly up-regulated of these genes were grouped into the largest Class VI (n = 1,422 genes), which contained genes upregulated by LPS stimulation regardless of (+)-JQ1 treatment, but where pre-treating with (+)-JQ1 resulted in a 2-fold decrease in expression compared with LPS stimulation alone on average. Pathway analysis using these genes showed enrichment for more than 20 pathways associated with inflammatory and viral response (Figure 2C) and included many pro-inflammatory genes such as Il12β, Il6, and Tnf (Figure 2D). Luminex analyses were used to determine if (+)-JQ1-mediated attenuation could also be observed at the protein level. Quantification of seven cytokines found in Class VI (IL-12p70, IL-27, IL-6, MCP-1, MCP-2, MIP-1α, and TNF- α) confirmed that (+)-JQ1 pre-treatment resulted decreased secretion relative to LPS stimulated cells alone (Figure 2E). Statistical analyses confirmed that expression of Il12β/IL-12p70 (Figure 2F) and Tnf/TNF-α (Figure 2G) at the mRNA (qRT-PCR, left) and protein levels (Luminex, right) was statistically significantly reduced with (+)-JQ1 pre-treatment as compared with LPS stimulation alone.

Class VII (n = 116 genes) and Class IX (n = 235 genes) consisted of genes induced by LPS stimulation alone, albeit less robustly than those in Class VI, but whose expression levels when treated with (+)-JQ1 treatment were lower than or similar to expression in unstimulated Il10-/- BMDMs, respectively (Figure 2B). Class VII genes were enriched for macrophage functions, such as antigen presentation and phagocytosis, while Class IX genes were involved in bacterial and viral response pathways as well as TNF signaling (Figure 2C). The remaining two classes (Class VIII, n = 385 genes; Class X, n = 563 genes) were down-regulated upon LPS stimulation regardless of (+)-JQ1 treatment, but whose expression was reduced to a greater extent when treated with (+)-JQ1 (Figure 2B). Neither class showed enrichment for immune pathways (Figure 2C).

Exogenous IL-10 Exerts Limited Anti-Inflammatory Effects on LPS-Induced Gene Expression

IL-10 negatively regulates expression of Il12β and other inflammatory genes (20, 22). We wanted to better understand how exogenous IL-10 affects the Il10-/- BMDM response to LPS and to compare this with (+)-JQ1 treatment. PCA revealed that IL-10 treatment prior to LPS stimulation resulted in transcriptional profiles markedly similar to untreated LPS stimulated Il10-/- BMDMs (Supplementary Figure 3B), including across genes in our ten previously defined LPS response classes (Supplementary Figure 3C). Expression fold-changes upon LPS stimulation of IL-10 treated samples varied across key inflammatory genes found in Class VI. Some genes, such as Il12β, did have attenuated responses comparable to (+)-JQ1 treated samples, but most other genes, such as Nod2, showed increased expression compared with (+)-JQ1 treatment (Supplementary Figure 3D). Similar changes associated with IL-10 were also observed at the protein level by Luminex (Supplementary Figure 3E). Quantitative analyses using qRT-PCR (left) and Luminex (right) for Il12β/IL-12p70 (Supplementary Figure 3F) and Tnf/TNF-α (Supplementary Figure 3G) confirmed that IL-10 pre-treatment significantly reduced LPS-induced expression, although (+)-JQ1 was slightly more effective in attenuating Il12β/IL-12p70. In contrast, IL-10 treatment had little effect on responses of Class VII and IX genes to LPS (Supplementary Figure 3D). Collectively, these results demonstrate that (+)-JQ1 treatment has a strong effect on genes induced by LPS stimulation, which were more extensive than with IL-10 treatment and resulted in down-regulation of genes implicated in the inflammatory response.

Effects of (+)-JQ1 on WT LPS-Stimulated Cells Has Less Pronounced Effects

To determine the effects of (+)-JQ1 pre-treatment on WT LPS-stimulated BMDMs, WT cells were cultured and stimulated in the exact same manner as Il10-/- BMDMs (Figure 1A) followed by qRT-PCR and luminex analyses to determine changes for select inflammatory genes. Like Il10-/- BMDMs, treatment of WT BMDMs with (+)-JQ1 did not affect cell viability (Supplementary Figure 4A) or cell cycle progression (Supplementary Figure 4B). Overall, WT BMDMs responded less robustly to LPS stimulation as demonstrated by the variation in the individual datapoints for Il12β/IL-12p70 (Supplementary Figure 4C) and Tnf/TNF-α (Supplementary Figure 4D). While (+)-JQ1 pre-treatment decreased expression at the gene expression and cytokine secretion level for both genes, these findings were not statistically significant (Supplementary Figures 4C, D). While IL-10 pre-treatment had a stronger effect at the transcriptional level as compared with cytokine secretion, none of these findings were statistically significant (Supplementary Figures 4C, D). Together, these data suggest that (+)-JQ1 and IL-10 exert anti-inflammatory effects on WT BMDMs, but the overall effects are less pronounced as compared with our findings in Il10-/- BMDMs.

(+)-JQ1 Restricts Chromatin Accessibility at Distal Regulatory Elements of Class VI Genes

Distal regulatory elements, such as enhancers, largely contribute to cell-specific gene expression programs defining cell identity (4, 5). Recent studies have suggested that BET protein inhibition alters use of cell-specific super-enhancers (5658). Therefore, we sought determine the relationship between effects of (+)-JQ1 on chromatin accessibility and altered expression at LPS-responsive genes. Across our ten gene expression classes, we found a similar number of distal accessible regions near member genes on average (+/- 25kb from the TSS; Supplementary Figures 5A, B). LPS stimulation overall increased distal chromatin accessibility in all classes, but the magnitude of this increase varied (Supplementary Figure 5C). Consistent with our genome-wide analyses, (+)-JQ1 treatment decreased the effect of LPS stimulation on chromatin accessibility changes for all 10 classes, most notably for Classes VI, VII, and IX that contain genes with attenuated expression (Figures 3A, B; Supplementary Figure 4C). In contrast, a similar analysis of IL-10 treated samples on Class VI, VII, and IX genes showed accessibility at promoter regions were even less altered on average than with (+)-JQ1 treatment (Supplementary Figure 5D). However, accessibility at distal regions associated with these classes more closely resembled those of untreated LPS stimulated samples, especially for Class VII, as compared with (+)-JQ1 treatment (Supplementary Figure 5E).

FIGURE 3
www.frontiersin.org

Figure 3 (+)-JQ1 prevents LPS-induced chromatin remodeling at putative AP-1 and IRF regulatory elements associated with class VI genes. Heatmaps of accessible chromatin during unstimulated (left), LPS stimulated (center) or LPS + (+)-JQ1 (right) conditions for Classes VI (top), VII (center), and IX (bottom) for ChIPSeeker identified (A) promoter regions and (B) distal regions. (C) Identification and quantification of sites containing known HOMER motifs for regions that increase in accessibility with LPS stimulation as compared with unstimulated, untreated BMDMs (FDR < 0.05 for differential analysis, FDR < 0.10 for motif enrichment). (D) Identification, motif clustering, and quantification of sites containing known HOMER motifs for regions that decrease in accessibility with (+)-JQ1 treatment prior to LPS stimulation as compared with LPS stimulation alone (FDR < 0.05 for differential analysis, FDR < 0.10 for motif enrichment). (E) Gene tracks of ATAC-seq peaks within 25kb of the transcription start site of Il12β for unstimulated (top), LPS + (+)-JQ1 (center) and LPS stimulated (bottom) conditions. Horizontal purple bars above the tracks are representative of peaks that are differentially increased during LPS stimulation compared with unstimulated BMDMs and are differentially decreased with (+)-JQ1 as compared with LPS alone (FDR < 0.05). Red line is representative of putative AP-1 motif identified by HOMER. Numbers to the right represent corresponding Il12β transcript levels. FDR adjusted P values for differential analysis determined using DESeq2 and for motifs using HOMER. N = 3 biological replicates.

Focusing more specifically on regions that significantly changed upon LPS stimulation, we found an enrichment for DARs near genes from Classes VI (n=504) and VII (n=47), but not Class IX (n=66) (hypergeometric test; p = 9.22x10-16, p = 7.56x10-4, p = 0.24, respectively). As Class VI was the largest and contains key inflammatory genes, we further investigated chromatin accessibility separately in promoter and distal regions. Both sets of regions showed increased chromatin accessibility upon LPS stimulation compared with unstimulated Il10-/- BMDMs (Figures 3A, B). Consistent with our genome-wide observations, (+)-JQ1 treatment showed little effect on promoter regions (Figure 3A), but distally we found an overall attenuation of accessibility changes and a reduced number of DARs in response to LPS stimulation (Figure 3B). Motif analysis of distal DARs associated with Class VI genes with increased accessibility in untreated LPS stimulated samples (n = 5,323 regions) showed enrichment for motifs of bZIP transcription factors, including many AP-1 family sub-units, ETS transcription factors, including PU.1, and IRF family members (Figure 3C). DARs with attenuated accessibility changes in (+)-JQ1 treated samples during LPS stimulation compared with LPS stimulation alone (n = 920 regions) revealed motifs for a subset of AP-1 family members as well as several IRF and RHD (NF-kB p65/p50) factors (Figure 3D). This included a putative AP-1 motif located downstream of the of the Il12β TSS (Figure 3E). DARs that decreased in accessibility during LPS stimulation alone or increased with (+)-JQ1 treatment with LPS were enriched for motifs from ETS family members (Supplementary Figures 5F, G).

Overall, these data suggest that (+)-JQ1 treatment restricts LPS-induced chromatin accessibility changes at distal regulatory elements near inflammatory genes. Interestingly, though, our data also suggests that exogenous IL-10 may also attenuate some LPS-induced chromatin changes, but primarily at promoters of LPS-inducible genes.

(+)-JQ1 Treatment Limits Onset of Experimental Colitis in Germ Free Il10-/- Mice

To determine whether (+)-JQ1 may mitigate colitis onset in vivo, we treated adult Il10-/- mice previously housed under germ free (GF) conditions with (+)-JQ1 (I.P. injections) followed by colonization with slurries made using fecal matter from WT C57BL/6J mice raised in specific pathogen free conditions (Figure 4A). Additional (+)-JQ1 treatments were administered on Days 3, 6, 9, and 12 and animals were sacrificed on Day 14 (Figure 4A). Treatment with (+)-JQ1 did not result in any significant changes in weight as compared with vehicle treated mice (Supplementary Figure 6A).

FIGURE 4
www.frontiersin.org

Figure 4 (+)-JQ1 treatment attenuates onset of microbiota-induced colitis in germ free Il10-/- mice. (A) Schematic of experimental workflow (n = 10+ mice/condition and represent at least 3 independent experiments). (B) Representative images of colon length and thickness in (+)-JQ1- and vehicle-treated mice. (C) Quantification of average colon length (n = 16 (+)-JQ1; n = 15 vehicle). (D) Representative images of H&E stained colons for proximal (top) and distal (bottom) portions of the colon in germ free (GF) mice (left), colonized mice treated with (+)-JQ1 harvested on Day 14 (center) and colonized mice treated with vehicle harvested on Day 14 (right). Scale represents 100 µM. (E) Quantification of histology scores based on H&E staining (n = 10 (+)-JQ1; n=8 vehicle). (F) Quantification of caecal Il12β expression by qPCR (n = 15/treatment). Significance values determined using a 2-tailed unpaired, non-parametric student’s t test. N.S. stands for Not Significant.

We isolated lamina propria mononuclear cells (LPMCs) and splenocytes from (+)-JQ1- and vehicle treated animals to identify changes in immune cell population sizes associated with treatment (Supplementary Figure 6B). Total CD45+ leukocytes constituted ~60% of singlets in LPMC populations and 75% of singlets in splenocyte populations in vehicle treated mice and did significantly change with (+)-JQ1 treatment (Supplementary Figure 6C). No significant differences in the population size of CD11b+CD11c-F4/80+ macrophages or total CD3+ T cells and CD19+ B cells isolated from the intestine or spleen were observed with (+)-JQ1 treatment (Supplementary Figures 6D, E). However, a reduction in the CD11bintCD11c- population was observed with (+)-JQ1 treatment (Supplementary Figure 6B). Based on these observations, we conclude that the dose and frequency of (+)-JQ1 given to these mice does not have any significant effect on bulk immune cell population size within the intestine or other peripheral organs.

Finally, we investigated the effect of (+)-JQ1 treatment on intestinal inflammation in these mice. Overall colon length and morphology remained unchanged with (+)-JQ1 treatment when compared with vehicle treated mice (Figures 4B, C). However, histologically, scored using a modified method from Berg et al. that compared colons with uninflamed uncolonized germ-free Il10-/- mice harvested on Day 0 (Figure 4D), (+)-JQ1 treated mice had statistically significantly lower scores compared with vehicle treated mice (Figure 4E). Additionally, we quantified caecal Il12β mRNA expression levels and found that (+)-JQ1 treatment also significantly decreased relative expression as compared with vehicle treated mice (Figure 4F), which correlates with our histological observations. Taken together, these data demonstrate that treatment with (+)-JQ1 attenuates severity of colitis onset in Il10-/- mice two weeks post colonization.

Discussion

The enhancer landscape is central to determining the gene expression program and overall identity of a cell (4, 5). Both the microenvironment and external stimuli influence macrophage enhancer utilization (4, 5, 811). Previously, we showed that putative enhancers regions are differentially accessible in a genotype- and stimulation-dependent manner in Il10-/- macrophages (9). In this study, we determined that BET protein inhibition using (+)-JQ1 is sufficient to attenuate LPS-induced chromatin remodeling at putative enhancers containing AP-1 and IRF motifs distal to the TSS of key inflammatory genes whose expression are induced by LPS and show concomitant attenuation with (+)-JQ1. BET proteins, particularly BRD4, increase in binding to chromatin in macrophages following LPS stimulation (30, 56). Furthermore, BRD4 is important for the formation of LPS-induced super-enhancers in close proximity to inflammatory genes such as Tnf and Nfkbia (56). Additional evidence has also suggested a potential role for BRD2 in controlling transcriptional activation by binding the chromatin insulator CTCF and the cohesion to promote cis-regulatory enhancer assembly (59). We propose that treatment with (+)-JQ1 inhibits BET protein-mediated nucleosome eviction thus preventing LPS-induced increases in chromatin accessibility at putative enhancers driving inflammatory gene expression (36). Our studies are limited by the absence of ChIP-seq datasets and functional assays that confirm the regions identified in this study are enhancers. Future studies should annotate these regions for H3K4me1 and H3K27ac marks as well as utilize luciferase-reporter assays and electrophoretic mobility shift assays to confirm enhancer activity and identify the exact transcription factors associated with promoting regulatory activity.

BET protein inhibition results in decreased posttranslational modification of transcription factors, including NF-kB, AP-1 and STAT5, thus limiting expression of genes that rely on these proteins (6062). Furthermore, BET protein recognition of K310ac on RelA is also central to NF-kB transcriptional activation (63). While we show (+)-JQ1 having a limited effect on attenuating LPS induced accessibility at putative p65 and p50 motifs, its impact on post translational modification of these subunits suggests different mechanisms through which BET protein inhibition attenuates inflammation. We also show loss of BET protein binding following LPS stimulation limits production of other transcription factors, such as IRF family members, which are responsible for additional chromatin modification to activate other inflammatory genes (64, 65).

To our knowledge, there are no studies that have determined if any of the chromatin-based effects of (+)-JQ1 are permanent following cessation in culture. Gibbons et al. recently demonstrated that (+)-JQ1 inhibits expression of IFN-γ in TH1 cells, but this blockade only resulted from the absence of BET protein binding at the promoter, not remodeling of chromatin associated with known enhancers for IFN-γ (34). Consequently, removal of (+)-JQ1 from culture resulted in the restoration of IFN-γ expression. In contrast, other studies have demonstrated blockade of chromatin modifiers is sufficient to have long-term consequences (66). Future studies are necessary to determine if (+)-JQ1-mediated inhibition of LPS-induced increased at putative AP-1 and IRF motifs are permanent following removal of (+)-JQ1 from culture and should include re-challenge with LPS to assess macrophage tolerization (67).

Previous publications have focused on pairwise differential gene expression analyses in macrophages between LPS and LPS + BET inhibitor or BET knock-out conditions, and consistently highlight the relative down regulation of inflammatory genes (30, 32, 37, 56). However, gene expression in unstimulated macrophages are seldom included in these studies despite the fact many inflammation-related surface receptors are expressed and central to macrophage functioning at baseline. Our study utilized an objective approach that presents a genome-wide annotation of LPS-stimulated gene expression changes and how (+)-JQ1 modulates this response relative to the expression observed in an unstimulated macrophage. While our findings regarding the attenuation of inflammatory genes with BET inhibition is consistent with prior studies, we demonstrate that the majority of genes induced by LPS stimulation still show increased expression with (+)-JQ1 treatment compared to unstimulated macrophages (Class VI). These data indicate that BET inhibition by (+)-JQ1 does not completely ablate the ability of macrophages to respond to a bacterial stimulus. However, our analytical approach identified additional gene expression dynamics associated with (+)-JQ1 treatment where a subset of LPS-induced genes had expression levels below that of unstimulated macrophages (Class VII) and another subset had expression levels that were comparable to unstimulated macrophages (Class IX). These include Fcgr1, Fcgr2b, Fcgr4 (Class VII), Tlr1, Tlr6, and Cd80 (Class IX), which are involved in antibody-mediated phagocytosis, bacterial signaling, and co-stimulation of T cells. Genomic data alone is not sufficient to determine the functions retained by macrophages in the presence of (+)-JQ1 and little data regarding BET inhibitor-mediated modulation of macrophage functioning exists (68). Therefore, it is necessary for future studies to determine if macrophages retain functional their phagocytic and functional capabilities upon and following BET inhibition.

We previously proposed that IL-10 signaling may promote the early establishment of accessible regulatory elements in the developing macrophage (9). Our findings here illustrate a novel role for IL-10 in facilitating chromatin remodeling, but it remains unclear if IL-10 interacts with chromatin directly or if these changes are prompted by downstream effectors in the IL-10 pathway. Co-culture of Il10-/- BMDMs with exogenous IL-10 is sufficient to alter the LPS-induced chromatin landscape but is distinct from the resulting chromatin landscape associated with BET inhibition. Chromatin profiles of samples treated with IL-10 or (+)-JQ1 prior to LPS stimulation more closely resemble unstimulated samples than untreated LPS stimulated samples, but these treatments generate distinct chromatin profiles. Interestingly, IL-10 was more potent in attenuating promoter accessibility changes, whereas (+)-JQ1 had a larger effect on distal regions. This observation, coupled with similar upregulation of LPS-induced genes as untreated LPS stimulated samples, suggests that IL-10 treatment is not sufficient to prevent chromatin changes that drive the upregulation of inflammatory genes. In contrast, (+)-JQ1 treatment predominantly resulted in the attenuation of LPS-induced gene expression, including genes involved in anti-inflammatory pathways promoted by IL-10, including Stat3 and Jak2 (69). Previously, it was shown the IL-10 serum levels were significantly lower in an LPS-induced endotoxic shock model when mice were treated with (+)-JQ1 (31). More recent findings have demonstrated that IL-10 producing B regulatory cells rely on BRD4 to promote expression of IL-10 (70). Together, these data, along with our own results, suggest that BET inhibition results in the repression of key anti-inflammatory pathways, in addition to inflammatory pathways. Future experiments should further interrogate the effects of (+)-JQ1 on macrophage LPS-induced gene expression with intact IL-10 signaling to determine the potential short and long-term consequences that may result due to the unintended inhibition of a key anti-inflammatory pathway.

The ability of BET protein inhibitors to attenuate inflammation by modulating the function of macrophages in colitis is largely unexplored. Although Cheung et al. demonstrated that treatment with BET inhibitors limited TH1 and TH17 differentiation and prevented adoptive T-cell transfer induced colitis, Weinerroither et al. revealed that the use of BET inhibitors at the onset of DSS-induced colitis exacerbated the phenotype as compared with DSS administration alone (37, 68). It is unclear whether changes in macrophage phenotype impacts T cell expansion and differentiation. Similarly, the DSS-model of colitis is an acute injury model that damages the intestinal epithelium, suggesting that (+)-JQ1 may prevent the expansion and differentiation of intestinal stem cells, resulting in worse tissue injury (71, 72). Our study presents the first evaluation of BET inhibitors as a method of mitigating onset of colitis severity in a genotype-driven, microbiota-dependent mouse model. There are limitations to our study. While our overall phenotype shows that (+)-JQ1 treatment results in modest reduction in colitis severity as compared with vehicle-treated mice, further studies are needed to identify the optimal dosing strategy that maximizes attenuation of inflammation without causing harm to the intestinal epithelium (71). The lack of differences in bulk T and B cells suggests (+)-JQ1 treatment may alter cell phenotype, as was previously observed in the adoptive T-cell transfer-induced colitis model when treated with a BET inhibitor (37). Understanding the effects of (+)-JQ1 on the genome wide regulatory landscape macrophages within the intestine is largely hindered by low yield of live macrophages an essential requirement for high quality/quantity RNA for sequencing studies. While our study focused on genomic changes associated with BMDMs, the decreased caecal expression of Il12β, a gene confirmed to be attenuated with (+)-JQ1 in BMDMs, suggests a similar attenuation of innate immune responses within the intestine.

In summary, our study provides evidence that BET inhibitors have the ability to overcome dysregulated inflammatory signaling in Il10-/- macrophages and accomplishes this through prevention of bacterial stimulation dependent chromatin remodeling. Due to the heterogeneous nature of both CD pathogenesis and presentation, it will be necessary to isolate LP MΦs from CD patients to identify underlying chromatin aberrancies driving inflammation and determine if our results can be recapitulated in ex vivo co-culture systems prior to the pursuit of clinical trials.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below:

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE183563

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE183564

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE183565.

Ethics Statement

The animal study was reviewed and approved by University of North Carolina at Chapel Hill IACUC (19-108.0).

Author Contributions

MHO, TSF, and SZS conceptualized the idea for the study and wrote the manuscript. MHO, RBS, TSF, and SZS designed the experiments. MHO, AB, ZJL, MRS, ECS, JH, and OKT performed the experiments. MHO, TF, MMKN, BPK, and JMS analyzed the sequencing data. WJ and JCA scored the histology. SZS and TSF supervised the project. SZS funded the project. All authors contributed to the article and approved the submitted version.

Funding

This work was funded in part through Helmsley Charitable Trust (SHARE Project 2), NIDDK P01DK094779, NIDDK 1R01DK104828, NIDDK P30-DK034987, NIDDK 1R01DK124617, NIH Ruth L. Kirschstein National Research Service Award Individual Predoctoral Fellowship (1F31DK122704), NIH T32 Genetics NIGMS Training Grant (T32-GM007092-43), NIH T32 Translational Medicine Training Grant (T32-GM122741), Crohn’s and Colitis Foundation Student Research Fellowship Award (528437), and Gnotobiotic Animal Facility.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We gratefully acknowledge the National Gnotobiotic Rodent Resource Center for providing all of the germ free mice used for this study (NIDDK P40OD010995, NIDDK P30DK034987, and Crohn’s and Colitis Foundation), the technical support from the UNC High Throughput Sequencing Facility (NCI P30-CA016086 & P30-ES010126), the Center for Gastrointestinal Biology and Disease Histology Core for processing and staining all histological specimens and Advanced Analytics Core for consulting and performance of the Luminex assay (NIDDK P30DK034987), and use of the LSRII by the UNC Flow Cytometry Core Facility (NCI P30-CA016086). We also thank members of the Sartor Lab, including Muyiwa Awoniyi, Bo Liu, and Akihiko Oka for generously providing SPF Il10-/- mice used for this study and for technical support for the in vivo experiments performed for this study, the Gulati lab for providing material for GF mouse colonization, Nina Nishiyama for assistance with the chi-square analysis, and Tamara Vital for providing technical support associated with the preliminary ATAC-seq analyses and supportive conversations throughout the progression of this project.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.856966/full#supplementary-material

References

1. van Furth R, Cohn ZA, Hirsch JG, Humphrey JH, Spector WG, Langevoort HL. The Mononuclear Phagocyte System: A New Classification of Macrophages, Monocytes, and Their Precursor Cells. Bull World Health Organ (1972) 46(6):845–52.

PubMed Abstract | Google Scholar

2. Lavin Y, Merad M. Macrophages: Gatekeepers of Tissue Integrity. Cancer Immunol Res (2013) 1(4):201–9. doi: 10.1158/2326-6066.CIR-13-0117

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Epelman S, Lavine KJ, Randolph GJ. Origin and Functions of Tissue Macrophages. Immunity (2014) 41(1):21–35. doi: 10.1016/j.immuni.2014.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lavin Y, Winter D, Blecher-Gonen R, David E, Keren-Shaul H, Merad M, et al. Tissue-Resident Macrophage Enhancer Landscapes Are Shaped by the Local Microenvironment. Cell (2014) 159(6):1312–26. doi: 10.1016/j.cell.2014.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Gosselin D, Link VM, Romanoski CE, Fonseca GJ, Eichenfield DZ, Spann NJ, et al. Environment Drives Selection and Function of Enhancers Controlling Tissue-Specific Macrophage Identities. Cell (2014) 159(6):1327–40. doi: 10.1016/j.cell.2014.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Smythies LE, Sellers M, Clements RH, Mosteller-Barnum M, Meng G, Benjamin WH, et al. Human Intestinal Macrophages Display Profound Inflammatory Anergy Despite Avid Phagocytic and Bacteriocidal Activity. J Clin Invest (2005) 115(1):66–75. doi: 10.1172/jci19229

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Smith PD, Smythies LE, Mosteller-Barnum M, Sibley DA, Russell MW, Merger M, et al. Intestinal Macrophages Lack CD14 and CD89 and Consequently Are Down-Regulated for LPS- and IgA-Mediated Activities. J Immunol (2001) 167(5):2651–6. doi: 10.4049/jimmunol.167.5.2651

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ostuni R, Piccolo V, Barozzi I, Polletti S, Termanini A, Bonifacio S, et al. Latent Enhancers Activated by Stimulation in Differentiated Cells. Cell (2013) 152157–71. doi: 10.1016/j.cell.2012.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Simon JM, Davis JP, Lee SE, Schaner MR, Gipson GR, Weiser M, et al. Alterations to Chromatin in Intestinal Macrophages Link IL-10 Deficiency to Inappropriate Inflammatory Responses. Eur J Immunol (2016) 46(8):1912–25. doi: 10.1002/eji.201546237

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Conaway EA, de Oliveira DC, McInnis CM, Snapper SB, Horwitz BH. Inhibition of Inflammatory Gene Transcription by IL-10 Is Associated With Rapid Suppression of Lipopolysaccharide-Induced Enhancer Activation. J Immunol (2017) 198(7):2906–15. doi: 10.4049/jimmunol.1601781

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Josefowicz SZ, Shimada M, Armache A, Li CH, Miller RM, Lin S, et al. Chromatin Kinases Act on Transcription Factors and Histone Tails in Regulation of Inducible Transcription. Mol Cell (2016) 64(2):347–61. doi: 10.1016/j.molcel.2016.09.026

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sartor RB. Mechanisms of Disease: Pathogenesis of Crohn’s Disease and Ulcerative Colitis. Nat Clin Pract Gastroenterol Hepatol (2006) 3(7):390–407. doi: 10.1038/ncpgasthep0528

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kamada N, Hisamatsu T, Okamoto S, Sato T, Matsuoka K, Arai K, et al. Abnormally Differentiated Subsets of Intestinal Macrophage Play a Key Role in Th1-Dominant Chronic Colitis Through Excess Production of IL-12 and IL-23 in Response to Bacteria. J Immunol (2005) 175(10):6900–8. doi: 10.4049/jimmunol.175.10.6900

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kobayashi T, Okamoto S, Hisamatsu T, Kamada N, Chinen H, Saito R, et al. IL23 Differentially Regulates the Th1/Th17 Balance in Ulcerative Colitis and Crohn’s Disease. Gut (2008) 57(12):1682–9. doi: 10.1136/gut.2007.135053

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kamada N, Hisamatsu T, Honda H, Kobayashi T, Chinen H, Kitazume MT, et al. Human CD14 + Macrophages in Intestinal Lamina Propria Exhibit Potent Antigen-Presenting Ability. J Immunol (2009) 183(3):1724–31. doi: 10.4049/jimmunol.0804369

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kamada N, Hisamatsu T, Okamoto S, Chinen H, Kobayashi T, Sato T, et al. Unique CD14+ Intestinal Macrophages Contribute to the Pathogenesis of Crohn Disease via IL-23/IFN-γ Axis. J Clin Invest (2008) 118(6):2269–80. doi: 10.1172/JCI34610

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Kobayashi T, Matsuoka K, Sheikh SZ, Elloumi HZ, Kamada N, Hisamatsu T, et al. NFIL3 Is a Regulator of IL-12 P40 in Macrophages and Mucosal Immunity. J Immunol (2011) 186(8):4649–55. doi: 10.4049/jimmunol.1003888

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Sheikh SZ, Kobayashi T, Matsuoka K, Onyiah JC, Plevy SE. Characterization of an Interferon-Stimulated Response Element (ISRE) in the Il23a Promoter. J Biol Chem (2011) 286(2):1174–80. doi: 10.1074/jbc.M110.147884

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Onyiah JC, Sheikh SZ, Maharshak N, Otterbein LE, Plevy SE. Heme Oxygenase-1 and Carbon Monoxide Regulate Intestinal Homeostasis and Mucosal Immune Responses to the Enteric Microbiota. Gut Microbes (2014) 5(2):220–4. doi: 10.4161/gmic.27290

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Sheikh SZ, Matsuoka K, Kobayashi T, Li F, Rubinas T, Plevy SE. Cutting Edge: IFN-γ Is a Negative Regulator of IL-23 in Murine Macrophages and Experimental Colitis. J Immunol (2010) 184(8):4069–73. doi: 10.4049/jimmunol.0903600

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Uno JK, Rao KN, Matsuoka K, Sheikh SZ, Kobayashi T, Li F, et al. Altered Macrophage Function Contributes to Colitis in Mice Defective in the Phosphoinositide-3 Kinase Subunit P110δ. Gastroenterology (2010) 139(5):1642–53. doi: 10.1053/j.gastro.2010.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kobayashi T, Matsuoka K, Sheikh SZ, Russo SM, Mishima Y, Collins C, et al. IL-10 Regulates Il12b Expression via Histone Deacetylation: Implications for Intestinal Macrophage Homeostasis. J Immunol (2012) 189(4):1792–9. doi: 10.4049/jimmunol.1200042

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Kobayashi T, Steinbach EC, Russo SM, Matsuoka K, Nochi T, Maharshak N, et al. NFIL3-Deficient Mice Develop Microbiota-Dependent, IL-12/23–Driven Spontaneous Colitis. J Immunol (2014) 192(4):1918–27. doi: 10.4049/jimmunol.1301819

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Steinbach EC, Kobayashi T, Russo SM, Sheikh SZ, Gipson GR, Kennedy ST, et al. Innate PI3K P110δ Regulates Th1/Th17 Development and Microbiota-Dependent Colitis. J Immunol (2014) 192(8):3958–68. doi: 10.4049/jimmunol.1301533

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Rennick DM, Fort MM, Davidson NJ. Studies With IL-10(-/-) Mice: An Overview. J Leukocyte Biol (1997) 61:389–96. doi: 10.1002/jlb.61.4.389

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Kühn R, Löhler J, Rennick D, Rajewsky K, Müller W. Interleukin-10-Deficient Mice Develop Chronic Enterocolitis. Cell (1993) 75:263–74. doi: 10.1016/0092-8674(93)80068-P

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sellon RK, Tonkonogy S, Schultz M, Dieleman LA, Grenther W, Balish E, et al. Resident Enteric Bacteria are Necessary for Development of Spontaneous Colitis and Immune System Activation in Interleukin-10-Deficient Mice. Infect Immun (1998) 66:5224–31. doi: 10.1128/iai.66.11.5224-5231.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Berg DJ, Davidson N, Kühn R, Müller W, Menon S, Holland G, et al. Enterocolitis and Colon Cancer in Interleukin-10-Deficient Mice are Associated With Aberrant Cytokine Production and CD4+ Th1-Like Responses. J Clin Invest (1996) 98:1010–20. doi: 10.1172/JCI118861

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Morgado-Pascual JL, Rayego-Mateos S, Tejedor L, Suarez-Alvarez B, Ruiz-Ortega M. Bromodomain and Extraterminal Proteins as Novel Epigenetic Targets for Renal Diseases. Front Pharmacol (2019) 10:1315. doi: 10.3389/fphar.2019.01315

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Nicodeme E, Jeffrey KL, Schaefer U, Beinke S, Dewell S, Chung CW, et al. Suppression of Inflammation by a Synthetic Histone Mimic. Nature (2010) 468(7327):1119–23. doi: 10.1038/nature09589

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Belkina AC, Nikolajczyk BS, Denis GV. BET Protein Function Is Required for Inflammation: Brd2 Genetic Disruption and BET Inhibitor JQ1 Impair Mouse Macrophage Inflammatory Responses. J Immunol (2013) 190(7):3670–8. doi: 10.4049/jimmunol.1202838

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Das A, Chai JC, Yang CS, Lee YS, Das ND, Jung KH, et al. Dual Transcriptome Sequencing Reveals Resistance of TLR4 Ligand-Activated Bone Marrow-Derived Macrophages to Inflammation Mediated by the BET Inhibitor JQ1. Sci Rep (2015) 5:16932. doi: 10.1038/srep16932

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Khan YM, Kirkham P, Barnes PJ, Adcock IM. Brd4 Is Essential for IL-1β-Induced Inflammation in Human Airway Epithelial Cells. PloS One (2014) 9(4):e95051. doi: 10.1371/journal.pone.0095051

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Gibbons HR, Mi DJ, Farley VM, Esmond T, Kaood MB, Aune TM. Bromodomain Inhibitor JQ1 Reversibly Blocks IFN-γ Production. Sci Rep (2019) 9(1):10280. doi: 10.1038/s41598-019-46516-x

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Conrad RJ, Fozouni P, Thomas S, Sy H, Zhang Q, Zhou MM, et al. The Short Isoform of BRD4 Promotes HIV-1 Latency by Engaging Repressive SWI/SNF Chromatin-Remodeling Complexes. Mol Cell (2017) 67(6):1001–12. doi: 10.1016/j.molcel.2017.07.025

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Devaiah BN, Case-Borden C, Gegonne A, Hsu CH, Chen Q, Meerzaman D, et al. BRD4 Is a Histone Acetyltransferase That Evicts Nucleosomes From Chromatin. Nat Struct Mol Biol (2016) 23(6):540–8. doi: 10.1038/nsmb.3228

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Cheung K, Lu G, Sharma R, Vincek A, Zhang R, Plotnikov AN, et al. BET N-Terminal Bromodomain Inhibition Selectively Blocks Th17 Cell Differentiation and Ameliorates Colitis in Mice. Proc Natl Acad Sci USA (2017) 114(11):2952–7. doi: 10.1073/pnas.1615601114

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Lassmann T, Hayashizaki Y, Daub CO. TagDust - A Program to Eliminate Artifacts From Next Generation Sequencing Data. Bioinformatics (2009) 25(21):2839–40. doi: 10.1093/bioinformatics/btp527

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics (2013) 29(1):15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Li B, Dewey CN. RSEM: Accurate Transcript Quantification From RNA-Seq Data With or Without a Reference Genome. BMC Bioinf (2011) 12:323. doi: 10.1186/1471-2105-12-323

CrossRef Full Text | Google Scholar

41. Love MI, Huber W, Anders S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol (2014) 15(12):550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: Interactive and Collaborative HTML5 Gene List Enrichment Analysis Tool. BMC Bioinf (2013) 14:128. doi: 10.1186/1471-2105-14-128

CrossRef Full Text | Google Scholar

43. Buenrostro JD, Wu B, Chang HY, Greenleaf WJ. ATAC-Seq: A Method for Assaying Chromatin Accessibility Genome-Wide. Curr Protoc Mol Biol (2015) 2015:21–9. doi: 10.1002/0471142727.mb2129s109

CrossRef Full Text | Google Scholar

44. Smith JP, Corces MR, Xu J, Reuter VP, Chang HY, Sheffield NC. PEPATAC: An Optimized ATAC-Seq Pipeline With Serial Alignments. bioRxiv (2021) 3(4):lqab101. doi: 10.1093/nargab/lqab101

CrossRef Full Text | Google Scholar

45. Langmead B, Salzberg SL. Fast Gapped-Read Alignment With Bowtie 2. Nat Methods (2012) 9(4):357–9. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Gaspar J. Improved Peak-Calling With MACS2. bioRxiv (2018) 496521. doi: 10.1101/496521

CrossRef Full Text | Google Scholar

47. Thomas R, Thomas S, Holloway AK, Pollard KS. Features That Define the Best ChIP-Seq Peak Calling Algorithms. Briefings Bioinf (2017) 18.(3):441–50 doi: 10.1093/bib/bbw035

CrossRef Full Text | Google Scholar

48. Yu G, Wang LG, He QY. ChIP Seeker: An R/Bioconductor Package for ChIP Peak Annotation, Comparison and Visualization. Bioinformatics (2015) 31(14):2382–3. doi: 10.1093/bioinformatics/btv145

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Quinlan AR, Hall IM. BEDTools: A Flexible Suite of Utilities for Comparing Genomic Features. Bioinformatics (2010) 26(6):841–2. doi: 10.1093/bioinformatics/btq033

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. DeepTools: A Flexible Platform for Exploring Deep-Sequencing Data. Nucleic Acids Res (2014) 42. doi: 10.1093/nar/gku365

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Chen H, Boutros PC. VennDiagram: A Package for the Generation of Highly-Customizable Venn and Euler Diagrams in R. BMC Bioinf (2011) 12:35. doi: 10.1186/1471-2105-12-35

CrossRef Full Text | Google Scholar

52. Benner C, Heinz S, Glass CK. HOMER - Software for Motif Discovery and Next Generation Sequencing Analysis. (2017).

Google Scholar

53. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime Cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell (2010) 38(4):576–89. doi: 10.1016/j.molcel.2010.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Castro-Mondragon JA, Jaeger S, Thieffry D, Thomas-Chollier M, van Helden J. RSAT Matrix-Clustering: Dynamic Exploration and Redundancy Reduction of Transcription Factor Binding Motif Collections. Nucleic Acids Res (2017) 45(13):119. doi: 10.1093/nar/gkx314

CrossRef Full Text | Google Scholar

55. Koelink PJ, Wildenberg ME, Stitt LW, Feagan BG, Koldijk M, van ‘t Wout AB, et al. Development of Reliable, Valid and Responsive Scoring Systems for Endoscopy and Histology in Animal Models for Inflammatory Bowel Disease. J Crohn’s Colitis (2018) 12(7):794–803. doi: 10.1093/ecco-jcc/jjy035

CrossRef Full Text | Google Scholar

56. Dey A, Yang W, Gegonne A, Nishiyama A, Pan R, Yagi R, et al. BRD4 Directs Hematopoietic Stem Cell Development and Modulates Macrophage Inflammatory Responses. EMBO J (2019) 38(7):e100293. doi: 10.15252/embj.2018100293

PubMed Abstract | CrossRef Full Text | Google Scholar

57. di Micco R, Fontanals-Cirera B, Low V, Ntziachristos P, Yuen SK, Lovell CD, et al. Control of Embryonic Stem Cell Identity by Brd4-Dependent Transcriptional Elongation of Super-Enhancer-Associated Pluripotency Genes. Cell Rep (2014) 9(1):234–47. doi: 10.1016/j.celrep.2014.08.055

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Brown JD, Feldman ZB, Doherty SP, Reyes JM, Rahl PB, Lin CY, et al. BET Bromodomain Proteins Regulate Enhancer Function During Adipogenesis. Proc Natl Acad Sci USA (2018) 115:2144–9. doi: 10.1073/pnas.1711155115

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Cheung KL, Zhang F, Jaganathan A, Sharma R, Zhang Q, Konuma T, et al. Distinct Roles of Brd2 and Brd4 in Potentiating the Transcriptional Program for Th17 Cell Differentiation. Mol Cell (2017) 65. doi: 10.1016/j.molcel.2016.12.022

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Toniolo PA, Liu S, Yeh JE, Moraes-Vieira PM, Walker SR, Vafaizadeh V, et al. Inhibiting STAT5 by the BET Bromodomain Inhibitor JQ1 Disrupts Human Dendritic Cell Maturation. J Immunol (2015) 194:3180–90. doi: 10.4049/jimmunol.1401635

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Tan YF, Wang M, Chen ZY, Wang L, Liu XH. Inhibition of BRD4 Prevents Proliferation and Epithelial–Mesenchymal Transition in Renal Cell Carcinoma via NLRP3 Inflammasome-Induced Pyroptosis. Cell Death Dis (2020) 11:239. doi: 10.1038/s41419-020-2431-2

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Zhou Z, Li X, Liu Z, Huang L, Yao Y, Li L, et al. A Bromodomain-Containing Protein 4 (BRD4) Inhibitor Suppresses Angiogenesis by Regulating AP-1 Expression. Front Pharmacol (2020) 11:1043. doi: 10.3389/fphar.2020.01043

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Huang B, Yang X-D, Zhou M-M, Ozato K, Chen L-F. Brd4 Coactivates Transcriptional Activation of NF-κb via Specific Binding to Acetylated RelA. Mol Cell Biol (2009) 29:1375–87. doi: 10.1128/mcb.01365-08

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Li P, Leonard WJ. Chromatin Accessibility and Interactions in the Transcriptional Regulation of T Cells. Front Immunol (2018) 9:2738. doi: 10.3389/fimmu.2018.02738

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Song R, Gao Y, Dozmorov I, Malladi V, Saha I, McDaniel MM, et al. IRF1 Governs the Differential Interferon-Stimulated Gene Responses in Human Monocytes and Macrophages by Regulating Chromatin Accessibility. Cell Rep (2021) 34(12):108891. doi: 10.1016/j.celrep.2021.108891

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Bintu L, Yong J, Antebi YE, McCue K, Kazuki Y, Uno N, et al. Dynamics of Epigenetic Regulation at the Single-Cell Level. Science (2016) 351(6274). doi: 10.1126/science.aab2956

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Foster SL, Hargreaves DC, Medzhitov R. Gene-Specific Control of Inflammation by TLR-Induced Chromatin Modifications. Nature (2007) 4477147):720–4. doi: 10.1038/nature05836

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Wienerroither S, Rauch I, Rosebrock F, Jamieson AM, Bradner J, Muhar M, et al. Regulation of NO Synthesis, Local Inflammation, and Innate Immunity to Pathogens by BET Family Proteins. Mol Cell Biol (2014) 34(3):415–27. doi: 10.1128/mcb.01353-13

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Verma R, Balakrishnan L, Sharma K, Khan AA, Advani J, Gowda H, et al. A Network Map of Interleukin-10 Signaling Pathway. J Cell Commun Signal (2016) 10(1):61–7. doi: 10.1007/s12079-015-0302-x

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Lee MB, Lee JH, Hong SH, You JS, Nam ST, Kim HW, et al. JQ1, a BET Inhibitor, Controls TLR4-Induced IL-10 Production in Regulatory B Cells by BRD4-NF-κb Axis. BMB Rep (2017) 50(12):640–6. doi: 10.5483/BMBRep.2017.50.12.194

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Bolden JE, Tasdemir N, Dow LE, van Es JH, Wilkinson JE, Zhao Z, et al. Inducible In Vivo Silencing of Brd4 Identifies Potential Toxicities of Sustained BET Protein Inhibition. Cell Rep (2014) 8(6):1919–29. doi: 10.1016/j.celrep.2014.08.025

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Nakagawa A, Adams CE, Huang Y, Hamarneh SR, Liu W, von Alt KN, et al. Selective and Reversible Suppression of Intestinal Stem Cell Differentiation by Pharmacological Inhibition of BET Bromodomains. Sci Rep (2016) 6:20390. doi: 10.1038/srep20390

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: macrophage, chromatin, inflammation, regulation, Crohn’s disease

Citation: Hoffner O’Connor M, Berglind A, Kennedy Ng MM, Keith BP, Lynch ZJ, Schaner MR, Steinbach EC, Herzog J, Trad OK, Jeck WR, Arthur JC, Simon JM, Sartor RB, Furey TS and Sheikh SZ (2022) BET Protein Inhibition Regulates Macrophage Chromatin Accessibility and Microbiota-Dependent Colitis. Front. Immunol. 13:856966. doi: 10.3389/fimmu.2022.856966

Received: 17 January 2022; Accepted: 16 February 2022;
Published: 24 March 2022.

Edited by:

Guan-Jun Yang, Ningbo University, China

Reviewed by:

Giuseppina Caretti, University of Milan, Italy
Yaoyao Zou, Guangdong Provincial People’s Hospital, China
Rebecca Furze, GlaxoSmithKline, United Kingdom

Copyright © 2022 Hoffner O’Connor, Berglind, Kennedy Ng, Keith, Lynch, Schaner, Steinbach, Herzog, Trad, Jeck, Arthur, Simon, Sartor, Furey and Sheikh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Terrence S. Furey, tsfurey@email.unc.edu; Shehzad Z. Sheikh, shehzad_sheikh@med.unc.edu

These authors have contributed equally to this work and share senior authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.