Selective Recognition of H3.1K36 Dimethylation/H4K16 Acetylation Facilitates the Regulation of All-trans-retinoic Acid (ATRA)-responsive Genes by Putative Chromatin Reader ZMYND8*

ZMYND8 (zinc finger MYND (Myeloid, Nervy and DEAF-1)-type containing 8), a newly identified component of the transcriptional coregulator network, was found to interact with the Nucleosome Remodeling and Deacetylase (NuRD) complex. Previous reports have shown that ZMYND8 is instrumental in recruiting the NuRD complex to damaged chromatin for repressing transcription and promoting double strand break repair by homologous recombination. However, the mode of transcription regulation by ZMYND8 has remained elusive. Here, we report that through its specific key residues present in its conserved chromatin-binding modules, ZMYND8 interacts with the selective epigenetic marks H3.1K36Me2/H4K16Ac. Furthermore, ZMYND8 shows a clear preference for canonical histone H3.1 over variant H3.3. Interestingly, ZMYND8 was found to be recruited to several developmental genes, including the all-trans-retinoic acid (ATRA)-responsive ones, through its modified histone-binding ability. Being itself inducible by ATRA, this zinc finger transcription factor is involved in modulating other ATRA-inducible genes. We found that ZMYND8 interacts with transcription initiation-competent RNA polymerase II phosphorylated at Ser-5 in a DNA template-dependent manner and can alter the global gene transcription. Overall, our study identifies that ZMYND8 has CHD4-independent functions in regulating gene expression through its modified histone-binding ability.

