- Open Access
Bromodomain and extra-terminal domain (BET) proteins regulate melanocyte differentiation
Epigenetics & Chromatin volume 13, Article number: 14 (2020)
Pharmacologic inhibition of bromodomain and extra-terminal (BET) proteins is currently being explored as a new therapeutic approach in cancer. Some studies have also implicated BET proteins as regulators of cell identity and differentiation through their interactions with lineage-specific factors. However, the role of BET proteins has not yet been investigated in melanocyte differentiation. Melanocyte inducing transcription factor (MITF) is the master regulator of melanocyte differentiation, essential for pigmentation and melanocyte survival. In this study, we tested the hypothesis that BET proteins regulate melanocyte differentiation through interactions with MITF.
Here we show that chemical inhibition of BET proteins prevents differentiation of unpigmented melanoblasts into pigmented melanocytes and results in de-pigmentation of differentiated melanocytes. BET inhibition also slowed cell growth, without causing cell death, increasing the number of cells in G1. Transcriptional profiling revealed that BET inhibition resulted in decreased expression of pigment-specific genes, including many MITF targets. The expression of pigment-specific genes was also down-regulated in melanoma cells, but to a lesser extent. We found that RNAi depletion of the BET family members, bromodomain-containing protein 4 (BRD4) and bromodomain-containing protein 2 (BRD2) inhibited expression of two melanin synthesis enzymes, TYR and TYRP1. Both BRD4 and BRD2 were detected on melanocyte promoters surrounding MITF-binding sites, were associated with open chromatin structure, and promoted MITF binding to these sites. Furthermore, BRD4 and BRD2 physically interacted with MITF.
These findings indicate a requirement for BET proteins in the regulation of pigmentation and melanocyte differentiation. We identified changes in pigmentation specific gene expression that occur upon BET inhibition in melanoblasts, melanocytes, and melanoma cells.
The bromodomain and extra-terminal (BET) family of proteins contains BRDT, BRD2, BRD3, and BRD4 . BRD2, BRD3, and BRD4 are ubiquitously expressed while BRDT is expressed mainly in testes. BET proteins have two bromodomains, defined as motifs of 110 amino acids that recognize acetylated lysines, as well as a conserved extra-terminal domain which mediates additional protein–protein interactions. BET family members play important roles in transcription by binding acetylated histones and interacting with components of the transcriptional machinery and chromatin remodeling enzymes. They can have overlapping, distinct, and even opposing roles in transcription depending on the promoter and pathway [2,3,4]. BRD2 and BRD4 have been shown to interact with the mediator complex while BRD4 can also interact with P-TEFb to promote transcription [5,6,7,8]. BRD4 contributes to the regulation of nucleosome and higher order chromatin structure [9, 10] but by a different mechanism than does BRD2, which plays a role unique from BRD4 by interacting with CTCF .
In recent years, BET proteins have emerged as drivers of tumorigenesis in diverse human cancers. Translocations that generate oncogenic BRD4 or BRD3 fusion proteins have been identified in NUT-midline carcinoma [12, 13]. In other malignancies, BRD4 associates with cancer-specific super enhancers, hematopoietic transcription factors in acute myeloid leukemia, YAP/TAZ, and other pro-tumorigenic transcription factors to promote proliferation and cancer progression [14, 15]. The discovery of the first BET inhibitor, JQ1 made it possible to target BET proteins as a therapeutic strategy against cancer . JQ1 is a stereo-specific, cell-permeable small molecule that binds competitively to acetyl-lysine recognition motifs present in BET proteins which removes them from chromatin and disrupts their tumorigenic functions [14, 17, 18]. Clinical and pre-clinical studies using JQ1 and next generation BET inhibitors are in progress for hematological cancers, and solid tumors including NUT-midline carcinoma and melanoma .
Although recent studies have elucidated BET proteins as potential therapeutic targets in melanoma, their function in normal melanocytes is not known. The melanocyte lineage arises from neural crest cell derived melanoblasts which differentiate into pigmented melanocytes and give skin and hair their characteristic color . The pigment, melanin, is synthesized in organelles called melanosomes and then transferred to surrounding keratinocytes in a process known as melanogenesis . This involves a series of chemical reactions catalyzed by tyrosinase (TYR) and tyrosinase related proteins (TYRP1). Additional proteins regulate melanosome structure, function, and distribution. The level of pigmentation is determined by the amount and type of melanin, as well as the size, number, and distribution of melanosomes. These are all determined by genetic factors and modulated by physiological and environmental conditions. The master regulator of melanocyte development, melanocyte inducing transcription factor (MITF, historically also called microphthalmia-associated transcription factor) coordinates with other melanocyte-specific transcription factors to regulate melanogenesis in normal melanocytes and play critical roles in melanoma . An increasing number of studies show that epigenetic regulators are also required for regulation of both normal melanocyte function and melanoma tumorigenicity . Furthermore, inherent features of the melanocyte differentiation program contribute to melanoma progression and resistance to therapeutics [24, 25]. Therefore, a better understanding of the function of BET proteins in melanocyte differentiation may provide additional insight into their function in melanoma.
In this study, we tested the hypothesis that BET proteins regulate melanocyte differentiation by interacting with MITF. We treated melanoblasts with the BET inhibitor, JQ1 and found that JQ1 prevented melanoblasts from differentiating into pigmented melanocytes. BET inhibition decreased expression of pigment-specific genes and slowed proliferation. Our experiments also revealed that both BRD4 and BRD2 interact with MITF and contribute to the regulation of melanocyte-specific gene expression.
Materials and methods
Melb-a cells (melanoblasts) were maintained in growth medium (RPMI, 10% fetal calf serum, 40 pM fibroblast growth factor and 10 ng/ml stem cell factor) and differentiated in differentiation medium [DMEM, 10% fetal calf serum, 2 nM alpha melanocyte stimulating hormone (NDP-alpha MSH) and 200 nM phorbol-myristate-acetate] as previously described . Neonatal human epidermal melanocytes (NHEMs) were isolated and cultured as previously described . Human embryonic kidney (HEK) 293T cells were obtained from ATCC. They were cultured in DMEM Media (Life Technologies, Grand Island, NY, USA) containing 10% fetal bovine serum. The pair of active (+) [(S)-enantiomer] and inactive (−) [(R)-enantiomer] JQ1 stereoisomers were purchased from BPS Biosciences (San Diego, CA, USA) and each dissolved in DMSO.
Cells were trypsinized and resuspended in media at the indicated times. Cell counts were taken using the Scepter 2.0 handheld automated cell counter (Millipore-Sigma, Burlington, MA, USA).
Cell cycle analysis
Approximately 1 × 106 cells were fixed with 100% ethanol for 1 h, stained with a propidium iodide (PI)-RNAse solution (Millipore-Sigma, Burlington, MA, USA) for 30 min and loaded on a FACS-Calibur (BD Biosciences, Franklin Lakes, NJ, USA) at the University of Toledo Flow Cytometry Core Facility. Data were analyzed using Cell Quest Pro (BD Biosciences).
Cells were stained with the Guava Nexin Annexin V reagent (Millipore-Sigma) as we previously described and analyzed on a FACS-Calibur (BD Biosciences) at the University of Toledo Flow Cytometry Core Facility .
The melanin assay was performed as previously described . Cells were counted and lysed in 0.1 M NaOH for 30 min. Absorbance was read using a spectrophotometer (Molecular Devices) at 475 nM and melanin content was calculated based on a standard curve using synthetic melanin. Melanin content was then normalized to cell number.
Total RNA was prepared from vehicle (DMSO) or (+)JQ1 treated Melb-a melanoblasts by Trizol extraction (Invitrogen, Carlsbad, CA, USA) as previously described . All experiments were performed in triplicate. Library preparation, quality assessment, and sequencing were performed by Bio Basics Inc. (Markham, Ontario, Canada). mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. First strand cDNA was synthesized using random hexamers and M-MuLV Reverse Transcriptase. Library fragments were purified and PCR amplified. Paired-end sequencing was performed on Illumina HiSeq 4000. The raw reads were filtered by removing reads containing N > 10% (where N represents bases that could not be determined). The Q score (quality value) of over 50% bases of the read was less than or equal to 5.
Raw RNA-seq data were converted to sequence reads using CASAVA. Raw reads were filtered to remove low-quality reads and reads with adaptors. The clean reads were mapped to the Ensembl mouse genome (grcm38/mm10) using Tophat2 (v2.0.9) . Differential gene expression between vehicle and (+)JQ1-treated samples was determined using the DESeq2 package 126.96.36.199 . The resulting p values were adjusted using Benjamini and Hochberg’s approach for controlling false discovery rate (FDR) and set at p < 0.05. Gene ontology (GO) analysis was performed using the GOseq R package . In addition, differential gene expression was analyzed with Kegg pathway analysis . Comparative analyses with previously published data sets were performed with the CLC Genomic Workbench (Qiagen, Hilden, Germany). For this comparative analysis, mouse and human gene homologues were matched using HOM Mouse Human Sequence.rpt (accessed from Jackson laboratories, http://www.informatics.jax.org; accessed 03/29/2019), two tab delimited files of Ensembl homog IDs and gene names (human and mouse, accessed from grch37.ensembl.org/biomart; accessed 03/29/2019), and when necessary, manual matching of homologues using the online directory GeneCards.
Published ChIP-seq datasets in GEO and Encode in bigWig format were visualized directly as Custom Tracks on the UCSC Genome Browser. To determine the number of overlapping peaks for MITF and BRD4, the GSM2527370 dataset was used for MITF, and the GSM2700494 dataset was used for BRD4. The BED files were downloaded from the Cistrome database: http://cistrome.org/db/#/ and imported into RStudio as GRanges. The overlapping peaks were identified using the peakPermTest function in the ChIPpeakAnno package available through Bioconductor. The peakPermTest function was used to perform a permutation test in order to determine the statistical significance of the overlap.
RNA isolation and quantitative real-time PCR
Total RNA was reverse transcribed into cDNA using the Quantitect Reverse Transcription kit (Qiagen). Quantitative PCR (qPCR) was performed in SYBR Green master mix (Qiagen) with an Applied Biosystems 7500 PCR and analyzed with the SDS software as described . Mouse Tyr, Tyrp1, Trpm1, Brd2, and Brd4 mRNA levels were normalized to mouse Rpl7. Human TYRP1 mRNA levels were normalized to human RPL19. All primers were obtained from Integrated DNA Technologies (IDT) (Coralville, IA, USA). Primer sequences for murine genes were: Tyr: forward: 5′-TTC AAA GGG GTG GAT GAC CG-3′ and reverse: 5′-GAC ACA TAG TAA TGC ATC C-3′, Tyrp1: forward: 5′-GCC CCA ACT CTG TCT TTT CTC AAT-3′ and reverse: 5′-GAT CGG CGT TAT ACC TCC TTA GC-3′, Trpm1: forward: 5′-CCT ACG ACA CCA AGC CAG AT-3′ and reverse: 5′-GAC GAC ACC AGT GCT CAC AC-3′, Brd2: forward: 5′-GCCCTTCTATAAGCCAGTGG-3′ and reverse: 5′-AAACTCCTGTGCATCCCG-3′, Brd4: forward: 5′-TAA AAA CTC CAA CCC CGA TG-3′ and reverse: 5′-TGC TCT CCG ACT CAG AGG AT-3′, Rrpl7: forward: 5′-GGA GGA AGC TCA TCT ATG AGA AGG-3′ and reverse: 5′-AAG ATC TGT GGA AGA GGA AGG AGC-3′. Primer sequences for human genes were: TYRP1: forward: 5′-TGG GAT CCA GAA GCA ACT TT-3′ and reverse: 5′-TGT GGT TCA GGA AGA CGT TG-3′, RPL19: forward: 5′-AAA CAA GCG GAT TCT CAT GG-3′ and reverse: 5′-TTG GTC TCT TCC TCC TTG GAT-3′.
The MITF (ab12039) and BRD4 (ab128874) antibodies were from Abcam (Cambridge, MA, USA). The BRD2 antibody (A302-583A-M) was from Bethyl (Montgomery County, Texas, USA). The BRD3 (sc-81202), TYR (sc-7833), and TYRP1 (sc-10443) antibodies were from Santa Cruz Technologies (Santa Cruz, CA, USA). The Tubulin (21485) antibody was from Cell Signaling Technology (Boston, MA, USA).
Cell extracts and immunoblots
Cell extracts were prepared and immunoblots were performed as described . Briefly, cell extracts were prepared in Triton Lysis Buffer [20 nM Tris–HCl, pH 7.4, 150 nM NaCl, 2 mM EDTA, 10% glycerol, 2 mM DTT, 1 mM PMSF, protease cocktail (Millipore-Sigma, Burlington, MA, USA), phosphatase inhibitor cocktail (Millipore-Sigma, Burlington, MA, USA)]. Protein transfer was performed in a buffer containing 20% methanol and the blots were probed with primary antibodies. Bands were visualized by enhanced chemiluminescence after incubation with species-matched secondary antibodies.
Control siRNA (211631978) and siRNAs targeting BRD2 (mm.Ri.Brd2.13.1), BRD3 (mm.Ri.Brd3.13.1), and BRD4 (mm.Ri.Brd4.13.1) were purchased from Integrated DNA Technology (Coralville, IA, USA). For some experiments, a pool of siRNAs targeting mouse Brd4 (M-041493), mouse BRD2 (M-043404), and a non-targeting siRNA (5′-UUCUCCGAACGUGUCACGU-3′) were obtained from Dharmacon (Lafayette, CO, USA). Transfection using Dharmafect 1 was performed on cells in growth media according to manufacturer’s instructions. After 48 to 72 h, growth media was replaced with differentiation media and cells were then harvested 48 h later.
HEK 293T cells were transfected with FLAG-tagged MITF construct using Lipofectamine LTX (Invitrogen, Carlsbad, CA, USA). Cells were harvested 48 h post-transfection. For endogenous co-immunoprecipitations, Melb-a cells were differentiated and harvested after 6 h. Co-immunoprecipitations were performed as previously described . Briefly, cell extracts were prepared in Triton lysis buffer and precleared for 30 min. 293T cells were transfected FLAG-tagged MITF were incubated with anti-FLAG M2 agarose (A2220) or normal mouse IgG agarose (AO919) (Millipore-Sigma, Burlington, MA, USA) Lysates from Melb-a cells were incubated overnight with an antibody to MITF, BRD2, BRD4, or a species-matched control antibody which was either normal mouse IgG (sc2025) or normal rabbit IgG (sc2027) (Santa Cruz Technologies). The lysates were then incubated for 4 h with Protein G or A Sepharose (GE Healthcare Life Sciences, Marlborough, MA, USA), and the beads washed three times with Triton Lysis Buffer. Immunocomplexes were eluted at 50 °C for 20 min in SDS loading buffer.
Chromatin immunoprecipitation studies (ChIPs)
Chromatin immunoprecipitation studies were performed with a control rabbit IgG, BRD2, or BRD4 antibody as previously described . Briefly, cells were cross-linked with formaldehyde (1% final concentration) for 10 min at room temperature in culture media. The reaction was then stopped by adding glycine to a final concentration of 125 mM. Cells were washed twice with phosphate buffered saline and nuclei were isolated in a buffer containing 50 mM HEPES–KOH, pH 7.5, 140 mM NaCl, 1 mM EDTA, 10% glycerol, 0.5% NP40, 0.25% TritonX, and protease inhibitors. Nuclei were washed with a buffer containing 10 mM Tris–HCl, pH 8.0, 200 mM NaCl, 1 mM EDTA, 0.5 mM EGTA and protease inhibitors. Nuclei were sonicated in 20 s pulses for 4 min in a buffer containing 10 mM Tris–HCl, pH 8.0. 100 mM NaCl, 9 mM EDTA, 0.5 mM EGTA, 0.1% Na-deoxycholate, 0.5% N-lauroylsarcosine and 1.1% Triton X-100. Complexes were incubated overnight with antibodies and then incubated with Protein A or G Sepharose (GE Healthcare Life Sciences, Marlborough, MA, USA). Complexes were washed five times in RIPA buffer (50 mM HEPES, 500 mM LiCl, 0.1 mM EDTA, 1.0% NP-40 and 0.7% Na-deoxycholate) and once with TE containing 50 mM NaCl. The immune complexes were eluted with 50 mM Tris [pH 8.0], 10 mM EDTA and 1.0% SDS and heated overnight at 65 °C to reverse crosslinks. DNA was then purified by digesting with proteinase K followed by phenol chloroform extraction and ethanol precipitation.
Primers used were as follows:
Murine Tyr promoter: Forward: 5′-AGT CAT GTG CTT TGC AGA AGA T-3′ and Reverse: 5′-CAG CCA AGA ACA TTT TCT CCT T-3′.
Murine Tyrp1 promoter: Forward: 5′-GCA AAA TCT CTT CAG CGT CTC-3′ and Reverse: 5′-AGC CAG ATT CCT CAC ACT GG-3′.
Murine IgH enhancer: Forward: 5′-GCC GAT CAG AAC CAG AAC ACC-3′ and Reverse: 5′- TGG TGG GGC TGG ACA GAG TGT TTC-3′.
Chromatin accessibility (FAIRE analysis)
The formaldehyde-assisted isolation of regulatory elements (FAIRE) assay was performed as previously described . Briefly, cells were cross-linked with 1% formaldehyde for 6 min at room temperature, then nuclei were isolated and sonicated as described for ChIPs. Chromatin was extracted twice with phenol/chloroform, back extracted with TE and then extracted with chloroform. The aqueous phase was then heated at 65 °C overnight to reverse crosslinks. DNA was purified as described for ChIPs. Primers used were the same as for ChIP studies.
Statistical significance was calculated by the Student’s t test using Graphpad Prism.
Treatment of Melb-a melanoblasts with (+)JQ1 inhibits visible pigmentation and melanin synthesis
Melb-a cells are unpigmented mouse melanoblasts that can be induced to differentiate into pigmented melanocytes over the course of several days . To determine if BET proteins regulate the process of melanogenesis, Melb-a melanoblasts were induced to differentiate in the presence or absence of the active stereoisomer of the BET protein inhibitor (+)JQ1. As previously reported, in the absence of (+)JQ1, Melb-a cells became progressively pigmented when they were induced to differentiate and synthesize melanin . However, treatment with (+)JQ1 inhibited visible pigmentation (Fig. 1a) and melanin synthesis (Fig. 1b). We found that melanogenesis in Melb-a cells was similarly inhibited by the BET inhibitor, PFI-1 (Additional file 1: Fig. S1).
The inhibitory effects of (+)JQ1 on pigmentation and melanin synthesis are dose dependent and reversible
We tested the dose-dependent effects of (+)JQ1 on melanogenesis by exposing Melb-a cells to (+)JQ1 concentrations ranging from 0 to 2 μM. Increasing concentrations of (+)JQ1 had progressively inhibitory effects on visible pigmentation (Fig. 2a). When we quantified melanin synthesis by normalizing to cell counts, we found that 125 nM (+)JQ1 was the lowest dose to significantly inhibit melanin synthesis and that the maximum inhibitory effect was achieved with 500 nM (+)JQ1.
Several studies indicate that the effects of BET inhibitors on biological processes are reversible once the drug is removed [34, 35]. To determine the durability of BET inhibition on melanogenesis, we induced differentiation of Melb-a cells in the presence of an inactive enantiomer of JQ1, (−)JQ1 or the active isomer (+)JQ1. After 72 h, the cells that were induced in the presence of (−)JQ1 were visibly pigmented while those induced in the presence of (+)JQ1 were not (Fig. 2c, top). We then split (+)JQ1-treated cells into two plates, one of which was cultured in the presence of (+)JQ1 and the other in the presence of (−)JQ1. We found that pigmentation was restored 3 days after replacement of (+)JQ1 with inactive (−)JQ1, while cells that were maintained in media with (+)JQ1 remained unpigmented (Fig. 2c, bottom). Thus, the inhibitory effect of BET inhibition on melanogenesis lasted several days, but was ultimately reversible after drug removal.
We next investigated whether BET inhibitors can suppress melanogenesis in normal human epidermal melanocytes (NHEMs). Unlike Melb-a melanoblasts which are initially unpigmented and can be induced to synthesize melanin, NHEMs are pigmented and constitutively produce melanin. We found that treatment of NHEMs with (+)JQ1 resulted in time-dependent de-pigmentation (Additional file 2: Fig. S2). This suggests that BET proteins are required for maintaining pigmentation of differentiated melanocytes.
(+)JQ1 inhibits Melb-a proliferation
We found treatment with (+)JQ1 decreased Melb-a proliferation during differentiation (Fig. 3a). To determine if the effects on proliferation were due to cell death or changes in the cell cycle, we differentiated Melb-a cells in the presence of (−)JQ1 or (+)JQ1 and evaluated both apoptosis and cell cycle progression at several time-points. We found no significant difference in apoptosis between (−)JQ1 and (+)JQ1 treated cells at any time point during the differentiation process (Fig. 3b). However, there was a significant increase in the number of cells in the G1 phase when they were induced to differentiate in the presence of (+)JQ1 (Fig. 3c). These results suggest that BET inhibition during melanocyte differentiation suppresses proliferation primarily by promoting G1 arrest.
BET inhibition alters the expression of melanocyte and melanoma-specific gene expression
To determine the effects of (+)JQ1 on gene expression, we performed RNA-seq of Melb-a cells that were differentiated in the presence or absence of (+)JQ1. After 48 h, we detected widespread effects on gene expression, with 4069 significantly up-regulated genes and 4430 significantly down-regulated genes (Additional file 3: Fig. S3A). GO-seq analysis indicated that nucleoside and protein binding, RNA processing and modification, and metabolic processes were enriched (Additional file 3: Fig. S3B, C). Kegg analysis indicated that there was significant down-regulation of metabolic pathways and up-regulation of ribosome and spliceosome-related gene expression (Additional file 3: Fig. S3D, E). Interestingly, Kegg analysis indicated up-regulation of cell cycle components (Additional file 3: Fig. S3) including Rb1, Ink4b (p15) and Cdkn1b (p27). We also noted a decrease in cyclin D1 expression. A change in the expression of these cell cycle regulators might explain how (+)JQ1 promotes G1 arrest (Fig. 3c).
Additional pathway analysis focused on melanocyte and melanoma relevant pathways revealed major cellular networks of melanogenesis, cellular differentiation, epigenetic regulation, and transcriptional regulation to be significantly enriched (Fig. 4a). Specifically, transcriptional networks controlled by SOX family members were down-regulated as were pathways associated with melanoma resistance and disease progression.
To understand the molecular basis for the observed effects of (+)JQ1 on pigmentation and to delineate a role for BET proteins in the regulation of melanogenesis, we interrogated the effects of (+)JQ1 on the expression of 49 curated “color” genes . This list includes curated “color” genes that either directly regulate melanin synthesis or proteins involved in melanosome construction and transport (excluding transcriptional and systemic regulators). We analyzed the transcriptional effects of BET inhibition on “color” genes in Melb-a cells and compared this data to a previously reported transcriptome analysis of (+)JQ1-treated 501Mel and SK-Mel147 melanoma cells . BET inhibition suppressed the expression of nearly half the “color” genes in Melb-a cells, and a smaller number (6 genes) in 501MEL and (6 genes) in SK-MEL147 melanoma cells (Fig. 4b, left). However, some “color” genes were up-regulated in BET inhibited Melb-a cells, a number similar to that in 501MEL and in SK-MEL147 melanoma cells (Fig. 4b, right). Interestingly, there was greater overlap in gene regulation by BET inhibition between Melb-a and either 501MEL or SK-MEL147 melanoma cells than between 501MEL and SK-MEL147. RAB27A was the only down-regulated gene amongst all three cell lines. Thus, BET inhibition has both overlapping and distinct effects on the expression of intrinsic regulators of pigmentation, which are markers of terminal melanocyte differentiation. These results suggest that (+)JQ1 potently inhibits differentiation of melanoblasts and mildly affects the differentiation status of melanoma cells.
Interrogation of a larger list of pigment-related genes (379 eGenes)  (including both positive and negative regulators of pigmentation) showed a similar trend with the expression of almost 50% of eGenes being differentially regulated as a result of BET inhibition in Melb-a cells and 20–30% of eGenes in melanoma cells (Additional file 4: Table S1). We noticed differential regulation of several transcriptional regulators and chromatin remodeling enzymes that impact pigmentation in all three cell lines (Fig. 4c). IRF4, a transcriptional regulator associated with freckles, blue eyes, and brown hair, was down-regulated in both Melb-a and 501MEL cells . SNAI2, a transcriptional repressor involved in melanocyte migration and associated with pigmentary disturbances in some types of Waardenburg syndrome, was down-regulated in all three cell lines . Interestingly, expression of the master regulator of melanocyte differentiation, Mitf, was modestly (< 1.5 fold) but significantly up-regulated in Melb-a cells. However, expression of PIAS3, a repressor of MITF activity , was strongly up-regulated (> 1.5 fold) in all three cell lines, suggesting that BET inhibition might suppress MITF activity. This could explain why many MITF target genes were down-regulated while Mitf expression was only slightly up-regulated.
Many of the pigment genes affected by BET inhibition are MITF target genes (Fig. 4b, denoted by asterisks). MITF also regulates genes involved in melanocyte proliferation and survival  as well as mast cell  and osteoclast function . We found that BET inhibition suppressed expression of 36% of validated MITF target genes  in Melb-a cells, half as many as in melanoma cells. BET inhibition increased expression of approximately 22% of the MITF target genes in all three cell lines, with overlapping effects on expression of genes involved in proliferation, survival, mast cell, and osteoclast function (Fig. 4d). There was also a high degree of overlap between genes affected by BET inhibition (Additional file 5: Table S2) and genes that were differentially regulated by MITF in melanoma cells (50% in Melb-a cells and 40% in 501Mel and SK-MEL147 cells) .
BET inhibition decreases TYR and TYRP1 expression in Melb-a and normal human epidermal melanocytes
To further investigate the role of BET proteins in the regulation of MITF target genes, we assayed the effects of (+)JQ1 treatment on gene expression during a time course of melanocyte differentiation. We confirmed earlier reports of a time-dependent increase in Tyr, Tyrp1, and Trpm1 mRNA levels during melanoblast differentiation (Fig. 5a) . However, treatment of melanoblasts with (+)JQ1 inhibited the expression of these genes at all time-points (Fig. 5a) and prevented accumulation of TYR and TYRP1 enzymes, which are required for melanin synthesis (Fig. 5b). In contrast, MITF protein levels were not noticeably affected by treatment with (+)JQ1 (Fig. 5c).
Similarly, the depigmenting effects of (+)JQ1 on human melanocytes (Additional file 2: Fig. S2) corresponded with decreased protein levels of TYRP1, occurring day 1 after treatment and decreased expression of TYR and MITF occurring over a period of 6 days (Fig. 5c). Treatment with (+)JQ1 suppressed TYRP1 at the mRNA level (Fig. 5d), but did not have a significant effect on TYR nor MITF mRNA levels during the time period tested (data not shown). Thus, BET inhibition of NHEMs suppressed expression of genes required for melanin synthesis, but the kinetics and gene-specific effects were different from that observed in differentiating Melb-a melanoblasts.
BRD4 and BRD2 are required for expression of melanocyte-specific genes
To determine if any or all of the BET protein family members are required for expression of melanocyte-specific genes, we transfected Melb-a cells with control siRNA or siRNAs that target BRD2, BRD3, or BRD4. We found that down-regulation of BRD2 or BRD4, but not BRD3 decreased TYR protein levels (Fig. 6a). Since BRD4 has been implicated in several cellular differentiation processes [34, 46,47,48,49], we investigated its role in melanocyte differentiation more thoroughly. In order to control for potential off-target effects, we used a second siRNA to deplete BRD4 from Melb-a cells. There was a decrease in TYR and TYRP1 at the protein level (Fig. 6b) and at the mRNA level (Fig. 6c). Because down-regulation of BRD2 also affected TYR protein levels, we used a second siRNA to deplete BRD2 and investigate its role in melanocyte investigation. Like BRD4, we found that BRD2 depletion inhibited TYR and TYRP1 expression at both the protein (Fig. 6b) and at the mRNA levels (Fig. 6d).
BRD4 and BRD2 occupy MITF-binding sites and are associated with open chromatin structure at TYR and TYRP1 loci
Since TYR and TYRP1 are MITF target genes, we interrogated previously published ChIP-seq datasets to determine if BRD4 and MITF co-occupy MITF-binding sites at the TYR and TYRP1 loci in melanocytes and melanoma cells. BRD4 and BRD2 peaks overlapped MITF and histone H3 lysine 27 acetylation (H3K27ac) peaks, at the TYR promoter region and at an upstream enhancer region and promoter of TYRP1 (Fig. 7a, b).
We also interrogated other representative MITF target genes that were either down-regulated (Additional file 6: Fig. S4) or up-regulated (Additional file 7: Fig. S5) by (+)JQ1. While BRD4 and BRD2 peaks as well as H3K27ac peaks overlapped MITF-binding sites at the majority of these representative genes, in both normal human melanocytes and melanoma cells, there were some MITF sites such as within the DCT, PIAS3, SLC45A2, and IRF4 loci where BRD4 binding was observed in SK-MEL-5 melanoma cells but not normal human melanocytes (Additional file 7: Fig. S5). Differences in BET protein binding in different cell lines could explain some of the cell line-specific effects of (+)JQ1 on gene expression. Interestingly, PIAS3, which was highly up-regulated by (+)JQ1, had several BRD2 peaks which did not overlap the MITF-binding site nor H3K27ac-rich regions which could potentially have an inhibitory effect on transcription (Additional file 7: Fig. S5).
While it was possible to obtain overlapping BRD4, BRD2, and MITF peaks on representative genes in different melanocytes and melanoma cells, we could not find BRD4, BRD2, and MITF ChIP-seq data for the same melanocyte or melanoma cell line to reliably quantify the overlap. However, when we compared Encode data for the K257 leukemia cell line using the GSM2527370 data for MITF and the GSM2700494 data for BRD4, we found 6937 MITF peaks and 8509 BRD4 peaks, with 1204 overlapping between MITF and BRD4. Using the peakPermTest function, we determined a p-value of 0.099 indicating highly significant overlap. This suggests that MITF and BRD4 are also likely to interact in melanocytes and melanoma cells. However, additional ChIP-seq studies will be needed to determine the extent of the MITF/BRD4/BRD2 overlap on melanocyte-specific loci.
Consistent with the ChIP-seq data showing overlapping peaks in human melanocytes and melanoma cells, we found that both BRD4 and BRD2 were bound to the Tyr and Tyrp1 promoters in Melb-a cells and that (+)JQ1 disrupted binding (Fig. 8a, b). Furthermore, FAIRE analysis indicated that (+)JQ1 also decreased chromatin accessibility at the MITF-binding sites in the Tyr and Tyrp1 promoters in Melb-a cells and in the TYRP1 promoter in human melanocytes (Fig. 8c, d).
We next determined that the decrease in chromatin accessibility at the Tyr and Tyrp1 promoters in (+)JQ1 treated cells was associated with a reduction in MITF binding on the Tyr and Tyrp1 promoters in (+)JQ1 treated cells (Fig. 9a). Furthermore, depletion of either BRD4 or BRD2 also decreased MITF binding on these promoters (Fig. 9b), suggesting that these BET proteins facilitate MITF binding to the promoters of these melanocyte-specific genes.
BRD4 and BRD2 interact with MITF
To assay for physical interactions between BRD4, BRD2, and MITF, HEK293T cells were transfected with epitope-tagged (FLAG) MITF and co-immunoprecipitation was performed with either an antibody to BRD4 or to BRD2. Endogenous BRD4 co-immunoprecipitated with transfected MITF in HEK293T cells (Fig. 10a). We also detected interactions between endogenous BRD4 and MITF in Melb-a cells (Fig. 10b). BRD2 also co-immunoprecipitated with MITF in HEK293T cells (Fig. 10c) and in Melb-a cells (Fig. 10d).
Small molecule inhibitors of the bromodomain and extra-terminal (BET) class of epigenetic regulators have been extensively studied for their anticancer activities, quickly advancing to human clinical trials. These inhibitors selectively bind to bromodomains of the BET family of proteins and disrupt binding to chromatin, thereby altering gene expression profiles that impact cancer cells. However, these inhibitors also affect normal biological processes. Therefore, it is critical to delineate the role of each BET protein in non-cancerous cells in order to more effectively utilize these drugs and to minimize toxicities .
Bromodomain and extra-terminal proteins regulate transcriptional networks required for normal organismal development, stem cell maintenance, cellular differentiation, and immune function in vivo. Mice with homozygous disruption of either Brd2 or Brd4 die during embryogenesis [51, 52]. BRD4 is required for maintaining the self-renewal and pluripotency of stem cells [53, 54] and BRD2 for lineage specification . Inducible knockdown of BRD4 by RNAi interferes with normal hematopoiesis and leads to epidermal hyperplasia and stem cell depletion . BRD4 disruption in myeloid cells compromises innate immunity . Furthermore, BRD4 is required for muscle and brown adipose tissue development . BET proteins also promote erythroid, osteoblast, muscle, and adipocyte differentiation [34, 46, 48, 49]. However, the function of BET proteins in melanocyte differentiation had not been investigated prior to this study.
Our findings reveal that inhibiting BET proteins with (+)JQ1 in melanoblasts suppresses differentiation, resulting in a striking decrease in visible pigmentation and melanin synthesis. Importantly, (+)JQ1 treatment resulted in reduced expression of Tyr mRNA and protein in differentiating Melb-a melanoblasts. Tyr is the gene encoding the rate-limiting enzyme in melanin synthesis, the loss of which results in a complete failure to synthesize melanin and is manifested as oculocutaneous albinism in individuals with TYR mutations . In addition to the quantity of melanin, pigmentation is determined by the type of melanin, size, number, and structure of melanosomes . We found that BET inhibition suppressed the expression of additional pigmentation genes regulating these melanocyte functions. The mRNA levels of many transcriptional regulators that affect the level of pigmentation, including Sox10 , Ets1 , and Irf4 , were also reduced while the expression of Tfap2 was strongly increased in both Melb-a melanoblasts and SK-MEL-147 melanoma cells. TFAP2, which promotes melanocyte differentiation by activating a subset of MITF target genes, is frequently lost in melanoma due to silencing by CpG methylation [60, 61]. Interestingly, the expression of Pias3, a potent inhibitor of MITF activity  was strongly up-regulated and the expression of an extensive number of MITF target genes was decreased. These findings suggest that BET inhibition prevented melanocyte differentiation downstream of Mitf expression. This mechanism could involve the reduction in the expression of other positive transcriptional regulators of MITF target genes or inhibition of MITF activity through the increase in the expression of PIAS3. Additionally, we considered the possibility that BET proteins are MITF coactivators and that (+)JQ1 inhibits this function.
Since pan-BET inhibitors such as (+)JQ1 disrupt the activities of all BET family members, regardless of the intended target, it is critical to delineate the role of each BET protein . Therefore, we further explored the regulation of melanogenic genes by depleting single BET proteins in our Melb-a model. Our studies determined that both BRD4 and BRD2 are required for expression of two critical melanogenic genes, Tyr and Tyrp1, indicating a role for both BRD2 and BRD4 in melanocyte differentiation. Similarly, BRD2 and BRD4 are required for erythroid gene activation , but play opposing roles during adipogenesis [46, 47, 62].
Although the gene expression changes that resulted from treatment with (+)JQ1 obtained in our RNA-seq data suggests an extensive role for BET proteins in melanocyte differentiation, a limitation of this study is that we cannot determine which specific BET proteins mediate the transcriptional changes. Moreover, the three somatic BET proteins can have overlapping or opposing functions in the regulation of gene expression. Interestingly, transcriptional profiling of the epithelial mesenchymal transition showed that gene expression changes elicited by (+)JQ1 treatment most closely paralleled changes elicited by BRD2 depletion rather than by depletion of BRD3 or BRD4 which elicited opposite transcriptional changes as (+)JQ1 . The opposing actions of different BET proteins could explain why some genes are activated while others repressed in our experiments. Future studies will employ an unbiased approach to determine which of the BET proteins play a prominent role in the activation or repression of the melanocyte differentiation program. However, since depletion of either BRD2 or BRD4 compromised tyrosinase expression, which is absolutely required for melanin synthesis, our data suggest that these two proteins play an important role in pigmentation. Therefore, we proceeded to investigate the mechanism by which they regulate melanocyte-specific gene expression.
We found that both BRD4 and BRD2 physically interact with MITF and that occupancy of both proteins overlap MITF-binding sites on MITF target promoters. Furthermore, depletion of either BRD4 or BRD2 compromised MITF recruitment to melanocyte-specific promoters. Similarly, BRD2, BRD3, and BRD4 interact with master regulators of adipocyte and erythroid differentiation [47, 63]. BET proteins interact with acetylated GATA 1 to promote erythroid differentiation and depletion of BET proteins can also compromise GATA1 recruitment to erythroid genes. Acetylation of the RelA subunit of the NF-κβ has also been demonstrated to promote interactions with BET proteins . Although there are no published reports of MITF being acetylated, there are several lysine residues upstream of the basic helix loop helix domain that are predicted to be acetylation sites . These are in regions of MITF that modulate its localization to the nucleus [65, 66]. Moreover, MITF like MYOD in myocytes, interacts with the histone acetyl transferase, CBP, which can affect DNA binding and activity [67, 68]. It is tempting to speculate that acetylation of these helix loop helix lineage-specific transcription factors also facilitates binding to bromodomain proteins such as BRD4. However, neither BRD4 nor BRD2 were identified in an unbiased screen of MITF interacting proteins . This may indicate that BRD4/BRD2/MITF complexes are less abundant than those of other MITF-binding partners that were identified in the interactome or that the BRD4/BRD2/MITF interaction is mediated by other proteins, weaker, or more difficult to detect. Interestingly, the BRG1 component of the SWI/SNF complex which is present in the MITF interactome is also known to physically interact with BRD4 and BRD2 [8, 54]. Therefore, BRG1 or a BRG1-associated factor (BAF) could potentially mediate BRD4 and BRD2 interactions with MITF. However, additional studies will be required to determine the mechanisms by which BRD4 and BRD2 interact with MITF.
There was a less striking effect of BET inhibition on the expression of differentiation markers in melanoma cells when we compared a published data set  with our results in melanoblasts. While this could be due to the variations in experimental conditions, differences in differentiation status may also be a cause. In our experiments, Melb-a cells were induced to differentiate and express previously silent genes, while in melanoma cells, these genes were constitutively on or off. Our experiments on normal human melanocytes showed that BET inhibition can cause de-pigmentation by reducing expression of TYR and TYRP1 in a time-dependent fashion. Interestingly, we observed that inhibition of TYRP1 expression occurred at the mRNA and protein levels, but that inhibition of TYR occurred only at the protein level, suggesting that BET inhibition may affect TYR at a post-transcriptional level. Mechanistically, the inhibition of TYR could result from decreased TYRP1, which is known to regulate TYR protein stability . The lack of a transcriptional effect on TYR expression by BET inhibition of differentiated melanocytes suggests that there are molecular distinctions regarding the induction of “differentiation” in melanoblasts compared with “de-differentiation” of melanocytes. This could imply that for some genes, BET proteins are required for induction but not for maintenance of gene expression. Alternatively, complete “de-differentiation” of melanocytes by BET inhibition may take longer than the duration of this study. Likewise, the observation that BET inhibition of melanoma cells had either no effect on TYR expression in SK-MEL-147 or increased its expression in 501MEL cells could be due to the distinct differentiation status of each cell line. In melanoma cells, the altered signaling pathways in SK-MEL-147 cells, which have an NRAS mutation and 501Mel cells, which have a BRAF mutation, could also be a reason for the differential response to BET inhibitors. Pharmacologic inhibition of mutant BRAF in melanoma cells activates melanocyte-specific gene expression . Therefore, it is possible that signaling pathways can impact upon the interaction between MITF and BET proteins, possibly through regulation of MITF posttranslational modifications, such as acetylation.
This study demonstrates that BET proteins play a role in melanocyte differentiation. We found that chemical inhibition of BET proteins by (+)JQ1 decreases melanin synthesis and de-regulates expression of melanocyte-specific gene expression in both normal melanocytes and in melanoma. We have elucidated a mechanism by which the BET protein, BRD4, interacts with MITF to regulate expression of genes important for melanin synthesis. Both MITF and the melanocyte differentiation program have been shown to confer aggressiveness and promote resistance to therapeutics [24, 25]. The observation that BET inhibition can alter this process can impinge upon therapeutic approaches.
Availability of data and materials
The RNA-seq dataset has been deposited into the NCBI Sequence Read Archive (SRA) database (https://www.ncbi.nlm.nih.gov/sra/), under the Accession number PRJNA560240.
Wu SY, Chiang CM. The double bromodomain-containing chromatin adaptor Brd4 and transcriptional regulation. J Biol Chem. 2007;282(18):13141–5.
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(6):1068–1080.e5.
Andrieu GP, Denis GV. BET proteins exhibit transcriptional and functional opposition in the epithelial-to-mesenchymal transition. Mol Cancer Res. 2018;16(4):580–6.
Deeney JT, Belkina AC, Shirihai OS, Corkey BE, Denis GV. BET bromodomain proteins Brd2, Brd3 and Brd4 selectively regulate metabolic pathways in the pancreatic beta-cell. PLoS ONE. 2016;11(3):e0151329.
Jiang YW, Veschambre P, Erdjument-Bromage H, Tempst P, Conaway JW, Conaway RC, et al. Mammalian mediator of transcriptional regulation and its possible role as an end-point of signal transduction pathways. Proc Natl Acad Sci USA. 1998;95(15):8538–43.
Patel MC, Debrosse M, Smith M, Dey A, Huynh W, Sarai N, et al. BRD4 coordinates recruitment of pause release factor P-TEFb and the pausing complex NELF/DSIF to regulate transcription elongation of interferon-stimulated genes. Mol Cell Biol. 2013;33(12):2497–507.
Wu T, Kamikawa YF, Donohoe ME. Brd4’s bromodomains mediate histone H3 acetylation and chromatin remodeling in pluripotent cells through P300 and Brg1. Cell Rep. 2018;25(7):1756–71.
Denis GV, McComb ME, Faller DV, Sinha A, Romesser PB, Costello CE. Identification of transcription complexes that contain the double bromodomain protein Brd2 and chromatin remodeling machines. J Proteome Res. 2006;5(3):502–11.
Wang R, Li Q, Helfer CM, Jiao J, You J. Bromodomain protein Brd4 associated with acetylated chromatin is important for maintenance of higher-order chromatin structure. J Biol Chem. 2012;287(14):10738–52.
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.
Hsu SC, Gilgenast TG, Bartman CR, Edwards CR, Stonestrom AJ, Huang P, et al. The BET protein BRD2 cooperates with CTCF to enforce transcriptional and architectural boundaries. Mol Cell. 2017;66(1):102–116.e7.
French CA, Miyoshi I, Kubonishi I, Grier HE, Perez-Atayde AR, Fletcher JA. BRD4-NUT fusion oncogene: a novel mechanism in aggressive carcinoma. Cancer Res. 2003;63(2):304–7.
French CA, Ramirez CL, Kolmakova J, Hickman TT, Cameron MJ, Thyne ME, et al. BRD-NUT oncoproteins: a family of closely related nuclear proteins that block epithelial differentiation and maintain the growth of carcinoma cells. Oncogene. 2008;27(15):2237–42.
Zanconato F, Battilana G, Forcato M, Filippi L, Azzolin L, Manfrin A, et al. Transcriptional addiction in cancer cells is mediated by YAP/TAZ through BRD4. Nat Med. 2018;24(10):1599–610.
Loven J, Hoke HA, Lin CY, Lau A, Orlando DA, Vakoc CR, et al. Selective inhibition of tumor oncogenes by disruption of super-enhancers. Cell. 2013;153(2):320–34.
Filippakopoulos P, Qi J, Picaud S, Shen Y, Smith WB, Fedorov O, et al. Selective inhibition of BET bromodomains. Nature. 2010;468(7327):1067–73.
Roe JS, Mercan F, Rivera K, Pappin DJ, Vakoc CR. BET bromodomain inhibition suppresses the function of hematopoietic transcription factors in acute myeloid leukemia. Mol Cell. 2015;58(6):1028–39.
Xu Y, Vakoc CR. Targeting cancer cells with BET bromodomain inhibitors. Cold Spring Harb Perspect Med. 2017;7(7):a026674.
Stathis A, Bertoni F. BET proteins as targets for anticancer treatment. Cancer Discov. 2018;8(1):24–36.
Vandamme N, Berx G. From neural crest cells to melanocytes: cellular plasticity during development and beyond. Cell Mol Life Sci. 2019;76(10):1919–34.
Costin GE, Hearing VJ. Human skin pigmentation: melanocytes modulate skin color in response to stress. FASEB J. 2007;21(4):976–94.
Levy C, Khaled M, Fisher DE. MITF: master regulator of melanocyte development and melanoma oncogene. Trends Mol Med. 2006;12(9):406–14.
Sarkar D, Leung EY, Baguley BC, Finlay GJ, Askarian-Amiri ME. Epigenetic regulation in human melanoma: past and future. Epigenetics. 2015;10(2):103–21.
Gupta PB, Kuperwasser C, Brunet JP, Ramaswamy S, Kuo WL, Gray JW, et al. The melanocyte differentiation program predisposes to metastasis after neoplastic transformation. Nat Genet. 2005;37(10):1047–54.
Johannessen CM, Johnson LA, Piccioni F, Townes A, Frederick DT, Donahue MK, et al. A melanocyte lineage program confers resistance to MAP kinase pathway inhibition. Nature. 2013;504(7478):138–42.
Marathe HG, Watkins-Chow DE, Weider M, Hoffmann A, Mehta G, Trivedi A, et al. BRG1 interacts with SOX10 to establish the melanocyte lineage and to promote differentiation. Nucleic Acids Res. 2017;45(11):6442–58.
Saladi SV, Wong PG, Trivedi AR, Marathe HG, Keenen B, Aras S, et al. BRG1 promotes survival of UV-irradiated melanoma cells by cooperating with MITF to activate the melanoma inhibitor of apoptosis gene. Pigment Cell Melanoma Res. 2013;26(3):377–91.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database issue):D480–4.
Keenen B, Qi H, Saladi SV, Yeung M, de la Serna IL. Heterogeneous SWI/SNF chromatin remodeling complexes promote expression of microphthalmia-associated transcription factor target genes in melanoma. Oncogene. 2010;29(1):81–92.
Sviderskaya EV, Wakeling WF, Bennett DC. A cloned, immortal line of murine melanoblasts inducible to differentiate to melanocytes. Development. 1995;121(5):1547–57.
Roberts TC, Etxaniz U, Dall’Agnese A, Wu SY, Chiang CM, Brennan PE, et al. BRD3 and BRD4 BET bromodomain proteins differentially regulate skeletal myogenesis. Sci Rep. 2017;7(1):6153.
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.
Montoliu L, Detting WS, Bennett DC. Color genes. European Society for Pigment Cell Research. 2011. http://www.espcr.org/micemut. Accessed 1 Jan 2020.
Fontanals-Cirera B, Hasson D, Vardabasso C, Di Micco R, Agrawal P, Chowdhury A, et al. Harnessing BET inhibitor sensitivity reveals AMIGO2 as a melanoma survival gene. Mol Cell. 2017;68(4):731–744.e9.
Zhang T, Choi J, Kovacs MA, Shi J, Xu M, Program NCS, et al. Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes. Genome Res. 2018;28(11):1621–35.
Praetorius C, Grill C, Stacey SN, Metcalf AM, Gorkin DU, Robinson KC, et al. A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway. Cell. 2013;155(5):1022–33.
Sanchez-Martin M, Rodriguez-Garcia A, Perez-Losada J, Sagrera A, Read AP, Sanchez-Garcia I. SLUG (SNAI2) deletions in patients with Waardenburg disease. Hum Mol Genet. 2002;11(25):3231–6.
Levy C, Nechushtan H, Razin E. A new role for the STAT3 inhibitor, PIAS3: a repressor of microphthalmia transcription factor. J Biol Chem. 2002;277(3):1962–6.
Strub T, Giuliano S, Ye T, Bonet C, Keime C, Kobi D, et al. Essential role of microphthalmia transcription factor for DNA replication, mitosis and genomic stability in melanoma. Oncogene. 2011;30(20):2319–32.
Ito A, Jippo T, Wakayama T, Morii E, Koma Y, Onda H, et al. SgIGSF: a new mast-cell adhesion molecule used for attachment to fibroblasts and transcriptionally regulated by MITF. Blood. 2003;101(7):2601–8.
Meadows NA, Sharma SM, Faulkner GJ, Ostrowski MC, Hume DA, Cassady AI. The expression of Clcn7 and Ostm1 in osteoclasts is coregulated by microphthalmia transcription factor. J Biol Chem. 2007;282(3):1891–904.
Hoek KS, Schlegel NC, Eichhoff OM, Widmer DS, Praetorius C, Einarsson SO, et al. Novel MITF targets identified using a two-step DNA microarray strategy. Pigment Cell Melanoma Res. 2008;21(6):665–76.
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(9):2144–9.
Lee JE, Park YK, Park S, Jang Y, Waring N, Dey A, et al. Brd4 binds to active enhancers to control cell identity gene induction in adipogenesis and myogenesis. Nat Commun. 2017;8(1):2217.
Stonestrom AJ, Hsu SC, Jahn KS, Huang P, Keller CA, Giardine BM, et al. Functions of BET proteins in erythroid gene expression. Blood. 2015;125(18):2825–34.
Najafova Z, Tirado-Magallanes R, Subramaniam M, Hossan T, Schmidt G, Nagarajan S, et al. BRD4 localization to lineage-specific enhancers is associated with a distinct transcription factor repertoire. Nucleic Acids Res. 2017;45(1):127–41.
Andrieu G, Belkina AC, Denis GV. Clinical trials for BET inhibitors run ahead of the science. Drug Discov Today Technol. 2016;19:45–50.
Houzelstein D, Bullock SL, Lynch DE, Grigorieva EF, Wilson VA, Beddington RS. Growth and early postimplantation defects in mice deficient for the bromodomain-containing protein Brd4. Mol Cell Biol. 2002;22(11):3794–802.
Shang E, Wang X, Wen D, Greenberg DA, Wolgemuth DJ. Double bromodomain-containing gene Brd2 is essential for embryonic development in mouse. Dev Dyn. 2009;238(4):908–17.
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.
Liu W, Stein P, Cheng X, Yang W, Shao NY, Morrisey EE, et al. BRD4 regulates Nanog expression in mouse embryonic stem cells and preimplantation embryos. Cell Death Differ. 2014;21(12):1950–60.
Fernandez-Alonso R, Davidson L, Hukelmann J, Zengerle M, Prescott AR, Lamond A, et al. Brd4-Brd2 isoform switching coordinates pluripotent exit and Smad2-dependent lineage specification. EMBO Rep. 2017;18(7):1108–22.
Bao Y, Wu X, Chen J, Hu X, Zeng F, Cheng J, et al. Brd4 modulates the innate immune response through Mnk2-eIF4E pathway-dependent translational control of IkappaBalpha. Proc Natl Acad Sci USA. 2017;114(20):E3993–4001.
Tomita Y, Takeda A, Okinaga S, Tagami H, Shibahara S. Human oculocutaneous albinism caused by single base insertion in the tyrosinase gene. Biochem Biophys Res Commun. 1989;164(3):990–6.
Mollaaghababa R, Pavan WJ. The importance of having your SOX on: role of SOX10 in the development of neural crest-derived melanocytes and glia. Oncogene. 2003;22(20):3024–34.
Saldana-Caboverde A, Perera EM, Watkins-Chow DE, Hansen NF, Vemulapalli M, Mullikin JC, et al. The transcription factors Ets1 and Sox10 interact during murine melanocyte development. Dev Biol. 2015;407(2):300–12.
Seberg HE, Van Otterloo E, Loftus SK, Liu H, Bonde G, Sompallae R, et al. TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF. PLoS Genet. 2017;13(3):e1006636.
Hallberg AR, Vorrink SU, Hudachek DR, Cramer-Morales K, Milhem MM, Cornell RA, et al. Aberrant CpG methylation of the TFAP2A gene constitutes a mechanism for loss of TFAP2A expression in human metastatic melanoma. Epigenetics. 2014;9(12):1641–7.
Zang K, Wang J, Dong M, Sun R, Wang Y, Huang Y, et al. Brd2 inhibits adipogenesis via the ERK1/2 signaling pathway in 3T3-L1 adipocytes. PLoS ONE. 2013;8(10):e78536.
Lamonica JM, Deng W, Kadauke S, Campbell AE, Gamsjaeger R, Wang H, et al. Bromodomain protein Brd3 associates with acetylated GATA1 to promote its chromatin occupancy at erythroid target genes. Proc Natl Acad Sci USA. 2011;108(22):E159–68.
Huang B, Yang XD, Zhou MM, Ozato K, Chen LF. Brd4 coactivates transcriptional activation of NF-kappaB via specific binding to acetylated RelA. Mol Cell Biol. 2009;29(5):1375–87.
Hou T, Zheng G, Zhang P, Jia J, Li J, Xie L, et al. LAceP: lysine acetylation site prediction using logistic regression classifiers. PLoS ONE. 2014;9(2):e89575.
Fock V, Gudmundsson SR, Gunnlaugsson HO, Stefansson JA, Ionasz V, Schepsky A, et al. Subcellular localization and stability of MITF are modulated by the bHLH-Zip domain. Pigment Cell Melanoma Res. 2018;32(1):41–54.
Sato S, Roberts K, Gambino G, Cook A, Kouzarides T, Goding CR. CBP/p300 as a co-factor for the Microphthalmia transcription factor. Oncogene. 1997;14(25):3083–92.
Polesskaya A, Harel-Bellan A. Acetylation of MyoD by p300 requires more than its histone acetyltransferase domain. J Biol Chem. 2001;276(48):44502–3.
Laurette P, Strub T, Koludrovic D, Keime C, Le Gras S, Seberg H, et al. Transcription factor MITF and remodeller BRG1 define chromatin organisation at regulatory elements in melanoma cells. eLife. 2015;4:e06857.
Kobayashi T, Imokawa G, Bennett DC, Hearing VJ. Tyrosinase stabilization by Tyrp1 (the brown locus protein). J Biol Chem. 1998;273(48):31801–5.
We thank Megan Dreier for help with immunoblots.
This research was supported by NIAMS R01AR059379, by a University of Toledo Biomedical Innovation grant, and a grant from the Melanoma Research Foundation to ILD. VS acknowledges support in part from VA Merit Award 1 I01BX002623 and NIH/NIAMS Skin Disease Research Center grant P30AR066524. FVF is grateful for the support by grants NCI CA154887, CRN-17-427258, and the Science Alliance on Precision Medicine and Cancer Prevention by the German Federal Foreign Office, implemented by the Goethe-Institute, Washington, DC, USA, and supported by the Federation of German Industries (BDI), Berlin, Germany.
Ethics approval and consent to participate
Consent for publication
ILD and AT hold a patent (Patent Number: 10,064,798 B1) on the use of BET bromodomain inhibitors for modulating pigmentation.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Fig. S1.
Melb-a cells were differentiated for the indicated number of days in vehicle (VC) or 500 nM PFI-1. Cells were pelleted and photographed.
Additional file 2: Fig. S2.
Neonatal human epidermal melanocytes (NHEMs) were cultured in the presence of vehicle (VC) or 500nM (+)JQ1 for the indicated number of days. Cells were pelleted and photographed.
Additional file 3: Fig. S3.
Melb-a cells were differentiated for 48 h in the presence of vehicle or 500nM +JQ1. RNA was isolated and subjected to RNA-seq. Differential gene expression ((+)JQ1 vs. Vehicle) was determined based on the average of three biological replicates (adj. p < 0.05). A. Volcano plot showing down-regulated gene expression (green) and up-regulated gene expression (red). B. Gene ontology analysis of up-regulated biological processes, cellular components, molecular function. C. Gene ontology analysis of down-regulated biological processes, cellular components, molecular function. D. Kegg Pathway analysis of up-regulated pathways. E. Kegg Pathway analysis of down-regulated pathways.
Additional file 4: Table S1.
The effect of BET inhibition on E-gene expression. E-genes that are significantly (p < 0.05) down-regulated (yellow) or up-regulated (blue) by BET inhibition of Melb-a, 501MEL, or SK-MEL147 cells. E-genes that are affected by BET inhibition in all three cell lines are shown in bold type.
Additional file 5: Table S2.
The effect of BET inhibition on MITF-target gene expression. MITF target genes that are significantly (p < 0.05) down-regulated (yellow) or up-regulated (blue) by BET inhibition of Melb-a, 501MEL, or SK-MEL147 cells. MITF target genes that are affected by BET inhibition in all three cell lines are shown in bold type.
Additional file 6: Fig. S4.
Publicly available ChIP-seq data in normal human melanocytes, SK-MEL-5 melanoma cells (GSM1968282), 501 melanoma cells (GSM1517751), and neonatal human epidermal melanocytes (NHEM) (GSM2842798) with the indicated antibodies at the MITF target genes that are down-regulated by (+)JQ1 treatment.
Additional file 7: Fig. S5.
Publicly available ChIP-seq data in normal human melanocytes, SK-MEL-5 melanoma cells (GSM1968282), 501 melanoma cells (GSM1517751), and neonatal human epidermal melanocytes (NHEM) (GSM2842798) with the indicated antibodies at the MITF target genes that are up-regulated by (+)JQ1 treatment.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Trivedi, A., Mehrotra, A., Baum, C.E. et al. Bromodomain and extra-terminal domain (BET) proteins regulate melanocyte differentiation. Epigenetics & Chromatin 13, 14 (2020). https://doi.org/10.1186/s13072-020-00333-z
- Bromodomain and extra-terminal domain
- Systems biology
- Melanocyte differentiation
- Transcriptional networks