Martijn R H Zwinderman1, Thamar Jessurun Lobo2, Petra E van der Wouden1, Diana C J Spierings2, Marcel A T M van Vugt3, Peter M Lansdorp2,4,5, Victor Guryev2, Frank J Dekker1. 1. Department of Chemical and Pharmaceutical Biology, Groningen Research Institute of Pharmacy, University of Groningen, 9713 AV Groningen, The Netherlands. 2. European Research Institute for the Biology of Ageing, University of Groningen, University Medical Center Groningen, 9713 AV Groningen, The Netherlands. 3. Department of Medical Oncology, University Medical Center Groningen, University of Groningen, 9713 GZ Groningen, The Netherlands. 4. Terry Fox Laboratory, British Columbia Cancer Agency, Vancouver, V5Z 1L3 British Columbia, Canada. 5. Department of Medical Genetics, University of British Columbia, Vancouver, V6T 1Z4 British Columbia, Canada.
Abstract
Following DNA replication, equal amounts of chromatin proteins are distributed over sister chromatids by re-deposition of parental chromatin proteins and deposition of newly synthesized chromatin proteins. Molecular mechanisms balancing the allocation of new and old chromatin proteins remain largely unknown. Here, we studied the genome-wide distribution of new chromatin proteins relative to parental DNA template strands and replication initiation zones using the double-click-seq. Under control conditions, new chromatin proteins were preferentially found on DNA replicated by the lagging strand machinery. Strikingly, replication stress induced by hydroxyurea or curaxin treatment and inhibition of ataxia telangiectasia and Rad3-related protein (ATR) or p53 inactivation inverted the observed chromatin protein deposition bias to the strand replicated by the leading strand polymerase in line with previously reported effects on replication protein A occupancy. We propose that asymmetric deposition of newly synthesized chromatin proteins onto sister chromatids reflects differences in the processivity of leading and lagging strand synthesis.
Following DNA replication, equal amounts of chromatin proteins are distributed over sister chromatids by re-deposition of parental chromatin proteins and deposition of newly synthesized chromatin proteins. Molecular mechanisms balancing the allocation of new and old chromatin proteins remain largely unknown. Here, we studied the genome-wide distribution of new chromatin proteins relative to parental DNA template strands and replication initiation zones using the double-click-seq. Under control conditions, new chromatin proteins were preferentially found on DNA replicated by the lagging strand machinery. Strikingly, replication stress induced by hydroxyurea or curaxin treatment and inhibition of ataxia telangiectasia and Rad3-related protein (ATR) or p53 inactivation inverted the observed chromatin protein deposition bias to the strand replicated by the leading strand polymerase in line with previously reported effects on replication protein A occupancy. We propose that asymmetric deposition of newly synthesized chromatin proteins onto sister chromatids reflects differences in the processivity of leading and lagging strand synthesis.
Post-translational
modifications (PTMs) of histones play important
roles in the regulation of nuclear organization and gene transcription.[1] Histone PTMs are relatively stable and heritable
during DNA replication.[2,3] Importantly, PTMs of parental
histones differ from those of de novo-synthesized
histones.[4] Therefore, near-symmetrical
deposition of old and new histones during DNA replication is crucial
to maintain similar chromatin states for the two sister chromatids
following cell division. The functional importance of balanced histone
inheritance was underscored by the observation that repressed chromatin
domains are preserved by local re-deposition of parental histones.[5] Clearly, mechanisms for accurate deposition of
parental and new histones following DNA replication are crucial.[6]Recent studies have provided insight into
histone deposition during
DNA replication by using immunoprecipitation of PTMs on histones that
are characteristic for either parental histones (H4K20me2) or new
histones (H4K5ac).[7,8] When this approach was combined
with labeling of new DNA using the thymidine analogue 5-ethynyl-2′-deoxyuridine
(EdU) to separate parental from newly synthesized DNA, a slight deposition
bias of new histones to the lagging strand during DNA replication
was revealed in mouse embryonic stem cells.[7] A similar method was applied to budding yeast cells treated with
hydroxyurea (HU), which revealed a slight bias of new histone deposition
onto the leading strand.[8] A clear explanation
for these opposing findings has not been reported. Of note, both studies
concluded that the structural integrity of the replisome is essential
to maintain near-symmetrical histone inheritance.The structural
integrity of the replication fork machinery is often
compromised in cancer due to replication stress.[9] Experimentally, replication stress can be induced by inhibition
of the enzyme ribonucleotide reductase (RNR) with HU.[10,11] RNR inhibition results in a reduction of nucleotides required for
DNA synthesis[12] and leads to uncoupling
of the replicative helicase from the leading strand DNA polymerase.[13] Helicase–polymerase uncoupling activates
a DNA damage response (DDR), which involves the cell cycle checkpoint
kinase ataxia telangiectasia and Rad3-related protein (ATR)[14,15] and p53-dependent transcriptional effects.[16] However, the connection between replication stress and DDR pathways
in relation to histone distribution following DNA replication is largely
unexplored.In this study, we have developed the double-click-seq
protocol
to enrich parental DNA bound to new chromatin proteins. In contrast
to previous methods involving immunoprecipitation of histone PTMs,
this approach is based on metabolic labeling of new chromatin proteins
(Figure A). The double-click-seq
provides stable labeling of de novo synthesized chromatin
proteins by co-translational incorporation of the methionine surrogate
azidohomoalanine (AHA), which allows enrichment of nucleosomes that
contain new chromatin proteins following a click reaction with a biotin
affinity tag. Using this methodology, we show that both replication
stress as induced by HU and inhibition of the DDR invert the generally
asymmetric distribution of new chromatin proteins onto replicated
DNA strands in human cells. We propose that asymmetric deposition
of newly synthesized chromatin proteins onto sister chromatids reflects
differences in the processivity of leading and lagging strand synthesis.
Figure 1
Double-click-seq
reveals a bias in the deposition of new chromatin
proteins toward the DNA strand replicated by lagging strand machinery.
(A) Schematic overview of the double-click-sequencing method. (B)
Effect of different culture methods on cell doubling time and the
extent of AHA incorporation in histone H4 determined by mass spectrometry.
Culturing hRPE-1 in 95% AHA medium supplemented with 5% DMEM for 20
h resulted in 3% AHA incorporation in histone H4 without a significant
impact on cell doubling times (black box). (C) Western blot and gel
electrophoresis analysis of labeled and unlabeled histones. Reactions
were performed on nuclei with or without the click reaction as indicated above. Positions of histones
are indicated to the right of the blot. (D) Average RFD at replication
initiation zones. A plot including confidence intervals can be found
in Supporting Information Figure S1B. (E)
Average bias of new histone deposition at replication initiation zones
under the untreated condition (red line) and with olaparib treatment
(yellow line). The black line shows the background signal in hRPE-1
cells (negative control and without second click reaction). Separate
plots of individual replicates including confidence intervals can
be found in Supporting Information Figure
S2. (F) Heatmap of new histone deposition bias at replication initiation
zones under the untreated condition.
Double-click-seq
reveals a bias in the deposition of new chromatin
proteins toward the DNA strand replicated by lagging strand machinery.
(A) Schematic overview of the double-click-sequencing method. (B)
Effect of different culture methods on cell doubling time and the
extent of AHA incorporation in histone H4 determined by mass spectrometry.
Culturing hRPE-1 in 95% AHA medium supplemented with 5% DMEM for 20
h resulted in 3% AHA incorporation in histone H4 without a significant
impact on cell doubling times (black box). (C) Western blot and gel
electrophoresis analysis of labeled and unlabeled histones. Reactions
were performed on nuclei with or without the click reaction as indicated above. Positions of histones
are indicated to the right of the blot. (D) Average RFD at replication
initiation zones. A plot including confidence intervals can be found
in Supporting Information Figure S1B. (E)
Average bias of new histone deposition at replication initiation zones
under the untreated condition (red line) and with olaparib treatment
(yellow line). The black line shows the background signal in hRPE-1
cells (negative control and without second click reaction). Separate
plots of individual replicates including confidence intervals can
be found in Supporting Information Figure
S2. (F) Heatmap of new histone deposition bias at replication initiation
zones under the untreated condition.
Results
Double-Click-Seq
Reveals that Bias of New Chromatin Proteins
is Skewed toward the Lagging Strand
We used immortalized
human retinal pigment epithelial (hRPE-1) cells to study chromatin
protein deposition onto replicated DNA. The hRPE-1 cell line has a
well-characterized genome and serves as a near-diploid model of normal
human cells. The initial step for development of the double-click-seq
protocol was optimization of the conditions for co-translational incorporation
of AHA as a replacement of methionine into newly synthesized methionine-containing
proteins over the time of one replication cycle in hRPE-1. We aimed
to find conditions with a minimal effect on the average cell doubling
time.[17] Cells were cultured under conditions
with various AHA and methionine concentrations, and histones were
extracted and analyzed with mass spectrometry in order to assess the
levels of AHA incorporation. Conditions were identified in which cells
were cultured in a medium containing a mixture of 95% AHA and 5% methionine
for 20 h to reach a 3% AHA incorporation in histone H4 (Figure B), which were subsequently
used throughout the study. The 20 h labeling period allowed us to
track histone deposition across a full cell cycle. Western blot analysis
showed that AHA was also incorporated into histones H2A/B and H3 (Figure C).After AHA
labeling of de novo synthesized proteins, nuclei
were isolated to remove cytoplasmic AHA-labeled proteins. The nuclei
were then exposed to an initial copper-catalyzed [3 + 2] cycloaddition
reaction between the azide group of AHA and an alkyne linked to biotin
(click), which enabled affinity enrichment of AHA-labeled histones
in complex with DNA (i.e., nucleosomes) and other
AHA-labeled proteins using streptavidin beads after MNAse digestion.
Subsequently, several washing steps were performed to rid all DNA
not bound to chromatin proteins. The remaining nucleosomal DNA was
then extracted from the enriched nascent nucleosomes by digestion
of chromatin proteins and streptavidin by proteinase K treatment,
after which the resulting enriched DNA was ligated to forked adaptors.
In a second click reaction, nascent DNA strands labeled with EdU[7] were coupled to an azide-linked biotin, thus
enabling a second streptavidin bead capture. Finally, the unlabeled
parental single stranded DNA was eluted from the beads with an alkaline
solution and amplified to construct a directional short read sequencing
library. The library was then sequenced with next generation sequencing
and the resulting paired-end data (summary of all libraries in Table S1) were processed to assign genome-wide
deposition of new chromatin proteins to either the forward or the
reverse strand around the center of replication initiation zones (termed
“origin centers” in all related figures). Autosomal
replication initiation zones were mapped using a publicly available
Okazaki-fragment sequencing (OK-seq) data set for hRPE-1 cells[18] (n = 9,608; Figure D and Supporting Information Figure S1A).The double-click-seq revealed
that the deposition bias of new chromatin
proteins was skewed toward the strand replicated by the lagging strand
machinery at mapped replication initiation zones (Figure E, red line and Figure F). This finding is in agreement
with the bias observed in mouse embryonic stem cells.[7] However, our experiments indicate a more pronounced partition
score as a measure of chromatin protein deposition bias of around
0.05 (calculated as the proportion of forward and reverse read counts,
see Figure S1C), whereas Petryk et al. reported a partition score of approximately 0.008
for differential histone PTMs. As we observe an average replication
fork directionality (RFD) of 0.42 (Supporting Information Figure S1D) where Petryk et al. report an average RFD of around 0.13, variability in the replication
initiation zone firing rate between the different cell types used
for both experiments can only partially explain the sixfold less pronounced
asymmetry observed through immunoprecipitation of distinctive PTMs.[7] The remaining difference might be the result
of signal dilution through the exchange of histone PTMs over time,
while our method stably labels new chromatin proteins. Nonetheless,
with the double-click-seq, we observed a clear bias of new chromatin
proteins around replication initiation zones toward the strand replicated
by lagging strand polymerases. It is worth noting that our estimate
of deposition bias might be conservative since turnover times range
from fast for histones H2A/H2B to slow for H4 (0.4% per hour).[6]We hypothesized that the deposition bias
of new chromatin proteins
toward the lagging strand is most likely the result of the relatively
less efficient capture of parental chromatin proteins by the lagging
strand, thus leaving the lagging strand with more new chromatin proteins
compared to the leading strand. Consequently, interfering with lagging
strand synthesis would result in a more pronounced deposition bias
of new chromatin proteins toward the lagging DNA strand. During lagging
strand synthesis, the 5′ flap of Okazaki fragments is excised
by the nuclease FEN1, upon which the fragments are ligated by DNA
ligase I (LIG1).[19] Under untreated conditions,
a fraction of Okazaki fragments escape LIG1-mediated ligation and
are processed by a poly(ADP-ribose) polymerase (PARP)-mediated back-up
route.[20] Thus, inhibition of PARP activity
interferes with the completion of lagging strand synthesis. In line
with our hypothesis, we observed that treatment with the PARP inhibitor
olaparib[21] increased the deposition bias
of new chromatin proteins onto the DNA strand replicated by lagging
strand machinery (Figure E, yellow line). This indicates that disturbance of lagging
strand synthesis increases the deposition bias of new chromatin proteins
to the lagging strand compared to the untreated condition.
Replication
Stress Inverts the General Asymmetric Deposition
of New Chromatin Proteins
Given the observed increase in
the deposition asymmetry of new chromatin proteins toward the lagging
strand upon interference with lagging strand synthesis, we wondered
what the effect of more general replication stress on new chromatin
protein deposition would be. Strikingly, HU-induced replication stress
completely inverted the bias of new chromatin proteins toward the
DNA strand replicated by leading strand synthesis (Figure A). An inversion of the deposition
bias of new chromatin proteins to the leading strand was also found
upon treatment with the DNA intercalator curaxin (Figure B), which induces replication
stress by triggering nucleosome unfolding.[22] We checked if changes in origin usage may underlie the inverted
chromatin protein deposition bias; however, we found that the OK-seq
profiles between untreated samples and HU treatment were highly correlated
(Spearman correlation 0.98). Therefore, we conclude that replication
stress, such as that induced by either HU or curaxin treatment, inverts
the asymmetry in new chromatin protein deposition relative to the
untreated samples. Analogously, the presence or absence of HU may
explain the apparent contradiction between the leading strand bias
of new histones found by Yu et al.(8) and the lagging strand bias of new histones found by Petryk et al.(7)
Figure 2
Replication stress inverts
the deposition asymmetry of new chromatin
proteins. (A) Average bias of new chromatin protein deposition at
replication initiation zones during replication stress induced with
HU. (B) Average bias of new chromatin protein deposition at replication
initiation zones during replication stress induced with curaxin. Separate
plots for individual replicates including confidence intervals can
be found in Supporting Information Figure
S2.
Replication stress inverts
the deposition asymmetry of new chromatin
proteins. (A) Average bias of new chromatin protein deposition at
replication initiation zones during replication stress induced with
HU. (B) Average bias of new chromatin protein deposition at replication
initiation zones during replication stress induced with curaxin. Separate
plots for individual replicates including confidence intervals can
be found in Supporting Information Figure
S2.
Inhibition of DDR Pathways
Also Inverts New Chromatin Protein
Deposition Asymmetry
As replication stress can lead to activation
of the DDR, we speculated that the effect of replication stress on
the deposition bias of new chromatin proteins might be increased by
inhibition of the DDR. The ATR DDR pathway plays an essential role
in suppressing replication stress.[15] Additional
DDR pathways involve the tumor suppressor p53.[16] We therefore investigated the involvement of the ATR- and
p53-mediated DDR pathways in the altered chromatin protein deposition
pattern during replication stress. Toward this aim, we employed p53-null
hRPE-1 cells. Lack of p53 expression was confirmed by Western blot
analysis (Figure A).
In untreated p53-null cells, a deposition bias of new histones to
the leading strand was observed (Figure B), indicating that loss of p53 already inverts
the deposition bias of new histones seen in untreated hRPE-1 cells.
The bias inversion increased further upon HU treatment of the p53-null
cells (Figure B).
Pharmacological inhibition of ATR in hRPE-1 cells also led to inversion
of the deposition bias of new histones toward the leading strand.
Again, HU treatment provided a stronger inversion of the deposition
bias of new chromatin proteins (Figure C). As an example, we included the heatmap of the deposition
bias of new chromatin proteins in cells treated with the ATR inhibitor
and HU, showing that chromatin protein bias inversion was observed
consistently at the replication initiation zones that were identified
(Figure D). We conclude
that both replication stress and inhibition of DDR pathways through
ATR or p53 inactivation disturb the asymmetry in the deposition of
new histones such that the bias is inverted from the lagging to the
leading strand.
Figure 3
Inhibition of DDR pathways also inverts the deposition
asymmetry
of new chromatin proteins (A). Western blot showing knock-out of p53
in hRPE-1 cells with gene editing. Equal protein loading was ensured
by measuring the protein concentration with the Bradford assay. (B)
Average bias of new histone deposition in p53-null hRPE-1 cells treated
with HU (light blue line) and without (pink line) at replication initiation
zones. (C) Average bias of new histone deposition at replication initiation
zones in hRPE-1 cells with either ataxia telangiectasia and Rad3-related
protein (ATR) inhibition (green line) or ATR inhibition and HU treatment
(dark blue line). (D) Heatmap of bias of new histone deposition at
replication initiation zones in cells treated with the ATR inhibitor
and HU. Separate bias plots of individual replicates including confidence
intervals can be found in Supporting Information Figure S2.
Inhibition of DDR pathways also inverts the deposition
asymmetry
of new chromatin proteins (A). Western blot showing knock-out of p53
in hRPE-1 cells with gene editing. Equal protein loading was ensured
by measuring the protein concentration with the Bradford assay. (B)
Average bias of new histone deposition in p53-null hRPE-1 cells treated
with HU (light blue line) and without (pink line) at replication initiation
zones. (C) Average bias of new histone deposition at replication initiation
zones in hRPE-1 cells with either ataxia telangiectasia and Rad3-related
protein (ATR) inhibition (green line) or ATR inhibition and HU treatment
(dark blue line). (D) Heatmap of bias of new histone deposition at
replication initiation zones in cells treated with the ATR inhibitor
and HU. Separate bias plots of individual replicates including confidence
intervals can be found in Supporting Information Figure S2.
Deposition Bias of New
Chromatin Proteins Is More Pronounced
in AT-Rich Regions
Finally, we investigated whether the chromatin
protein deposition bias depends on common characteristics of genomic
regions. We categorized replication initiation zones into actively
transcribed (active) or not actively transcribed (inactive) (Figure C), early replicated
or late replicated (Figure F), and AT-rich or GC-rich zones (Figure I). We then determined the average deposition
bias of new chromatin proteins around each category of initiation
zones for both untreated hRPE-1 cells (“untreated”, Figure A,D,G) and cells
treated with the ATR inhibitor and HU (“replication stress”, Figure B,E,H). We conclude
that the deposition bias tends to be stronger in late-replicated and
AT-rich zones. Slight differences between replication initiation zones
(Figure C,F,I) were
observed too, but these differences were not large enough to account
for the observed differences in chromatin protein deposition bias
(Figure A,B).
Figure 4
Deposition
bias of new chromatin proteins is more pronounced in
AT-rich regions. (A,D,G) Average bias of new chromatin protein deposition
at replication initiation zones under untreated conditions. (B,E,H)
Average bias of new chromatin protein deposition at replication initiation
zones with ataxia telangiectasia and Rad3-related protein (ATR) inhibition
and HU treatment. (C,F,I) Average RFD at replication initiation zones.
For (A,B,C), the replication initiation zones were categorized into
actively transcribed (red) and not actively transcribed (orange) regions.
For (D,E,F), the zones were categorized into actively transcribed
(light blue) and not actively transcribed (dark blue) regions. For
(G,H,I), the zones were categorized into AT-rich (black) and GC-rich
(gray) regions. Separate plots of individual replicates split into
quartiles and including confidence intervals can be found in Supporting Information Figure S3.
Deposition
bias of new chromatin proteins is more pronounced in
AT-rich regions. (A,D,G) Average bias of new chromatin protein deposition
at replication initiation zones under untreated conditions. (B,E,H)
Average bias of new chromatin protein deposition at replication initiation
zones with ataxia telangiectasia and Rad3-related protein (ATR) inhibition
and HU treatment. (C,F,I) Average RFD at replication initiation zones.
For (A,B,C), the replication initiation zones were categorized into
actively transcribed (red) and not actively transcribed (orange) regions.
For (D,E,F), the zones were categorized into actively transcribed
(light blue) and not actively transcribed (dark blue) regions. For
(G,H,I), the zones were categorized into AT-rich (black) and GC-rich
(gray) regions. Separate plots of individual replicates split into
quartiles and including confidence intervals can be found in Supporting Information Figure S3.Genes that are actively transcribed are generally replicated
in
the early S phase[23] and DNA regions replicated
in the early S phase generally have higher GC-content.[24] In order to estimate the relative contribution
of these genomic features to histone deposition asymmetry, we performed
a multiple regression analysis taking all three factors into account
(Table S2). The GC-content proved to be
the most consistently significant predictor of histone deposition
bias. It was significant (uncorrected p-value <
0.01) for all used conditions, with the exception of the p53 knock-out
cells. The GC-content beta values indeed signified a bias increasing
effect in AT-rich zones across all conditions. Replication timing
was a significant predictor for the bias in cells treated with the
ATR inhibitor, p53 knock-out cells, and for one out of two replicates
for cells treated with curaxin and cells treated with both the ATR
inhibitor and HU. Transcriptional status did not have a significant
effect on the deposition bias of new chromatin proteins (Table S2).
Discussion
We
developed the double-click-seq method to study the incorporation
of new chromatin proteins into nascent chromatin. We tracked de novo synthesized chromatin proteins in cultured human
cells and showed that they are enriched under untreated conditions
in DNA replicated by the lagging strand machinery, especially in AT-rich
regions. We found a more pronounced lagging strand bias in the deposition
of new chromatin proteins than reported previously[7] and we found this to be the case around all mapped replication
initiation zones.Most likely, redistribution of parental chromatin
proteins onto
the DNA strand replicated by either leading or lagging strand polymerases
favors the DNA strand that finishes replication first. The leading
strand polymerase replicates DNA in a continuous fashion and directly
follows unwinding of the parental DNA strands by the replicative helicase,
such that the leading strand has the highest chance to capture a majority
of the parental chromatin proteins that are evicted ahead of the replication
fork. The discontinuously replicated lagging strand by default would
receive a majority of new chromatin proteins to fill the gaps between
re-deposited parental chromatin proteins. In support of this model,
we show that interfering with the completion of lagging strand synthesis
with the PARP inhibitor olaparib increases the asymmetry in new chromatin
protein deposition toward the lagging strand. Thus, with olaparib,
the replicated lagging strand becomes less able to capture displaced
parental histones and receives even more new histones than would normally
be the case.The asymmetry of new chromatin protein deposition
toward the lagging
strand underwent a striking inversion toward the leading strand upon
treatment with the replication stress inducers HU and curaxin. Replication
stress is known to induce uncoupling between the helicase and the
leading strand polymerase, which leads to single-stranded DNA (ssDNA)
that requires stabilization by binding of replication protein A (RPA).[25] Interestingly, the asymmetry in RPA occupancy,
as an indirect measure of ssDNA, toward the lagging strand has been
found to invert upon replication stress.[25] Similarly, another study has shown that replication stress induces
strand switching of the DNA clamp Proliferating Cell Nuclear Antigen
(PCNA), such that more PCNA is found on the leading strand under stress
conditions compared to untreated conditions.[26] Both observations indicate that replication stress increases the
proportion of ssDNA on the leading strand to a greater extent than
on the lagging strand. The independently operating lagging strand
polymerase is apparently better equipped to handle replication stress
than its leading strand counterpart, perhaps in part due to repeated
relocation of the lagging strand polymerase following completion of
an Okazaki fragment.[27] Therefore, when
replication stressors target the replisome as a whole, the lagging
strand polymerase is able to continue synthesis more effectively than
the uncoupled leading strand polymerase, which is reflected in the
bias inversions of RPA, PCNA, and new histones. Notably, RPA at the
ssDNA is an important activator for ATRIP-ATR signaling,[28,29] involving multiple downstream effects, including on chromatin composition,[30] and could explain local differences in histone
deposition. In essence, the DNA strand that is replicated by the lagging
strand polymerase “leads” during replication stress.Accumulation of RPA on the leading strand following replication
stress-induced helicase-polymerase uncoupling triggers a DDR that
attempts to resolve the aberrant fork structure.[13] Clearly, the DDR does not prevent inversion of new chromatin
protein deposition asymmetry upon treatment with HU or curaxin. Previously,
it has been demonstrated that chronic treatment of chicken DT40 cells
with a low-dose of HU for 7 days stochastically induced a reduction
in the expression of certain active genes, which was connected to
a loss of the chromatin marks H3K4me3 and H3K9/14ac.[11] It was postulated that the observed loss of histone marks
was the result of uncoupling of histone recycling from DNA synthesis
due to replication stress-induced helicase–polymerase uncoupling.
Our results provide more direct evidence for this hypothesis. Moreover,
the results of Papadopoulou et al. illustrate that
asymmetric deposition of new histones may accumulate over multiple
cell divisions and cause epigenetic instability.[11]Additionally, we observed that functional DDR pathways
are important
in maintaining the characteristic pattern of chromatin protein deposition
since pharmacological inhibition of ATR and p53 knockout inverted
the histone deposition pattern already without HU-induced replication
stress. ATR is constitutively active in the S phase to sense ongoing
DNA replication and repress FOXM1 activity.[31] After completion of the S phase, a drop in ATR activity releases
the brakes on FOXM1 activity and allows cells to enter mitosis. Apparently,
inhibiting this intrinsic S/G2 checkpoint alters asymmetric
histone deposition. Similarly, inactivation of TP53 not only potentiated the observed bias inversion of new chromatin
proteins in HU-induced replication stress but also caused inversion
by itself. This alludes to a constitutive function of p53 in maintaining
replication fork integrity.While OK-seq profiles for hRPE-1
cells treated with HU for 3 h
or without HU were strongly correlated, we did not determine the effect
of a 20 h treatment with HU, the ATR inhibitor, olaparib, or p53 knock-out
on OK-seq profiles. We can therefore not exclude the possibility that
longer HU treatment and/or other treatments result in firing of more
dormant origins, causing the asymmetry to seem weaker than it truly
is. Additionally, some secondary cell divisions cannot be excluded
because of the long labeling time. Such events might again weaken
the asymmetry signal but not enhance it or change its direction.Furthermore, we observe that chromatin protein deposition asymmetry
is increased close to replication initiation zones located in AT-rich
DNA. AT-rich DNA is less stable than GC-rich DNA due to a difference
in hydrogen bonding.[32] This difference
in DNA stability leads to a higher rate of DNA unwinding in AT-rich
compared to GC-rich regions, not taking the chromatic context into
account.[33] Such an increase in DNA unwinding
rate may increase leading strand synthesis speed in untreated conditions
and helicase–polymerase uncoupling during replication stress.
In both cases, this would lead to a greater difference in leading
and lagging strand processivity in AT-rich compared to GC-rich regions
and thus result in an increase in the asymmetry of new chromatin protein
deposition.We note that EdU labeling affects the cell cycle
by an increase
of the G2/M phase (Supporting Information Figure S4) in line with literature describing that EdU delays the
cell cycle but does not affect the rate of elongation of replication
forks during synthesis.[34] Since application
of EdU is currently the only technology to separate new and old DNA,
it is not possible to include a control experiment for the effect
of EdU itself. This demonstrates, on one hand, the need for a Blanc
experiment as a reference for further experimental variation (Figure E, black line) and,
on the other hand, the need to advance this type of technology further
to enable analysis of histone deposition bias under less perturbing
conditions.In our hands, metabolic labeling of new histones
with AHA is limited
to approximately 3% of cellular histone H4 over the course of a 20
h incubation period. All histone isoforms contain a methionine that
can be replaced by AHA during de novo synthesis.
Therefore, a labeling period that covers a full cell cycle leads to
AHA incorporation in histones H3.3 and H2A/B, which have a replication-independent
turnover.[35] Histone H4 may also turnover
independently from replication, but this seems mostly confined to
centromeric regions.[36] Regardless, replication-independent
turnover of histones could decrease the chromatin protein deposition
asymmetry but not enhance it or change its direction. Also, we considered
the possibility that the asymmetric distribution of newly synthesized
proteins onto replicated DNA could reflect other proteins than histones,
for example, proteins that transiently associate with DNA during replication
such as RPA and PCNA. We dismissed this possibility because (1) histones
are by far the most abundant proteins associated with DNA at any time,
(2) in our double-click-seq method, we sequence DNA typically long
after replication (by pulling down recently synthesized proteins and
sequencing associated parental DNA template strands), and (3) we enrich
for DNA fragments <200 bp and exclude fragments with size <145
bp from analysis. We also considered the fact that some old and newly
synthesized histones will be exchanged after DNA replication before
being pulled down in our method. Such post replication histone exchange
events are expected to diminish the observed asymmetry in the distribution
of new versus old chromatin proteins over DNA replicated
by leading versus lagging strand synthesis. Inversion
of this asymmetry upon replication stress is furthermore difficult
to reconcile with histone exchanges after DNA replication. For a more
general analysis of strand-specific protein binding at replication
forks, we refer to the eSPAN method.[26]In conclusion, our observations indicate that the deposition of
newly synthesized proteins onto nascent chromatin is inherently biased
toward the lagging strand and that this bias increases upon PARP inhibition.
Replication stress induced by HU or curaxin, genetic inactivation
of p53, or pharmacological ATR kinase inhibition inverts the deposition
bias of new chromatin toward the leading strand. These findings can
be united by a model that includes differences in processivity of
leading and lagging strand synthesis. Normally, polymerase ε
on the leading strand is tightly coupled to the replicative helicase,
whereas polymerase δ on the lagging strand is operating independently.[37] As a result, re-deposition of parental chromatin
is biased toward the faster completed leading strand and deposition
of new chromatin is biased toward the lagging strand (Figure A). Helicase–polymerase
uncoupling, due to a lack of polymerase ε processivity during
replication stress, results in stretches of ssDNA on the leading strand.[37,38] The independently operating polymerase δ maintains processivity
and therefore completes replication faster than polymerase ε.[38] Together, this results in increased re-deposition
of parental chromatin proteins and decreased deposition of new chromatin
proteins on DNA replicated by polymerase δ, thus explaining
the inversion of new chromating protein deposition asymmetry toward
the strand replicated by polymerase ε during replication stress
(Figure B).
Figure 5
Helicase–polymerase
uncoupling during replication stress
inverts deposition of chromatin proteins exemplified by histones.
(A) Model of lagging strand (green arrow heads) bias of new histones
(white spheres) under untreated conditions following the leading strand
(pink arrow) bias of old histones (red spheres), with polymerase ε
(orange) coupled to the helicase (yellow). (B) Model of leading strand
bias of new histone deposition following helicase uncoupling of polymerase
ε. Helicase uncoupling results in lagging of leading strand
synthesis. This results in a leading strand bias for new histones
resulting from a biased re-deposition of old histones onto the faster
completed lagging strand.
Helicase–polymerase
uncoupling during replication stress
inverts deposition of chromatin proteins exemplified by histones.
(A) Model of lagging strand (green arrow heads) bias of new histones
(white spheres) under untreated conditions following the leading strand
(pink arrow) bias of old histones (red spheres), with polymerase ε
(orange) coupled to the helicase (yellow). (B) Model of leading strand
bias of new histone deposition following helicase uncoupling of polymerase
ε. Helicase uncoupling results in lagging of leading strand
synthesis. This results in a leading strand bias for new histones
resulting from a biased re-deposition of old histones onto the faster
completed lagging strand.
Data
Availability
The data sets generated during this study are
available at EBI
ArrayExpress under accession number E-MTAB-8624. The code generated
during this study is available at GitHub; https://github.com/thamarlobo/histone_deposition_analysis.git.
Authors: Pauline N C B Audergon; Sandra Catania; Alexander Kagansky; Pin Tong; Manu Shukla; Alison L Pidoux; Robin C Allshire Journal: Science Date: 2015-04-03 Impact factor: 47.728
Authors: Keith A Menear; Claire Adcock; Robert Boulter; Xiao-ling Cockcroft; Louise Copsey; Aaron Cranston; Krystyna J Dillon; Jan Drzewiecki; Sheila Garman; Sylvie Gomez; Hashim Javaid; Frank Kerrigan; Charlotte Knights; Alan Lau; Vincent M Loh; Ian T W Matthews; Stephen Moore; Mark J O'Connor; Graeme C M Smith; Niall M B Martin Journal: J Med Chem Date: 2008-09-19 Impact factor: 7.446
Authors: Hana Hanzlikova; Ilona Kalasova; Annie A Demin; Lewis E Pennicott; Zuzana Cihlarova; Keith W Caldecott Journal: Mol Cell Date: 2018-07-05 Impact factor: 17.970
Authors: Joshua C Saldivar; Stephan Hamperl; Michael J Bocek; Mingyu Chung; Thomas E Bass; Fernanda Cisneros-Soberanis; Kumiko Samejima; Linfeng Xie; James R Paulson; William C Earnshaw; David Cortez; Tobias Meyer; Karlene A Cimprich Journal: Science Date: 2018-08-24 Impact factor: 47.728