ZMYND8 (zinc finger MYND (Myeloid, Nervy and DEAF-1)-type containing 8) is a putative chromatin reader/effector harboring a PWWP domain, a bromodomain, and a PHD type zinc finger. It associates with the CHD4 protein of NuRD chromatin remodeling complex implicated in gene transcription, cell cycle progression, and genome integrity (1,2). Reader proteins specifically recognize histone post-translational modifications (PTMs) 5 and translate such recognition(s) into meaningful biological outcomes by virtue of either their intrinsic activities or those of their interacting partners (3). These readers may interpret histone PTMs via "monovalent" or "multivalent" recognition (4). In monovalent recognition, a chromatin reader recognizes one histone PTM. In contrast, in multivalent recognition, readers recognize multiple histone modifications intra/inter-nucleosomally. An accessible surface is provided by these readers (such as a cavity or surface groove), which accommodates the modified histone residues, and determines the modification (e.g. acetylation versus methylation) or state specificity (such as monoversus trimethylation of lysine) (3).
ZMYND8 is a well known component of the transcription coregulator complex and is associated with several demethylase machinery components, including KDM5A, KDM5C, or LSD1 (5,6). Through its ability to interact with Xenopus RCoR2, ZMYND8 plays a significant role in embryonic neural differentiation (7). Apart from this, ZMYND8 is also involved in T-cell lymphoma and breast and cervical cancer (8 -10). Another interesting feature is that ZMYND8 is significantly involved in transcription activation (5). ATRA, a vitamin A metabolite, is a well known inducer of transcription of several genes. It modulates RA receptor (retinoic acid receptor/retinoid X receptor) dimerization at the retinoic acid-response element (RARE) (11). We show for the first time that ZMYND8 is ATRAresponsive and a substantial number of genes regulated by ZMYND8 harbor the RARE sequence. Interestingly, we found that ZMYND8 is involved in regulating transcription initiation through its interaction with the RNA polymerase II complex in a DNA-mediated manner. Furthermore, ZMYND8 can regulate global gene expression in a CHD4-independent manner. In an attempt to understand its mode of recruitment to chromatin, we found that it has a selective interaction with H3.1K36Me2/ Chromatosome Isolation-Chromatosome was prepared from freshly harvested HeLa cells as described elsewhere (13). Briefly, nuclear pellet of HeLa cells was digested with micrococcal nuclease (Sigma 0.2 units/l) and extracted with TE buffer for 1 h. Concentrated chromatin was separated by sucrose gradient (5-40%) ultracentrifugation using Sorvall WXUltra100 (Thermo Scientific) in an AH650 rotor for 16 h at 35,000 rpm. Chromatin fractions were analyzed in agarose gel; chromatosomal fractions were pooled and dialyzed in TE buffer for further analysis.
CelluSpots Peptide Array Analysis-MODified Histone Peptide Array (CelluSpots array) was procured from Active Motif (catalogue no. 13005). CelluSpots array analysis was done as per standard protocol (14,15). Briefly, the array was blocked with the Blocking Buffer overnight at 4°C, washed three times with Wash Buffer, and incubated with GST-fused protein at a concentration of 1 M. The array was probed with anti-GST antibody.
Peptide pulldown assays were done as per standard protocol (16). Briefly, biotinylated histone peptides were incubated with GST-tagged protein in IP Buffer (50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 0.05% Nonidet P-40, 1 mM DTT) at 4°C, overnight. The reactions were subjected to Streptavidin beads (GE Healthcare) binding for 2 h at 4°C and washed with IP Buffer, and the complex was pulled down and analyzed by Western blotting.
Protein-Protein Interaction Studies-The core histone and chromatosome interactions with GST and GST fusion proteins were performed as per standard protocol (5,14). Briefly, proteins were incubated in an equal molar ratio in binding buffer (50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 0.05% Nonidet P-40, 1 mM DTT), pulled down by pre-blocked glutathione-Sepharose beads, and washed in binding buffer, and the protein complex was analyzed by Western blotting.
Fluorescence Spectroscopy and Analysis of the Binding Data-The fluorescence measurements were performed in an LS 55 luminescence spectrometer (PerkinElmer Life Sciences) at 25°C. In the case of H3K36Me3, H3K36Me2 , H3K36Me1, and H3K36 unmodified peptides, tryptophan fluorescence measurements were carried out on samples containing 0.45 M ZMYND8 (His-PHD/Bromo/PWWP) and dissolved in 20 mM HEPES (pH 6.5), 200 mM NaCl, 300 mM ammonium sulfate, 1 mM ␤-mercaptoethanol, while progressively increasing the peptide concentration in each case until there was no further change in the emission intensity with an increase in the peptide concentration. Samples were excited at 295 nm, and emission was recorded at 320 nm with a scan speed of 200 nm/s, and averaged over three scans.
In case of FAM-conjugated H4K16-acetylated and H4K16unmodified peptides, fluorescence measurements were carried out using 0.6 M of the peptides dissolved in 20 mM Tris-HCl (pH 6.5), 150 mM NaCl, 1 mM ␤-mercaptoethanol at 25°C with increasing addition of ZMYND8 (His-PBP). Samples were excited at 494 nm, and emission was recorded at 520 nm with scan speed 200 nm/s and averaged over three scans.
All experimental points for binding isotherms were fitted by least squares analysis as shown in Equation 1.
where ⌬F is the change in fluorescence emission intensity at 340 nm (excitation of tryptophan residue(s) at 295 nm) for each point of the titration curve. ⌬F max is the same parameter when the ligand (ZMYND8 His-PBP) is totally bound to histone. C p is the concentration of histone, and C 0 is the initial concentration of the ligand. A double-reciprocal plot was used for determination of ⌬F max using Equation 2, In the case of titration of ZMYND8 (His-PBP) with H4K16 unmodified and H4K16Ac histone peptides, the data analysis was done following the same theory as reported above, and the non-linear curve fitting was done following Equations 1 and 2. The only difference being that in this case, C 0 is the initial concentration of the FAM-conjugated histone H4K16 unmodified and C p is the concentration of the ligand ZMYND8 His-PBP.
Molecular Modeling Studies-Bromo-PWWP module of ZMYND8 was generated by MODELLER (18) package using chain A of 4N4G as a reference structure.
Coimmunofluorescence and Confocal Microscopy-Coimmunofluorescence staining was performed as per standard protocol (19). Briefly, the cells were fixed with 4% paraformaldehyde, permeabilized with 1% Triton X-100, and blocked with 3% BSA. Cells were incubated with anti-FLAG antibody (Sigma) and either anti-H4K16Ac (Active Motif) or H3K36Me2 (Abcam) antibody for 1 h. Following washes with PBST, the cells were incubated with Alexa Fluor 488-and Alexa Fluor 594-conjugated secondary antibodies for 1 h at room temperature. The coverslips were washed with PBST and mounted after DAPI staining. OLYMPUS FLUOVIEW FV1000 confocal scanning microscope was used for imaging.
Co-IP-Cross-linked cells were subjected to co-IP as delineated elsewhere (20). In brief, after cross-linking the cells were lysed with 50 mM HEPES (pH 7.5), 150 mM NaCl, 1.5 mM MgCl 2 , 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 0.5% sodium deoxycholate, 5% glycerol, 1 mM DTT, 2 mM PMSF, 1 mM trichostatin A, 5 mM sodium butyrate, and complete protease inhibitor mixture and incubated on ice for 1 h followed by centrifugation at 13,000 rpm for 10 min at 4°C. After pre-clearing the lysates, immunoprecipitation was done with the antibodies, followed by washes with the same buffer, and the complexes were analyzed by Western blotting.
Chromatin Immunoprecipitation (ChIP)-ChIP assays were performed as per standard protocol (22). Briefly, cross-linking was done with 1% formaldehyde and then the reaction was stopped by 0.125 M glycine, followed by cell lysis in Cell Lysis buffer (5 mM PIPES (pH 8.0), 85 mM KCl, 0.5% Nonidet P-40 (with fresh protease inhibitor)). The isolated nuclei were resuspended in nuclei lysis buffer (50 mM Tris-HCl (pH 8.0), 10 mM EDTA, 1% SDS (with fresh protease inhibitor)). After sonication of the chromatin, pre-clearing was done. Subsequently, chromatin was immunoprecipitated with anti-ZMYND8 (Sigma Prestige Group), anti-H4K16Ac (Active Motif), anti-H3K36Me2 (Abcam), anti-H4 (Abcam), and anti-H3 (Abcam) antibodies overnight at 4°C. Dynabeads were blocked with 0.3 mg/ml salmon sperm DNA and set for interaction. Beads were washed with RIPA buffer, high salt buffer, LiCl buffer, and TE buffer, consecutively. Following RNase A and proteinase K treatment, the beads were kept for de-cross-linking at 65°C. Phenol/chloroform extraction followed by ethanol precipitation was performed. The DNA pellet was dissolved in H 2 O and then qPCR analysis was done using gene-specific primers.
ChIP assay in SH-SY5Y cells was performed after treating with 10 M ATRA or DMSO for 4 days. The primers used in ChIP assays were designed inside the gene body or promoter elements using the UCSC Genome Browser. The detailed primer information is available on request.
Western Blot Analysis-Whole cell extracts were prepared with Laemmli Buffer (4% SDS, 20% glycerol, and 120 mM Tris-HCl (pH 6.8)) and sonicated followed by boiling at 100°C. The samples were electrophoresed on 7.5, 11, or 15% SDS-PAGE and transferred on nitrocellulose membrane followed by blocking with nonfat skimmed milk and probed for specific antibodies. Primary antibodies used are listed in Table 1.
Reporter Construct and Luciferase Assay-The 527-bp upstream promoter region of mouse Zmynd8 gene was amplified by PCR using mouse genomic DNA as template and then cloned into pGL3 basic vector (Promega) at KpnI/HindIII site. Similarly, the 873-bp upstream region of human ZMYND8 gene was cloned into pGL3-basic vector at KpnI/SacI site. The primer sequences used to clone the promoter regions are mentioned in Table 2 where the restriction enzyme sites are under-lined. All constructs were sequenced to confirm their identity. The mutation within the RARE in both human and mouse promoters was done by using QuikChange II XL site-directed mutagenesis kit (Agilent).
For reporter assay, 5 ϫ 10 4 cells were seeded on 12-well culture plates. After 24 h, the pGL3-ZMYND8 construct (0.4 g) was transiently transfected with Lipofectamine 2000 (Invitrogen). After 4 h of transfection, the medium was replaced with a fresh incomplete one supplemented with either RA or DMSO for the next 16 h. Each transfection was normalized with pRL-CMV vector (0.04 g). The following day, cells were harvested, and firefly/Renilla luciferase activity was determined using the manufacturer's protocol. Each transfection was performed in triplicate, and the experiments were repeated three times.
Quantitative Real Time PCR (qRT-PCR)-Total RNA was extracted using the RNeasy mini kit (Qiagen). Complementary DNA was synthesized from 2 g of total RNA using RevertAid fast strand cDNA synthesis kit (Thermo Scientific) according to the manufacturer's protocol followed by qRT-PCR using ABI-SYBR GREEN mix (ABI). qRT-PCR was performed using 7500 real time PCR machine. Primers used are listed in Table 2.
Microarray Analysis-Total cellular RNA was extracted, and RNA quality was determined using a Nanodrop spectrophotometer (NanoDrop Technologies, Wilmington, DE) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). Samples with RNA Integrity number of greater than 7.0 were taken for Affymetrix Human Whole Genome Gene Chip (Human HT Prime View Gene Chip)-based microarray hybridization as per the manufacturer's recommended protocol (Affymetrix, Santa Clara, CA). Raw data were RMA (Robust Multiarray Average)-normalized, and baseline transformation was done to median of all samples using Gene Spring GX 12.5 software (Agilent Technologies).
Statistical Analysis of Differentially Expressed Genes from Microarray-Differentially expressed probe sets (genes) upon ZMYND8 siRNA treatment in comparison with control siRNA-treated cells were identified by applying a fold change threshold of absolute fold change greater than or equal to 1.1, and a statistically significant t test p value threshold was adjusted for the false discovery rate of less than or equal to 0.05. The same was visualized using a volcano plot to understand the differential expression pattern. Unsupervised hierarchical clustering of differentially expressed genes upon treatment in comparison with untreated were done using Euclidian algorithm with Centroid linkage rule to identify gene clusters whose expression levels are significantly reproduced across the replicates.
Biological Pathways and Gene Ontology Enrichment Analysis of Microarray Data-Differentially expressed gene list was subjected to biological significance analysis by GOElite tool. 21,887 protein coding genes were used as the background, and differentially expressed gene lists were used as query. The database of Gene Ontology categories, Wikipathways, KEGG Pathways, Pathway Commons, Pheno Ontology, Diseases, Protein Domains, transcription factor targets, and tissue expression was configured for significance analysis. Each query list was subjected to over-representation analysis against each of the above databases. A Z score and permutation or Fisher's Exact Test p value were calculated to assess over-representation of enriched biological categories.

