Next Article in Journal
Identification of Reference Genes for Quantitative Gene Expression Studies in Three Tissues of Japanese Quail
Next Article in Special Issue
Origin Firing Regulations to Control Genome Replication Timing
Previous Article in Journal
Cytological Evaluations of Advanced Generations of Interspecific Hybrids between Allium cepa and Allium fistulosum Showing Resistance to Stemphylium vesicarium
Previous Article in Special Issue
On the Interplay of the DNA Replication Program and the Intra-S Phase Checkpoint Pathway
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 3
10.3390/genes10030196
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Cohesin-Mediated Genome Architecture Does Not Define DNA Replication Timing Domains

by
Phoebe Oldach
and
Conrad A. Nieduszynski
*
Sir William Dunn School of Pathology, University of Oxford, South Parks Road, Oxford OX1 3RE, UK
*
Author to whom correspondence should be addressed.
Genes 2019, 10(3), 196; https://doi.org/10.3390/genes10030196
Submission received: 14 January 2019 / Revised: 7 February 2019 / Accepted: 19 February 2019 / Published: 4 March 2019
(This article belongs to the Special Issue DNA Replication Timing: Where, When, How and Why?)
Browse Figures
Versions Notes

Abstract

:
3D genome organization is strongly predictive of DNA replication timing in mammalian cells. This work tested the extent to which loop-based genome architecture acts as a regulatory unit of replication timing by using an auxin-inducible system for acute cohesin ablation. Cohesin ablation in a population of cells in asynchronous culture was shown not to disrupt patterns of replication timing as assayed by replication sequencing (RepliSeq) or BrdU-focus microscopy. Furthermore, cohesin ablation prior to S phase entry in synchronized cells was similarly shown to not impact replication timing patterns. These results suggest that cohesin-mediated genome architecture is not required for the execution of replication timing patterns in S phase, nor for the establishment of replication timing domains in G1.

1. Introduction

The double helix structure of DNA immediately lends itself to a mechanism for genome replication so simple and so elegant that the question of faithful transmission of hereditary information, to some, seemed solved as soon as it was posed [1]. However, behind the elegance of the underlying mechanism, the process of genome replication is in fact a complex, high-stakes, and intricately regulated cellular undertaking. A significant factor necessitating regulation is the sheer scale of the process: in a human cell, 12 billion bases are replicated in about eight hours, with replication initiating from somewhere between 12,000 and 250,000 replication origins spread across 23 chromosome pairs [2]. Intriguingly, not all of these origins fire in any given S phase. Of those that do fire, not all fire at the same time. Instead, clusters of contiguous origins appear to fire in synchrony, leading to zones of hundreds of kilobases to megabases of DNA making up distinct replication timing domains (RD) [3]. The patterns of replication timing domains, driven by selective origin firings, lead to nuanced genome-wide replication timing profiles distinct to cell type and differentiation status [4,5,6,7].
Work to characterize the cellular determinants of replication timing has uncovered a strong correlation to genome organization. This has been assayed through both cytological studies of replication foci over the course of S phase [8,9,10,11] and chromosome conformation capture studies [12,13,14,15]. Genome structure appears primarily to emerge from two fundamental organizing principles: (1) a loop-based organization dependent upon the proteins cohesin and CTCF results in domains of ≤1 Mb that are defined by enriched internal interactions and known as topologically associating domains (TADs) [16,17,18], and (2) compartmentalization into like chromatin states results in alternating 1–10 Mb domains of active or inactive chromatin [16,19,20]. Of these different classes of genome structure, there is evidence that replication timing domains align well with compartments: early-replicating loci are significantly enriched for marks of open and active chromatin, or compartment A, whereas late-replicating loci share marks with compartment B [3,21].
Beyond the concordance between replication timing domains and genome compartments, the smaller TADs offer a compelling possibility for a unit of replication timing domain regulation. Their boundaries were found to map almost one-to-one to replication timing domain boundaries [6], and cohesin was found to be enriched at origins as well as at TAD boundaries [22]. A link between TAD-based genome structure and replication origin function was further supported by the finding that both are established concomitantly in early G1 [23,24] and change correspondingly over the course of cellular differentiation [25]. The loop-based mechanism by which cohesin is hypothesized to drive genome organization is furthermore in agreement with work showing a correlation between chromatin loop size and interorigin distance: the relationship held as reprogramming of differentiated cells to embryonic cells resulted in both a shortening of interorigin distances and chromatin loop domains [26]. Together, this abundance of correlative evidence has made a convincing case for a “replication domain model,” in which the stable and spatially distinct structural unit of TADs gives rise to the stable functional unit of replication timing domains [6,27].
Due to cohesin’s critical role in chromosome segregation, studies of the impact of cohesin on genome architecture and replication timing have until recently been limited by the inability to carry out a clean cohesin knock-out or complete depletion over the course of a single cell cycle. A study which used cohesin knockdown by small interfering RNA (siRNA) showed a reduction in active origins and an increase in interorigin distance in response to depleted cohesin levels [22]. However, this work may have been limited by the slower approach of regulating cohesin at the transcriptional level, and furthermore did not explore an effect on replication timing domains. The development of the auxin-inducible degron system [28,29] has made this sort of research much more tractable, allowing for a rapid and near-complete depletion of cohesin. Chromosome conformation studies of cohesin ablation over a single cell cycle via a degron tag on cohesin’s Scc1 subunit have shown a rapid loss of TADs, while compartment structure remained intact [30].
In this study, an auxin-inducible degron tag on the Scc1 subunit of cohesin is exploited to test whether the observation of a strong correlation between cohesin-mediated genome architecture and replication timing is substantiated as a causative relationship.

2. Materials and Methods

2.1. Cell Culture

HCT116 cells [29] were grown in MyCoys 5A media (Sigma M9309, St. Louis, MO, USA) with 10% FBS (Gibco 10500, Thermo Fisher Scientific, Waltham, MA, USA), 100 U/mL penicillin, and 100 µg/mL streptomycin (Gibco 15140-122, USA) at 37 °C with 5% CO2. Cells were treated with 500 µM indole-3-acetic acid (IAA, Sigma I5148, USA) to induce loss of mAID-tagged Scc1.

