Construction of histone H4-HaloTag expression plasmids
The human histone H4 (H4C1) coding region was amplified from the total cDNA library of RPE-1 cells and displaced H2B in PB-EF1a-H2B-HaloTag-IRES-Puro previously made in our laboratory. The plasmid for tailless histone H4-HaloTag, H4(∆1–19)-HaloTag, was obtained from PB-EF1a-H4-HaloTag-IRES-Neo by incomplete amplification of the H4 coding region, excluding amino acid residues 1–19. Later, PB-EF1a-H4(∆1–19)-HaloTag-IRES-Neo was used to make the constructs to express histone H4 quadruple mutants, H4K(5,8,12,16)mut-HaloTag, by reconstruction of the H4 N-terminal tail with PCR primers including four identical pre-synthesized mutations: AAG (Lys) -> AGG (Arg), CAG (Gln), or GCG (Ala). To achieve stable and comparable ectopic expression levels of the fusion protein in HeLa cells, the coding regions for H4-HaloTag, H4(∆1–19)-HaloTag, and H4K(5,8,12,16)mut-HaloTag were all transferred to the Flp-In system vector (V602020; Invitrogen). The resulting plasmids were named pEF5-H4-HaloTag-FRT, pEF5-H4(∆1–19)-HaloTag-FRT, and pEF5-H4K(5,8,12,16)mut-HaloTag-FRT, respectively.
H4K16 point substitution mutants, H4K16R- and H4K16A-HaloTag, were obtained through PCR amplification of pEF5-H4-HaloTag-FRT with appropriate pairs of mutagenic PCR primers. H4K16Q-HaloTag was first obtained from PCR amplification of PB-EF1a-H4-HaloTag-IRES-Puro with mutagenic PCR primers. Then the H4K16Q-HaloTag fragment was transferred to the Flp-In system vector to obtain pEF5-H4K16Q-HaloTag-FRT.
HeLa cells containing the flippase recognition target (FRT) site in their genome, HeLa FRT-bla, were used as the model cell line. HeLa cells (Maeshima et al. 2006) were cultured in Dulbecco’s Modified Eagle Medium (DMEM, D5796-500ML; Sigma-Aldrich) supplemented with 10% fetal bovine serum (FBS, ABB213138; HyClone) at 37°C in 5% CO2.
To establish HeLa cells with stable ectopic expression of H4-HaloTag or mutant H4-HaloTag, the corresponding plasmids were co-transfected with pOG44 (V600520; Thermo Fischer Scientific) using the Effectene Transfection Reagent kit (301425; QIAGEN) in a 12-well plate. Two days following the transfection, the cell populations were expanded on 100-mm dishes with the medium containing 200 µg/mL hygromycin B (10687010; Thermo Fischer Scientific) for positive selection. The medium was refreshed every 3–4 days. After about 2 weeks, the clones were isolated. For stable ectopic expression of H4(Δ1–19)-HaloTag with PiggyBac system, PB-EF1a-H4(∆1–19)-HaloTag-IRES-Neo was co-transfected with pCMV-hyPBase (Yusa et al. 2011) using the Effectene Transfection Reagent kit in a 12-well plate. The selection was performed with 800 µg/mL G418 (08973-14; Nacalai Tesque).
Expression and localization of H4-HaloTag in cells
To assess the localization of H4-HaloTag as well as mutant H4-HaloTag in HeLa cells, the cells were seeded on 1 mg/mL poly-L-lysine-coated (P1524-500MG; Sigma-Aldrich) coverslips (C018001; Matsunami). The HaloTag was labeled with an excess amount of ligand, 5 nM tetramethylrhodamine (TMR, 8251; Promega), for 30 min at 37°C in 5% CO2. The rest of the procedures were conducted at room temperature.
After briefly rinsing with 1×PBS, the cells on coverslips were fixed with 1.85% formaldehyde (064–00406; Wako) in 1×PBS for 15 min. Excess molecular crosslinking and fixation were quenched with 50 mM glycine (077–00735; Wako) in 1×PBS for 5 min. The cells were then permeabilized with 0.5% Triton X-100 (T-9284; Sigma-Aldrich) in 1×PBS for 5 min. After briefly rinsing with 1×PBS, the DNA in cells was stained with 0.5 µg/mL 4’,6-diamidino-2-phenylindole (DAPI, 10236276001; Roche) in 1×PBS for 5 min. Finally, after rinsing with 1×PBS twice, the coverslips were mounted onto micro slide glasses (S011120; Matsunami) in p-phenylenediamine (PPDI; 20 mM Hepes [pH 7.4], 1 mM MgCl2, 100 mM KCl, 78% glycerol, and 1 mg/mL paraphenylene diamine [695106-1G; Sigma-Aldrich]) and sealed with nail polish (T and B; Shiseido).
Z-stack images (every 0.4 µm along the Z axis, 20 sections in total) of the cells were taken using a DeltaVision Elite (or DeltaVision Ultra) epifluorescence microscope (Applied Precision) equipped with an Olympus PlanApoN 60× NA 1.42 objective, sCMOS camera, InsightSSI light (~ 50 mW), and the four-color standard filter set. DeltaVision image acquisition software, SoftWorx 7.X, was used to capture and perform deconvolution of acquired images in all cases. Image fusion, projection over the entire nucleus, and quantification were done using Fiji (ImageJ).
H4-Halo labeled nucleosome pull-down, purification of nucleosomal DNA, and library preparation.
Cell nuclei were isolated as previously described (Lewis and Laemmli 1982; Shimamoto et al. 2017) with minor modifications. Briefly, HeLa cells were maintained at 37°C and 5% CO2 in DMEM supplemented with 10% FBS on a cell culture plate. Collected cells were washed with nuclei isolation buffer composed of 3.75 mM Tris-HCl (pH 7.5), 20 mM KCl, 0.5 mM EDTA, 0.05 mM spermine, 0.125 mM spermidine, 0.1% Trasylol (14716; Cayman Chemical), and 0.1 mM phenylmethylsulfonyl fluoride (PMSF; P7626-1G; Millipore Sigma) by two cycles of centrifugal spins at 193×g for 7 min at 23°C. The pellets were then resuspended in nuclei isolation buffer containing 0.05% Empigen BB detergent (45165-50ML; Sigma-Aldrich) (nuclei isolation buffer+) and immediately homogenized with 10 downward strokes using a tight Dounce pestle. After a 5-min centrifugation of the cell lysates at 440×g, the pellets were washed with nuclei isolation buffer+. After centrifugation at 440×g at 4°C for 5 min, the supernatant was removed and resuspended in 10 ml of nuclei isolation buffer+. This step was repeated four times before the samples were ready for nucleosome purification.
Chromatin purification was carried out as described by (Ura and Kaneda 2001; Iida et al. 2023), with some modifications. The nuclei (equivalent to ~ 2 mg of DNA) in nuclei isolation buffer (10 mM Tris–HCl, pH 7.5, 1.5 mM MgCl2, 1.0 mM CaCl2, 0.25 M sucrose, and 0.1 mM PMSF) were digested with 50 U of micrococcal nuclease (LS004797; Worthington) at 30°C for 1 h. In this digestion step, 1 µM HaloTag PEG-Biotin Ligand (G859A; Promega) or an equivalent amount of DMSO was added. After centrifugation at 440×g at 4°C for 5 min, the nuclei were lysed with lysis buffer (10 mM Tris–HCl, pH 8.0, 5 mM EDTA, and 0.1 mM PMSF) on ice for 5 min. The lysate was dialyzed against dialysis buffer (10 mM HEPES-NaOH, pH 7.5, 0.1 mM EDTA, and 0.1 mM PMSF) at 4°C overnight using Slide-A-Lyzer (66380; Thermo Fisher Scientific). The dialyzed lysate was centrifuged at 20,400×g at 4°C for 10 min. The supernatant was recovered and used as the purified nucleosome fraction (mainly mono nucleosomes). To verify complete MNase digestion, DNA was purified, electrophoresed on a 1.5% agarose gel, and visualized by staining with ethidium bromide (Fig. 2b).
The biotinylated or untreated nucleosome fraction (28 µg) was diluted with an equal volume of 2× binding buffer (10mM HEPES-NaOH pH 7.5, 200 mM NaCl, 2 mM DTT, and cOmplete™ EDTA-free protease inhibitor cocktail (11873580001; Roche). As shown previously (Ide et al. 2021), the chromatin solution was added to streptavidin-FG beads (TAS8848 N1170, Tamagawa-Seiki) pre-equilibrated with coating buffer (PBS, 2.5% BSA, 0.05% Tween 20) for 1 h. The mixtures were incubated at 5 ºC for 2 h on a shaker. The beads were collected with a magnetic rack and subjected to the washing procedure as described by Gatto et al. (Gatto et al. 2022). Briefly, the beads were washed once with 1 mL of cold 1× binding buffer (10 mM HEPES-NaOH pH 7.5, 100 mM NaCl, 1 mM DTT,), twice with 1 mL of cold wash buffer 1 (HEPES-NaOH pH 7.5, 140 mM NaCl, 1% Triton X-100, 0.5% NP40, 0.1% SDS), twice with 1 mL of cold wash buffer 2 (HEPES-NaOH pH 7.5, 360 mM NaCl, 1% Triton X-100, 0.5% NP40, 0.1% SDS), twice with 1 mL of cold wash buffer 3 (HEPES-NaOH pH 7.5, 250 mM LiCl, 0.5% Triton X-100, 0.5% NP40) and once with TE buffer (10 mM Tris–HCl pH 8.0, 1 mM EDTA). Beads were resuspended in 50 µL TE buffer.
For DNA extraction, 40 µL of beads underwent RNase digestion (4 µL RNase A, 1 mg/mL, 30 min at 37°C) and subsequent proteinase K digestion in the presence of SDS (4 µL of 20 mg/mL proteinase K, and 10 µL of 10% SDS, at 1,200 rpm and 50 ºC on shaker for 1 h). Beads were removed by magnetic separation and DNA extracted using AMPure XP beads (A63880, Beckman Coulter) according to the manufacturer’s instructions, and eluted in Milli-Q water. Total DNA content of each sample was measured using Qubit (Q32851, Thermo), and the quality of DNA samples was assessed by an Agilent 2100 Bioanalyzer using an Agilent High Sensitivity DNA kit (5067 − 4626, Agilent). cDNA libraries were synthesized by the ThuPLEX DNA-seq Kit (R400675, Takara). The size distributions of the libraries were checked by an Agilent 2100 Bioanalyzer using an Agilent High Sensitivity DNA kit. Pooled amplicon library was sequenced with paired-end 2 × 100 bp reads on the Illumina NovaSeq 6000 platform.
For protein analysis of the HaloTag pull-down fraction (Fig. 2C), the remaining 10 µL of beads were resuspended with an equal volume of 2× Laemmli sample buffer (Laemmli 1970) supplemented with 10% 2-mercaptoethanol (133–1457; Wako) and incubated at 99 ºC for 10 min. Beads were removed by centrifugation at 10,000×g for 3 min. The Input sample and the pull-down samples from the biotinylated and untreated nucleosome (negative control) were separated using SDS–polyacrylamide gel electrophoresis and were subjected to western blotting by anti-HaloTag antibody (G9211, Promega).
For data analysis of purified H4-Halo labeled nucleosomal DNA, the nf-core (Ewels et al. 2020, 2022) ChIP-seq pipeline (nfcore/chipseq: version 2.0.0) was utilized, employing a Docker configuration profile and executed using the default parameters. The human genome hg19, retrieved from Illumina’s iGenomes, was used for the reference genome. Peaks were detected using the broad mode of MACS2 (Zhang et al. 2008). To evaluate the overlap between the H4-HaloTag labeled regions and contact domains (i.e., Hi-C A-compartment and histone modification patterns), we compared the H4-HaloTag regions with previously published data as described in (Nozaki et al. 2023). The contact domain data from HeLa defined by (Rao et al. 2014) were downloaded from Gene Expression Omnibus, accession number GSE63525 (GSE63525_HeLa_Arrowhead_domainlist.txt). Compartment annotations on the HeLa cell genome are SNIPER predictions obtained from (Xiong and Ma 2019). Peak lists on hg19 genome defined by the ChIP-seq analysis for the histone modifications were obtained from the ENCODE portal (Consortium 2012) (Davis et al. 2018) with the following identifiers: H3K4me2 (ENCFF108DAJ), H3K4me3 (ENCFF447CLK), H3K9ac (ENCFF723WDR), H3K27ac (ENCFF144EOJ), H3K79me2 (ENCFF916VLX), H3K36me3 (ENCFF001SVY), H3K4me1 (ENCFF162RSB), and H3K27me3 (ENCFF252BLX).
Single nucleosome tracking and analysis
Image processing and preparation for single nucleosome tracking were performed using ImageJ Fiji. Sequential 16-bit images of single H4-HaloTag molecules were kept (for Spot-On analysis) or reduced to an 8-bit grayscale (for MSD analysis). The background signal was subtracted using the rolling ball method (radius = 50 pixels). The region corresponding to the cell nucleus was manually extracted using the “Polygon selections” method. The centroids of detected fluorescent dots in each image were determined, and their trajectories were calculated using u-track (MATLAB package (Jaqaman et al. 2008); for MSD analysis) or TrackMate (Fiji package (Tinevez et al. 2017); for Spot-On analysis). The u-track trajectories were then used to calculate the displacement distributions and the mean squared displacements (MSD) of fluorescent dots in Python. The originally calculated 2D MSD values were multiplied by 1.5 (4 to 6 Dt) to obtain 3D MSD values. MSD plots with appropriate statistical analyses between various phenotypes were made on Python.
For angle distribution analysis (Iida et al. 2022; Izeddin et al. 2014), first, the consecutive tracked positions {(x0, y0), (x1, y1), ⋯, (xn, yn), ⋯} of single nucleosomes on the XY plane were converted into a set of displacement vectors ∆rn = (xn + 1 − xn, yn + 1 − yn)t. Then, the angle between two vectors ∆rn and ∆rn + 1 was calculated for every trajectory in our experiments, and the normalized polar histogram was plotted in Python. The angle distribution was normalized by 2, and the magnitudes correspond to the probability density. The color spectrum indicates each angle in a circular range from 0॰ to 360॰.
The TrackMate trajectories were submitted for the 2-state Spot-On analysis as described in (Hansen et al. 2018). The frame rate was set at 10 ms, the number of time points was set at 10, and the maximum diffusion coefficient of the bound fraction, Dbound (max), was set at 0.03 µm2/s.