Results
ZMYND8 Is a Histone H3/H4-interacting Protein-The human transcription factor ZMYND8 is a component of the transcription coregulator network having a wide range of interacting partners, including chromatin remodeling complexes and histone demethylase/deacetylase enzymes (5). It has been reported as a novel DNA damage response factor that recruits NuRD complex to damaged chromatin (23). Its involvement in chromatin template-dependent processes, viz. transcription and repair, prompted us to speculate that ZMYND8 could be a chromatin-associated protein. Indeed, we found that ZMYND8 had a tight association to chromatin by fractionation experiment, where even 0.5% Nonidet P-40 treatment could not release the protein to the non-chromatin fraction in HEK293 cells (Fig. 1A). Histone H3 and tubulin were used as nuclear and cytosolic markers, respectively (Fig. 1A). Next, we investigated the ability of ZMYND8 to interact with core histones in HEK293 cells. By co-IP with anti-ZMYND8 antibody, we observed a strong interaction with histone H3 and H4 ( Fig. 1B) but not with histone H2A and H2B. However, the ability of ZMYND8 to interact with core histones could also be an indirect effect of preferential binding of ZMYND8 to other chromatin protein(s). To elucidate that, we performed an in vitro interaction of the PBP module of ZMYND8 with recombinant histones and core histones. Although ZMYND8 did not interact with individual histone H3 and H4 (Fig. 1C, panel I and panel II), it showed a strong interaction with histone H3 and H4 of core histone (Fig. 1D). These results indicate that the interaction with core histones is direct, although there is a requirement of a histone post-translational modification. Similar interaction of the PBP module of ZMYND8 with purified chromatosome from HeLa cells also showed its association to histone H3 and H4 (Fig. 1E). These results clearly indicate that ZMYND8 is a core histone H3/H4-interacting protein having a tight association with chromatin.
The selective interaction of ZMYND8 with H3.1K36-methylated and H4K16-acetylated peptides was further quantified by fluorescence spectroscopy. The PBP module was titrated with H3.1K36Me0/Me1/Me2/Me3 peptides, and quenching of the fluorescence intensity of the tryptophan residue at the peptide binding pocket of the protein was measured. No appreciable binding was seen with H3.1K36Me0 peptide. H3.1K36Me2 peptide (K d ϳ0.23 M) showed 2-fold higher affinity compared with H3.1K36Me3 (K d ϳ0.45 M) (Fig. 2F, panel I, and Table 3). In comparison, binding affinity of H3.1K36Me1 (K d ϳ19.7 M) was 85-fold weaker compared with H3.1K36Me2 peptide (Fig.  2F, panel I, and Table 3). Because the bromodomain, which is crucial for the interaction with H4K16Ac peptide, is devoid of tryptophan residues, we performed extrinsic fluorophore conjugation to H4 or H4K16Ac peptides with FAM. Subsequently, the FAM-conjugated peptides were titrated with the PBP module in increasing concentrations to score the interaction. Although no appreciable binding was seen for histone H4 peptide with the PBP module, H4K16Ac peptide showed significant interaction (K d ϳ0.75 M) (Fig. 2F, panel II, and Table 3). Thus, the specificity in the interaction of ZMYND8 with modified histone H3.1K36Me2/H4K16Ac was reconfirmed.
As ZMYND8 shows similar domain organization (Bromo-PWWP) as well as the amino acid sequence identity to ZMYND11 (33% for the Bromo-PWWP module) (Fig. 2G), Bromo-PWWP domain of ZMYND8 was modeled with reference to ZMYND11 structure (Fig. 2H). The modeled bromodomain structure consists of a four-helix bundle interspersed by ZA and BC loops, which is a conserved canonical structural feature (Fig. 2I). The asparagine residue (Asn-228), which is involved in the acetylated lysine recognition, has also been highlighted in stick representation (Fig. 2I). The PWWP domain adopts a five-bladed ␤-barrel fold with an ␣-helix at the C terminus (Fig. 2J). The H3.1K36Me2-binding aromatic residues (Phe-288, Trp-291, Phe-307, and Asp-312) are in stick representation in the PWWP structure (Fig. 2J). We mutated the active site residues of the bromodomain (i.e. Tyr-227 and Asn-228) and the aromatic cage of PWWP (Phe-307) to alanine. The Y227A/N228A double mutation in PBP module of ZMYND8 showed compromised binding to H4K16Ac peptide (Fig. 2K, panel I, compare lanes 5 and 6). Similarly, PWWP F307A mutation in the PBP module showed minimal interaction to H3.1K36Me2 (Fig. 2L, panel I, compare lanes 5 and 6). Y227A/N228A double mutant, however, showed no alteration in H3.1K36Me2 binding as compared with the wild type protein (Fig. 2K, panel II, compare lanes 5 and 6). In parallel, wild type and F307A mutant proteins showed significant H4K16Ac binding (Fig. 2L, panel II, compare lanes 5 and 6). This indicates that Tyr-227 and Asn-228 residues are critical for H4K16Ac binding, and the aromatic cage in PWWP is critical for H3.1K36Me2 binding. These results confirm that the bromodomain and the PWWP motif are responsible for interacting with H4K16Ac and H3.1K36Me2 respectively, although the comprehensive organization of PBP, harboring certain critical residues in the modified histone binding pocket, is essential for mediating this interaction.
ZMYND8 Interacts with H3.1K36Me2 and H4K16Ac Histones ex Vivo-We subsequently wanted to score the binding preference of ZMYND8 for canonical or variant histone in an ex vivo context for which FLAG-H3.1-or FLAG-H3.3-transfected HEK293 cells were subjected to M2-agarose pull down and immunoblotted with anti-ZMYND8 antibody. We observed a preferential interaction of histone H3.1 with ZMYND8 and histone H3.3 with ZMYND11 (Fig. 3A, panel I) and subsequently quantified each of these interactions (Fig. 3A). A sequence alignment of ZMYND8 and ZMYND11 (Fig. 3A, panel II) and H3.1 and H3.3 (Fig. 3A, panel III) was performed to understand these binding preferences.
We subsequently performed co-IP assays to understand its modified histone binding preference ex vivo. Endogenous ZMYND8 showed robust interaction with H3K36Me2, compared with H3K36Me3 (Fig. 3B). Here, H4K16Ac was also found to be associated with endogenous ZMYND8 (Fig. 3B). Genome-wide distribution of H3K36Me2 and H4K16Ac epigenetic marks has been previously reported (27)(28)(29). The NSD2 (Nuclear receptor SET Domain-containing family of histone lysine methyltransferase) target genes, for example TGFA, MET, and RRAS2, were found to have high H3K36Me2 occupancy (Fig. 3C, panel I). However, several key developmental genes, including HOXA9, WNT1, and TP53, showed a high abundance of H4K16Ac (Fig. 3D, panel  I). Interestingly, we observed a significant occupancy of ZMYND8 at these H3K36Me2-and H4K16Ac-enriched gene loci in HEK293 cells (Fig. 3, C, panel II, and D, panel II). In contrast, ZMYND8 had meager binding to the region of the genome with a low abundance of H3K36Me2 and H4K16Ac levels (marked as -ve Control in panels I and II of  3, C and D), implicating that the recruitment of the protein to the genome is through its ability to interact with these modified histones.