2.2. Cell Sorting

Cell sorting was performed as described previously [31]. In brief, asynchronous cells were treated or not treated with IAA (500 µM) for two hours and then treated with BrdU (100 µM) for two hours. Cells were washed, trypsinized and harvested, and fixed with 75% ethanol. Fixed cells were prepared for fluorescence activated cell sorting (FACS) with RNase A (250 µg/mL) and propidium iodide (PI, 50 µg/mL). Cells were sorted on the basis of PI stain into four bins spanning S phase to collect a minimum of 2.5 × 105 cells per bin.

2.3. Synchronization

Cells were treated with 2.5 mM thymidine for 24 h and released by washing into fresh media. Three hours after release, cells were treated with 100 ng/mL nocodazole for eight hours, at which point mitotic cells were collected via shake-off. Cells were plated into fresh media without nocodazole, and 30 min later were treated with 500 µM auxin or left untreated. Cells were collected at four, seven, or ten hours post release from nocodazole, and BrdU was added to 100 µM two hours before the collection of each cell population (i.e., at two, five, and eight hours).

2.4. Flow Cytometry

To confirm the loss of Scc1 upon auxin treatment, cells were fixed with 75% ethanol and prepared for flow cytometry assessment of mClover signal.
To confirm the cell cycle staging of BrdU-pulsed synchronized samples, cells were fixed and permeabilized with 75% ethanol. DNA was denatured via incubation in 0.2 mg/mL pepsin in 2 M HCl for 20 min. Cells were washed twice with wash solution composed of: 1% bovine serum albumin (BSA)/0.5% (v/v) Tween 20 in Dulbecco’s phosphate buffered saline (dPBS). Cells were then incubated in wash solution with anti-BrdU antibody (BD347580) at 1:200 for one hour. Cells were washed twice, then incubated with AlexaFluor 647 Goat anti-Mouse (Invitrogen A32728, Waltham, MA, USA) antibody at 1:100 in wash solution for one hour, followed by two washes. Cells were finally treated with RNase A (250 µg/mL) and PI (50 µg/mL) before acquisition on a BD LSRFortessa.

2.5. RepliSeq

DNA was extracted from each sample using phenol/chloroform followed by an isopropanol precipitation. The extracted DNA was sonicated to an average fragment size of 300 bases. DNA was purified using SPRI beads, repaired with the NEBNext End Prep kit (NEB) followed by ligation of sequencing adaptors using the NEBNext UltraII ligation kit. DNA was purified using AMPure XP beads (Beckman Coulter). BrdU pulldown was carried out as per Peace et al. (2016) [32]. Briefly, DNA was heat-denatured and snap-cooled, then incubated overnight with anti-BrdU antibody (BD347580). The sample was then incubated for one hour with Protein G Dynabeads (ThermoFisher), washed three times with IP buffer (PBS/0.0625% (v/v) Triton X-100), and once with Tris-EDTA (TE), and then eluted into TE/1% (w/v) sodium dodecyl sulfate (SDS) at 65 °C. Eluted DNA was cleaned using AMPure XP beads, and then amplified with Illumina indexing primers using 16 cycles of the NEBNext Ultra II kit. Amplified libraries were cleaned, quantified using NEBNext Library Quant kit, and checked for fragment length using Tapestation.
Libraries were sequenced using 75 cycles on an Illumina NextSeq 500 for a minimum of ten million reads per library.
Sequenced HCT116 raw fastq files and bigwig files reporting processed log2 replication timing scores are available from the NCBI GEO database (accession number GSE124025). The processed replication timing scores can also be visualized via a UCSC genome browser hub: https://ln1.path.ox.ac.uk/groups/nieduszynski/Oldach2019/Oldach2019_hub.txt.

2.6. Immunofluorescence

