the term replication stress describes the slowing or stalling of  replication forks by endogenously or exogenously derived impediments to dna polymerization ( zeman and cimprich , 2014 ) . replication stressors can be local factors , such as dna damage or secondary structures that affect forks randomly as they are encountered , or global ones , such as nucleotide pool depletion or imbalance that simultaneously slows all forks ( poli et  al . , 2012 , it is now recognized that replication stress induced by nucleotide pool imbalance is an important consequence of the activation of some oncogenes , which drive cells into s phase without upregulation of nucleotide supply ( bester et  al . , 2011 ) . the resulting loss of polymerase processivity is thought to lead to localized uncoupling of the replicative helicase and polymerase and formation of tracts of single - stranded dna ( byun et  al . , 2005 , pacek and walter , 2004 ) . while this normally induces checkpoint activation and senescence ( bartkova et  al . , 2006 , , in cells that can bypass the checkpoint , such replication stress provides a fertile source of genetic instability , particularly in the vicinity of fragile sites and sites capable of forming secondary structures ( de and michor , 2011 , tsantoulis et  al . , 2008 ) . in addition to the extensive genetic changes that have been well documented in many types of cancer , there are also extensive local and global alterations in histone and dna modifications . the consequent changes in chromatin structure are accompanied by significant dysregulation of gene expression ( timp and feinberg , 2013 , berdasco and esteller , 2010 ) , which , since it is not accompanied by changes in the dna sequence , may be considered epigenetic ( berger et  al . , 2009 ) . these epigenetic changes could act alongside genetic instability to produce clonal variation within a tumor , upon which selective pressure can act , and so may contribute to tumor evolution . mutations in histone and dna - modifying enzymes , and even histone proteins themselves , have been found in several cancers and are likely to explain at least some of the observed epigenetic instability ( timp and feinberg , 2013 ) . however , it is not clear that mutations in histone - modifying enzymes account for all the alterations observed in different cancer types . we recently provided evidence that deficiencies in enzymes responsible for replicating g quadruplex ( g4 ) structures , such as the specialized dna polymerase rev1 and helicases fancj , wrn , and blm , can lead to localized changes in histone modifications and gene expression ( sarkies et  al . , 2010 , sarkies et  al . , 2012 , schiavone et  al . , 2014 ) . g4s can form within motifs comprising four short runs of dg bases , separated by linker sequences . the dg bases in the motif form planar hoogsteen - bonded quartet structures that can stack on top of each other , resulting in an often highly thermodynamically stable secondary structure , the g4 ( reviewed in maizels and gray , 2013 ) . we proposed that persistent replication fork stalling at g4s in mutants such as rev1 or fancj leads to pathologically long daughter strand gap formation , resulting in local uncoupling of dna synthesis from parental histone recycling . this , in turn , leads to loss of the histone modifications present on the parental chromatin , which , if in the vicinity of a gene promoter , results in changes in transcription ( sarkies et  al . , 2010 , sarkies et  al . , 2012 , schiavone et  al . , a prediction of this model is that global replication stressors that lead to loss of processive dna polymerization with uncoupling of the replicative helicase and polymerase also should promote epigenetic instability by dissociating dna synthesis from histone recycling . here we test this hypothesis by examining the effect of hydroxyurea ( hu)-induced nucleotide pool depletion on the epigenetic stability of a sensitive reporter locus , bu-1 , in chicken dt40 cells ( sarkies et  al . we show that chronic treatment with low - dose hu induces stochastic instability of bu-1 expression , characterized by loss of the chromatin marks h3k4me3 and h3k9/14ac seen in the normally active locus . this instability depends significantly on the presence of a g4 motif 3 of the promoter , oriented to stall the leading strand of a fork heading toward the transcription start site ( tss ) . the presence of this g4 motif not only increases the rate at which bu-1 expression is lost , but is additionally associated with phosphorylation of h2ax and appearance of the heterochromatic mark h3k9me3 . this is consistent with the g4 acting to focus dna damage induced by the global replication stress imposed by hu , with the damage leading to repression of the locus . further , we show that , across the genome , chronic exposure to hu results in an altered pattern of gene expression similar to that seen in cells lacking the g4-unwinding helicases fancj , wrn , and blm , and that affected genes are enriched in g4 motifs . together , these observations indicate that nucleotide depletion can combine with naturally occurring dna secondary structures to promote epigenetic instability . we first sought conditions in which we could culture dt40 cells in low - dose hu such that replication is slowed but completed ( alvino et  al . , 2007 ) . we therefore exposed wild - type dt40 cells to a range of hu concentrations and monitored their doubling time . the cells were able to proliferate for over a week in up to 150  m hu ( figure  1a ) . at this dose , their doubling time increased from 12.3 to 32.7  hr , recovering when the hu was washed out ( figure  1a ) . to determine the effect of low - dose hu on replication dynamics , we performed dna molecular combing after pulse labeling the cells with halogenated nucleotides ( figure  s1a ) 3  days after initiating culture in hu . average fork velocity decreased from 1.26 to 0.71 kb / min ( figure  1b ) , with a compensatory decrease in average interorigin distance from 72 to 40 kb ( figure  1b ) . consistent with these perturbed replication dynamics , cell - cycle profiles revealed a significant accumulation of cells in s phase while in hu ( figure  s1b ) . we have reported previously that replication - dependent transcriptional instability associated with g4 motifs can be monitored by following expression of a surface marker , bu-1a , in dt40 cells ( sarkies et  al . , 2012 , schiavone et  al . , 2014 ) . the bu-1 locus contains prominent g4 motifs 3.5 kb downstream of the tss and 3 kb upstream . both are orientated to be g - rich on the feature strand with respect to the bu-1 transcript ( figure  1c ) . epigenetic instability of bu-1 in rev1 cells is entirely dependent on the  + 3.5 g4 motif , and it requires the motif to be orientated such that its g - rich strand forms on the leading strand of a replication fork entering the locus from the 3 end ( figure  1c ; schiavone et  al . , 2014 ) . we have reported previously that the bu-1 locus is bidirectionally replicated , meaning that during any given s phase there is a 50% probability of the  + 3.5 g4 being replicated on the leading strand template ( schiavone et  al . , 2014 ) . growth of wild - type dt40 cells in 150  m hu resulted in the appearance of a bu-1a population as cells divided over the course of 7  days ( figure  1d ) . surface expression of bu-1a correlates closely with transcript abundance ( sarkies et  al . , 2012 ) , and this held true for bu-1a clones recovered after hu treatment ( figure  s1c ) . to estimate the rate at which bu-1a variants are formed in hu , we performed a fluctuation analysis by expanding multiple parallel populations of 10 bu-1a cells in hu for 7  days , after which we monitored the appearance of bu-1a variants . this revealed a striking degree of expression instability despite the small number of cell cycles through which the cells had passed ( figure  1e ) . using our previously described monte carlo simulation of bu-1a loss as a replication - dependent phenomenon ( schiavone et  al . , 2014 ) , we estimated a per - division probability of generating of a bu-1a state during culture in hu of c. 0.15 . to obtain additional evidence that this induced transcriptional instability of bu-1 reflected decreased dna polymerase processivity , we asked whether bu-1 variants could be induced by aphidicolin . aphidicolin slows replication by directly inhibiting dna polymerases , particularly pol ( oguro et  al . , 1979 ) , and a low dose induces replication stress ( pacek et  al . , 2006 ) . dt40 cells were able to proliferate in up to 150  m aphidicolin for 10  days and , as with low - dose hu , this resulted in substantial instability of bu-1a expression ( figure  1f ) . we have shown previously that removal of the  + 3.5 g4 motif from both alleles of bu-1 in rev1-deficient cells results in complete stabilization of expression of the locus ( schiavone et  al . , 2014 ) . we therefore examined the extent to which this motif also accounted for the observed hu - induced instability of bu-1a expression in wild - type dt40 . we grew wild - type cells lacking the  + 3.5 g4 on both alleles , bu-1 ( schiavone et  al . , 2014 ) , in hu and assessed the frequency of bu-1a variants after 7  days by fluctuation analysis ( see figure  e3 in schiavone et  al this revealed that removal of the  + 3.5 g4 motif resulted in a significant reduction in the rate at which bu-1a variants were generated . however , it did not result in complete stabilization of the locus ( figures 2a and 2b , i and ii ) . we considered the possibility that the residual instability could be due to the 3.0 g4 upstream of the tss . however , deleting this motif ( figure  s2 ) had no impact on hu - induced instability ( figure  2b , iii ) . to confirm the contribution of the  + 3.5 g4 motif , this resulted in the return of the high - level hu - induced instability of bu-1 expression observed in wild - type cells ( figure  2c , i ) . however , this was not seen if the motif was mutated to render it incapable of forming an intramolecular g4 ( figure  2c , ii ) or when it was inverted so that the g4 structure would form on the lagging strand template ( figure  2c , iii ) . thus , hu treatment alone can induce instability of bu-1a expression , but its effect is significantly potentiated by the presence of a g4-forming sequence orientated to stall the leading strand replication of a fork heading toward the tss . we next investigated the basis for the hu - induced generation of  bu-1a variants . to test whether the bu-1 state is permanent , we isolated five bu-1 clones at the end of 1-week growth in hu and cultured them for a further 3  weeks in hu - free medium . the clones remained stably bu-1 with no evidence of reversion to bu-1 , suggesting that this was a permanent change . we considered the possibility that  bu-1 cells resulted from genetic changes in the locus , although the observed rate of mutation would be extraordinarily high for this to be the case . we therefore sequenced the region around the  + 3.5g4 to look for mutation of the motif and used pcr with restriction digestion to detect larger deletions ( figure  s3 ) . neither assay revealed any evidence of genetic instability consistent with the formation of bu-1 variants being an epigenetic event . we therefore examined the pattern of histone modification at the bu-1 promoter by chromatin immunoprecipitation ( chip ) from bulk populations of cells exposed to hu . bu-1 is a transcriptionally active locus characterized by high levels of h3k4me3 around its tss . after 48-hr treatment with hu , we observed a small but not significant loss of h3k4me3 at the bu-1 promoter ( figure  3a ) , consistent with the size of the population of bu-1a cells generated by this time point ( figure  1d ) . however , after 7-day treatment , we observed a more significant loss of h3k4me3 correlating with the much larger population of bu-1a cells at this time point ( figures 1d and 2a ) . the loss of  h3k4me3 was accompanied by a reduction in h3k9/14 acetylation ( figure  s4 ) . seven days of hu treatment also induced a marked increase in h3k9me3 at the bu-1 promoter ( figure  3a ) . it  has been proposed previously that hu - induced displacement of parental h3/4 and its buffering by the histone chaperone asf1 may lead to unscheduled heterochromatinization by ectopic deposition of pre - marked histones upon their release from asf1 ( jasencakova et  al . , 2010 , schwab et  al . , 2013 ) . alternatively , the appearance of h3k9me3 may result from dna damage - induced heterochromatinization , which has been observed following double - strand breaks ( ayrapetov et  al . breaks can arise from fork collapse in hu ( petermann et  al . , 2010 ) , and thus , if unscheduled incorporation of h3 with k9 methylation was responsible , then an increase in h3k9me3 would be observed irrespective of whether the  + 3.5 g4 motif was present . however , if localized g4-induced dna damage was responsible , then the appearance of h3k9me3 would be dependent on the  + 3.5 g4 motif . we therefore examined h3k4me3 and h3k9me3 at the promoter of bu-1 in cells lacking the  + 3.5 g4 motif . after 7  days in hu , h3k4me3 was reduced ( figure  3b ) , but to a lesser extent than in wild - type cells ( figure  3a ) , consistent with the reduced rate at which bu-1a variants are generated in cells lacking the  + 3.5 g4 motif ( figure  2b ) . however , we observed no associated increase in  h3k9me3 ( figure  3b ) . to monitor the extent to which hu induced dna damage in the two situations , we performed chip for phosphorylated h2ax ( h2ax ) ( rogakou et  al . , 1998 ) at the bu-1 promoter . h2ax was enriched 2.5-fold at the bu-1 promoter in wild - type cells after 7  days in hu , but not enriched in cells lacking the  + 3.5 g4 motif grown under the same conditions ( figure  3c ) . this favors a model in which heterochromatinization of bu-1 in hu is promoted by dna damage , likely from fork collapse associated with the  + 3.5 g4 motif . we next examined the extent to which hu - induced changes in histone modifications were permanent by performing chip at the bu-1 promoter in the five stable bu-1a clones discussed above . this revealed that promoter h3k4me3 remained low , showing that loss of this mark was permanent ( figure  3d ) . thus , while h3k9me3 is induced by hu in cells containing the  + 3.5 g4 motif , this mark is not essential to maintain the bu-1a state . this may be because it is installed only transiently during repair of hu - induced dna damage in the locus , or because cells in which h3k9me3 persists are growth disadvantaged and are lost from the population . the data thus far were consistent with a working hypothesis that reduced polymerase processivity increases the probability of g4 formation at the  + 3.5 g4 motif through exposure of more single - stranded dna within the replisome , which , in turn , focuses replication stalling at this site . implicit in this model is the idea that the  + 3.5 g4 can form during replication but that it is usually rapidly resolved to maintain fork progression . we therefore reasoned that trapping the g4 structure using a g4-binding ligand also might induce instability of bu-1a expression in otherwise wild - type cells further , we predicted that g4 ligands and hu would act synergistically to destabilize expression of the locus . to test these ideas , we treated cells with the g4 ligand n - methyl mesoporphyrin ix ( nmm ) ( nicoludis et  al . , 2012 ) . we first identified the maximum dose at which the cells retained normal viability and global replication dynamics . at 2  m nmm , the fork rate , as assessed by molecular combing , was 1.13  kb / min compared with 1.26 kb / min in wild - type cells , with no significant change in the inter - origin distance ( figure  s5 ) . nonetheless , fluctuation analysis for bu-1a loss in wild - type and bu-1 cells cultured for 7  days in 2  m nmm revealed instability in bu-1a expression in wild - type cells , but not cells lacking the  + 3.5 g4 motif ( figure  4a ) . since 2  m nmm does not in itself significantly reduce global fork rates and the agent will only interact with the formed g4 structure , not with just the  linear dna sequence ( ren and chaires , 1999 ) , this observation is consistent with transient formation of g4s during normal replication . we next asked whether combining nmm - induced g4 stabilization with hu - induced reduction in polymerase processivity led to a further destabilization of bu-1 expression . interestingly , use of both drugs together resulted in significant toxicity , meaning that we had to reduce the dose of each drug by 50% in order to carry out the fluctuation analysis . as expected , growth of cells in hu at 75  m or nmm at 1  m individually had little effect on stability of bu-1a expression ( figure  4b ) . however , the combination of nmm and hu at these doses resulted in a significant increase in bu-1a instability , revealing a marked synergy between hu - induced replication stress and g4 stabilization . we replaced the natural  + 3.5 g4 motif with a series of four g4 motifs of varying in  vitro thermal stabilities ( schiavone et  al . , 2014 ) . all four motifs ( g4 14 ) potentiated the formation of bu-1a variants upon treatment with hu ( figure  5a ) . interestingly , we observed no correlation between the degree of potentiation by the motifs and the in  vitro melting temperature of the equivalent oligonucleotides ( figure  5b ) . however , there was a significant trend toward greater potentiation of bu-1a loss being associated with longer non - g loops in the range of 1 to 9  bp ( figure  5c , solid line ) . to explore this further , we also tested a single repeat of the g4 motif containing human ceb1  mini - satellite ( piazza et  al . , 2012 ) , which has 18  bp between its first three and last run of dgs . this g4 motif , but not a mutated form that is incapable of forming a g4 structure in  vitro ( piazza et  al . , 2012 ) , also potentiated bu-1a instability after treatment with hu . however , this was not to a greater extent than the natural  + 3.5 g4 dna with its central 9-bp loop , suggesting that there may be a limit after which lengthening the loop has no further effect . finally , we asked whether we could detect genome - wide evidence of an interaction between hu and g4s . we therefore performed affymetrix expression microarray analysis on cells before and after culture in 150  m hu . three parallel cultures of dt40 were treated with 150  m hu for 7  days , or mock treated , and then recovered into normal medium for 7  days , after which rna was prepared for array hybridization . despite this relatively short treatment , a total of 2,937 of 12,920 unique genes exhibited a change in expression of > 0.25 log2 units with p  < 0.05 , with an approximately equal number of genes being upregulated and downregulated ( figure  6a ) . we previously have observed a similarly large number of dysregulated genes in cells deficient in the 5-3 g4-unwinding helicase fancj and in double mutants for the 3-5 helicases wrn and blm ( sarkies et  al . , further , we found a highly significant overlap in the identity of dysregulated genes in the two sets , the direction in which their expression changed , and the association of the dysregulated genes with g4 motifs ( sarkies et  al . we anticipated that if transcriptional dysregulation by hu was linked with g4s that there might be a significant similarity in the gene set  altered by hu and the sets altered by loss of fancj and wrn / blm . indeed , the overlap in the identities of the genes dysregulated in all three conditions was highly significant ( figure  6b ) , as were the pairwise correlations in the direction of the change in expression ( figure  6c ) . nearly 68% of the 6,061 genes within the venn diagram in figure  6b have a g4 motif within 1 kb upstream of the tss and the end of the body of the gene in comparison with 59% of the 6,859 genes in the remainder of the array ( p  < 1   10 ) ( table s1 ) . to ascertain whether the overlaps in the identity of dysregulated genes reflected the perturbation of common pathways in the three datasets , we analyzed the functional annotation terms associated with the genes in each set and in the overlap sets using david ( https://david.ncifcrf.gov ; huang et  al . while treatment with hu resulted in dysregulation of genes with gene ontology ( go ) terms associated with cellular stress and nucleotide metabolism , a large number of miscellaneous go terms also were enriched to a similar degree ( table s2 ) . significantly , despite the large number of genes overlapping in the three datasets , there was no evidence of their being members of common pathways ( figure  s6 ) . this is consistent with much of the dysregulation of expression resulting from processes that are not related to a coordinated physiological response either to treatment with hu or ablation of fancj or wrn and blm helicases . nonetheless , the degree of overlap in the dysregulated transcriptomes in these three conditions suggest that cells treated with hu and cells lacking fancj and wrn / blm face similar challenges . however , the enrichment of g4s in affected genes , while statistically significant , is still relatively modest , suggesting that other factors , such as secondary effects or other dna secondary structures , may be contributing as well . numerous lines of evidence have linked replication stress with  genetic instability ( halazonetis et  al . , 2008 , zeman and cimprich , 2014 ) . imbalanced or depleted nucleotide pools during replication are an important cause of such stress and can arise from the expression of oncogenes uncoupling entry into s phase from upregulation of nucleotide supply ( bester et  al . , 2011 ) . importantly , the dna damage resulting from replication stressors like hu or aphidicolin that act directly on the replicative  dna polymerases is not randomly distributed across the genome , but is instead focused on sites that often have features that make them potentially problematic to replicate even under ideal conditions ( tsantoulis et  al . , 2008 ) . many of these hotspots also correspond to classical fragile sites , in which chromosome breaks are observed after replication stress . such sites have been linked to regions depleted in replication origins , meaning that single forks have to traverse long distances ( letessier et  al . , 2011 ) . thus , a combination of regions of low fork density and problematic structures may focus sites of fork collapse under conditions of global replication stress ( wickramasinghe et  al . , 2015 ) . the mechanisms by which replication stress leads to epigenetic  changes are less well explored . alterations in chromatin composition and structure are common features of cancer cells ( berdasco and esteller , 2010 , hansen et  al . , 2011 , timp and feinberg , 2013 ) and are particularly associated with g4-dense breakpoint hotspots ( de and michor , 2011 ) . although the cell line we used in this study , dt40 , is itself transformed , we found no evidence of significant stress to the dna replication program under normal growth conditions . however , growth of the cells in  low - dose hu recapitulated the key features of the acutely stressed replication observed in oncogene - expressing primary cells ( bester et  al . we have been able to explore directly the interaction between global replication stress induced by nucleotide depletion and a dna secondary structure to demonstrate how they conspire to exacerbate replication - dependent epigenetic instability . we have provided evidence that two parallel epigenetic perturbations contribute to permanent and transient epigenetic changes following an episode of replication stress . the first mechanism relates to the uncoupling of the activity of replicative helicase and polymerase ( byun et  al . , 2005 , pacek and walter , 2004 ) during hu treatment . this has been shown to lead to interruption of the normal flow of histones from the parental to the nascent daughter strands , with the histone chaperone asf1 buffering the displaced h3/h4 ( jasencakova et  al . , 2010 ) . groth and colleagues suggested that their release from asf1 might lead to local alterations in epigenetic state of chromatin due to unscheduled incorporation of inappropriately marked histones ( jasencakova and groth , 2010 , jasencakova et  al . , 2010 ) . ( 2013 ) invoked this model to explain an increase in heterochromatin formation in cells deficient in fancj , suggesting that failure to unwind lagging - strand template g4 structures in fancj - deficient cells led to unscheduled deposition of histones bearing marks that would lead to h3k9me3 and heterochromatin formation . however , this model does not adequately explain the bidirectional changes in gene expression changes seen either in fancj cells or in wild - type cells  treated with hu ( sarkies et  al . in contrast , the model we have developed previously , in which loss of processive replication at g4s leads to localized loss of parental histone mark recycling , could explain both derepression of loci , such as -globin ( sarkies et  al . , 2010 ) , and loss of activation ,  as can be observed in the bu-1 locus ( sarkies et  al . however , since the mechanism by which h3k4me3 is maintained during replication is poorly understood , the precise mechanisms by which replication impediments disturb the maintenance of this mark remain to be fully elucidated this mark is induced alongside the loss of h3k4me3 when the  + 3.5 g4 motif is present and cells are exposed to hu . importantly , the appearance of this mark of heterochromatin is accompanied by h2ax phosphorylation , a marker of dna damage ( rogakou et  al . , 1998 ) . g4 motifs have been linked to hotspots of genetic instability and translocation ( de and michor , 2011 ) . further , the g4 ligand pyridostatin , which acts similarly to nmm , leads to localized h2ax accumulation at g4 motifs across the genome , suggesting the formation of dna breaks ( rodriguez et  al . , dna breaks have been shown to induce transcriptional repression ( shanbhag et  al . , 2010 ) and to induce h3k9me3 even in a normally euchromatic locus ( ayrapetov et  al . , 2014 ) . thus , the appearance of h2ax and h3k9me3 in the bu-1 locus only in hu - treated cells containing the  + 3.5 g4 is consistent with collapse or incision of replication forks , already stressed by nucleotide depletion , that have stalled at the g4 . we therefore propose that nucleotide depletion can give rise to loss of parental h3k4me3 and the appearance of h3k9me3 by distinct mechanisms . h3k4me3 is lost stochastically as a result of interruption of parental histone recycling , a mechanism that is locally exacerbated by the presence of a g4 motif . in contrast , we propose that h3k9me3 may reflect protective transient heterochromatinization of the locus during repair of breaks resulting from hu - induced fork collapse at g4 structures ( figure  7 ) . in the case of loci like bu-1 , in which hu - induced epigenetic instability is linked to g4 formation , an interesting question is whether hu results in a greater opportunity for g4 formation during replication or diminished g4 resolution . it is not currently possible to formally distinguish these possibilities and indeed it is likely that elements of both are true . notably , the observation that the g4-binding ligand nmm can induce bu-1a expression instability at a dose that does not significantly impact on global replication dynamics provides strong evidence that g4 structures form readily during normal replication and that they are usually promptly resolved . the synergy between replication stress caused by nucleotide depletion and structured dna is of considerable potential importance to understanding the development of cancer . epigenetic changes are prevalent in many cancer types , although their origin is unclear , and likely complex ( berdasco and esteller , 2010 , timp and feinberg , 2013 ) . recently , several instances of epigenetic instability in cancer have been linked to mutations in histone or dna - modifying enzymes . however , the widespread and often apparently random nature of epigenetic changes in tumors suggests that other processes also may be at work . we have suggested previously that delayed replication of g  quadruplex structures could contribute to the epigenetic diversity of cancer ( sarkies and sale , 2012 ) . however , mutations in enzymes that may cause this form of epigenetic instability , for example rev1 , fancj , wrn , and blm ( sarkies et  al . , 2010 , sarkies et  al . , 2012 ) , are rarely observed in sporadic cancers . replication stress , on the other hand , is emerging as an important feature of cancer cells , particularly in the early stages of their evolution ( bester et  al . , 2011 , di micco et  al . , 2006 , thus , we suggest that some of the epigenetic changes seen in tumors may be explained by problems managing replication blocks . consistent with this idea , both copy number variations and changes in dna methylation patterns in cancer have been linked to g4 motifs ( de and michor , 2011 ) . finally , it is worth noting that hu is used extensively in treatment of hemoglobinopathies , such as sickle cell disease and thalassaemia , as it can re - induce expression of the fetal -globin gene , ameliorating the effects of the defective adult globins found in these disorders ( platt et  al . , 1984 ) . importantly , the effect of hu on -globin expression is unlikely to be specific , since chronic exposure to the drug leads to quite widespread changes in erythroid gene expression ( flanagan et  al . , 2012 ) . although the -globin locus in humans has no g4 motifs in the immediate vicinity of its promoter , its key transcriptional regulator , bcl11a ( bauer et  al . , 2013 ) , has a high density of g4 motifs on both sides of its tss . it will , therefore , be interesting to explore whether the mechanisms we propose here could help explain the action of hu on fetal globin expression . dt40 cells culture and the strategy for removing and replacing the  + 3.5 g4 motif in the bu-1 locus have been described previously ( schiavone et  al . , 2014 ) . genetic manipulation of the  + 3.5 g4 motif was performed in the bu-1a allele of cells in which the motif had been removed from the bu-1b allele to avoid the transvection - like effect between the alleles ( schiavone et  al . , 2014 ) . oligonucleotides are listed in table s3 . for fluctuation analysis , 150  m hu ( sigma - aldrich , h8627 ) after 7  days , cells at a concentration between 0.2 and 1   10 were stained for 20  min at 37c with anti - bu-1a - phycoerythrin ( 1:100 , santa cruz biotechnology clone 5k98 , 70447 ) . bu-1a expression was assessed by flow cytometry using an lsrii cytometer ( becton dickinson ) . experiments with aphidicolin ( sigma - aldrich , a0781 ) and nmm ( frontier scientific , nmm580 ) were conducted in 96-well plates starting with ten cells expanded for 10  days . bu-1a cells were isolated after hu treatment using a moflo sorting cytometer ( dako - cytomation ) . chip was performed as described previously ( nelson et  al . , 2006 ) with modifications . following a 10-min incubation at room temperature with 1% ( v / v ) formaldehyde , glycine was added to 0.2  m for 5  min . the extracted nuclei were sonicated at 4c using a bioruptor water bath sonicator ( diagenode ) with 30 cycles of 30  s separated by 30-s intervals . sheared chromatin samples were resuspended in dilution buffer ( 1.1% triton x-100 , 1.2  mm edta , 16.7  mm tris [ ph 8.0 ] , 167  mm nacl supplemented with pmsf , and a protease inhibitor cocktail ) . for immunoprecipitation , lysates were incubated overnight with the following antibodies at 4c : histone h3 ( 1:100 , cell signaling technology , 2650 ) , h3k4me3 ( 1:100 , cell signaling technology , 9727 ) , h3k9/14ac ( 1:200 , millipore , 17 - 615 ) , h3k9me3 ( 1:200 , abcam , ab8898 ) , h2ax ( 1:50 , abcam , ab2893 ) , and the negative control normal rabbit igg ( millipore ) . following overnight incubation at 4c with tumbling and four washing steps , the qpcr was performed with power sybr green master mix ( applied biosystems , 4367659 ) on an abi prism real - time cycler with the following cycle times : 50c for 2  min , 90c for 10  min , 45 cycles of 90c for 15  s plus 60c for 1  min . the cdna was made from 5  mg mrna with super rt ( ht biotechnology ) and oligodt primer in a final volume of 40  l . dna molecular combing was conducted 3  days into culture with 150  m hu . it was performed and analyzed as previously described ( guilbaud et  al . , 2011 ) . rna was extracted from three independent wild - type cell populations treated for 7  days with hu and allowed to recover for another 7  days , as well as from three untreated parallel controls . microarray analysis was performed using r ( http://www.r-project.org/ ) and its bioconductor packages ( gentleman et  al . , 2004 ) . raw cel files were processed using the robust multichip average ( rma ) algorithm available in the affy package ( gautier et  al . , 2004 ) . genes that showed a change of > 0.25 log2 units relative to the mean wild - type intensity , with a p value of  < 0.05 ( t test ) , were identified as exhibiting statistically significant transcriptional dysregulation . custom written r  scripts were used to identify and plot genes co - dysregulated between different mutants . venn diagrams were generated with the limma package ( smyth , 2004 ) , and significance for the overlaps was calculated using fisher s hypergeometric distribution .
summarynucleotide pool imbalance has been proposed to drive genetic instability in cancer . here , we show that slowing replication forks by depleting nucleotide pools with hydroxyurea ( hu ) can also give rise to both transient and permanent epigenetic instability of a reporter locus , bu-1 , in dt40 cells . 
 hu induces stochastic formation of bu-1low variants in dividing cells , which have lost the h3k4me3 present in untreated cells . 
 this instability is potentiated by an intragenic g quadruplex , which also promotes local h2ax phosphorylation and transient heterochromatinization 
 . genome - wide , gene expression changes induced by hu significantly overlap with those resulting from loss of the g4-helicases fancj , wrn , and blm . 
 thus , the effects of global replication stress induced by nucleotide pool depletion can be focused by local replication impediments caused by g quadruplex formation to induce epigenetic instability and changes in gene expression , a mechanism that may contribute to selectable transcriptional changes in cancer .