Minimal PBP Module of ZMYND8 Is Critical for Its Recruitment to the Modified Histones of Chromatin-
We hypothesized that the PBP module of ZMYND8 is the minimal requirement for its chromatin-binding ability. To confirm this, FLAG-conjugated constructs of wild type and PBP-deleted ZMYND8 were transfected into HeLa cells, and immunofluorescence staining with anti-FLAG and anti-H3K36Me2/-H4K16Ac was performed. The extent of co-localization of H3K36Me2 with wild type (Pearson's coefficient 0.78) was more significant than with PBP-deleted (Pearson's coefficient 0.47) (Fig. 4A, compare panels I and II). Similarly, H4K16Ac and wild type showed better co-localization (Pearson's coefficient 0.68) than with the mutant (Pearson's coefficient 0.29) (Fig. 4A, compare panels III and IV). Co-IP assays with anti-FLAG antibody also established that unlike wild type, PBP deleted-ZMYND8 showed a poor interaction with H3K36Me2 and H4K16Ac marks (Fig. 4B). However, both wild type and PBP-deleted ZMYND8 showed comparable interaction with HDAC1 and CHD4 (Fig. 4B). To further confirm whether the MYND domain is also involved in mediating the interaction with histones, we expressed and purified the GST-MYND domain of ZMYND8. In a pulldown assay with H4K16Ac and H3.1K36Me2 peptides, specific interaction could be detected with GST-PBP but not GST-MYND (Fig. 4C). Finally, the requirement of the PBP module of the protein for being targeted to its binding sites was assessed. Here, we selected TP53 and RRAS2 (24 kb) as they showed significant enrichment of H4K16Ac and H3K36Me2, respectively (Fig. 4D). The genes harboring H3K36Me2 and H4K16Ac had high occupancy of wild type but not the PBP-deleted construct as seen by ChIP assays (Fig. 4D). These results clearly indicate that the minimum PBP module of ZMYND8 is critical for its recruitment to chromatin through its modified histone interacting partners.  ). B, interaction of ZMYND8 with modified histones. Endogenous ZMYND8 was co-immunoprecipitated from HEK293 cells and immunoblotted for ZMYND8, H3K36Me2, H3K36Me3, and H4K16Ac. IgG pull down served as a negative control. C, ChIP assays were performed in HEK293 cells with ␣-H3K36Me2 and ␣-ZMYND8 antibodies to check for enrichment of H3K36Me2 (panel I) and ZMYND8 (panel II) at RRAS2, TGFA, and MET gene loci and a negative control region (depleted of H3K36Me2 and H4K16Ac). Two separate regions of each of the genes (5 and 24 kb for RRAS2, 22 and 57 kb for TGFA, and 3 and 26 kb for MET, downstream from start site) were selected for designing primers. D, similarly, ChIP assays were performed in HEK293 cells with ␣-H4K16Ac and ␣-ZMYND8 antibodies to show enrichment of H4K16Ac (panel I) and ZMYND8 (panel II) at HOXA9, WNT1, and TP53 gene loci and a negative control region (depleted of H3K36Me2 and H4K16Ac). Two separate regions of 646 and 1890 bp downstream from the start site were selected for designing primers from HOXA9 gene. WNT1 primers were for promoter region, whereas the TP53 primer was designed 11 kb downstream. Quantitative PCR was done, and relative fold was plotted normalizing by IgG. For H3K36Me2 and H4K16Ac, ChIP relative fold was further normalized to endogenous H3 and H4, respectively. At least three separate experiments were performed. Error bars show standard deviation.