BrdU Replication Foci: Cells were seeded onto coverslips (thickness #1.5) and grown to 80% confluency. Media was replaced and cells were treated with IAA (500 µM) or left untreated for 3.5 h at the end of which all cells were pulsed with BrdU (50 µM) for 30 min.
Upon completion of the BrdU pulse, coverslips were washed with PBS and fixed with 3% (w/v) formaldehyde/dPBS, washed, and then permeabilized with 0.2% (v/v) Triton X-100 in dPBS. DNA was denatured with 4 M HCl treatment for ten minutes, after which coverslips were washed and blocked with 5% (w/v) BSA in PBS with 0.1% Tween 20 (PBST). Cells were incubated for one hour with anti-BrdU antibody (BD347580) 1:200 in block solution. After three washes, coverslips were incubated for one hour in the dark with AlexaFluor Rabbit anti Mouse 568 (Invitrogen A11061) at 1:1000 in 5% BSA. Coverslips were washed and nuclei were stained with DAPI. Coverslips were finally mounted with VectaShield H1000. DAPI signal was used to identify cells for image acquisition and BrdU-positive cells were classified by eye into five S phase stages by researchers blinded to the treatment status of the images [9].
Scc1 IF: Asynchronous cultures of HCT116-Scc1-mAID cells were seeded on coverslips and treated for two hours with IAA (500 µM) or left as controls. At the end of the treatment, coverslips were washed with PBS then fixed with 3% formaldehyde. Cells were permeabilized with 0.2% Triton X-100 in PBS, washed with PBST, then blocked for one hour in 5% BSA in PBST. Coverslips were incubated for one hour in anti-Scc1 (Abcam 154769) at 1:500 in block. After three washes, coverslips were incubated for one hour in Alexafluor 647 Goat anti Rabbit (Invitrogen A21245) 1:1000 in block. Coverslips were washed, stained with DAPI, and mounted in Vectashield.

2.7. RepliSeq Computational Analysis

Reads were aligned to the human genome (hg38) with STAR [33]. Uniquely mapping reads were binned into 1 kb windows and these windows were filtered to remove bins with signal spikes (Z score > 2.5). Reads were then binned into 50 kb windows. Early over late replication timing (E/L RT) was calculated as the log2(S phase bin 1/S phase bin 3) reads for the asynchronous sample and log2(4h/10h) for synchronized replicates. Quantile normalization and Loess smoothing with a span of 300 kb were applied to E/L RT values. RepliSeq data for wildtype HCT116 cells (Accessions: 4DNFI5BZJXDE, 4DNFIC4VUF86) and IMR90 cells (Accessions: 4DNFILOYZWEM, 4DNFIKQDQCNB) were gathered from the 4D Nucleome data portal [34], and a Pearson correlation was used to quantify the degree of similarity between the profiles.
A Z score was calculated based on the difference between treated and control replication timings to identify significantly changed bins. A false discovery rate (FDR) adjustment was applied to p values to correct for multiple testing.
Weighted RT values were calculated as the weighted average of the four S phase bins for asynchronous samples, or three S phase bins for synchronized samples. Publicly accessible data for histone marks, CTCF binding sites, and other genome characteristics in HCT116 cells was extracted from ENCODE [35] and the UCSC Table Browser [36] (Supplemental Table S1). ProSeq data from Rao et al. (2017) was used to identify genes with differential expression in response to cohesin loss upon auxin treatment [30]. Briefly, significantly changed genes were classified as those with an absolute fold change greater than two, an adjusted p value less than 0.05, and at least 0.5 reads per kilobase million (RPKM). These cutoffs yielded 74 genes that were induced upon auxin treatment and 16 that were downregulated. The genomic coordinates for these genes was extracted using BioMart [37]. BEDOPS was used to calculate the mean weighted RT for auxin-treated or control samples across annotated genomic ranges [38].
A mappability score from Karimzadeh et al. [39] was used to mask unmappable regions of the genome. Analysis was carried out in R [40] and data were visualized using ggplot2 [41] and the UCSC Genome Browser [42].
The 3D Genome Browser [43] was used to visualize and extract TAD calls from control and auxin-treated HCT116-Scc1-mAID cell HiC data from Rao et al. (2017), yielding a list of 2307 domains in untreated cells and 883 domains in auxin-treated cells. Juicer [44] was used to calculate compartment boundaries from untreated HCT116-Scc1-mAID cell HiC data. TAD boundaries that fell within 500 kb of a compartment boundary were considered separately to control for the effect of compartment status on replication timing. Compartment-independent TADs were classed as “Scc1-dependent” if they were not present (within 500 kb) in the contact domains of auxin-treated HCT116-Scc1-mAID HiC data, or “Scc1-invariant” if they remained in the contact domains of auxin-treated cells.

3. Results

3.1. Efficient and Rapid Degradation of Endogenous Scc1

It has been previously shown that degradation of the Scc1 subunit of the cohesin complex leads to an inability of the complex to associate with DNA [30]. Here, we used an HCT116 cell line with the endogenous Scc1 tagged with mini-AID and mClover for rapid, auxin-inducible degradation of the cohesin subunit and thus destabilization of the cohesin complex [29]. To confirm loss of the Scc1 fusion protein after two hours of auxin treatment, we analyzed mClover signal in treated and control cells using microscopy (Figure 1A) and flow cytometry (Figure 1B). These two methods yielded comparable values: after two hours of auxin treatment, 73% of cells were mClover-negative using flow cytometry (n = 40,649), and 78% were mClover-negative using microscopy (n = 54) (Figure 1B,C). Cells that retained residual mClover signal after auxin treatment had a notably reduced intensity of signal relative to untreated levels (Figure 1A). To further confirm that loss of mClover related to loss of Scc1, immunofluorescence (IF) was carried out for Scc1. A total of 81% of cells imaged were negative for the Scc1 stain after two hours of auxin treatment and those cells that retained any Scc1 signal after treatment were markedly reduced in signal intensity (Figure 1A,C). These results agree with previous work that has demonstrated a loss of Scc1 within an hour of auxin treatment in the same cell line [29,30].

3.2. Acute Loss of Cohesin in S Phase Does Not Perturb Replication Patterns

To analyze the impact of cohesin loss on DNA replication timing, we have used RepliSeq on asynchronously growing auxin-treated and untreated control cells [31,32]. Asynchronously growing cells were treated with auxin for two hours or left as untreated controls, pulsed with BrdU for a further two hours, and then FACS enriched on DNA content into four S phase fractions. DNA was extracted from each fraction, then nascent DNA was immunoprecipitated and subjected to high-throughput sequencing. The proportion of nascent DNA in each fraction served as a measure of the relative replication timing across the genome. In the mock-treated samples, replication timing was found to be in good agreement with previously reported HCT116 cell replication timing data (Supplemental Figure S1) [45]. The overall profile of replication timing did not change in response to auxin-induced Scc1-ablation (Figure 2A). Quantitatively, the correlation between the control and treated profiles was comparable to the correlation between the control and a previously-reported HCT116 replication timing profile (Figure 2C). To facilitate identification of regions with changed replication timing, the difference of replication timing in control versus auxin-treated samples was calculated for each genomic locus. Replication timing was only significantly changed (FDR-adjusted p < 0.05) in 0.34% of mappable genomic bins (189 bins out of 56,034 bins with a mappability score > 0.5) [39]. Significant differences were enriched in bins with low mappability and dispersed across the genome (median distance between significant bins: 925,000 bases), suggesting that these differences were a consequence of experimental noise rather than domains of biological relevance (Figure 2B). Furthermore, loci with significantly-changed replication timing were not enriched at TAD boundaries, as might be expected if these changes resulted from loss of TAD architecture (Figure 2D).
As genome-wide analysis did not demonstrate widespread perturbations to the replication timing profile in response to the loss of cohesin, genomic regions marked by features hypothesized to show a change were specifically assessed. These included genomic regions with features characteristic of early replication (super enhancers, DNase hotspots, CpG islands, H3K4 methylated regions, H3K9 and H3K27 acetylated regions), late replication (H3K9me3, polycomb binding sites), CTCF binding sites, or genes with expression changes in response to auxin-induced loss of Scc1. None of these sets of sites showed a differential response (Figure 2E, Supplemental Figure S6).
To investigate the impact of cohesin ablation on the cytological patterns of replication timing, auxin-treated or control BrdU-pulsed Scc1-mAID-mClover cells were classed into five S phase categories based on their BrdU foci patterning (Figure 3A) [9]. The treatment with auxin and resultant loss of cohesin did not alter the abundance of cells across each of the five classes (Figure 3B). In combination with the finding that auxin-induced loss of Scc1 does not change progression through S phase (Supplemental Figure S3), this suggests that auxin treatment did not change the patterns of replication characteristic to different stages of S phase.

3.3. Loss of Cohesin from Early G1 Does Not Perturb Replication Timing

To assess whether cohesin is required for the establishment, rather than execution, of replication timing, it was necessary to ablate cohesin from cells before the onset of S phase. This also diminished the potential for artifacts arising from the impact of cohesin ablation on chromosome segregation. Cells were synchronized with a 24-h thymidine block and released for three hours, followed by nocodazole treatment for eight hours. Cells were released from nocodazole by mitotic shake-off, and 30 min after the shake-off, auxin was introduced for treated cells. Both auxin-treated and untreated cells showed similar cell cycle progression following release (Supplemental Figure S3).
Synchronized auxin-treated or untreated control cells were pulsed with BrdU at two, five, and eight hours post-release. Two hours after the addition of BrdU to each synchronized timepoint, the cells were harvested and DNA was extracted. BrdU-labeled nascent DNA was immunoprecipitated and subjected to high-throughput sequencing. The overall replication timing profiles for synchronized cells with cohesin depleted from early G1 showed the same pattern to those of untreated synchronized cells (Figure 4A). The variability between treated and control cells was comparable to that between control cells and an independent control HCT116 dataset (Figure 4C). A total of 219 bins showed significant changes in response to cohesin ablation (out of 56,034 50 kb bins with a map score > 0.5), and these were spread across the genome (median distance between bins 700,000 bases) and enriched in unmappable regions such as centromeres (Figure 4B). Just 14 of these bins were shared with the significant bins from the asynchronous sample, and in four of these bins, the direction of the change was different, which is suggestive of noisy or unmappable loci rather than a consistent biological effect (Table 1).
As with the replication timing data from asynchronous cells, no specific genomic features showed a particular change in replication timing (Figure 4E).

3.4. Impact of Genome Organization on Replication Timing

The lack of response of replication timing domains to cohesin ablation led to the hypothesis that compartment boundaries may be a more robust delineator of replication timing than TAD boundaries. To assess this possibility, TAD calls from the control and auxin-treated cells of Rao et al. (2017) [30] were collected from the 3D Genome Browser [43]. Upon cohesin ablation, Rao et al. saw a pervasive loss of TAD structure while higher level compartment structure remained (Figure 5A) [30]. In keeping with this, of the 2026 compartment-independent TAD boundaries (i.e., not within 500 kb of a compartment boundary) mapped in the untreated HCT116 cells, only 887 remained following auxin treatment. TAD boundaries that disappeared in the auxin-treated dataset were classed as “Scc1-dependent boundaries” (n = 1139), while those that were shared (within 500 kb) between control and auxin-treated datasets were classed as “invariant” (n = 887), and any that fell within 500 kb of an A/B compartment domain boundary were analyzed separately (n = 1378).
Replication timing values for auxin-treated and control cells were assessed along 10 Mb windows to either side of Scc1-invariant and dependent boundaries. Both control and auxin-treated datasets showed a similar relationship between replication timing domains and TAD boundaries: the loss of TAD organization did not lead to changes in replication timing across TAD boundaries (Figure 5B). Furthermore, Scc1-invariant and compartment boundaries were found to show a stronger relationship to replication timing domains than Scc1-dependent boundaries, as evidenced by a sharper change in replication timing at the boundary (Figure 5B).

4. Discussion & Conclusions

The hypothesis that TADs act as a fundamental unit of organization for replication timing has gained widespread support due to the near one-to-one mapping of replication timing domains and TADs [6], the concurrent establishment of both domains in early G1 [24], and their coordinated changes over the course of cellular differentiation [13,25]. Using the Scc1-degron system, we were able to explicitly test whether the loss of cohesin resulted in a perturbation to replication timing patterns. We found that neither acute loss of cohesin in S phase nor prolonged loss of cohesin from early G1 resulted in significant changes to replication timing patterns. Thus, cohesin is not required for the execution nor establishment of replication timing domains. This is in good agreement with recent work showing that depletion of TAD boundary protein CTCF also does not result in significant genome-wide perturbations to replication timing patterns [46].
The RepliSeq technique specifically informs about replication timing domains and does not shed light on specific origin usage. Thus, it remains feasible that cohesin-dependent structure does in fact play a role in the specific choice of sites of origin firing within a domain, while the timing of higher-order domains remains independent of TAD structure. This could reconcile the apparent lack of domain timing change to the changes in interorigin distance seen previously in cohesin knockdown experimentation [22].
The lack of response of replication timing domains to cohesin loss and the stronger relationship seen between cohesin-invariant and compartment boundaries than cohesin-dependent boundaries to changes in replication timing suggests that replication timing domains may be regulated by the same organizing principles that drive A/B compartmentalization rather than TADs. While this could in part be explained by the fact that neighboring TADs can share similar chromatin states, and thus may be likely to replicate at similar times regardless of being members of distinct regulatory modules, there is increasing evidence to support a model in which compartmentalization-based structure rather than loop-based structure drives replication timing. Recently-discovered early replication control elements (ERCE) have been shown to drive robust CTCF-independent interactions in a manner akin to compartmentalization, and their deletion leads to changes in replication timing [46]. Further work characterizing the cellular mechanisms driving genome compartmentalization will open pathways to perturb this system and elucidate the causal relationship between replication timing and genome topology.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/3/196/s1, Supplemental Figure S1. Newly generated RepliSeq data agrees with previously published RepliSeq data for HCT116 cells; Supplemental Figure S2. Sample preparation for asynchronous RepliSeq experiment; Supplemental Figure S3. Auxin treatment post release from nocodazole by mitotic shake-off does not perturb synchronous progression through S phase; Supplemental Figure S4. Sample preparation for synchronous RepliSeq experiment; Supplemental Figure S5. Additional plots of control versus auxin-treated replication timing at genomic loci with chromatin states or protein binding sites associated with characteristically early or late replication timing; Supplemental Figure S6. Loci with significant changes to replication timing, in asynchronous, synchronized, or both experiments; Supplemental Table S1. Source for genomic features annotated in HCT116 cells.

Author Contributions

Conceptualization, P.O. and C.A.N.; investigation, P.O.; writing—original draft preparation, P.O.; writing—review and editing, C.A.N.; supervision, C.A.N.

Funding

This work was supported by Biotechnology and Biological Sciences Research Council grants BB/N016858/1 and Wellcome Trust Investigator Award 110064/Z/15/Z. PFO is generously supported by a Wellcome Trust studentship award 203727/Z/16/Z.

Acknowledgments

The authors are very grateful to Masato Kanemaki for sharing the Scc1-AC Tir1 HCT116 cell line. We thank David Gilbert for sharing unpublished results. We thank Line Eriksen and Michal Maj of the Dunn School Flow Cytometry Suite for facilitating FACS. We thank Amanda Williams and Rebecca Busby for their support at the Zoology Sequencing Facility. We would in particular like to thank the Nieduszynski lab for helpful discussions and substantial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Watson, J.D.; Crick, F.H.C. Molecular Structure of Nucleic Acids: A Structure for Deoxyribose Nucleic Acid. Nature 1953, 171, 737–738. [Google Scholar] [CrossRef] [PubMed]
  2. Petryk, N.; Kahli, M.; D’Aubenton-Carafa, Y.; Jaszczyszyn, Y.; Shen, Y.; Silvain, M.; Thermes, C.; Chen, C.-L.; Hyrien, O. Replication landscape of the human genome. Nat. Commun. 2015, 7, 10208. [Google Scholar] [CrossRef] [PubMed]
  3. Fu, H.; Baris, A.; Aladjem, M.I. Replication timing and nuclear structure. Curr. Opin. Cell Biol. 2018, 52, 43–50. [Google Scholar] [CrossRef] [PubMed]
  4. Hadjadj, D.; Denecker, T.; Maric, C.; Fauchereau, F.; Baldacci, G.; Cadoret, J.-C. Characterization of the replication timing program of 6 human model cell lines. Genomics Data 2016, 9, 113–117. [Google Scholar] [CrossRef] [PubMed]
  5. Ryba, T.; Hiratani, I.; Sasaki, T.; Battaglia, D.; Kulik, M.; Zhang, J.; Dalton, S.; Gilbert, D.M. Replication Timing: A Fingerprint for Cell Identity and Pluripotency. PLoS Comput. Biol. 2011, 7, e1002225. [Google Scholar] [CrossRef] [PubMed]
  6. Pope, B.D.; Ryba, T.; Dileep, V.; Yue, F.; Wu, W.; Denas, O.; Vera, D.L.; Wang, Y.; Hansen, R.S.; Canfield, T.K.; et al. Topologically associating domains are stable units of replication-timing regulation. Nature 2014, 515, 402–405. [Google Scholar] [CrossRef] [PubMed]
  7. Hansen, R.S.; Thomas, S.; Sandstrom, R.; Canfield, T.K.; Thurman, R.E.; Weaver, M.; Dorschner, M.O.; Gartler, S.M.; Stamatoyannopoulos, J.A. Sequencing newly replicated DNA reveals widespread plasticity in human replication timing. Proc. Natl. Acad. Sci. USA 2010, 107, 139–144. [Google Scholar] [CrossRef] [PubMed]
  8. Visser, A.E.; Eils, R.; Jauch, A.; Little, G.; Bakker, P.J.M.; Cremer, T.; Aten, J.A. Spatial Distributions of Early and Late Replicating Chromatin in Interphase Chromosome Territories. Exp. Cell Res. 1998, 243, 398–407. [Google Scholar] [CrossRef] [PubMed]
  9. Dimitrova, D.S.; Berezney, R. The spatio-temporal organization of DNA replication sites is identical in primary, immortalized and transformed mammalian cells. J. Cell Sci. 2002, 115, 4037–4051. [Google Scholar] [CrossRef] [PubMed]
  10. Jackson, D.A.; Pombo, A. Replicon Clusters Are Stable Units of Chromosome Structure: Evidence that Nuclear Organization Contributes to the Efficient Activation and Propagation of S Phase in Human Cells. J. Cell Biol. 1998, 140, 1285–1295. [Google Scholar] [CrossRef] [PubMed]
  11. Nakayasu, H.; Berezney, R. Mapping replicational sites in the eucaryotic cell nucleus. J. Cell Biol. 1989, 108, 1–11. [Google Scholar] [CrossRef] [PubMed]
  12. Lieberman-Aiden, E.; van Berkum, N.L.; Williams, L.; Imakaev, M.; Ragoczy, T.; Telling, A.; Amit, I.; Lajoie, B.R.; Sabo, P.J.; Dorschner, M.O.; et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science 2009, 326, 289–293. [Google Scholar] [CrossRef] [PubMed]
  13. Ryba, T.; Hiratani, I.; Lu, J.; Itoh, M.; Kulik, M.; Zhang, J.; Schulz, T.C.; Robins, A.J.; Dalton, S.; Gilbert, D.M. Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types. Genome Res. 2010, 20, 761–770. [Google Scholar] [CrossRef] [PubMed]
  14. Moindrot, B.; Audit, B.; Klous, P.; Baker, A.; Thermes, C.; De Laat, W.; Bouvet, P.; Mongelard, F.; Arneodo, A. 3D chromatin conformation correlates with replication timing and is conserved in resting cells. Nucleic Acids Res. 2012, 40, 9470–9481. [Google Scholar] [CrossRef] [PubMed]
  15. Yaffe, E.; Farkash-Amar, S.; Polten, A.; Yakhini, Z.; Tanay, A.; Simon, I. Comparative Analysis of DNA Replication Timing Reveals Conserved Large-Scale Chromosomal Architecture. PLoS Genet. 2010, 6, e1001011. [Google Scholar] [CrossRef] [PubMed]
  16. Fudenberg, G.; Imakaev, M.; Lu, C.; Goloborodko, A.; Abdennur, N.; Mirny, L.A. Formation of Chromosomal Domains by Loop Extrusion. Cell Rep. 2016, 15, 2038–2049. [Google Scholar] [CrossRef] [PubMed]
  17. Nasmyth, K. Disseminating the Genome: Joining, Resolving, and Separating Sister Chromatids During Mitosis and Meiosis. Annu. Rev. Genet. 2001, 35, 673–745. [Google Scholar] [CrossRef] [PubMed]
  18. Wutz, G.; Várnai, C.; Nagasaka, K.; Cisneros, D.A.; Stocsits, R.R.; Tang, W.; Schoenfelder, S.; Jessberger, G.; Muhar, M.; Hossain, M.J.; et al. Topologically associating domains and chromatin loops depend on cohesin and are regulated by CTCF, WAPL, and PDS5 proteins. EMBO J. 2017, 36, 3573–3599. [Google Scholar] [CrossRef] [PubMed]
  19. Nuebler, J.; Fudenberg, G.; Imakaev, M.; Abdennur, N.; Mirny, L.A. Chromatin organization by an interplay of loop extrusion and compartmental segregation. Proc. Natl. Acad. Sci. USA 2018, 115, E6697–E6706. [Google Scholar] [CrossRef] [PubMed]
  20. Schwarzer, W.; Abdennur, N.; Goloborodko, A.; Pekowska, A.; Fudenberg, G.; Loe-Mie, Y.; Fonseca, N.A.; Huber, W.; Haering, C.; Mirny, L.; et al. Two independent modes of chromatin organization revealed by cohesin removal. Nature 2017, 551, 51–56. [Google Scholar] [CrossRef] [PubMed]
  21. Lubelsky, Y.; Prinz, J.A.; DeNapoli, L.; Li, Y.; Belsky, J.A.; MacAlpine, D.M. DNA replication and transcription programs respond to the same chromatin cues. Genome Res. 2014, 24, 1102–1114. [Google Scholar] [CrossRef] [PubMed]
  22. Guillou, E.; Ibarra, A.; Coulon, V.; Casado-Vela, J.; Rico, D.; Casal, I.; Schwob, E.; Losada, A.; Méndez, J. Cohesin organizes chromatin loops at DNA replication factories. Genes Dev. 2010, 24, 2812–2822. [Google Scholar] [CrossRef] [PubMed]
  23. Dimitrova, D.S.; Gilbert, D.M. The Spatial Position and Replication Timing of Chromosomal Domains are Both Established in Early G1 Phase. Mol. Cell 1999, 4, 983–993. [Google Scholar] [CrossRef]
  24. Dileep, V.; Ay, F.; Sima, J.; Vera, D.L.; Noble, W.S.; Gilbert, D.M. Topologically-associating domains and their long-range contacts are established during early G1 coincident with the establishment of the replication timing program. Genome Res. 2015, 25, 1104–1113. [Google Scholar] [CrossRef] [PubMed]
  25. Dixon, J.R.; Jung, I.; Selvaraj, S.; Shen, Y.; Antosiewicz-Bourget, J.E.; Lee, A.Y.; Ye, Z.; Kim, A.; Rajagopal, N.; Xie, W.; et al. Chromatin architecture reorganization during stem cell differentiation. Nature 2015, 518, 331–336. [Google Scholar] [CrossRef] [PubMed]
  26. Lemaitre, J.M.; Danis, E.; Pasero, P.; Vassetzky, Y.; Méchali, M. Mitotic remodeling of the replicon and chromosome structure. Cell 2005, 123, 787–801. [Google Scholar] [CrossRef] [PubMed]
  27. Dileep, V.; Rivera-Mulia, J.C.; Sima, J.; Gilbert, D.M. Large-Scale Chromatin Structure-Function Relationships during the Cell Cycle and Development: Insights from Replication Timing. Cold Spring Harb. Symp. Quant. Biol. 2015, 80, 53–63. [Google Scholar] [CrossRef] [PubMed]
  28. Nishimura, K.; Fukagawa, T.; Takisawa, H.; Kakimoto, T.; Kanemaki, M. An auxin-based degron system for the rapid depletion of proteins in nonplant cells. Nat. Methods 2009, 6, 917–922. [Google Scholar] [CrossRef] [PubMed]
  29. Natsume, T.; Kiyomitsu, T.; Saga, Y.; Kanemaki, M.T. Rapid Protein Depletion in Human Cells by Auxin-Inducible Degron Tagging with Short Homology Donors. Cell Rep. 2016, 15, 210–218. [Google Scholar] [CrossRef] [PubMed]
  30. Rao, S.S.P.; Huang, S.-C.; Glenn St Hilaire, B.; Engreitz, J.M.; Perez, E.M.; Kieffer-Kwon, K.-R.; Sanborn, A.L.; Johnstone, S.E.; Bascom, G.D.; Bochkov, I.D.; et al. Cohesin Loss Eliminates All Loop Domains. Cell 2017, 171, 305.e24–320.e24. [Google Scholar] [CrossRef] [PubMed]
  31. Marchal, C.; Sasaki, T.; Vera, D.; Wilson, K.; Sima, J.; Rivera-Mulia, J.C.; Trevilla-García, C.; Nogues, C.; Nafie, E.; Gilbert, D.M. Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq. Nat. Protoc. 2018, 13, 819–839. [Google Scholar] [CrossRef] [PubMed]
  32. Peace, J.M.; Villwock, S.K.; Zeytounian, J.L.; Gan, Y.; Aparicio, O.M. Quantitative BrdU immunoprecipitation method demonstrates that Fkh1 and Fkh2 are rate-limiting activators of replication origins that reprogram replication timing in G1 phase. Genome Res. 2016, 26, 365–375. [Google Scholar] [CrossRef] [PubMed]
  33. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  34. Dekker, J.; Belmont, A.S.; Guttman, M.; Leshyk, V.O.; Lis, J.T.; Lomvardas, S.; Mirny, L.A.; O’Shea, C.C.; Park, P.J.; Ren, B.; et al. The 4D nucleome project. Nature 2017, 549, 219–226. [Google Scholar] [CrossRef] [PubMed]
  35. Consortium, T.E.P. An integrated encyclopedia of DNA elements in the human genome. Nature 2012, 489, 57–74. [Google Scholar] [CrossRef] [PubMed]
  36. Karolchik, D.; Hinrichs, A.S.; Furey, T.S.; Roskin, K.M.; Sugnet, C.W.; Haussler, D.; Kent, W.J. The UCSC Table Browser data retrieval tool. Nucleic Acids Res. 2004, 32, D493–D496. [Google Scholar] [CrossRef] [PubMed]
  37. Zerbino, D.R.; Achuthan, P.; Akanni, W.; Amode, M.R.; Barrell, D.; Bhai, J.; Billis, K.; Cummins, C.; Gall, A.; Girón, C.G.; et al. Ensembl 2018. Nucleic Acids Res. 2018, 46, D754–D761. [Google Scholar] [CrossRef] [PubMed]
  38. Neph, S.; Kuehn, M.S.; Reynolds, A.P.; Haugen, E.; Thurman, R.E.; Johnson, A.K.; Rynes, E.; Maurano, M.T.; Vierstra, J.; Thomas, S.; et al. BEDOPS: High-performance genomic feature operations. Bioinformatics 2012, 28, 1919–1920. [Google Scholar] [CrossRef] [PubMed]
  39. Karimzadeh, M.; Ernst, C.; Kundaje, A.; Hoffman, M.M. Umap and Bismap: Quantifying genome and methylome mappability. Nucleic Acids Res. 2018, 46, e120. [Google Scholar] [CrossRef] [PubMed]
  40. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  41. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  42. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, Y.; Song, F.; Zhang, B.; Zhang, L.; Xu, J.; Kuang, D.; Li, D.; Choudhary, M.N.K.; Li, Y.; Hu, M.; et al. The 3D Genome Browser: A web-based browser for visualizing 3D genome organization and long-range chromatin interactions. Genome Biol. 2018, 19, 151. [Google Scholar] [CrossRef] [PubMed]
  44. Durand, N.C.; Shamim, M.S.; Machol, I.; Rao, S.S.P.; Huntley, M.H.; Lander, E.S.; Aiden, E.L. Juicer Provides a One-Click System for Analyzing Loop-Resolution Hi-C Experiments. Cell Syst. 2016, 3, 95–98. [Google Scholar] [CrossRef] [PubMed]
  45. Weddington, N.; Stuy, A.; Hiratani, I.; Ryba, T.; Yokochi, T.; Gilbert, D.M. ReplicationDomain: A visualization tool and comparative database for genome-wide replication timing data. BMC Bioinforma. 2008, 9, 530. [Google Scholar] [CrossRef] [PubMed]
  46. Sima, J.; Chakraborty, A.; Dileep, V.; Michalski, M.; Rivera-Mulia, J.C.; Trevilla-Garcia, C.; Klein, K.N.; Bartlett, D.; Washburn, B.K.; Paulsen, M.T.; et al. Identification of cis elements for spatio-temporal control of DNA replication. Cell 2018, 176, 1–15. [Google Scholar] [CrossRef]
Genes 10 00196 g001 550
Figure 1. Auxin treatment results in loss of Scc1. (A) Immunofluorescence (IF) for Scc1 and mClover in HCT116 Scc1-mAID-mClover cells under control conditions or two hours of auxin treatment. (B) Flow cytometry analysis for mClover fluorescence using both HCT116 Scc1-mAID-mClover cell line and parental HCT116 cell line as a negative control. The mClover-positive gate was set based on the parental cell line at 102 fluorescence units. (C) Comparison of the direct IF read-out of loss of Scc1 and indirect read out of loss of mClover through both flow cytometry and IF.
Figure 1. Auxin treatment results in loss of Scc1. (A) Immunofluorescence (IF) for Scc1 and mClover in HCT116 Scc1-mAID-mClover cells under control conditions or two hours of auxin treatment. (B) Flow cytometry analysis for mClover fluorescence using both HCT116 Scc1-mAID-mClover cell line and parental HCT116 cell line as a negative control. The mClover-positive gate was set based on the parental cell line at 102 fluorescence units. (C) Comparison of the direct IF read-out of loss of Scc1 and indirect read out of loss of mClover through both flow cytometry and IF.
Genes 10 00196 g001
Genes 10 00196 g002 550
Figure 2. Acute loss of cohesin in S phase does not perturb replication timing. (A) Example locus overlay of replication timing (RT, log-scaled ratio of read counts in early over late S phase bin) for auxin-treated versus control cells. (B) RT plots for four example chromosomes with significantly-changed loci marked with black dots. The gray box highlights the region shown in (A). (C) Pearson correlation calculated for RT values across the genome compared to previously-published 4D Nucleome RepliSeq data. (D) Violin plot showing the distribution of Z scores for control versus auxin-treated RT difference across the genome, separating out genomic bins covering a TAD boundary. (E) RT for control versus auxin-treated cells across the whole genome, and across genomic coordinates annotated for super enhancers, H3K9me3, CTCF, and genes that showed induction or repression upon auxin treatment in the Scc1-mAID cell line. The gray topology plot in each subfigure shows the whole-genome behavior and genome features with dense annotations (CTCF and H3K9me3) have a topology map in dark blue.
Figure 2. Acute loss of cohesin in S phase does not perturb replication timing. (A) Example locus overlay of replication timing (RT, log-scaled ratio of read counts in early over late S phase bin) for auxin-treated versus control cells. (B) RT plots for four example chromosomes with significantly-changed loci marked with black dots. The gray box highlights the region shown in (A). (C) Pearson correlation calculated for RT values across the genome compared to previously-published 4D Nucleome RepliSeq data. (D) Violin plot showing the distribution of Z scores for control versus auxin-treated RT difference across the genome, separating out genomic bins covering a TAD boundary. (E) RT for control versus auxin-treated cells across the whole genome, and across genomic coordinates annotated for super enhancers, H3K9me3, CTCF, and genes that showed induction or repression upon auxin treatment in the Scc1-mAID cell line. The gray topology plot in each subfigure shows the whole-genome behavior and genome features with dense annotations (CTCF and H3K9me3) have a topology map in dark blue.
Genes 10 00196 g002
Genes 10 00196 g003 550
Figure 3. Acute loss of cohesin in S phase does not perturb patterns of replication foci. (A) Replication foci patterns were assayed using BrdU IF in asynchronous cells after auxin or control treatment. (B) Number of cells in each S phase stage were counted for two biological repeats of auxin-treated and untreated samples.
Figure 3. Acute loss of cohesin in S phase does not perturb patterns of replication foci. (A) Replication foci patterns were assayed using BrdU IF in asynchronous cells after auxin or control treatment. (B) Number of cells in each S phase stage were counted for two biological repeats of auxin-treated and untreated samples.
Genes 10 00196 g003
Genes 10 00196 g004 550
Figure 4. Loss of cohesin from early G1 does not perturb replication timing. (A) Example locus overlay of replication timing (RT, log-scaled ratio of four hour timepoint read counts over ten hour timepoint counts) for auxin-treated versus control cells. (B) RT plots for four example chromosomes with significantly-changed loci marked with black dots. The gray box highlights the region shown in (A). (C) Pearson correlation calculated for RT values across the genome compared to previously-published 4D Nucleome RepliSeq data. (D) Violin plot showing the distribution of Z scores for control versus auxin-treated RT difference across the genome, separating out genomic bins covering a TAD boundary. (E) RT for control versus auxin-treated cells across the whole genome and across genomic coordinates annotated for super enhancers, H3K9me3, CTCF, and genes that showed induction or repression upon auxin treatment in the Scc1-mAID cell line. The gray topology plot in each subfigure shows the whole-genome behavior and genome features with dense annotations (CTCF and H3K9me3) have a topology map in dark blue.
Figure 4. Loss of cohesin from early G1 does not perturb replication timing. (A) Example locus overlay of replication timing (RT, log-scaled ratio of four hour timepoint read counts over ten hour timepoint counts) for auxin-treated versus control cells. (B) RT plots for four example chromosomes with significantly-changed loci marked with black dots. The gray box highlights the region shown in (A). (C) Pearson correlation calculated for RT values across the genome compared to previously-published 4D Nucleome RepliSeq data. (D) Violin plot showing the distribution of Z scores for control versus auxin-treated RT difference across the genome, separating out genomic bins covering a TAD boundary. (E) RT for control versus auxin-treated cells across the whole genome and across genomic coordinates annotated for super enhancers, H3K9me3, CTCF, and genes that showed induction or repression upon auxin treatment in the Scc1-mAID cell line. The gray topology plot in each subfigure shows the whole-genome behavior and genome features with dense annotations (CTCF and H3K9me3) have a topology map in dark blue.
Genes 10 00196 g004
Genes 10 00196 g005 550
Figure 5. Replication timing across TAD boundaries versus compartment boundaries. (A) Visualization of HiC data from Rao et al. (2017) [30] in untreated cells versus auxin-treated cells. The left panel spans 50 Mb to highlight compartments and the right zooms in to 6 Mb from within that locus (gray box) to highlight changes to TAD structure. (B) Replication timing metaplots for windows 5 Mb on either side of TAD boundaries. Invariant boundaries are those called in both control and auxin-treated cells, Scc1-dependent boundaries are those that were lost upon auxin treatment.
Figure 5. Replication timing across TAD boundaries versus compartment boundaries. (A) Visualization of HiC data from Rao et al. (2017) [30] in untreated cells versus auxin-treated cells. The left panel spans 50 Mb to highlight compartments and the right zooms in to 6 Mb from within that locus (gray box) to highlight changes to TAD structure. (B) Replication timing metaplots for windows 5 Mb on either side of TAD boundaries. Invariant boundaries are those called in both control and auxin-treated cells, Scc1-dependent boundaries are those that were lost upon auxin treatment.
Genes 10 00196 g005
Table
Table 1. Genomic loci showing significant difference in both asynchronous and synchronized experiments.
Table 1. Genomic loci showing significant difference in both asynchronous and synchronized experiments.
ChrStartEndRT Control—RT IAA
AsynchronousSynchronous
chr3198100000198150000−3.4403373.242949
chr449300000493500002.7230622.296154
chr41194000001194500003.486406−1.410933
chr41900000001900500004.479181−1.897292
chr419005000019010000012.678844−5.136778
chr419010000019015000025.976415−10.545368
chr54965000049700000−3.2582034.272367
chr51812500001813000002.921178−1.623559
chr756850000569000002.7792651.355865
chr91500002000003.343623−1.447844
chr131143000001143500003.4607911.561016
chr131143500001143643287.8563823.064328
chr1690100000901500004.566427−3.221925
chr21635000064000003.521525−2.910266

Share and Cite

MDPI and ACS Style

Oldach, P.; Nieduszynski, C.A. Cohesin-Mediated Genome Architecture Does Not Define DNA Replication Timing Domains. Genes 2019, 10, 196. https://doi.org/10.3390/genes10030196

AMA Style

Oldach P, Nieduszynski CA. Cohesin-Mediated Genome Architecture Does Not Define DNA Replication Timing Domains. Genes. 2019; 10(3):196. https://doi.org/10.3390/genes10030196

Chicago/Turabian Style

Oldach, Phoebe, and Conrad A. Nieduszynski. 2019. "Cohesin-Mediated Genome Architecture Does Not Define DNA Replication Timing Domains" Genes 10, no. 3: 196. https://doi.org/10.3390/genes10030196

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop