Cell lines and establishment of stable cell lines
Mouse 10T1/2 (Reznikoff et al. 1973) and CIRAS-3 cells (Egan et al. 1987) (a gift from Dr. James R. Davie at the University of Manitoba) were cultured at 37°C with 5% CO2 in Dulbecco’s Modified Eagle’s medium (DMEM) (D5796-500ML, Sigma-Aldrich) supplemented with 10% fetal bovine serum (FBS) (FB-1061/500, Biosera).
To stably express H2B-Halo in 10T1/2 and CIRAS-3 cells, the transposon system was used. The constructed plasmid pPB-EF1a-IB-H2B-HaloTag (Nagashima et al. 2019) was cotransfected with pCMV-hyPBase (provided from the Sanger Institute with a materials transfer agreement) to 10T1/2 and CIRAS-3 cells with the Effectene Transfection Reagent kit (301425; QIAGEN). Transfected 10T1/2 cells were selected with blasticidin S (10 µg/ml) (029-18701, Wako). Transfected CIRAS-3 cells were selected with blasticidin S (10 µg/ml) and G418 (400 µg/ml) (16512-81, Nakalai).
Western blot
To check the expression level of H2B-HaloTag, transfected 10T1/2 and CIRAS-3 cells were lysed in Laemmli sample buffer (Laemmli 1970) supplemented with 10% 2-mercaptoethanol (133–1457; Wako) and incubated at 95°C for 10 min to denature their proteins. The cell lysates, equivalent 2x105 10T1/2 or CIRAS-3 cells per well, were subjected to SDS–polyacrylamide gel electrophoresis (PAGE) (14%) and subsequent Western blotting. For Western blot, the fractionated proteins in the gel were transferred to a polyvinylidene difluoride membrane (PVDF) (IPVH00010, Millipore) by a semidry blotter (BE-320, BIO CRAFT). After blocking with 10% or 3% skim milk (190-12865, FujiFilm), the membrane-bound proteins were probed by the anti-H2B rabbit (1:10,000 dilution; ab1790, Abcam), anti-HaloTag mouse (1:1000; G9211, Promega), anti-H3K9me3 mouse (1:1,000 dilution; provided by Prof. Hiroshi Kimura), or anti-H3K27me3 mouse (1:1000; provided by Prof. Hiroshi Kimura) antibody, followed by the appropriate secondary antibody, anti-rabbit (1:10,000 dilution; 170–6515, Bio-Rad) or anti-mouse (1:5000 dilution; 170–6516, Bio-Rad) horseradish peroxidase (HRP)–conjugated goat antibody. Chemiluminescence reactions were used (WBKLS0100, Millipore) and detected by EZ-Capture MG (AE-9300H-CSP, ATTO).
Imaging and quantification of HaloTag-labeled cells
To examine H2B-HaloTag localization in 10T1/2 and CIRAS-3 cells, cells grown on poly-l-lysine-coated (P1524-500MG, Sigma-Aldrich) coverslips (C018001, Matsunami) were treated with 5 nM HaloTag TMR ligand (8251, Promega) overnight at 37°C in 5% CO2. The cells were fixed with 1.85% FA (064–00406, Wako) on coverslips at room temperature for 15 min followed by quenching with 50 mM glycine (077–00735, Wako) for 5 min, permeabilized with 0.5% Triton X-100 (T-9284, Sigma-Aldrich) for 5 min, washed with 1x HMK (20 mM HEPES pH 7.5, 1 mM MgCl2, 100 mM KCl), and stained with 4′,6-diamidino-2- phenylindole (DAPI) (0.5 µg/ml) (10236276001, Roche) for 5 min, before being embedded in PPDI (20 mM HEPES (pH 7.4), 1 mM MgCl2, 100 mM KCl, 78% glycerol, and paraphenylene diamine (1 mg/ml) (695106-1G, Sigma-Aldrich)). Optical sectioning images were recorded with a 0.2-µm step size using a DeltaVision Personal Microscope (Applied Precision) with an Olympus PlanApo N 60x objective (NA, 1.42) and a scientific complementary metal-oxide semiconductor (sCMOS) camera. The four-color standard filter set was also equipped. DeltaVision acquisition software, Softworx, was used to project deconvolved Z-stacks to cover the whole nucleus.
Chromosome spreads
To observe and count the chromosome numbers of 10T1/2 and CIRAS-3 cells, mitotic chromosome spreads were made. Cells were treated with colcemid (0.5 µg/ml) (045-16963, Wako) for 1 h at 37°C. The following steps were performed at room temperature except for the final fixation step. After washing with phosphate-buffered saline (PBS), cells were trypsinized and resuspended in hypotonic buffer (75 mM KCl) for 8 min. Cells were then gently fixed by repeating the following treatment three times: fixation buffer (methanol: acetic acid, 3:1) for 5 min, centrifuged at 850×g for 4 min, and then resuspended in new fixation buffer. Lastly, the cells were completely fixed by treatment with 200 µl of fixation buffer at − 30°C for at least 30 min. For imaging, 10 µl of the fixed cell suspension was dropped onto coverslips after mixing and completely dried by incubating at 60°C for 30 min. Dried cells were stained with DAPI (0.5 µg/ml) for 5 min, followed by PPDI mounting. Images of mitotic chromosome spreads were obtained using DeltaVision Personal Microscopy with an Olympus PlanApoN 60x objective (NA 1.42) and an sCMOS camera. The number of chromosomes in each cell was counted manually.
Nuclear volume analysis
Nuclear volumes were measured (Iida et al. 2022) as follows. The cells were grown on poly-l-lysine-coated (P1524-500MG, Sigma-Aldrich) coverslips (C018001, Matsunami). Subsequent processes were performed at room temperature. The cells were washed with 1x PBS, fixed with 1.85% FA for 15 min followed by quenching with 50 mM glycine (077–00735, Wako) for 5 min, permeabilized with 0.5% Triton X-100 for 5 min, washed with 1x HMK, stained with DAPI (0.5 µg/ml) for 5 min, and washed with 1x HMK for 3 min three times before being embedded in PPDI. Z-stack images (every 0.4 µm, 20 to 40 sections in total) of the cells were obtained using a FLUOVIEW FV1000 confocal laser scanning microscope (OLYMPUS) equipped with an Olympus UPLANSAPO 60x W objective (NA 1.20) at room temperature. Obtained Z-stack images were loaded to Imaris (Bitplane AG, Zurich, Switzerland) and converted to Imaris 3D image files (.ims). To calculate the nuclear volume, the tool “Surface” was used. A threshold value was automatically determined to include all DAPI signals. Note that only well-isolated nuclei were recorded and analyzed.
Cell growth and cell cycle
Cells were seeded in 6-well plates (1.0×104 cells/ml). The number of proliferating cells was examined microscopically at several time points. The growth rates of 10T1/2 and CIRAS-3 cells were calculated from exponentially fitted growth curves.
Flow cytometry was performed to check the cell cycle profiles of 10T1/2 and CIRAS-3 cells. Collected cells were fixed in 70% ethanol at − 30°C for over 30 min. After fixation, cells were centrifuged at 603×g for 1 min, and the cell pellets were washed with PBS containing 1% BSA (bovine serum albumin; A9647-100G, Sigma-Aldrich). After centrifugation, cell pellets were resuspended in 600 µl of PBS containing 1% BSA, ribonuclease A (50 µg/ml) (10109169001, Sigma-Aldrich), and propidium iodide (40 µg/ml) (P4170-10MG, Sigma-Aldrich) and then further incubated for 30 min at 37°C with light protection. Lastly, cells were analyzed by BD Accuri C6 Plus (BD FACS, San Jose, CA). At least 20,000 cells were used for each analysis, and the obtained results were displayed as histograms. The percentage of cell cycle stage distributions in the G0/G1, S, and G2/M phases was analyzed by the ModFit LT software (Verity Software House, Topsham, ME).
Single-nucleosome imaging
CIRAS-3 and 10T1/2 cells stably expressing H2B-Halo were cultured on 35-mm dishes (153066, Nunc). H2B-Halo incorporated in nucleosomes was fluorescently labeled with 110 pM HaloTag TMR ligand (for 10T1/2 cells) and 100 pM HaloTag TMR ligand (for CIRAS-3 cells) for 20 min at 37°C in 5% CO2 and then washed once with 1x PBS. Cells were trypsinized and resuspended in phenol red-free DMEM (21063-029, Thermo Fisher Scientific) supplemented with 10% FBS, and then seeded on fibronectin (354008, Corning) coated glass-base dishes (3970-035, IWAKI). A live-cell chamber (INU-TIZ-F1, Tokai Hit) and digital gas mixer (GM-8000, Tokai Hit) were used to maintain cell culture conditions (37°C, 5% CO2, and humidity) during microscopy. Single nucleosomes were observed using an inverted Nikon Eclipse Ti microscope with a 100-mW Sapphire 561-nm laser (Coherent) and sCMOS ORCA-Flash 4.0 camera or ORCA-Fusion BT (Hamamatsu Photonics). Live cells labeled with H2B-Halo-TMR were excited by the 561-nm laser through an objective lens (100x PlanApo TIRF, NA 1.49; Nikon) and detected at 575 to 710 nm. An oblique illumination system with the TIRF unit (Nikon) was used to excite fluorescent nucleosome molecules within a limited thin area in the cell nucleus and reduce background noise (Fig. 2A). Sequential image frames were acquired using MetaMorph software (Molecular Devices) or NIS-Elements AR (Nikon) at a frame rate of 50 ms under continuous illumination.
Single-nucleosome tracking analysis
Image processing, single-nucleosome tracking, and single-nucleosome movement analysis were performed as previously described (Iida et al. 2022; Nozaki et al. 2023). Briefly, sequential images were converted to a 16-bit grayscale, and the background noise signals were subtracted with the rolling ball background subtraction (radius, 50 pixels) of Fiji/ImageJ (Schindelin et al. 2012). The nuclear regions in the images were manually extracted. Following this step, the position of each fluorescent dot in each image was determined and its trajectory was tracked with u-track ( (Jaqaman et al. 2008), MATLAB package). We previously calculated the SD of the 2D movement of immobilized nucleosomes per 50 ms in FA-fixed RPE-1 cells (n = 10 nucleosomes) to determine the accuracy of the position of H2B-Halo nucleosomes. The localization accuracy was 15.55 nm (Nagashima et al. 2019). For single-nucleosome movement analysis, the displacement distribution and the mean square displacement (MSD) of the fluorescent dots were calculated based on their trajectory using lab-made Python programs. The originally calculated MSD was in 2D. To obtain the 3D value, the 2D value was multiplied by 1.5 (4 to 6 Dt). Graphs and statistical analyses of the obtained single-nucleosome MSD between various conditions were performed using R. For chemical fixation, cells grown on poly-l-lysine-coated glass-based dishes were incubated with 2% FA in 1x HBSS at 37°C for 15 min and washed with 1x PBS at 4°C for at least 3 h.
Angle distribution analysis
Angle distribution analysis was performed as previously described (Iida et al. 2022; Izeddin et al. 2014). For the tracked consecutive points [(x0, y0), (x1, y1), ⋯, (xn, yn), ⋯] of a single nucleosome on the xy plane, we converted the data into a set of displacement vectors, Δrn = (xn + 1 – xn, yn + 1 – yn)t. Then, we calculated the angle between the two vectors Δrn and Δrn + 1. We carried out this procedure for all the points of each trajectory in our experiments. Lastly, we plotted the normalized polar histogram with our Python program. The angle distribution was normalized by 2π, and the values corresponded to the probability density.
Anchorage-independent cell growth analysis
Cell culture in a soft agar layer was performed as previously described (Xu et al. 2021; Zhao et al. 2015) with a 0.5% agar (010-08725, FujiFilm) base layer containing 10% FBS in DMEM with penicillin/streptomycin and a 0.35% agar growth layer containing 10% FBS in DMEM. Colonies that had three or more cells were scored at 21 days after the cell seeding.
In situ Hi-C library preparation
In situ Hi-C libraries in this study were prepared following the protocols described in (Kawaguchi and Tanaka 2023; Schloissnig et al. 2021), with some modifications. One million of 10T1/2 or CIRAS-3 cells were collected into 50 ml conical tubes with 10% FCS/DMEM and then fixed with 1% formaldehyde for 10 min at room temperature. The fixation was halted by adding ice-cold glycine (125 mM final concentration) and tubes were kept on ice for 15 min with occasional inversion for mixing. Cells were collected by centrifugation (350xg for 10 min at 4 ℃) and washed twice with ice-cold PBS (350xg for 10 min at 4 ℃). Cell pellets were then resuspended in ice-cold lysis buffer (10 mM Tris-HCl pH 7.5, 10 mM NaCl, 0.2% NP-40, protease inhibitor-EDTA free (P8340, SIGMA)) and incubated on ice for 30 min with occasional agitation. Afterward, the nuclei were collected by centrifugation (300xg for 10 min at 4 ℃) and washed with 500 µl of 1.25x NEBuffer3 at the same speed. The nuclei were resuspended in 1 ml of 1.25x NEBuffer3 and transferred to DNA low-binding tubes (500 µl each). SDS (0.4% final concentration) was added and the tubes were incubated with agitation (950 rpm) for 2 h at 37 ℃. Triton X-100 (3.3% final concentration) was added to quench SDS, and the nuclei were incubated with agitation (950 rpm) for 2 h at 37 ℃. Nuclei were pelleted with centrifugation (100xg for 1 min at room temperature) and the pellets were resuspended in 450 µl of 1x DpnII digestion buffer (R0543M, NEB). DpnII digestion (NEB; 100 U per 0.5 million cells) was performed with agitation (950 rpm) for 2 h at 37 ℃. The enzyme was then inactivated by heating the tubes to 65 ℃ for 20 min. Overhang filling was carried out using biotin-14-dATP, dCTP, dGTP, and dTTP, along with Klenow (NEB; 50 U per 0.5 million cells) for 1 h at 37 ℃ with repeated agitation (repeated cycles of 700 rpm for 10 sec and rest for 30 sec for 1 h in a thermal cycler). Ligation was performed overnight at 18 ℃ with 2000 U of T4 DNA ligase. De-crosslinking with Protease K in SDS buffer was done overnight at 65°C, followed by RNase A treatment for 2 h at 37 ℃. Subsequently, DNA extraction was performed with phenol/chloroform and EtOH precipitation steps, and the ligated gDNA products were eluted in 20 µl of TE buffer. The yield was determined using a Qubit3 device, and 5 µg of biotinylated DNA was used for the library preparation. To remove excess biotin, non-ligated fragment ends were then treated with T4 DNA polymerase for 30 min at 37°C, and the reaction was stopped by adding EDTA (10 mM final concentration). DNA was sonicated using Covaris S220 system to generate DNA fragments with the size peak around 400 bp (Covaris settings were duty factor: 10%; peak incident power: 5; cycles per burst: 200; time: 65 sec). After end repair with T4 DNA polymerase, Large (Klenow) Fragment, and T4 DNA polynucleotide kinase (all from NEB) in the presence of dNTPs in T4 DNA ligation buffer (for 30 min at room temperature), the DNA was purified with QIAGEN mini purification kit. A double-sized selection using AMPure XP beads (Beckman) was performed (see (Kawaguchi and Tanaka 2023; Schloissnig et al. 2021)). Firstly, the ratio of the AMPure XP bead volume to DNA sample volume was adjusted to 0.6:1, and then the second step was performed with the final ratio of 0.9:1. Following two washes with 80% EtOH, the DNA was eluted in Elution buffer (QIAGEN). From DNA pools, Hi-C ligation products were isolated using pre-washed MyOne Streptavidin C1 Dynabeads (# 65001, Life Technologies) on a magnet stand in binding buffer (5 mM Tris pH 8.0, 0.5 mM EDTA, 1 M NaCl) for 30 min at room temperature. dA-tailing was carried out on the beads. dATP was added with Klenow exo- (for 1 h at 37°C), then the enzyme was heat-inactivated (for 20 min at 65°C). After two washes in binding buffer and one wash in T4 DNA ligation buffer, Illumina adapters were ligated onto Hi-C ligation products bound to streptavidin beads in T4 DNA ligase while slowly rotating (for 2 h at room temperature). After washing twice with wash buffer (5 mM Tris pH 8.0, 0.5 mM EDTA, 1 M NaCl, 0.05% Tween-20) and once with binding buffer, the Hi-C product bound beads were resuspended in a final volume of 20 µl 1x NEBuffer2. Captured Hi-C DNA was amplified with PCR amplification cycles using NEBNext® High-Fidelity 2X PCR Master Mix (# M0541S, NEB) for 7 cycles. The final Hi-C libraries were purified with 1x volume of AMPure XP beads. The concentration of the Hi-C library was determined using Tapestation (Agilent), and the Hi-C libraries were sequenced with the PE150 setting on NovaSeq 6000. We repeated the above experiment for the reproducibility and Hi-C sequence data are available in DDBJ database as DRA017367.
Hi-C data analysis
Hi-C data were analyzed using the nf-core (Ewels et al. 2020) Hi-C pipeline (nf-core/hic: 2.0.0) (Servant et al. 2023), executed within a Docker configuration profile. The Mus musculus genome mm10 (http://www.ncbi.nlm.nih.gov/assembly/GCF_000001635.20), retrieved from Illumina’s iGenomes, was used for the reference genome. HiC-Pro (Servant et al. 2015) was employed, which performs sequence reads mapping with Bowtie2 (Langmead and Salzberg 2012) and identifies valid chromosomal interactions. Subsequently, Cooler (Abdennur and Mirny 2020) was used for constructing and normalizing contact maps through its balancing algorithm, and compartment callings were conducted with cooltools (Venev et al. 2022). To detect chromosomal translocations and breakpoints from the Hi-C data, HiNT (Wang et al. 2020) software (version 2.1.9) was run in the HiNT-TL mode with the default parameter settings. The visualization of contact maps and A/B compartment tracks in Fig. 6 were generated using CoolBox (v0.3.9) (Xu et al. 2021). To visualize the strength of interactions between compartments, the saddle plot analysis was performed using the R-package GENOVA (v1.0) (van der Weide et al. 2021).