ZMYND8 Is Itself an ATRA-responsive Gene and Regulates
Other RA-inducible Genes-To understand the role of ZMYND8 in regulating transcription in chromatin context, we used ATRA, which is a well studied transcription inducer (11). At first, we attempted to check whether there was an alteration in expression of ZMYND8 upon ATRA treatment. For this, we selected SH-SY5Y cells that showed elevated expression of TAU, TUJ1, NAV1.2, and SNAP25 (30 -32) and down-regulation of REST, a pluripotency marker leading to neurite formation (33) in response to RA. A similar alteration of gene expres-sion upon RA stimulation was observed (Fig. 5A). ZMYND8 gene showed a progressive hyper-expression after 2 or 4 days of ATRA treatment (Fig. 5A). Similarly, ZMYND8 increased in protein abundance with progressive days of RA treatment (Fig.  5B) where TAU was used as a positive control.
The increase in ZMYND8 in both RNA and protein levels led us to perform its promoter analysis. We observed the presence of RARE sequence at the upstream region of ZMYND8 in both human and mouse genome and subsequently cloned them in pGL3-basic vector to check their promoter activity. Indeed, we FIGURE 4. Chromatin binding module of ZMYND8 is essential for its histone interaction. A, HeLa cells were transiently transfected with FLAG-WT ZMYND8 or FLAG-⌬PBP ZMYND8, co-stained with ␣-FLAG and either ␣-H3K36Me2 or ␣-H4K16Ac antibodies. The Pearson's coefficient for FLAG-WT ZMYND8 and H3K36Me 2 /H4K16Ac was Ͼ0.5, whereas that for FLAG-⌬PBP ZMYND8 and H3K36Me 2 /H4K16Ac was Ͻ0.5. B, HEK293 cells were transiently transfected with FLAG-WT or FLAG-⌬PBP ZMYND8, and whole cell extracts were subjected to M2-agarose FLAG pull down and immunoblotted with ␣-FLAG, ␣-H3K36Me2, ␣-H4K16Ac, ␣-HDAC1, and ␣-CHD4 antibodies. C, interaction of GST-PBP or GST-MYND with biotinylated H4K16Ac and H3.1K36Me2 peptides. Western blot analysis was done with ␣-GST antibody. D, ChIP assay was performed after transiently transfecting FLAG-WT or FLAG-⌬PBP ZMYND8 in HEK293 cells with ␣-FLAG antibody. Relative fold enrichment at TP53 (11 kb) and RRAS2 (24 kb) gene loci was scored by quantitative PCR normalized to IgG. At least three separate experiments were performed. Error bars show standard deviation. FEBRUARY 5, 2016 • VOLUME 291 • NUMBER 6 found an ATRA-induced firing of the luciferase fusion gene constructs in both mouse (Neuro2A) and human (SH-SY5Y) cells (Fig. 5C, panels III and IV). Mutating the RARE sequence led to a significant reduction in ZMYND8 promoter activity in both the cells (Fig. 5C, panels III and IV). These results clearly indicate that ZMYND8 is itself an ATRA-responsive gene.

