Abstract
Natural killer (NK) cells are innate lymphocytes that play a pivotal role in the immune surveillance and elimination of transformed or virally infected cells. Using a chemo-genetic approach, we identify BET bromodomain containing proteins BRD2 and BRD4 as central regulators of NK cell functions, including direct cytokine secretion, NK cell contact-dependent inflammatory cytokine secretion from monocytes as well as NK cell cytolytic functions. We show that both BRD2 and BRD4 control inflammatory cytokine production in NK cells isolated from healthy volunteers and from rheumatoid arthritis patients. In contrast, knockdown of BRD4 but not of BRD2 impairs NK cell cytolytic responses, suggesting BRD4 as critical regulator of NK cell mediated tumor cell elimination. This is supported by pharmacological targeting where the first-generation pan-BET bromodomain inhibitor JQ1(+) displays anti-inflammatory effects and inhibit tumor cell eradication, while the novel bivalent BET bromodomain inhibitor AZD5153, which shows differential activity towards BET family members, does not. Given the important role of both cytokine-mediated inflammatory microenvironment and cytolytic NK cell activities in immune-oncology therapies, our findings present a compelling argument for further clinical investigation.
Introduction
Natural killer cells are cytolytic lymphocytes belonging to the innate immune system and are involved in anti-viral and anti-tumor responses () and are recognized as major players in immune-mediated anti-tumor therapies (). Their function is regulated by the activation of a number of activating and inhibitory receptors that bind to specific ligands expressed on the surface of target cells. In particular, NK cells mediate their cytolytic function through the engagement of activating receptors, such as NKG2D, DNAM-1, NKp30, NKp46, and NKp44 (, ), or following pro-inflammatory cytokine stimulation (). Following activation, NK cells mediate killing of target cells through two major pathways that require direct contact between NK cells and their target cells (). The first mechanism involves killing mediated by cytotoxic molecules (e.g. perforin and granzymes) that are stored in secretory granules of lysosomal origin (). The second pathway involves the engagement of death receptors with their ligands (e.g. Fas/FasL) that results in caspase-dependent apoptosis. Moreover, NK cells are poised to release cytokines such as IFN-γ, TNF-α and growth factors that can initiate inflammatory responses mediated by both the innate and the adaptive arm of the immune system.
Bromodomains are proteins that contain modules of ~110 amino acids that recognize and bind acetylated lysine residues in histones and other proteins. Recognition of acetylated chromatin marks by BRDs enables the regulation of gene expression through a wide range of activities. BRDs can act as scaffolds that enable the recruitment of large protein complexes or they can act as transcription factors themselves. Additionally, BRDs contain several catalytic domains that enable them to act as methyltransferases, ATP-dependent re-modellers or histone acetyltransferases and helicases () Bromodomain and extra-terminal domain (BET) proteins are a family of transcriptional mediators that regulate gene expression (, ). The BET family consists of BRD2, BRD3, BRD4, and BRDT which bind histone acetylated lysine residues via two highly conserved amino-terminal bromodomains, BD1 and BD2 found in each family member () (Figure 1). A number of BET bromodomain inhibitors have been developed to understand the utility of BETs in oncology, with each having different specificities for BD1 and BD2 (, , ). However, no BET bromodomain inhibitor can reliably distinguish between different BET family members ().
Figure 1
Despite extensive work undertaken to explore the function and mechanisms regulating BET bromodomains in myeloid and T cells (
Materials and Methods
Reagents
IL-15 (R&D systems), IL-2(PeproTech), Pam3CSK (InvivoGen), MCSF (PeproTech), GMCSF (PeproTech). CD3/28 activation beads were purchased from Invitrogen. JQ1(+) and JQ1(−) were used at a concentration of 1 µM, while AZD5153 was used at a concentration of 0.1 µM for all experiments, unless otherwise stated. Antibodies for flow cytometry were purchased from BioLegend. Media and sera were tested for endotoxin before being used in experiments.
Cell Isolation and Cell Culture Experiments
NK cells were isolated from either venous bloody obtained from healthy volunteers, or from platelet pheresis residues obtained from the Oxford National Blood Transfusion Service. Peripheral blood was obtained from Rheumatoid Arthritis patients attending the early RA clinic at Northwick Park Hospital, London. The study was approved by the London Riverside Research Ethics Committee (REC) 07/H0706/81 and the Oxford Research Ethics Committee 06/Q1606/139. All samples were obtained ethically, and the material was used in accordance with the terms of the informed consent. All patients were diagnosed according to the American College of Rheumatology (ACR) Eular 2010 criteria.
Human peripheral blood mononuclear cells were isolated by Ficol density gradient centrifugation, and CD56+ cells were isolated using the Dynabeads® UntouchedTM Human NK cell kit (Invitrogen), as per manufacturer’s instructions. Isolated NK cells were cultured in IMDM (Gibco) supplemented with 5% heat inactivated fetal calf serum.
Monocytes were isolated using the Pan Monocyte Isolation kit (Miltenyi) and were cultured in the presence of MCSF (10 ng/ml) and GMCSF (10 ng/ml) for 24 h before being stimulated with Pam3CSK (30 ng/ml). CD4+ T cells were isolated using Dynabeads® Untouched™ Human CD4 kit (Invitrogen) and cultured in the presence of IL-2 (50 ng/ml) for 24 h before being stimulated with CD3/28 activation beads (1:1 ratio, Invitrogen), as per manufacturer’s instructions. CD8+ T cells were isolated using Dynabeads® Untouched™ Human CD8 kit (Invitrogen) and then cultured in the presence of IL-2 for 24 h before being stimulated with CD3/28 activation beads (Invitrogen), as per manufacturer’s instructions.
IC50 ELISA Experiments
ELISA (eBioscience) was used to measure the concentration of IFN-γ cytokines within the cell culture supernatant, which was performed according to manufacturer’s conditions. All standards and samples were measured in triplicate.
Flow Cytometry
Cells were analyzed using a BD LSR Fortessa™ flow cytometer following staining with fluorochrome conjugated antibodies. Prior to staining, cells were fixed using the fix buffer set (Biolegend), according to the manufacturer’s instructions. Antibodies used within this study were conjugated anti-CD3 (UCHT1, Biolegend), anti-CD56 (N901, Beckman Coulter), anti-NKp30 (P30-15, Biolegend), anti-NKp44 (P44-8.1, BD), anti-NKp40 (29A1.4, Biolegend), anti-TNF-α (Mab11, eBioscience), anti-IFN-γ (485.B3, Biolegend). Annexin V/PI staining was performed according to the manufacturer’s protocol (Biolegend). The data was analyzed using FlowJo 10.7.1 software.
Quantitative RT-PCR
Total RNA was extracted using TRIzol reagent (Invitrogen) and a Direct-zol RNA miniprep kit (Zymo Research). Complimentary DNA was generated using a SuperScript II RT kit (Invitrogen), according to the manufacturer’s instructions. Reverse transcription PCR using specific primers was used to determine gene expression. Gene expression levels were normalized relative to the expression of β-actin. All samples were measured in triplicate. Specific primers used within this study are:
| Gene | Forward primer | Reverse primer |
|---|---|---|
| IFNG | TCGGTAACTGACTTGAATGTCCA | TCGCTTCCCTGTTTTAGCTGC |
| TNF | CCTCTCTCTAATCAGCCCTCTG | GAGGACCTGGGAGTAGATGAG |
| NCR1 | TGGACCCGAAGTGATCTCG | TCCTTGAGCAGTAAGAACATGC |
| NCR2 | GGCTCTCAGGCACAATCCAAG | GCTGAAGCCTCCTTACACCA |
| NCR3 | CCCCTGAGATTCGTACCCTG | CTCCACTCTGCACACGTAGAT |
Locked Nucleic-Acid Knockdown Experiments
Knockdown experiments were performed using locked nucleic-acids (LNA). LNAs are nucleic acids with a methylene bridge connecting the 2′ oxygen and 4′ carbon of the ribose ring. This results in increased stability and target affinity, which results in efficient knockdown without the requirement for transfection reagents or transduction techniques (
NK Cell Killing Experiments
NK cells killing experiments were performed as previously described (
In parallel, NK cell degranulation assay was performed as previously described (
RNA Isolation and Sequencing Library Preparation
Total RNA was extracted using TRIzol reagent (Invitrogen) and a Direct-zol RNA miniprep kit (Zymo Research) according to the manufacturer’s protocol. The quantity of RNA was assessed using a high sensitivity Tapestation (Agilent) to determine the RNA Integrity Number (RIN). All samples had a RIN value within the range of 8.0–9.5.
RNA libraries were prepared using a NEBNext® Ultra™ RNA library prep kit for Illumina® using TruSeq indexes. Libraries were prepared according to the manufacturer’s protocols. The final libraries were pooled and sequenced on a NextSeq 500 instrument (Illumina) using a paired-end run 2 × 41bp, to a minimum depth of 20 million paired-end reads/sample.
RNA Sequencing Analysis
Reads sequenced following RNA-seq library preparation were mapped to the reference genome GRCh38 (hg38) using hisat v0.1.6 (
ChIP-Seq and Library Preparation
Chromatin immunoprecipitation (ChIP) was performed as previously described (
ChIP Sequencing Analysis
In order to analyze the ChIP-seq datasets, a computational pipeline was generated using scripts from the CGAT toolkit and CGAT core workflow manager (
Single-Cell Library Preparation
Single-cell libraries were generated using the drop-seq protocol (
Single-Cell RNA Sequencing Analysis
Raw sequencing data (FASTQ files) were processed using a CGAT-core computational pipeline, which is available at https://github.com/Acribbs/scflow. Reads were demultiplexed and then aligned to the GRCh38 (hg38) assembly reference genome using Salmon alevin v1.3.0 (
Full description of the clustering performed is described in the Seurat documentation. Briefly, following quality control and filtering, gene expression for each cell was normalized and transformed. Highly variable genes that account for cellular heterogeneity were identified and cells were aligned using Harmony (
Statistical Analyses
All values are presented as means ± S.D. Mann–Whitney U test and Kruskal–Wallis with Dunn’s test were used for multiple comparisons. Wilcoxon matched-pairs test was used for paired analyses. All calculations were performed using R or GraphPad prism software.
Data and Software Availability
RNA-seq, ChIP-seq and single-cell RNA-seq datasets are deposited within GEO under the accession number GSE156423. A computational pipeline was generated using the CGAT-core workflow manager and the CGAT toolkit (
Results
BET Bromodomain Targeting Inhibits NK Cell Mediated Inflammatory Function
The hallmark of an NK cell is its ability to rapidly secrete a number of pro-inflammatory cytokines in response to activation. In a previous small molecule epigenetic compound screen, we identified a number of pan-BET bromodomain inhibitors, including JQ1(+) and PFi-1 that can reduce the production of pro-inflammatory IFN-γ cytokine (
Figure 2

JQ1(+) bromodomain inhibitor inhibits NK cell mediated inflammatory function. (A) The measurement of IFN-γ in the culture supernatant following stimulation of NK cells with IL-15 (10 ng/ml) and 24 h of treatment with JQ1(+) and the negative control JQ1(−). (B) Flow cytometry analysis of 7-AAD and Annexin-V expression as markers for cell death in NK cells, following 48 h of treatment with DMSO, JQ1(+) or JQ1(−). (C) Flow cytometry analysis of IFN-γ staining for NK cells stimulated with IL-15 (10 ng/ml) and treated for 24 h with DMSO, JQ1(+) or JQ1(−). The frequency of IFN-γ positive cells is shown for CD56brightCD57−, CD56dimCD57− and CD56dimCD57+ NK cell populations. (D) Flow cytometry analysis of CD14+ monocytes following co-culture with NK cells stimulated with IL-15 (10 ng/mL) and pre-treated with DMSO, JQ1(+) or JQ1(−) for 24 h. Cells were stained with CD14 and TNF-α and then analyzed by flow cytometry. (E) The frequency of IFN-γ positive cells in NK cells isolated from the peripheral blood of drug-naïve RA patients. (F) The frequency of IFN-γ positive cells in total NK cell populations isolated from RA synovial tissue. For (E, F), P values were calculated using a Wilcoxon matched-pairs test. For (C, D), P values were calculated using a Mann–Whitney U test. *P < 0.05. Error bars show mean ± SD (n = 3).
Knockdown of BRD2 and BRD4 Suppresses the Inflammatory Function of NK Cells
To evaluate the specificity of JQ1(+) compound on NK cell function, we used locked nucleic acid (LNA) knockdown of BRD2 and BRD4, the previously validated targets of JQ1(+). We show efficient knockdown of both BRD2 and BRD4 when compared to a scrambled control LNA, as measured using Rt-PCR (Figure 3A). Both BRD2 and BRD4 were found to significantly reduce the expression of IFN-γ mRNA following IL-15 stimulation (Figure 3B). Next we measured the expression of IFN-γ protein expression by flow cytometry. These data confirm a significant reduction in the IFN-γ protein expression following BRD2 and BRD4 knockdown (Figure 3C). The reduction in IFN-γ was apparent in all NK cell subsets (CD56bright, CD56dimCD57+ and CD56dimCD57−), with the knockdown effect being greatest in the CD56bright NK cell population, which was an expected observation (Figure 3D). Taken together, these data demonstrate that both BRD2 and BRD4 BET proteins are key regulators of the inflammatory response within NK cell populations.
Figure 3

BRD2 and BRD4 regulate NK cell inflammatory function. NK cells were cultured in the presence of IL-15 and either scrambled control (SC), BRD2 or BRD4 LNA oligonucleotide (n = 5). (A) Knockdown efficiency was assessed using qPCR. (B) The expression of IFNG gene expression was measured using qPCR. (C) The expression of IFN-γ protein in total NK cells was measured using flow cytometry. (D) NK cells were further gated based on the expression of CD56 and CD57 and then intracellular IFN-γ was measured. A representative figure shows the median fluorescent intensity of IFN-γ among CD56brightCD57−, CD56dimCD57− and CD56dimCD57+ NK cell populations. The right panel shows a representative flow cytometry histogram of IFN-γ expression in each NK cell subpopulation. The data are from three independent experiments. The data represents means ± S.D. in (A–D). P values were calculated using Kruskal–Wallis with Dunn’s multiple comparison test *P < 0.05. MFI, median fluorescence intensity.
BET Targeting in NK Cells Imparts an Anti-Inflammatory Phenotype and a Reduction in the Expression of Cytolytic Markers
The observed anti-inflammatory phenotype following treatment with JQ1(+) prompted us to investigate the global transcriptional changes in NK cells following BET bromodomain inhibition. Given the non-selective property of JQ1(+) between BET bromodomain members, we included the small molecule inhibitor AZD5153 to specifically delineate the effect of BET bromodomain inhibition in NK cells. Unlike the JQ1(+) monovalent BET inhibitor, AZD5153 interacts with both bromodomains simultaneously and has greater specificity for BRD2 than JQ1(+), which has a preferential affinity for BRD4 (
Figure 4

JQ1(+) and AZD5153 impact on NK cell inflammatory function. NK cells were cultured in the presence of IL-15 (10 ng/ml) and treated with either DMSO, JQ1(+) or JQ1(−) for 24 h. (A) A volcano plot showing the gene expression differences between JQ1(+) and DMSO treated NK cells. (B) A volcano plot showing the gene expression differences between AZD5153 and DMSO treated NK cells. (C) A Bar chart showing the number of upregulated and downregulated genes following JQ1(+) or AZD5153 treated NK cells. (D) Box plots showing the normalized counts of BRD2 and BRD4 following treatment with DMSO, JQ1(+) or AZD5153. (E) Pathway enrichment analysis showing the metacore process pathways for the differences in gene expression between JQ1(+) and AZD5153 treated NK cells. (F) Heatmap of gene expression changes in three individual NK cell donors upon treatment with either DMSO, JQ1(+) or JQ1(–). The percent expression of NKp30 (G), NKp44 (H) and NKp46 (I) as measured by flow cytometry. (J) The relative quantification of NCR3 following culture with scrambled control (SC), BRD2 or BRD4 LNA, as measured by qPCR. Data in (G–I) represents mean ± S.D. and is measured in three independent donors. P values were calculated using Kruskal–Wallis with Dunn’s multiple comparison test *P < 0.05, **P < 0.01.
Identification of Commonly Regulated Gene Expression Signatures Across Major Inflammatory Cell Populations
BET bromodomain inhibition has been shown to elicit a general anti-inflammatory effect in many immune cell types (
Figure 5

BET bromodomain inhibitors have a broad anti-inflammatory effect on the immune system. (A) A schematic showing the experimental setup. NK, CD8, CD4 and monocytes were isolated from healthy peripheral blood and then cultured in the presence of cytokines for 24 h before being stimulated with either IL-15 (10 ng/ml), CD3/28 activation beads (1:1 ratio) or MCSF/GMCSF (10 ng/ml). (B) A Venn diagram showing the overlapping differentially regulated genes for NK, CD4, CD8 and monocytes following 24 h of JQ1(+) treatment. (C) A Venn diagram showing the overlapping differentially regulated genes for each NK, CD4, CD8 and monocytes following 24 h of AZD5153 treatment. (D) A heatmap showing the common genes regulated between JQ1(+) and AZD5153 across all inflammatory cells. (E–G) Box plots showing the normalized counts of BATF, BATF3 and IKZF4 transcription factors in NK, CD4, CD8 and monocytes following JQ1(+) or AZD5153.
NK Cell Treatment With JQ1(+) but not AZD5153 Results in Compromised NK Cell Mediated Cytotoxicity
Having identified that JQ1(+) but not AZD5153 treatment inhibits a number of key NK cell mediated cytotoxicity genes, we next investigated whether this could impact on the cytolytic function of NK cells. NK cells were pre-treated with either DMSO control or BET inhibitors JQ1(+) or AZD5153 for 24 h and then washed twice before measuring cytolytic activity using CD107a surface expression and cell death of K562 target cells by Annexin-V and Propidium iodide uptake (
Figure 6

BRD2 inhibition impairs NK cell killing function. (A) NK cells were stimulated with IL-15 (10 ng/ml) and then treated with DMSO, JQ1(+) or AZD5153 for 24 h. NK cells were then co-cultured with K562 cells for 4 h and CD107a was measured by flow cytometry (B). Flow cytometry was also used to measure cell death in CD56- K562 cells by measuring the expression of Annexin-V and propidium iodide uptake. (C) An LDH release assay was used to measure NK cell mediated cell lysis of K562 cells following 24 h of treatment with DMSO, JQ1(−), JQ1(+) or AZD5153. (D) An LDH release assay showing the % cytotoxicity of Jurkat cells following 24 h of NK cell treatment with DMSO, JQ1(−), JQ1(+) or AZD5153. (E) NK cells were cultured in the presence of scrambled control (SC), BRD2 or BRD4 LNA and then CD107a expression was measured following co-culture with K562 cells. (F) The level of killing of K562 cells was determined by measuring propidium iodide uptake. (G) Blocking antibodies against IgG1, NKp30, NKp44, NKp46 were added to NK cells stimulated with IL-15 (10 ng/ml). An LDH release assay was then performed and the % cytotoxicity of K562 cells was measured using an LDH release assay. (A, B, E, F) are representative of three independent experiments. (C, D, G) show mean ± S.D. of three independent experiments. P values were calculated using Kruskal–Wallis with Dunn’s multiple comparison test *P < 0.05, **P < 0.01. n.s., not significant.
Single-Cell Sequencing Reveals a Heterogeneous Response to BET Bromodomain Inhibition in NK Cells
Next we investigated the effect of BET bromodomain inhibitors on cytolytic killing using single-cell killing assay. To examine the transcriptional response to BET inhibitors in single-cells, we pre-treated NK cells with either DMSO control, JQ1(+), or AZD5153 for 24 h and then co-cultured them in the presence of K562 target cells for 4 h. We performed droplet-based single-cell sequencing (11,970 total cells) across all three conditions (
Figure 7

Single-cell transcriptomics reveals a heterogeneous response to BET bromodomain inhibition. NK cells were stimulated with IL-15 (10ng/mL) and treated with DMSO, JQ1(+) or AZD5153 for 24 h. NK cells were then co-cultured with K562 cells for 2 h and single-cell RNA sequencing was performed. (A) A Uniform Manifold Approximation and Projection (UMAP) plot showing clustering and differences between DMSO, JQ1(+) and AZD5153 treated NK cells in co-culture with K562 cells (B) The average gene expression in highly inflammatory NK clusters 1 and 5 between DMSO, JQ1(+) and AZD5153 treated NK cells. (C) Bar charts showing the relative proportion for each cluster in all treatment groups. The right bar chart shows the relative cell proportions of NK cells to K562 cells. (D) UMAP plots showing the expression of NCR3 between DMSO, JQ1(+) and AZD5153 treated NK cells. (E) UMAP plots showing the expression of NCR1 between DMSO, JQ1(+), and AZD5153 treated NK cells.
We next compared the relative proportion of K562 to NK cell clusters to determine the level of NK cell mediated cytolytic killing (Figure 7C). In line with our killing experiments we identified a greater number of K562 cells in JQ1(+) pre-treated NK cells (27%) when compared to both AZD5153 (15%) and DMSO control (8%). Given that our previous flow cytometry experiments showed a clear reduction in the expression of NK cell activation markers following JQ1(+) treatment, we next evaluated the expression of NCR1 and NCR3. Despite generally low expression levels for these transcripts, the expression of both NCR1 and NCR3 was absent within the JQ1(+) pre-treated NK cells and thus confirms our previous bulk RNA-seq results (Figures 7D, E). Although the expression of both genes was reduced in AZD5153 treated cells, a proportion of cells still expressed both NCR1 and NCR3. Overall, these results confirm our previous bulk RNA-seq and NK cell killing experiments showing that there is a significantly attenuated cytolytic response when NK cells are treated with JQ1(+), but the effect on cytolytic function is not as apparent when cultured with AZD5153.
BET Bromodomain Inhibition Reduces Both BRD2 and BRD4 Binding to Promoters
In order to understand the mechanisms of how BRD4 and BRD2 regulate NK cell inflammatory and cytolytic function, we performed genome-wide occupancy studies using chromatin immunoprecipitation for BRD4 and BRD2, followed by sequencing. These studies revealed a significantly reduced occupancy of BRD4 and a modest but distinct reduction in BRD2 occupancy throughout the genome following. JQ1(+) treatment (Figures 8A, B). The reduction of BRD4 following JQ1(+) treatment is most striking when compared to BRD2 (Figure 8C). However, both NCR3 and IFNG show a marked decrease in both BRD4 and BRD2 occupancy following treatment with JQ1(+) (Figures 8D, E). Altogether, these results are consistent with a positive role of gene transcription for both BRD4 and BRD2 within NK cells. Moreover, the profile of BRD2 occupancy is more restricted than that of BRD4, suggesting that BRD2 may regulate fewer genes than that of BRD4.
Figure 8

Genome-wide BRD2 and BRD4 occupancy following stimulation of NK cells with IL-15 (10 ng/ml) and 24 h of JQ1(+) treatment of NK cells. (A) Coverage plot showing average read counts centered around the transcriptional start site (TSS) from BRD4 ChIP-seq following DMSO or JQ1(+) treatment of NK cells. (B) Coverage plot showing average read counts centered around the TSS from BRD2 ChIP-seq following DMSO or JQ1(+) treatment of NK cells. (C) Heatmap of BRD4 and BRD2 following treatment of NK cells with DMSO or JQ1(+) covering 1 kb upstream and downstream of the center of called peaks. (D) Representative ChIP-seq coverage profiles for BRD2 and BRD4 across the NCR3 locus. (E) Representative ChIP-seq coverage profiles for BRD2 and BRD4 across the IFNG locus.
Discussion
Our work provides novel insight into the epigenetic regulation of NK cells within an inflammation and cancer context. A handful of studies identified DNA methylation (
We were also able to show that a second BET bromodomain inhibitor, AZD5153 also imparts an anti-inflammatory effect on NK cells. Our single-cell experiments were able to show a significant response of specific NK cell subsets to both JQ1(+) and AZD5153 treatment. In these experiments, highly inflammatory NK cells show the highest level of gene transcriptional changes, with a significant reduction in a broad range of inflammatory genes. Our bulk RNA-seq data revealed that JQ1(+) treatment led to a downregulated NK cell cytolytic response, which impacted NK cell killing of cancer cells, in part through the downregulation of NK cell activating receptors. Consistent with these findings, JQ1 has been shown to downregulate target ligands in neuroblastoma cells and downregulate NK cell activation receptors NKG2D and DNAM-1 on NK cells (
JQ1(+) is a monovalent BET inhibitor that does not engage simultaneously the two bromodomains found in the BET family members. It has been shown to be more effective in hematopoietic cancers, but is far less effective in other solid tumors, such as breast cancer and cervical cancers (
Much more is known about the regulation of gene expression by BRD4 than BRD2 and very little is known about both of their functions in NK cell biology. BRD4 recruits PTEF-b to sites of active transcription, which in turn phosphorylates RNA polymerase II Ser2 during elongation (
Funding
The study was supported through funding from the Kennedy Trust for Rheumatology Research, the National Institute for Health Research Oxford Biomedical Research Unit (UO, MF), Cancer Research UK (CRUK, UO), Arthritis Research UK (20522) (UO), the Bone Cancer Research Trust (BCRT) (AC and UO), and a Leducq Epigenetics of Atherosclerosis Network program grant from the Leducq Foundation (UO). UO was supported by the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007–2013) under REA grant agreement n° [609305]. AC was supported by the Medical Research Council (MRC) CGAT program (G1000902), MRC Career Development Fellowship (MR/V010182/1) and a CRUK Oxford Centre Development Fund Award (CRUKDF-0318-AC[AZ]), with support from AstraZeneca. PF was supported by the Medical Research Council (MR/N010051/1).
Statements
Data availability statement
RNA-seq, ChIP-seq and single-cell RNA-seq datasets are deposited within GEO under the accession number GSE156423. A computational pipeline was generated using the CGAT-core workflow manager and the CGAT toolkit (
Ethics statement
The studies involving human participants were reviewed and approved by the Oxford Research Ethics Committee and the London Riverside Research Ethics Committee (REC numbers: 07/H0706/81 and 06/Q1606/139). The patients/participants provided their written informed consent to participate in this study.
Author contributions
AC, MF, and UO: conceptualization and designed research. AC, HO, VV-A, and UO contributed new reagents/analytical tools. AC, PF, MP, and GW: performed research. HP: provided patient samples and metadata. AC, PF, and UO analyzed data. AC and UO: writing-original manuscript draft. AC and UO: writing-review and editing. All authors contributed to the article and approved the submitted version.
Acknowledgments
We would like to thank all patients who donated blood for this study and Jingwen Zhang for her comments during the drafting of this manuscript.
Conflict of interest
VV-A was employed by AstraZeneca. HO was employed by Roche Innovation Center Copenhagen A/S.
The remaining 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.
The authors declare that this study received funding from AstraZeneca. The funder was involved in the interpretation of the data and reviewed the article before the decision to submit it for publication.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.626255/full#supplementary-material
References
1
TrinchieriG. Biology of natural killer cells. Adv Immunol (1989) 47:187–376. doi: 10.1016/S0065-2776(08)60664-1
2
NicholsonSEKeatingNBelzGT. Natural killer cells and anti-tumor immunity. Mol Immunol (2019) 110:40–7. doi: 10.1016/j.molimm.2017.12.002
3
MorettaLBottinoCPendeDCastriconiRMingariMCMorettaA. Surface NK receptors and their ligands on tumor cells. Semin Immunol (2006) 18:151–8. doi: 10.1016/j.smim.2006.03.002
4
LanierLL. NK cell recognition. Annu Rev Immunol (2005) 23:225–74. doi: 10.1146/annurev.immunol.23.021704.115526
5
ZwirnerNWDomaicaCI. Cytokine regulation of natural killer cell effector functions. BioFactors (2010) 36:274–88. doi: 10.1002/biof.107
6
SmythMJCretneyEKellyJMWestwoodJAStreetSEYagitaHet al. Activation of NK cell cytotoxicity. Mol Immunol (2005) 42:501–10. doi: 10.1016/j.molimm.2004.07.034
7
TophamNJHewittEW. Natural killer cell cytotoxicity: how do they pull the trigger? Immunology (2009) 128:7–15. doi: 10.1111/j.1365-2567.2009.03123.x
8
FujisawaTFilippakopoulosP. Functions of bromodomain-containing proteins and their roles in homeostasis and cancer. Nat Rev Mol Cell Biol (2017) 18:246–62. doi: 10.1038/nrm.2016.143
9
CochranAGConeryARSimsRJ. 3rd, Bromodomains: a new target class for drug development. Nat Rev Drug Discovery (2019) 18:609–28. doi: 10.1038/s41573-019-0030-7
10
ShiJVakocCR. The mechanisms behind the therapeutic activity of BET bromodomain inhibition. Mol Cell (2014) 54:728–36. doi: 10.1016/j.molcel.2014.05.016
11
RenCZhangGHanFFuSCaoYZhangFet al. Spatially constrained tandem bromodomain inhibition bolsters sustained repression of BRD4 transcriptional activity for TNBC cell growth. Proc Natl Acad Sci USA (2018) 115:7949–54. doi: 10.1073/pnas.1720000115
12
FilippakopoulosPQiJPicaudSShenYSmithWBFedorovOet al. Selective inhibition of BET bromodomains. Nature (2010) 468:1067–73. doi: 10.1038/nature09504
13
BradburyRHCallisRCarrGRChenHClarkEFeronLet al. Optimization of a Series of Bivalent Triazolopyridazine Based Bromodomain and Extraterminal Inhibitors: The Discovery of (3R)-4-[2-[4-[1-(3-Methoxy-[1,2,4]triazolo[4,3-b]pyridazin-6-yl)-4-piperidyl]phen oxy]ethyl]-1,3-dimethyl-piperazin-2-one (AZD5153). J Med Chem (2016) 59:7801–17. doi: 10.1021/acs.jmedchem.6b00070
14
NicodemeEJeffreyKLSchaeferUBeinkeSDewellSChungCWet al. Suppression of inflammation by a synthetic histone mimic. Nature (2010) 468:1119–23. doi: 10.1038/nature09589
15
RhyasenGWHattersleyMMYaoYDulakAWangWPetterutiPet al. AZD5153: A Novel Bivalent BET Bromodomain Inhibitor Highly Active against Hematologic Malignancies. Mol Cancer Ther (2016) 15:2563–74. doi: 10.1158/1535-7163.MCT-16-0141
16
FilippakopoulosPKnappS. Targeting bromodomains: epigenetic readers of lysine acetylation. Nat Rev Drug Discovery (2014) 13:337–56. doi: 10.1038/nrd4286
17
KagoyaYNakatsugawaMYamashitaYOchiTGuoTAnczurowskiMet al. BET bromodomain inhibition enhances T cell persistence and function in adoptive immunotherapy models. J Clin Invest (2016) 126:3479–94. doi: 10.1172/JCI86437
18
GibbonsHRMiDJFarleyVMEsmondTKaoodMBAuneTM. Bromodomain inhibitor JQ1 reversibly blocks IFN-gamma production. Sci Rep (2019) 9:10280. doi: 10.1038/s41598-019-46516-x
19
SteinCAHansenJBLaiJWuSVoskresenskiyAHogAet al. Efficient gene silencing by delivery of locked nucleic acid antisense oligonucleotides, unassisted by transfection reagents. Nucleic Acids Res (2010) 38:e3. doi: 10.1093/nar/gkp841
20
CribbsAHookwayESWellsGLindowMObadSOerumHet al. Inhibition of histone H3K27 demethylases selectively modulates inflammatory phenotypes of natural killer cells. J Biol Chem (2018) 293:2422–37. doi: 10.1074/jbc.RA117.000698
21
AlterGMalenfantJMAltfeldM. CD107a as a functional marker for the identification of natural killer cell activity. J Immunol Methods (2004) 294:15–22. doi: 10.1016/j.jim.2004.08.008
22
KimDLangmeadBSalzbergSL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods (2015) 12:357–60. doi: 10.1038/nmeth.3317
23
LiaoYSmythGKShiW. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics (2014) 30:923–30. doi: 10.1093/bioinformatics/btt656
24
OrlandoDAChenMWBrownVESolankiSChoiYJOlsonERet al. Quantitative ChIP-Seq normalization reveals global modulation of the epigenome. Cell Rep (2014) 9:1163–70. doi: 10.1016/j.celrep.2014.10.018
25
SimsDIlottNESansomSNSudberyIMJohnsonJSFawcettKAet al. CGAT: computational genomics analysis toolkit. Bioinformatics (2014) 30:1290–1. doi: 10.1093/bioinformatics/btt756
26
CribbsALuna-ValeroSGeorgeCSudberyIBerlanga-TaylorASansomSet al. CGAT-core: a python framework for building scalable, reproducible computational biology workflows version 2; peer review: 1 approved, 1 approved with reservations]. F1000Research (2019) 8:377–90. doi: 10.12688/f1000research.18674.1
27
LangmeadBTrapnellCPopMSalzbergSL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol (2009) 10:R25. doi: 10.1186/gb-2009-10-3-r25
28
KentWJZweigASBarberGHinrichsASKarolchikD. BigWig and BigBed: enabling browsing of large distributed datasets. Bioinformatics (2010) 26:2204–7. doi: 10.1093/bioinformatics/btq351
29
MacoskoEZBasuASatijaRNemeshJShekharKGoldmanMet al. Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell (2015) 161:1202–14. doi: 10.1016/j.cell.2015.05.002
30
SrivastavaAMalikLSmithTSudberyIPatroR. Alevin efficiently estimates accurate gene abundances from dscRNA-seq data. Genome Biol (2019) 20:65. doi: 10.1186/s13059-019-1670-y
31
StuartTButlerAHoffmanPHafemeisterCPapalexiEMauckWM,3et al. Comprehensive Integration of Single-Cell Data. Cell (2019) 177:1888–902 e21. doi: 10.1016/j.cell.2019.05.031
32
KorsunskyIMillardNFanJSlowikowskiKZhangFWeiKet al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods (2019) 16:1289–96. doi: 10.1038/s41592-019-0619-0
33
RouillardADGundersenGWFernandezNFWangZMonteiroCDMcDermottMGet al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database (Oxford) (2016) 2016. doi: 10.1093/database/baw100
34
KerscherBBarlowJLRanaBMJolinHEGogoiMBartholomewMAet al. BET Bromodomain Inhibitor iBET151 Impedes Human ILC2 Activation and Prevents Experimental Allergic Lung Inflammation. Front Immunol (2019) 10:678. doi: 10.3389/fimmu.2019.00678
35
KleinK. Bromodomain protein inhibition: a novel therapeutic strategy in rheumatic diseases. RMD Open (2018) 4:e000744. doi: 10.1136/rmdopen-2018-000744
36
McInnesLHealyJUMAPMJ,. Uniform Manifold Approximation and Projection for Dimension Reduction. arXiv (2018) 1:63. doi: 10.21105/joss.00861
37
Luetke-EverslohMCicekBBSiracusaFThomJTHamannAFrischbutterSet al. NK cells gain higher IFN-gamma competence during terminal differentiation. Eur J Immunol (2014) 44:2074–84. doi: 10.1002/eji.201344072
38
RossiLEAvilaDESpallanzaniRGZiblatAFuertesMBLapyckyjLet al. Histone deacetylase inhibitors impair NK cell viability and effector functions through inhibition of activation and receptor expression. J Leukoc Biol (2012) 91:321–31. doi: 10.1189/jlb.0711339
39
LeeSParkHYKangSYKimSJHwangJLeeSet al. Genetic alterations of JAK/STAT cascade and histone modification in extranodal NK/T-cell lymphoma nasal type. Oncotarget (2015) 6:17764–76. doi: 10.18632/oncotarget.3776
40
VenezianiIFruciDCompagnoneMPistoiaVRossiPCifaldiL. The BET-bromodomain inhibitor JQ1 renders neuroblastoma cells more resistant to NK cell-mediated recognition and killing by downregulating ligands for NKG2D and DNAM-1 receptors. Oncotarget (2019) 10:2151–60. doi: 10.18632/oncotarget.26736
41
ZuberJShiJWangERappaportARHerrmannHSisonEAet al. RNAi screen identifies Brd4 as a therapeutic target in acute myeloid leukaemia. Nature (2011) 478:524–8. doi: 10.1038/nature10334
42
MertzJAConeryARBryantBMSandyPBalasubramanianSMeleDAet al. 3rd, Targeting MYC dependence in cancer by inhibiting BET bromodomains. Proc Natl Acad Sci USA (2011) 108:16669–74. doi: 10.1073/pnas.1108190108
43
DawsonMAPrinjhaRKDittmannAGiotopoulosGBantscheffMChanWIet al. Inhibition of BET recruitment to chromatin as an effective treatment for MLL-fusion leukaemia. Nature (2011) 478:529–33. doi: 10.1038/nature10509
44
ChaidosACaputoVGouvedenouKLiuBMarigoIChaudhryMSet al. Potent antimyeloma activity of the novel bromodomain inhibitors I-BET151 and I-BET762. Blood (2014) 123:697–705. doi: 10.1182/blood-2013-01-478420
45
AsanganiIADommetiVLWangXMalikRCieslikMYangRet al. Therapeutic targeting of BET bromodomain proteins in castration-resistant prostate cancer. Nature (2014) 510:278–82. doi: 10.1038/nature13229
46
ShiJWangYZengLWuYDengJZhangQet al. Disrupting the interaction of BRD4 with diacetylated Twist suppresses tumorigenesis in basal-like breast cancer. Cancer Cell (2014) 25:210–25. doi: 10.1016/j.ccr.2014.01.028
47
LeeMBLeeJHHongSHYouJSNamSTKimHWet al. JQ1, a BET inhibitor, controls TLR4-induced IL-10 production in regulatory B cells by BRD4-NF-kappaB axis. BMB Rep (2017) 50:640–6. doi: 10.5483/BMBRep.2017.50.12.194
48
BandukwalaHSGagnonJTogherSGreenbaumJALampertiEDParrNJet al. Selective inhibition of CD4+ T-cell cytokine production and autoimmunity by BET protein and c-Myc inhibitors. Proc Natl Acad Sci USA (2012) 109:14532–7. doi: 10.1073/pnas.1212264109
49
ItzenFGreifenbergAKBoskenCAGeyerM. Brd4 activates P-TEFb for RNA polymerase II CTD phosphorylation. Nucleic Acids Res (2014) 42:7577–90. doi: 10.1093/nar/gku449
50
YoshidaHBansalKSchaeferUChapmanTRiojaIProektIet al. Brd4 bridges the transcriptional regulators, Aire and P-TEFb, to promote elongation of peripheral-tissue antigen transcripts in thymic stromal cells. Proc Natl Acad Sci USA (2015) 112:E4448–57. doi: 10.1073/pnas.1512081112
51
RahmanSSowaMEOttingerMSmithJAShiYHarperJWet al. The Brd4 extraterminal domain confers transcription activation independent of pTEFb by recruiting multiple proteins, including NSD3. Mol Cell Biol (2011) 31:2641–52. doi: 10.1128/MCB.01341-10
52
DenisGVMcCombMEFallerDVSinhaARomesserPBCostelloCE. Identification of transcription complexes that contain the double bromodomain protein Brd2 and chromatin remodeling machines. J Proteome Res (2006) 5:502–11. doi: 10.1021/pr050430u
53
PengJDongWChenLZouTQiYLiuY. Brd2 is a TBP-associated protein and recruits TBP into E2F-1 transcriptional complex in response to serum stimulation. Mol Cell Biochem (2007) 294:45–54. doi: 10.1007/s11010-006-9223-6
54
GyurisADonovanDJSeymourKALovascoLASmilowitzNRHalperinALet al. The chromatin-targeting protein Brd2 is required for neural tube closure and embryogenesis. Biochim Biophys Acta (2009) 1789:413–21. doi: 10.1016/j.bbagrm.2009.03.005
55
VelisekLShangEVeliskovaJChachuaTMacchiaruloSMaglakelidzeGet al. GABAergic neuron deficit as an idiopathic generalized epilepsy mechanism: the role of BRD2 haploinsufficiency in juvenile myoclonic epilepsy. PloS One (2011) 6:e23656. doi: 10.1371/journal.pone.0023656
56
BelkinaACNikolajczykBSDenisGV. BET protein function is required for inflammation: Brd2 genetic disruption and BET inhibitor JQ1 impair mouse macrophage inflammatory responses. J Immunol (2013) 190:3670–8. doi: 10.4049/jimmunol.1202838
57
CheungKLZhangFJaganathanASharmaRZhangQKonumaTet al. Distinct Roles of Brd2 and Brd4 in Potentiating the Transcriptional Program for Th17 Cell Differentiation. Mol Cell (2017) 65:1068–80.e5. doi: 10.1016/j.molcel.2016.12.022
58
HildnerKEdelsonBTPurthaWEDiamondMMatsushitaHKohyamaMet al. Batf3 deficiency reveals a critical role for CD8alpha+ dendritic cells in cytotoxic T cell immunity. Science (2008) 322:1097–100. doi: 10.1126/science.1164206
59
SopelNGraserAMoussetSFinottoS. The transcription factor BATF modulates cytokine-mediated responses in T cells. Cytokine Growth Factor Rev (2016) 30:39–45. doi: 10.1016/j.cytogfr.2016.03.004
60
BetzBCJordan-WilliamsKLWangCKangSGLiaoJLoganMRet al. Batf coordinates multiple aspects of B and T cell function required for normal antibody responses. J Exp Med (2010) 207:933–42. doi: 10.1084/jem.20091548
61
CheungKLuGSharmaRVincekAZhangRPlotnikovANet al. BET N-terminal bromodomain inhibition selectively blocks Th17 cell differentiation and ameliorates colitis in mice. Proc Natl Acad Sci USA (2017) 114:2952–7. doi: 10.1073/pnas.1615601114
Summary
Keywords
BET bromodomain, BRD2, BRD4, NK cell, epigenetics (DNA methylation histone modifications)
Citation
Cribbs AP, Filippakopoulos P, Philpott M, Wells G, Penn H, Oerum H, Valge-Archer V, Feldmann M and Oppermann U (2021) Dissecting the Role of BET Bromodomain Proteins BRD2 and BRD4 in Human NK Cell Function. Front. Immunol. 12:626255. doi: 10.3389/fimmu.2021.626255
Received
05 November 2020
Accepted
13 January 2021
Published
26 February 2021
Volume
12 - 2021
Edited by
Tanapat Palaga, Chulalongkorn University, Thailand
Reviewed by
Michaela Semeraro, Assistance Publique Hopitaux De Paris, France; Marco Di Gioia, Boston Children’s Hospital and Harvard Medical School, United States
Updates

Check for updates
Copyright
© 2021 Cribbs, Filippakopoulos, Philpott, Wells, Penn, Oerum, Valge-Archer, Feldmann and Oppermann.
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: Adam P. Cribbs, adam.cribbs@ndorms.ox.ac.uk; Udo Oppermann, udo.oppermann@ndorms.ox.ac.uk
This article was submitted to Molecular Innate Immunity, a section of the journal Frontiers in Immunology
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.