ZMYND8, Regulates ATRA-responsive Genes
Next, we looked for the occupancy of ZMYND8 on the different ATRA-responsive genes, including TUJ1, TAU, DRD2, NAV1.2 and SNAP25 (34) by ChIP assays (Fig. 5, D and E). We have selected primers for two different regions of the gene, promoter (harboring the RARE element), and gene body. Upon ATRA treatment, the occupancy of ZMYND8 at the promoter as well as gene body of TUJ1, TAU, DRD2, NAV1.2, and SNAP25 showed significant increase (Fig. 5, D and E). The enrichment of ZMYND8 at the promoter and gene body of the ATRA-responsive genes could be through its ability to interact with transcription machinery components or modified histones. We at first investigated the cause of selective enrichment of ZMYND8 at the gene body. H3K36Me2 and H4K16Ac have been shown to be significantly enriched at the gene body by ChIP-Seq studies in different cell lines (27)(28)(29). To understand whether the preferential recruitment of ZMYND8 to the gene body of the ATRA-inducible genes in SH-SY5Y cells is through its modified histone-binding ability, we analyzed the histone PTM on these genes (Fig. 5, F and G). We observed an increased occupancy of H3K36Me2 and H4K16Ac by RA in SH-SY5Y cells in TUJ1, TAU, DRD2, NAV1.2 and SNAP25 (Fig. 5, F and G). The reason of recruit-ment of ZMYND8 at the promoter elements of these genes has been investigated later.
In addition, we monitored the altered level of different epigenetic marks in ATRA-treated SH-SY5Y cells. Among H3K36Me1, H3K36Me2, and H3K36Me3 signatures, H3K36Me2 mark showed a significant elevated level (Fig. 5H). However, RA treatment increased H4K16Ac appreciably, although the H3K9Ac signature did not show significant change (Fig. 5H). Thus, ATRA modulates specific histone modifications not only in gene-specific manner but also in a global context. Our next attempt was to examine the involvement of ZMYND8 in the regulation of these ATRA-responsive genes. We knocked down ZMYND8 in DMSO or RA-treated SH-SY5Y cells. The reduction of ZMYND8 expression upon siRNA transfection was scored by qRT-PCR and Western blotting analysis (Fig. 5, I and J). Knocking down ZMYND8 compromised the expression of different ATRA-responsive genes, including TUJ1, TAU, DRD2, NAV1.2, and SNAP25 significantly (Fig. 5K). However, the expression of REST was marginally increased upon silencing of ZMYND8 (Fig. 5K). The expression of glial markers, including GFAP and VIM (36,37), remained almost unaltered following ZMYND8 silencing (Fig.  5K). These results clearly show that ZMYND8 gets recruited to the ATRA-responsive genes and regulates their expression.
ZMYND8 Shows Selective Interaction with RNA Polymerase II Phosphorylated at Ser-5 in a DNA Template-dependent Manner-We have previously shown that upon ATRA treatment, ZMYND8 recruitment at the promoter element (harboring RARE sequence) was significantly high (Fig. 5D). Hence, we hypothesized its possible association with RNA polymerase II transcription machinery. Interestingly, by performing coimmunoprecipitation assays, we found that ZMYND8 interacts with transcription-initiation-competent RNA pol II S5P but not with transcription elongation/termination-competent RNA pol II S2P complex (Fig. 6A, compare lanes 1 and 3). Reciprocal interactions indicated that RNA pol II S5P but not S2P interacted with ZMYND8 (Fig. 6, C and D). Here, interac-tion with RNA pol II non-C-terminal domain was used as control (Fig. 6, A and B). We subsequently wanted to understand whether there is any alteration in its preferential interaction with RNA pol II S5P in presence or absence of ATRA. By immunoprecipitating RNA pol II S5P, the interaction with ZMYND8 did not alter in DMSO-or ATRA-treated SH-SY5Y cells (Fig.  6E). The interaction of ZMYND8 and RNA pol II machinery could be DNA template-dependent or -independent. Because the interaction between the two did not alter the presence or absence of ATRA, we predicted this complex formation onto the DNA template. We set immunoprecipitation with RNA pol II S5P after DNase I digestion of the lysates. In accordance with our speculation, we indeed see a sharp decrease in its association with ZMYND8 upon DNase I digestion of lysates (Fig. 6F). To further confirm the role of ZMYND8 in regulating transcription of ATRA-responsive genes, we performed ChIP assays with RNA pol II S5P following silencing of ZMYND8 gene expression. Indeed, we observed a significant reduction of RNA pol II S5P recruitment to chromatin in the absence of ZMYND8, even after ATRA induction (Fig. 6G). These results clearly indicate that ZMYND8 is a key component of transcription machinery complex having a significant role in mediating transcription initiation.
ZMYND8 Regulates Global Gene Transcription in CHD4-independent Manner-ZMYND8 is a component of the coregulator complex network having an extensive interaction with transcription machinery (5,38). Recently ZMYND8 has been established to be a component of the NuRD chromatin remodeling complex (39) and can influence gene transcription (5). In an attempt to comprehend whether ZMYND8 can regulate global gene expression only when it remains associated with the NuRD complex, we performed microarray analysis in ZMYND8 knocked down HeLa cells. Overall, 331 genes were up-regulated and 438 genes were down-regulated upon ZMYND8 silencing. The cluster diagram is represented in Fig.  7A. The microarray data has been validated by qRT-PCR analysis where candidate up-regulated (E2F3, TFEB, and BMP8A) FIGURE 5. ZMYND8 is ATRA-responsive and regulates other RA-inducible genes. A, ATRA treatment of SH-SY5Y cells for 2 and 4 days, respectively, leads to an increase in ZMYND8 RNA level. TAU, TUJ1, NAV1.2, SNAP25, and REST were used as controls. Relative mRNA level was plotted for ATRA-treated differentiated cells over DMSO-treated undifferentiated cells and normalized to GAPDH. At least three individual experiments were performed. Error bars show standard deviation. B, alteration in the expression of ZMYND8 protein abundance in SH-SY5Y cells upon ATRA treatment for 2 and 4 days. Western blotting was done by probing with ␣-ZMYND8, ␣-TAU, and ␣-GAPDH antibodies. GAPDH was used as loading control. At least three individual experiments were performed. Error bars show standard deviation. Blots were quantified using ImageJ software from National Institutes of Health. C, schematic diagram shows the upstream regulatory region of mouse (panel I) and human (panel II) ZMYND8 genes with the arbitrary position of RARE and the primer pairs used to clone those regions into pGL3-basic vector. The respective reporter constructs were transiently transfected into Neuro2A (panel III) and SH-SY5Y (panel IV) cells, respectively, with/without ATRA treatment followed by luciferase assay. The activities are shown as mean fold enhancement compared with the empty vector after normalization with Renilla luciferase activity. Here, the constructs containing mutated RARE sequence are named as RARE-mut. Each transfection was performed in triplicate, and the experiments were repeated at least three times. Error bar shows standard deviation. D-G, relative enrichment of ZMYND8 on the RARE harboring sequence of promoter of RA-responsive genes (TUJ1, TAU, DRD2, NAV1.2, and SNAP25) was observed (D). Similar enrichment of ZMYND8 (E), H3K36Me2 (F), and H4K16Ac (G) onto gene body region of TUJ1, TAU, DRD2, NAV1.2, and SNAP25 genes was monitored. ChIP assays were performed after 4 days of ATRA treatment in SH-SY5Y cells with ␣-ZMYND8 antibody. Relative fold was calculated by normalizing ZMYND8 with IgG for both DMSO (control) and ATRA-treated cells. At least three separate experiments were performed. Error bars show standard deviation. H, alteration in the expression of different epigenetic marks upon ATRA treatment for 2 and 4 days. Western blotting was done by probing with ␣-H3K36Me1, ␣-H3K36Me2, ␣-H3K36Me3, ␣-H4K16Ac, and ␣-H3K9Ac antibodies from whole cell extracts and quantified by ImageJ software from National Institutes of Health. ␣-H3 and ␣-H4 are used as loading controls for respective modifications. At least three separate experiments were done. Error bars show standard deviation. I, ZMYND8 silencing was done in 2 days of ATRA-or DMSO-treated SH-SY5Y cells. Similar experiments were performed with non-targeting siRNA, as negative control. Total mRNA level was analyzed by quantitative PCR, and GAPDH was used for normalization. Reference level was considered as 1. J, ZMYND8 knockdown in SH-SY5Y cells was scored by Western blotting after ATRA treatment for 2 days, with ␣-ZMYND8 antibody. GAPDH was used for normalization. Blots were quantified using ImageJ software from National Institutes of Health. K, alterations in the expression of neuronal markers (TUJ1, TAU, DRD2, NAV1.2, and SNAP25), glial markers (GFAP and VIM), and pluripotency marker REST were measured after ZMYND8 siRNA (or a non-targeting siRNA as negative control) transfection during 2 days of ATRA or DMSO treatment in SH-SY5Y cells. Total mRNA level was scored by quantitative PCR normalizing to GAPDH. Relative mRNA level was plotted as ATRA-treated differentiated cells over DMSO-treated undifferentiated cells. Reference level considered as 1. At least three separate experiments were done. Error bars show standard deviation. FEBRUARY 5, 2016 • VOLUME 291 • NUMBER 6

JOURNAL OF BIOLOGICAL CHEMISTRY 2675
and down-regulated (HIPK3, BARD1, and ZNF687) genes showed differential expression upon ZMYND8 silencing (Fig.  7B). Previous studies have shown that silencing CHD4 (an important chromatin remodeling subunit of NuRD) can alter global gene expression in HeLa cells (40). To understand the NuRD-associated functions of ZMYND8, we analyzed the differentially expressed gene sets upon silencing ZMYND8 with that of CHD4 (40). On careful examination, we found that only 40 up-and 71 down-regulated genes were commonly regulated by these factors (Fig. 7C). A detailed comparison of the clusters (Fig. 7D) and overlapping Venn diagram (Fig. 7E) of up/downregulated genes upon ZMYND8 or CHD4 silencing has been performed. Biological Analysis Network for the common gene regulatory pathways is represented in Fig. 7F. These results clearly indicate that ZMYND8 has a significant CHD4-independent role on a global platform.
ZMYND8 Is a Key Regulator of a Subset of ATRA-inducible Genes-While analyzing the microarray results upon silencing ZMYND8, we found that a number of genes, which were regulated by ZMYND8, also harbored the RARE motif and hence FIGURE 6. ZMYND8 interacts with RNA polymerase II phospho-Ser5 in a DNA-dependent manner. A, ZMYND8 was co-immunoprecipitated from SH-SY5Y cells and immunoblotted with ␣-RNA polymerase II non-CTD (epitope at N terminus), ␣-RNA polymerase II phospho-Ser-5, and ␣-RNA polymerase II phospho-Ser-2 antibodies. IgG serves as negative control. B-D, similarly, RNA polymerase II non-CTD (B), RNA pol II phospho-Ser-5 (C), and RNA pol II phospho-Ser-2 (D) was co-immunoprecipitated from SH-SY5Y cells and immunoblotted with ␣-ZMYND8 antibody. E, co-IP of RNA pol II phospho-Ser-5 from 4 days ATRA-or DMSO-treated SH-SY5Y cells followed by immunoblotting with ␣-ZMYND8 antibody. IgG serves as negative control. F, RNA pol II phospho-Ser-5 was coimmunoprecipitated from SH-SY5Y cells after DNase I treatment of lysates and immunoblotted with ␣-ZMYND8 antibody. IgG serves as negative control. G, ChIP assays were performed after ZMYND8 siRNA or a non-targeting siRNA transfection during 2 days of ATRA or DMSO treatment in SH-SY5Y cells with ␣-RNA polymerase II phospho-Ser-5 antibody. Relative fold was calculated by normalizing RNA pol II phospho-Ser-5 with IgG and represented as ATRA-treated differentiated cells over DMSO-treated undifferentiated cells.
can be considered as ATRA-responsive. A detailed list of differentially expressed genes as a sequel to ATRA treatment has been previously reported in HeLa cells (41). Upon ZMYND8 silencing, we found an overlap of at least 90 up-and 83 downregulated candidates that are ATRA-inducible (Fig. 8A). Detailed cluster analyses of the up/down-regulated genes, which are simultaneously ATRA-responsive, have been represented in Fig. 8B. We further validated candidate up-and down-regulated genes upon silencing ZMYND8, which also harbors RARE motifs by qRT-PCR (Fig. 8C). Biological network for differentially expressed genes upon silencing of ZMYND8 as well as harboring the RARE motif are represented in Fig. 8D. These results clearly indicate that ZMYND8 is indeed a regulator of a fraction of ATRA-inducible genes in a global context.

Discussion
Multifunctional transcription factors harbor conserved chromatin-binding modules and show differential affinity for chromatin (42,43). Essentially, these transcription factors are histone-interacting proteins (14), like ZMYND8, which has a distinct preference for modified histone H3 and H4 (Fig. 3B). Furthermore, ZMYND8 shows a preference for canonical histone H3.1, unlike ZMYND11 (Fig. 3A, panel I). A sequence alignment between these two proteins shows that certain key residues, including Asn-263 and Arg-265 (marked in blue), which were involved in mediating H3.3 binding through the H-bonding interaction in ZMYND11, is replaced by Cys-263 and Lys-265 in ZMYND8 (Fig. 3A, panel II). This possibly makes a difference in its choice for canonical histone H3.1 over its variant H3.3. However, among the five amino acids that are different between H3.1 and H3.3, Ala-31 of canonical H3.1 is substituted by Ser-31 in variant H3.3 and is probably an important determining factor for the preference of ZMYND8 (Fig.  3A, panel III).
The transcription factors with conserved chromatin-binding modules preferentially bind to modified histones (14). For HeLa cells. E, overlapping Venn diagram for common differentials between ZMYND8 and CHD4. The pink circle represents up-regulated genes in siCHD4, and the green circle is for down-regulated genes in siCHD4; the yellow circle is for up-regulated genes in siZMYND8, and the blue circle represents down-regulated genes in siZMYND8. As they have common differential genes, so the overlapping areas are in different colors. F, networks for ZMYND8 knocked down and CHD4 knocked down HeLa cells based on the common biology. The genes whose regulation are opposite in the two networks are bordered in black. The genes are colored according to their fold change and sized based on their p value which is Ͻ0.05. Green and red ellipses denote down-regulation and up-regulation, respectively, and the processes are denoted by a blue rectangle. FEBRUARY 5, 2016 • VOLUME 291 • NUMBER 6 example, the CREB-binding protein bromodomain binds to H3K36Ac/H3K56Ac/H4K20Ac (44,45) in a monovalent manner, and bromodomain PHD finger transcription factor/PHD finger binds to H4K16Ac/H3K4Me3 (15) in a multivalent fashion. In our study, ZMYND8 preferentially binds to H3.1K36Me2/Me3 (Fig. 2, C, panel I, and F, panel I) and

ZMYND8, Regulates ATRA-responsive Genes
H4K16Ac (Fig. 2, D, panel I, and F, panel II). Although isolated bromodomain of ZMYND8 could bind to H4K16Ac (Fig. 2D, panel II), PWWP domain could interact with H3.1K36Me2/ Me3 only when it was integrated along with bromodomain and PHD finger (Fig. 2C, panels I-III). For stabilization of conformation of the aromatic pocket in PWWP domain, the presence FIGURE 8. ZMYND8 modulates RARE-harboring genes. A, Venn diagram showing the overlap between up-and down-regulated genes in ZMYND8 knocked down and total RARE-harboring genes in HeLa cells. B, clustering and heat maps of up-and down-regulated genes in ZMYND8 knocked down and total RARE-harboring genes in HeLa cells. C, validation of microarray analysis. Bars show candidate RARE-harboring, up-and down-regulated genes after knocking down ZMYND8 in HeLa cells. 18S rRNA was used for normalization. Reference level was considered as 1. At least three separate experiments were done. Error bars show standard deviation. D, gene regulatory network for differentially expressed genes in ZMYND8 knocked down and total RARE-harboring genes in HeLa cells. Genes are in ellipses, and processes are in rectangles, and miRNAs are triangular in shape. The genes whose regulation are opposite in the two networks are bordered in black. Green denotes down-regulation, and red denotes up-regulation. The color shade of nodes differs due to their fold change. The size of the nodes is as per their p value which is Ͻ0.05. of other supporting domains, i.e. bromodomain and PHD finger, is essential. This comprehensive reader function can also be seen in case of ZMYND11 (14). To understand the binding preference of ZMYND8 for H4K16Ac in comparison with ZMYND11, we performed sequence alignment between the two proteins (Fig. 2G) and subsequently molecular modeling studies (Fig. 2, H-J). Several bromodomain proteins harbor YN residues in their active site (46). We found that ZMYND8 has Tyr-Y227 and Asn-228 amino acids in its bromodomain, which are absent in ZMYND11. The site-directed mutagenesis approach confirmed the Y227A/N228A mutant has reduced H4K16Ac binding ability (Fig. 2K). To find out the preference for modified histone recognition by ZMYND8 in cells, we performed the co-IP experiment that supported the better interaction between ZMYND8 and H3K36Me2 as compared with H3K36Me3 (Fig. 3B), in conformity with in vitro fluorescence titration assays (Fig. 2F). In parallel, the co-IP assay also confirmed the interaction of ZMYND8 with H4K16Ac ex vivo (Fig.  3B).
H3K36 methylation by Set2 recruits the Rpd3S histone deacetylase complex, which deacetylates the neighboring nucleosomes preventing cryptic transcription initiation within the gene body (47). Because ZMYND8 binds to H3K36Me2, it might have additional roles in modulating nucleosomal integrity over the coding regions preventing cryptic transcription initiation. Remarkably, the human orthologue of Drosophila melanogaster MOF (males-absent-on-the-first, a key member of the MYST family (hMOF)-mediated acetylation of H4K16) causes a recruitment of BRD4 and pTEFB to chromatin leading to a release of RNA polymerase II pausing (48,49). In this context ZMYND8 through its chromatin reader ability might be instrumental in recruiting these factors onto chromatin thereby promoting processive transcription. The di-and trimethylated H3K36 at the gene body shows a distinct distribution pattern, where H3K36Me2 is adjacent to promoters and H3K36Me3 accumulates at the 3Ј end of the gene (47,50). Although H3K36Me2 leads to elevated H4K16Ac, H3K36Me3 compensates it (50). The preferential interaction of ZMYND8 to H3K36Me2/H4K16Ac indicates its ability to modulate higher order chromatin structure and hence transcription. This could be possible through its ability to act as a chromatin reader, which further helps in recruiting a chromatin remodeling machinery complex steering transcription activation. Interestingly, as ZMYND8-enriched genes (RRAS2, TGFA, MET, HOXA9, TP53, and WNT1) (Fig. 3, E and F, panel II) are intricately involved in pathways regulating cell proliferation, differentiation, and development (51)(52)(53)(54)(55)(56), our results indicate a possible involvement of ZMYND8 in these processes. We elucidated that a higher occupancy of ZMYND8 to these developmentally regulated genes is due to its chromatin reader function, as there is also a significant enrichment of H3K36Me2/ H4K16Ac marks in these genes (Fig. 3, C and D, panel I). Furthermore, we observed an increased recruitment of ZMYND8 at the gene body of ATRA-inducible genes (TUJ1, TAU, DRD2, NAV1.2 and SNAP25) through its modified histone-binding ability (Fig. 5, E-G). However, we also observed robust enrichment of ZMYND8 at the promoter harboring the RARE sequence of these RA-responsive genes (Fig. 5D) and speculated its association with the core transcription machinery.
As RNA polymerase II transcribes a gene, there is a distinct alteration in its C-terminal domain phosphorylation status (57), which is instrumental in a dynamic association of factors that interact with the enzyme. We observed a strong interaction of ZMYND8 with RNA pol II S5P in a DNA-mediated manner (Fig. 6, A, C, and F). This interaction could be steered by some other factor(s) that promotes transcription initiation and brings ZMYND8 and RNA pol II complex to the promoter site. ZMYND8 might also be assisting in the recruitment of the RNA pol II complex through its selective binding to the H3K36Me2enriched region of the genome located at the promoter proximal site. ZMYND11, however, modulates the activity of RNA pol II S2P at the 3Ј end of the gene through its H3.3K36Me3 recognizing ability (14). Interestingly, the association of ZMYND8 with RNA pol II S5P is retained in absence of ATRA (Fig. 6E), further confirming that a stoichiometric complex formation between the two requires a DNA template. Indeed, we observed that the interaction is significantly reduced in the presence of DNase I (Fig. 6F). Furthermore, in the absence of ZMYND8, there was a significant reduction in RNA pol II S5P occupancy even upon ATRA induction, indicating that ZMYND8 is an indispensable component of initiating RNA pol II machinery (Fig. 6G).
Finally, we wanted to understand the role of ZMYND8, an integral component of transcription machinery, in regulating global gene expression. Previous studies have reported that ZMYND8 is a functional component of the transcription repression complex NuRD (58,59) and targets it to the DNA damage response site (23). CHD3/4 is responsible for the ATPdependent chromatin remodeling function of NuRD (35,60). In this study, we tried to understand the CHD4-independent function of ZMYND8 in perspective of gene expression. Microarray analysis clearly indicated that ZMYND8 regulates a large class of genes that are unique from those regulated by CHD4 (Fig. 7). Interestingly, subsets of genes that are differentially expressed upon silencing of ZMYND8 also harbor the RARE motif and hence are ATRA-responsive (Fig. 8). ZMYND8 itself harbors the RARE sequence and is regulated by RA (Fig. 5, B and C). Because ZMYND8 has a preference for canonical histone H3.1 over variant H3.3, we also analyzed whether the gene targets of ZMYND8 have a depleted occupancy of H3.3. However, a significant number of ZMYND8 target genes is found to harbor H3.3 (data not shown). It would be interesting to understand the mechanism of transcription regulation by ZMYND8 in H3.1 versus H3.3-enriched genes.
In summary, we found that ZMYND8, a transcription regulator, is a putative chromatin reader. Through its specific motifs as well as amino acid residues, it exerts its chromatin-binding ability. Furthermore, it shows its preferential interaction with canonical histone H3.1K36Me2 and H4K16Ac. Interestingly, ZMYND8 is itself an RA-inducible gene, and through its modified histone-binding ability it regulates the expression of other ATRA-responsive genes. Furthermore, we show that ZMYND8 is indeed associated with transcription initiation-competent RNA pol II S5P in a DNA template-dependent manner. Although it is a transcription factor, the mechanism underlying ZMYND8, Regulates ATRA-responsive Genes FEBRUARY 5, 2016 • VOLUME 291 • NUMBER 6 the gene expression regulation by ZMYND8 has remained elusive. Hence, this study is an attempt to decipher the so far uncharacterized ability of ZMYND8 to modulate gene expression. Further study may uncover its possible role in lineage commitment during development.
Author Contributions-C. D. and S. R. conceived the study and designed the experiments. S. A., S. S., M. B., I. S., S. S., and D. K. S. performed the experiments and analyzed the data. C. D., S. R., M. B., S. A., and S. S. wrote the paper. All authors discussed the paper and commented on the manuscript.