- Research article
- Open Access
The essential Schizosaccharomyces pombe Pfh1 DNA helicase promotes fork movement past G-quadruplex motifs to prevent DNA damage
BMC Biology volume 12, Article number: 101 (2014)
G-quadruplexes (G4s) are stable non-canonical DNA secondary structures consisting of stacked arrays of four guanines, each held together by Hoogsteen hydrogen bonds. Sequences with the ability to form these structures in vitro, G4 motifs, are found throughout bacterial and eukaryotic genomes. The budding yeast Pif1 DNA helicase, as well as several bacterial Pif1 family helicases, unwind G4 structures robustly in vitro and suppress G4-induced DNA damage in S. cerevisiae in vivo.
We determined the genomic distribution and evolutionary conservation of G4 motifs in four fission yeast species and investigated the relationship between G4 motifs and Pfh1, the sole S. pombe Pif1 family helicase. Using chromatin immunoprecipitation combined with deep sequencing, we found that many G4 motifs in the S. pombe genome were associated with Pfh1. Cells depleted of Pfh1 had increased fork pausing and DNA damage near G4 motifs, as indicated by high DNA polymerase occupancy and phosphorylated histone H2A, respectively. In general, G4 motifs were underrepresented in genes. However, Pfh1-associated G4 motifs were located on the transcribed strand of highly transcribed genes significantly more often than expected, suggesting that Pfh1 has a function in replication or transcription at these sites.
In the absence of functional Pfh1, unresolved G4 structures cause fork pausing and DNA damage of the sort associated with human tumors.
DNA helicases are essential for genome stability. They have critical roles in DNA replication, repair and recombination. Multiple human hereditary disorders are linked to mutations in helicase genes. For example, mutations in three of the five human RecQ helicases are associated with increased cancer risk and/or premature aging . A point mutation in the human PIF1 (hPIF1) DNA helicase (hPIF1 L319P) is present in certain families with increased risk of breast cancer and not detected in unaffected controls . This mutation changes a conserved residue within the 21-amino acid Pif1 signature motif that characterizes this family of DNA helicases ,.
Pif1 family helicases are found in the genomes of organisms from all three kingdoms ,. Most eukaryotes, including Schizosaccharomyces pombe and humans, encode a single Pif1 family helicase, while Saccharomyces cerevisiae encodes two, ScPif1 and ScRrm3. Pfh1, the S. pombe Pif1 helicase, is essential for maintenance of both the nuclear and mitochondrial genomes . In nuclear DNA, it facilitates replication fork progression through many sites that impede fork progression, such as highly transcribed RNA polymerase II and III genes, replication fork barriers within both the ribosomal DNA (rDNA) and the mating-type locus, and converged replication forks ,. At the mating type locus and rDNA, Pfh1 helps forks move past stable protein complexes. In the absence of Pfh1, double strand breaks (DSBs) occur specifically at these natural fork impediments .
Pfh1’s role at a different class of hard to replicate sites, stable non-canonical DNA secondary structures such as the G-quadruplex (G4), has not been systematically explored. G4 structures are stable DNA secondary structures held together by multiple stacked guanine quartets . G4 structures can form within a single DNA molecule (intra-strand) or between different DNA molecules (inter-strand). In virtually all genomes examined so far, DNA sequences that are capable of forming intra-strand G4 structures in vitro (G4 motifs) are observed. G4 motifs are highly enriched in G-rich telomeric DNA, where they affect telomerase action and end protection ,. In addition to telomeres, there are more than 300,000 sites in the human genome with the potential to form G4 DNA ,, and G4 structures can be detected in human cultured cells with G4 specific antibodies ,. Moreover, G4 structures were more frequent in FANCJ helicase-deficient human cultured cells . In bacteria, budding yeast, and humans, G4 motifs are common in rDNA and promoter regions -.
Although G4 motifs are not frequent sites of DNA damage in wild type (WT) S. cerevisiae, in pif1 mutant cells, replication forks slow and often break at these sites . In WT cells, ScPif1 binds a subset of G4 motifs, and this subset is more likely to be associated with fork slowing and DNA breakage in its absence. Moreover, in cells lacking ScPif1, G4 motifs induce gross-chromosomal rearrangements (GCRs) ,. Although G4 motifs do not induce genome instability in rrm3 cells, G4-induced GCR events are particularly elevated in pif1 rrm3 cells, suggesting that ScRrm3 acts as a backup for ScPif1 in suppressing G4-induced DNA damage. Consistent with a role for ScPif1 at G4 motifs in vivo, ScPif1 and four of four tested bacterial Pif1 helicases are particularly robust unwinders of G4 structures, even under single-cycle conditions . In human cells, some of the binding sites of the G4 stabilizing agent pyridostatin, co-localize with the binding of hPIF1, suggesting a role of hPIF1 in resolving G4 DNA .
Here, we investigated the relationship between G4 motifs and the S. pombe Pfh1 helicase. Compared to S. cerevisiae, S. pombe has a chromosome structure more similar to that of higher eukaryotes and, like human cells, encodes a single Pif1 family helicase. Thus, it is a good model for understanding the functions of hPIF1. We used computational methods to map intra-strand G4 motifs within the genome of S. pombe and three other sequenced Schizosaccaromyces yeasts to determine their association with genomic features. G4 motifs were significantly enriched in rDNA, telomeres, meiotic DSB hot spots, gene promoters, nucleosome-depleted regions (NDRs), untranslated regions (UTRs), and dubious open reading frames (ORFs), but depleted in ORFs. Using chromatin immunoprecipitation in combination with deep sequencing (ChIP-seq), we found that Pfh1 was bound near approximately 20% of the G4 motifs in the assembled S. pombe nuclear genome and that fork slowing and DNA damage, as indicated by association with Cdc20, the leading strand DNA polymerase, and phosphorylated H2A (γ-H2A), respectively, were associated with G4 motifs in Pfh1-depleted cells. Together, our data suggest that Pfh1 is needed to unwind G4 structures; when this unwinding does not occur, forks slow and often break. This increased genome instability in Pfh1-depleted cells could explain the association of hPIF1 mutations with cancer.
Identifying G4 motifs in fission yeasts
We performed a genome-wide search for DNA sequences with the potential to form G4 structures on the genomes of the four available Schizosaccaromyces species . We identified all sequences that contain four runs of three or more guanine base pairs (bp), ‘G-islands’, separated by ‘loop’ regions of no more than 25 bp (G≥3 N1–25)3 G≥3 . Hereafter, sequences matching this pattern are called ‘G4 motifs’. Regions with more than four G-islands separated by ≤25 bp were counted as a single G4 motif. Excluding repetitive DNA (see below), the S. pombe genome contained 446 G4 motifs that match this query pattern with a density of 0.036 G4 motifs/kilobase (kb) (Figure 1A; Table 1). The density of G4 motifs was similar across all three S. pombe chromosomes (Figure 1A). The S. octosporus and S. cryophilus genomes contained a similar number and density of G4 motifs as S. pombe (Table 1; Additional file 1). The more distantly related S. japonicus had roughly four times the number (1,757) and density (0.16 G4 motifs/kb) of G4 motifs as the other fission yeast species, most likely due to the higher GC content of its genome (44% versus 36% in S. pombe).
The above estimates for the number of G4 motifs do not include G4 motifs within repetitive telomeric DNA and most of the rDNA. Telomeric DNA was not included in this analysis because it is not in the S. pombe genome assembly. However, analysis of S. pombe telomeric DNA in a library containing 18 kb of sequenced telomeric DNA (Webb and Zakian, submitted) revealed that it had two orders of magnitude more G4 motifs per kb (4.5 G4 motifs/kb) than bulk nuclear DNA [see Additional file 2]. The presence of G4 motifs in telomeric DNA is not surprising as the S. pombe telomeric sequence, G2-8TTAC, is GC-rich . However, the density is nearly three times higher than expected in random sequences of the same GC content (1.7 G4 motifs per kb). Assuming 300 bps of telomeric DNA at each of the six chromosome ends, we estimate 1 to 2 G4 motifs per telomere or approximately 10 telomeric G4 motifs per haploid genome.
The S. pombe genome assembly contains only three full rDNA repeats, a small fraction of the estimated approximately 300 rDNA repeats in nuclear DNA . Each of the three copies of the 10.9 kb rDNA repeat had five G4 motifs (Figure 1B). Assuming that these G4 motifs are present in all of the rDNA repeats, there are approximately 1,500 G4 motifs in S. pombe rDNA, accounting for almost 80% of the G4 motifs in nuclear DNA. The density of G4 motifs in the rDNA (0.45 G4/kb) was significantly greater than expected from random sequences of the same GC content (0.08 G4/kb; P = 0.003) and more than ten times higher in rDNA than the average in the nuclear genome (0.45 G4/kb versus 0.036 G4/kb), even though the GC content of the two is similar (38% versus 36% GC content in rDNA versus nuclear genome). Each of the five rDNA G4 motifs was on the non-transcribed strand (Figure 1B): three were in transcribed regions, two in the 28S and one near the start of the 18S rRNA. Of the three G4 motifs in transcribed regions, the one in 18S and one of the two in 28S were conserved in sequence and position between S. pombe and S. japonicus  (sequence data for the rDNA in the other Schizosaccaromyces species are not available). The high density of G4 motifs is also true for S. cerevisiae  and human rDNA . In both S. pombe and S. cerevisiae, G4 motifs are found only on the non-transcribed strand , suggesting a potential role for G4 structures in the rRNA. However, the most highly conserved G4 motif, which is in the 18S rRNA, does not form a G4 structure in an existing crystal structure of the S. cerevisiae ribosome .
We also determined the G4 motif content of the mitochondrial (mt) DNA of the four fission yeast species, which range in size from 19.4 kb in S. pombe to 80.1 kb in S. japonicus (Table 1). No G4 motifs were present in S. pombe or S. octosporus mtDNA, while S. cryophilus and S. japonicus contained one and five motifs, respectively. These results are in sharp contrast to the very high density of G4 motifs in the AT-rich S. cerevisiae mtDNA compared to its nuclear genome (0.37 G4 motifs/kb mtDNA versus 0.055 motifs/kb nuclear DNA) . Hereafter, we will address the events that occur in the 446 G4 motifs found in the published S. pombe genome assembly.
Evolutionary conservation of G4 motifs across fission yeasts
There was little evolutionary conservation of G4 motif locations among the four fission yeast species. In each pairwise combination, only a small number of S. pombe G4 motifs (20 to 31) overlapped a G4 motif in the aligned homologous location in another species (Table 1). Moreover, different motifs were maintained between different pairs of species. Excluding the rDNA repeats, only five motifs were conserved between S. pombe, S. octosporus and S. cryophilius, and only one location, the promoter of cdc13 +, had a G4 motif in all four species.
This low level of evolutionary conservation is not surprising due to the evolutionary divergence of the available fission yeast genomes. S. octosporus, S. cryophilus and S. pombe diverged from their last common ancestor more than 100 million years ago (mya), while S. japonicus diverged more than 200 mya . Studies of the evolutionary turnover of G4 motifs  and other regulatory elements - in yeasts find that most regulatory elements are not conserved over these timescales.
The genomic distribution of G4 motifs
To investigate potential functions for S. pombe G4 motifs, we analyzed the distribution of the G4 motifs with respect to multiple genomic features, such as highly expressed genes, NDRs, meiotic DSB hotspots and so on. (See the Methods for the full list). For these analyses, we first computed the number of overlaps between the G4 motifs and a given feature of interest. Then, to evaluate if the observed association was more or less than expected by chance, we created 1,000 sets of ‘control’ regions. Each of these 1,000 control sets contained 446 random genomic regions—one for each of the actual G4 motifs. Each of the 446 regions in a single control set matched the chromosome, length and GC content of a different observed G4 motif, so the average region length and GC content for each of the 1,000 sets matched that of the actual G4 motifs. Then, for each of these 1,000 sets of control regions, we computed the number of overlaps with the feature of interest. By comparing the number of observed G4 motif overlaps to the 1,000 overlap counts from the control sets, we obtained an empirical estimate of the likelihood of the observed association by chance (that is, a P-value). To account for the testing of multiple hypotheses, for each enrichment test, we report q-values, which are the false discovery rate (FDR) analogue of P-values and correspond to the FDR if a particular test is called significant . See Methods for more information.
Using this approach, we found that G4 motifs were significantly associated with several genomic features (Table 2; Additional file 3). G4 motifs were more likely to occur in the promoters of RNA polymerase II-transcribed genes (q <0.003), NDRs (q <0.003), meiotic DSB hot spots (q <0.003), 3′ and 5′ UTRs (q = 0.013 and q = 0.003), and within dubious genes (q <0.003) than expected by chance. In contrast, G4 motifs were significantly depleted from ORFs of protein-coding genes (q <0.004), including essential and highly transcribed genes. However, when G4 motifs were found within ORFs, they were significantly more likely to occur on the transcribed strand than expected by chance (219/303, 72%; P <1E-12, binomial test). In these cases, the G4 motif would not be present in the mRNA. Thus, any function of these G4 motifs would likely be carried out in DNA. G4 motifs were not significantly associated with long terminal repeats (LTRs), tRNA genes, 5S rRNA genes, origins of replication or centromeres (q >0.05 for all; Additional file 3). Similar association patterns were found when decreasing the loop length from 25 bp to 12 bp (data not shown).
Many G4 motifs are Pfh1 associated
Pfh1 is a replisome component that moves with the leading strand polymerase ɛ (Sabouri et al., in preparation). Thus, we anticipated that if some (or all) G4 motifs slow DNA replication, even in wild type cells, they would have higher binding by both Pfh1 and Cdc20, the catalytic subunit of DNA polymerase ɛ. If Pfh1 promotes replication fork progression past G4 motifs, the Cdc20 association at these sites should be even higher in Pfh1-depleted cells. To test these hypotheses, we used a S. pombe strain expressing epitope-tagged Pfh1, isolated Pfh1-associated DNA by immunoprecipitation (ChIP), and then sequenced the associated DNA (ChIP-seq). The input DNA for the ChIP was sequenced as a control. Sites with significant Pfh1 occupancy were identified with the Model-based Analysis of ChIP-seq peak calling software (MACS; Zhang et al. ), using a stringent cutoff for both ChIP and input DNA (P <10-5). The same strategy was used in all of the ChIP-seq analyses in this paper (Methods section). With these methods, we identified 621 high confidence Pfh1-associated sites in DNA from asynchronously growing cells. Two of these peaks mapped to the rDNA, although not to the G4 motifs in the rDNA (Figure 1B). The assembled S. pombe genome lacks telomeric DNA, so this analysis did not assess Pfh1 association with telomeres. However, ChIP-qPCR shows that telomeres were also Pfh1 associated .
Of the 621 Pfh1 peaks, 76 (12%) were ≤300 bp, the shearing size of the ChIP DNA, from a G4 motif. Several Pfh1 peaks were associated with more than one G4 motif, so in total, 90 (20%) of the 446 G4 motifs in the assembled nuclear genome were Pfh1 associated (Figure 1A). The observed association between Pfh1 and G4 motifs was significantly greater than expected by chance (P = 0.002). However, it was significantly lower than expected when taking the GC-content of the G4 motifs into account using our control region sets (P = 0.016). We also validated the ChIP-seq peaks by quantitative PCR (ChIP-qPCR). We compared the association of Pfh1-13Myc or an untagged otherwise isogenic control strain to a tRNA gene (tRNA glu.05), a previously known Pfh1-binding site , three GC-rich sites and three G4 motifs (see Methods for details). Pfh1 was significantly associated with all these sites compared to the control strain (Figure 1C). Together, these results suggest that Pfh1 is present not only at many G4 motifs but also at many other sites, especially at other GC-rich sequences, consistent with its being a multi-functional DNA helicase . This finding is also consistent with the behavior of ScPif1, which binds preferentially to G-rich regions, even those unable to form G4 structures, in vivo and in vitro .
Replication forks pause near G4 motifs in Pfh1-depleted cells
To monitor fork progression at G4 motifs in the presence and absence of Pfh1, we epitope tagged Cdc20 and performed ChIP-seq in WT and Pfh1-depleted cells. Although all sites in the genome are Cdc20-associated at their time of replication, sites where replication forks move slowly are expected to have elevated Cdc20 binding, as seen with DNA Pol2, the catalytic subunit of DNA polymerase ε in S. cerevisiae ,. To deplete cells of Pfh1, an essential protein, we used the thiamine repressible nmt81 promoter . Growth of nmt-Pfh1-GFP cells in thiamine for 12 hours reduces Pfh1 expression, so that Pfh1 is no longer detected by western blot analysis ,. Hereafter, cells treated in this way are referred to as ‘Pfh1-depleted cells’ and untreated cells as ‘WT’/Pfh1-expressing’ cells.
Using these methods, we identified 485 sites of high Cdc20 binding in WT cells, including 50 G4 motifs (Figures 1A and 2A). This number of G4 motifs was not significantly different from the expected association based on the 1,000 control sets (P = 0.529; Figure 2B, blue arrow). Thus, in WT cells, G4 motifs were not enriched among sites of replication fork pausing compared to other GC-rich regions.
There were 517 high confidence Cdc20 peaks in Pfh1-depleted cells, modestly more than the number in WT cells. However, in contrast to WT cells, G4 motifs were highly enriched among the high Cdc20 binding sites in Pfh1-depleted cells (77 of the G4 motifs were associated with a Cdc20 site; Figure 2B, red arrow; P = 0.003). Moreover, G4 motifs associated with Pfh1 in WT cells were much more likely to show high Cdc20 association in Pfh1-depleted cells compared to G4 motifs not associated with Pfh1 (P = 1.9E-11; Table 3). For example, 43% of the Pfh1-associated G4 motifs were also high Cdc20 binding sites in Pfh1-depleted cells. In contrast, only 11% of the G4 motifs not associated with Pfh1 had high Cdc20 binding in Pfh1-depleted cells. This pattern was also true for G4 motifs and Cdc20 peaks from Pfh1-expressing cells (P = 1.5E-17; Table 3). Thus, while G4 motifs as a group were not significantly associated with high Cdc20 in WT cells, those G4 motifs that had high Pfh1 binding were also significantly associated with high Cdc20 binding.
To validate the high Cdc20-binding sites, we performed ChIP-qPCR in WT and Pfh1-depleted cells with the same primer pairs used to confirm Pfh1 association (Figure 2C). At the tRNA gene, and all three G4 motifs, we found significantly higher Cdc20 occupancy in Pfh1-depleted cells compared to WT cells (Figure 2C). Cdc20 occupancy in WT cells and Pfh1-depleted cells did not increase at the three GC-rich regions (Figure 2C).
DNA damage occurs near G4 motifs
Phosphorylation of histone H2A is one of the earliest responses to a DSB. In S. pombe, H2A phosphorylation occurs in an area of approximately 25 kb on either side of a DSB, with the highest peaks around 5 kb from the break site . To determine if fork pausing near G4 motifs resulted in DNA damage, we performed ChIP-seq using antibodies that recognize phosphorylated H2A (γ-H2A) ,. As in the Cdc20 experiments, nmt-Pfh1-GFP cells were grown with or without thiamine for 12 hours and then processed for ChIP-seq and ChIP-qPCR.
We identified 179 γ-H2A peaks in Pfh1-expressing cells and 582 γ-H2A peaks in Pfh1-depleted cells. These peaks were associated with 77 and 177 G4 motifs, respectively (Figure 3A). Even in the presence of Pfh1, G4 motifs were significantly enriched within 5 kb of γ-H2A peaks (Figure 3B, blue arrow and histogram; P = 0.021; 77 associations). This association was even stronger in Pfh1-depleted cells, with 177 G4 motifs with high γ-H2A levels (Figure 3B, red arrow and histogram; P = 0.014). Thus, G4 motifs were near sites of DNA damage in both Pfh1-expressing and Pfh1-depleted cells, but the number of damage-associated G4 motifs was much higher in the absence of Pfh1.
As for Cdc20 occupancy, we validated the γ-H2A peaks with ChIP-qPCR using anti-γ-H2A antibodies in WT and Pfh1-depleted cells (Figure 3C). Both the tRNA gene and the three G4 motifs had significantly increased γ-H2A levels in Pfh1-depleted cells compared to WT cells (Figure 3C). We did not detect a significant increase for the three investigated GC-rich regions (Figure 3C).
To determine if DNA damage was more likely to occur at individual G4 sites in Pfh1-depleted versus WT cells, we compared the P-values for the γ-H2A peaks at a given G4 motif in Pfh1-expressing and Pfh1-depleted cells. We used the same number of ChIP-seq reads from both contexts so that P-values from the two conditions could be compared. G4 motifs without an overlapping peak in a given context were assigned a P-value of 1. The mean P-value decreased from 0.003 in Pfh1-expressing cells to 1E-14 in Pfh1-depleted cells, and the peak P-values in Pfh1-depleted cells were consistently more significant [see Additional file 4; P approximately 0, Wilcoxon signed-rank test]. These findings suggest that the probability of DNA damage at a given G4 motif is higher in Pfh1-depleted versus Pfh1-expressing cells.
In the absence of Pfh1, G4 motifs associated with replication fork stalling are more likely to result in double strand breaks
To combine the observations on G4 motifs, Pfh1 presence, replication fork slowing and DNA damage in Pfh1-depleted cells, we analyzed the association between high Cdc20 occupancy and γ-H2A levels at G4 motifs in Pfh1-expressing and Pfh1-depleted cells. When Pfh1 was expressed, G4 motifs with high Cdc20 occupancy were no more or less likely to be near a site of DNA damage, as marked by γ-H2A, than G4 motifs without Cdc20 (P = 0.23, Fisher’s exact test; Table 4). However, when the same test was performed on data from Pfh1-depleted cells, there was a highly significant association between high Cdc20 occupancy at G4 motifs and nearby damage (P = 8.7E-6; Table 4). These data suggest that in WT cells Pfh1 prevents breakage of forks that pause at G4 sites.
Features of G4 motifs that are Pfh1- and DNA damage-associated
In total, 251 of the 446 G4 motifs in the S. pombe genome were associated with Pfh1, fork pausing (Cdc20) and/or DNA damage (γ-H2A). Based on their relationship with Pfh1, we defined three classes of G4 motifs. The first class contained 90 G4 motifs that were Pfh1 associated in WT cells (Class I). The second class consisted of 100 Pfh1-sensitive G4 motifs; that is, these motifs were sites of fork slowing and/or DNA breakage only in Pfh1-depleted cells (Class II). The third class consisted of 106 G4 motifs that were sites of fork slowing and/or DNA damage in both WT and Pfh1-depleted cells (Class III). By definition, there was no overlap between Class II and Class III, but Class I motifs could also be in either Class II or III.
Nearly 40% of Class I motifs were not sites of genome instability in either WT or Pfh1-depleted cells. This finding suggests that these G4 motifs do not form G4 structures, at least during S phase, or that these motifs are resolved by a different helicase in Pfh1-depleted cells. Only eight (9%) of the Pfh1-sensitive (Class II) motifs were Pfh1-associated (Class I) in WT cells, a surprising finding (see Discussion). In contrast, Pfh1 was detected at 40% of Class III sites, but this binding was not sufficient to prevent damage in WT cells. However, pausing, as monitored by levels of Cdc20 occupancy, was higher at 70% of these sites in Pfh1-depleted cells compared to WT cells (P = 0.015, Additional file 5); that is, the presence of Pfh1 did facilitate replication at many Class III sites.
We explored genomic features associated with the three categories of G4 motifs to see if any attributes distinguished them from each other and those seen when all G4 motifs were considered together. The 90 Class I G4 motifs lacked the associations seen when considering G4 motifs overall, except for being enriched at meiotic DSB hotspots. The only other significant association for this class was overrepresentation within the 500 most highly transcribed RNA polymerase II genes (36 of the 90 motifs; q <0.008). Remarkably, in almost all (88%) of these cases, the G4 motif was on the transcribed strand of the highly transcribed gene. This strand bias was significantly stronger than the bias observed for G4 motifs overall (87% versus 72%; P = 0.0002, Fisher’s exact test). This enrichment was particularly striking given the significant depletion of G4 motifs in ORFs when all G4 motifs were considered (q <0.004; Table 2). Genomic attributes for Class II and III G4 motifs were similar to the patterns observed for G4 motifs overall (Additional file 3).
Although Pif1 family helicases are found in almost all eukaryotes, virtually all in vivo evidence for their role at G4 motifs comes from budding yeast. To determine if the deleterious impact of unresolved G4 structures and the positive role of Pif1 family helicases at these structures holds true in other organisms, we used an integrated computational and experimental approach in S. pombe, an organism that diverged from S. cerevisiae more than a billion years ago. These studies are particularly important because budding yeast is unusual in encoding two Pif1 helicases, ScPif1 and ScRrm3, while most eukaryotes encode only one (reviewed in ). The two S. cerevisiae helicases, ScPif1 and ScRrm3, have multiple, often conflicting, roles in genome integrity, so it is not clear how findings on ScPif1 and ScRrm3 translate to organisms with a single Pif1 helicase. In addition, by multiple criteria, S. pombe chromosomes are more similar to mammalian chromosomes than S. cerevisiae chromosomes . Thus, the sole S. pombe Pif1 family helicase, Pfh1, is a more apt model for the mammalian enzyme.
Consistent patterns in the genomic distribution of G4 motifs across three diverse species
The genomic features associated with G4 motifs in S. pombe were strikingly similar to those seen in the S. cerevisiae and human nuclear genomes ,,, supporting the conservation of G4 motif biology across more than one billion years of evolution. In particular, three functional regions are G4-rich and a fourth is G4-poor in the three organisms. First, promoters of RNA polymerase II transcribed genes in S. pombe contained more G4 motifs than expected (Table 2). Likewise in budding yeast and humans, G4 motifs are enriched within 850 bp and 1 kb, respectively, of the transcriptional start site, suggesting a common regulatory function for G4 motifs ,. This correlation also agrees with our finding that S. pombe G4 motifs were enriched in NDRs, regions that are found in the majority of promoters . Second, as in budding yeast and humans, the S. pombe rDNA was enriched with G4 motifs ,,,. In S. pombe and S. cerevisiae, they were only present on the non-transcribed strand ,. The fact that the G4 motifs are over-represented in rDNA in evolutionarily diverse organisms argues that this arrangement has functional importance. For example, formation of G4 structures in the non-transcribed strand could facilitate high transcription rates by sequestering the transcribed template to prevent re-annealing to the G4-rich complementary strand. Third, telomeres in yeast, human and many other species contain G4 motifs. S. pombe telomeres are no exception; they have the highest density of G4 motifs of any region in the S. pombe genome: 4.5 G4 motifs per kb, which is even higher than expected from the high GC content of telomeric DNA. G4 structures form in vivo in a cell cycle dependent manner in ciliate and human telomeric DNA ,,, and their presence is proposed to protect ends from nuclease degradation, and/or affect telomerase recruitment -. The fact that telomeres bear constitutive single-stranded G-tails makes them strong candidates for G4 formation, unless proteins or other structures (for example, t-loops) prevent their formation. The fourth conserved association of G4 motifs in all three organisms is their depletion from ORFs (Table 2) ,.
Taken together, these common patterns in the genomic distribution of G4 motifs add to the increasing evidence that G4 structures have regulatory functions that are maintained by selection and that likely counterbalance their negative effects on genome stability ,,,. However, in certain contexts, such as ORFs, the price of their negative effects may be too steep, leading to selection against G4 motifs in these regions.
Pfh1 suppresses G4-induced genomic instability in S. pombe
We provide multiple lines of evidence that support the importance of Pfh1 in suppressing G4-induced genomic instability in S. pombe. First, 20% of G4 motifs had high Pfh1 occupancy, consistent with the possibility that Pfh1 acts at a subset of G4 motifs in vivo. This number is likely an underestimate. For example, only 9% of the Class II Pfh1-sensitive sites had significant Pfh1 binding; yet their dependence on Pfh1 argues that all of these sites may bind Pfh1. We attribute the lack of detectable binding at most Class II sites to the speed with which Pif1 family helicases unwind G4 structures , the stringent criteria used to identify binding sites, and technical difficulties detecting DNA helicases by ChIP . In contrast to Class II sites, Pfh1 was detected at 40% of Class III sites. This finding suggests that G4 structures persist longer at Class III than at Class II sites. For example, at Class III G4 motifs, G4 structures may reform after Pfh1 action, leading to repeated cycles of Pfh1 binding and unwinding at the site. This hypothesis could also explain why Class III G4 motifs were sites of replication fork slowing and/or DNA damage even in WT cells. Our data also likely underestimate the fraction of G4 motifs that form G4 structures: other helicases might act at other G4 motifs, while some G4 structures may form only outside of S phase or only in specific growth conditions. Of course, some G4 motifs may rarely or never form G4 structures.
The second result supporting a role for Pfh1 at G4 motifs is that replication pausing at G4 motifs, as monitored by high Cdc20 occupancy, was much higher in Pfh1-depleted than in WT cells (Figure 2). Moreover, G4 motifs that were associated with Pfh1 in WT cells were much more likely to be associated with replication fork pausing than those that were not bound by Pfh1 (Table 3). Third, G4 motifs were significantly associated with DNA damage, as indicated by the presence of γ-H2A. Some G4 motifs were sites of damage even in the presence of Pfh1, but this association and γ-H2A levels were considerably stronger in its absence (Figure 3; Additional file 4). Moreover, when Pfh1 was depleted, G4 motifs with replication fork pausing (high Cdc20) were much more likely to result in DNA damage than G4 motifs without pausing (62% versus 35%; Table 1).
As in budding yeast and humans, G4 motifs in S. pombe were underrepresented within ORFs (Table 2). However, when S. pombe G4 motifs did occur in ORFs, they were enriched on the transcribed strand (72%; P <1E-12). This enrichment was particularly marked for the Pfh1-associated G4 motifs (87%). Furthermore, this enrichment was especially high among the top 500 most highly transcribed RNA polymerase II genes (94%; P <0.001; Additional file 6). Although a recent report demonstrated that highly expressed genes may be biased for high ChIP-seq signals , fork slowing in highly expressed S. pombe genes and its increase in Pfh1-depleted cells, are detected by two-dimensional gels as well as by ChIP . The G4 motifs in ORFs tended to fall within the first half of the ORF (62%), but G4 motifs were observed near the ends of genes as well [see Additional file 7]. A similar strand bias was seen when only considering G4 motifs that overlapped 5′ and 3′ UTRs, 74 % (35/47) and 77% (37/48), respectively, were found on the transcribed strand. A G4 motif on the transcribed strand is expected to inhibit RNA polymerase progression. One possibility is that this class of G4 motifs regulates transcription elongation: transcription could start but then pause at the G4 motif until the regulated recruitment of Pfh1 allows G4 unwinding and resumption of transcription. This type of regulation is particularly appealing in multi-cellular organisms where developmentally regulated genes are often controlled by activation of a paused RNA polymerase . Alternatively, Pfh1 bound to G4 motifs within highly transcribed genes might facilitate RNA removal and inhibit R-loop formation. This speculation is based on the unusual property of budding yeast Pif1, which has higher unwinding activity on RNA/DNA compared to DNA/DNA hybrids .
The evolutionary conservation of G4 motif enrichment in promoters, rDNA and telomeres in two distantly related yeasts and humans argues that they have functions in each of these regions. Their depletion from the protein-coding ORFs in the three organisms, as well as their association with fork stalling and DNA breakage in the two yeasts, demonstrates that their positive roles come with a negative impact on genome integrity. These negative effects are mitigated in both S. cerevisiae and S. pombe by the action of Pif1 family helicases. We propose that processing of G4 structures by Pif1 family helicases is a common and evolutionarily conserved mechanism.
ChIP and preparation of libraries for DNA sequencing
All yeast strains used in this study are listed in Additional file 8. ChIP experiments were performed as described previously . Briefly, cells were cross-linked in 1% formaldehyde at 25°C for five minutes. The chromatin was sheared to an average of approximately 300 bps length with a Covaris E220 system and immunoprecipitated with anti-Myc antibody (Clontech Laboraties, Mountain View, California, USA Cat. nr 631206), γ-H2A antibody (a kind gift from C. Redon) or anti-HA antibody (Santa Cruz biotechnologies, Dallas, Texas, USA Cat nr. sc.7392x). Both input and immunoprecipitated DNA were purified and quantified by real-time PCR using primers for ade6 and STE. A sample of 20 ng DNA was used to prepare a sequencing library with the TruSeq DNA sample preparation kit v2 (Illumina, San Diego, California, USA Cat nr. FC-121-2001) as described by the manufacturer. The libraries were sequenced with the Illumina HiSeq 2000 sequencing platform at Princeton University (Princeton, NJ, USA). Two biological replicates were sequenced for each experiment.
ChIP combined with quantitative PCR
ChIP experiments were performed as described above. Both input and immunoprecipitated DNA were purified and quantified by quantitative PCR (Roche Diagnostics, Indianapolis, Indiana, USA Lightcycler® 96 instrument) with either primer pairs for tRNA glu.05, GC_1, GC_2, GC_3, G4_1, G4_2, or G4_3 [see Additional file 9 for primer sequences]. For γ-H2A ChIP-qPCR, the primer pairs were designed within a 5 kb region of the tRNA gene, the GC-rich region, and the G4-motifs [see Additional file 9]. For strains expressing Pfh1-13Myc, an untagged strain was used as a control (YSP3; Additional file 8). The IP/input ratio was calculated by dividing immunoprecipitated DNA by input DNA. The data presented are an average of three independent cultures.
Sequence analysis and peak calling
The 101 bp single-end sequencing reads were analyzed using tools available in the Galaxy platform . Reads were mapped to the S. pombe genome (Pombase_09052011)  using Bowtie , permitting two mismatches in the seed (seed length 28 bp). Reads were mapped to the genome for both input and ChIP samples. Only uniquely mapped reads were included in the alignments. Peaks were identified with MACS1.4  using the following settings: bandwidth of 300 (DNA shearing size), P-value cutoff of 10-5 and tag size of 101. Input DNA was used as the control. The average peak size was 1,838 ± 983 bp (± standard deviation) for the 621 Pfh1-associated regions, 1,526 ± 717 bp for the 485 high Cdc20-associated regions in Pfh1-expressing cells and 1,666 ± 860 bp for the 517 high Cdc20-associated regions in Pfh1-depleted cells. The 179 high γ-H2A peaks in Pfh1-expressing cells were 1,486 ± 742 bp, and the 582 γ-H2A peaks in Pfh1-depleted cells were 1,974 ± 1071 bp. All sequence data are deposited with the GEO accession number GSE59178.
G4 motif identification
The genome sequences and alignments for S. pombe, S. octosporus, S. cryophilus and S. japonicus were downloaded from the supporting web site for Rhind et al.  on May 3, 2012. We identified the location of all DNA sequences in the nuclear and mt genomes with the potential to form G4 structures (‘G4 motifs’) using a previously described regular expression search program with a minimum G-island length of 3 and a maximum loop length of 25 bp . Because the current S. pombe genome assembly does not include telomeres, we also scanned the sequence of DNA in a recently sequenced telomere library (Webb and VAZ, submitted). Using alignments of the nuclear genomes of the four fission yeasts, we analyzed G4 motif evolutionary conservation. In this analysis, any G4 motif for which the aligned portion of another genome contained a G4 motif was considered conserved between the two species.
Analysis of G4 motif associations with genomic features
To explore potential functions for G4 motifs, we compared their genomic distribution in S. pombe with known functional genomic regions. We considered the following genome annotations: ORFs, dubious genes, essential genes, highly expressed genes , origins of replication from OriDB , telomeres (Webb and VAZ, submitted), LTRs, 5S rRNA, tRNAs, centromeres, NDRs [see Additional file 10]  called using Podbat , rDNA repeats and meiotic DSB hot spots . Unless otherwise noted, these annotations were taken from PomBase . The highly transcribed genes were determined from the RNA-seq data collected by . We took the 500 genes with the highest expression in fragments per kilobase per million reads (FPKM) from their combined dataset, which summarized expression patterns from cells in log phase, glucose depletion, early stationary phase, and heat shock [see Additional file 6]. Results were similar when we considered only the 100 most highly expressed genes.
Regions that overlapped a G4 motif were considered to be associated with the G4. We also computed the association of G4 motifs with the Pfh1, Cdc20 and γ-H2A peaks determined here by ChIP-seq. For the γ-H2A analysis, a 5 kb window on both sides was used to determine G4 motif associations, as DNA damage results in maximal phosphorylation in a roughly 5 kb region on either side of the break . For Pfh1 and Cdc20 peaks, a window of 300 bp, the DNA shearing size, was used, but our main results held with windows of 0 and 500 bp.
To evaluate the significance of observed associations between G4 motifs and genomic annotations, we performed simulations to obtain an empirical P-value. The number of associations expected at random if there were no functional relationship between two sets of genomic regions (for example, G4 motifs and essential genes) depends on many factors, including their lengths, distribution across the chromosomes, and nucleotide contents. To account for these factors, we generated 1,000 random ‘control’ sets of genomic regions. Each of these sets consisted of 446 individual genomic regions—one for each of the observed G4 motifs. Each of these individual regions was randomly placed on the genome with three constraints. First, it had to be on the same chromosome as its corresponding G4 motif; second, it had to have the same length as the G4 motif; and third, it had to have the same GC content as the G4 motif. As a result, each of the 1,000 random control region sets had the same length, chromosome and GC content distribution as the actual G4 motifs. To obtain a P-value for an observed association between the G4 motifs and a genomic annotation of interest, we compared it to the number of associations for each of the 1,000 control region sets with the same annotation. The P-value was the number of the random control sets in which a more extreme association was observed. When indicated, we also considered control region sets that were not constrained to match the GC content of the G4 motifs. To account for the testing of multiple hypotheses in this analysis, we computed q-values from the P-values using the qvalue program . The q-value for a test is the FDR that results if that test is called significant. All associations with q less than 0.05 were considered significant.
Bernstein KA, Gangloff S, Rothstein R: The RecQ DNA helicases in DNA repair. Annu Rev Genet. 2010, 44: 393-417. 10.1146/annurev-genet-102209-163602.
Chisholm KM, Aubert SD, Freese KP, Zakian VA, King MC, Welcsh PL: A genomewide screen for suppressors of Alu-mediated rearrangements reveals a role for PIF1. PLoS One. 2012, 7: e30748-10.1371/journal.pone.0030748.
Bochman ML, Judge CP, Zakian VA: The Pif1 family in prokaryotes: what are our helicases doing in your bacteria?. Mol Biol Cell. 2011, 22: 1955-1959. 10.1091/mbc.E11-01-0045.
Bochman ML, Sabouri N, Zakian VA: Unwinding the functions of the Pif1 family helicases. DNA Repair (Amst). 2010, 9: 237-249. 10.1016/j.dnarep.2010.01.008.
Pinter SF, Aubert SD, Zakian VA: The Schizosaccharomyces pombe Pfh1p DNA helicase is essential for the maintenance of nuclear and mitochondrial DNA. Mol Cell Biol. 2008, 28: 6594-6608. 10.1128/MCB.00191-08.
Sabouri N, McDonald KR, Webb CJ, Cristea IM, Zakian VA: DNA replication through hard-to-replicate sites, including both highly transcribed RNA Pol II and Pol III genes, requires the S. pombe Pfh1 helicase. Genes Dev. 2012, 26: 581-593. 10.1101/gad.184697.111.
Steinacher R, Osman F, Dalgaard JZ, Lorenz A, Whitby MC: The DNA helicase Pfh1 promotes fork merging at replication termination sites to ensure genome stability. Genes Dev. 2012, 26: 594-602. 10.1101/gad.184663.111.
Bochman ML, Paeschke K, Zakian VA: DNA secondary structures: stability and function of G-quadruplex structures. Nat Rev Genet. 2012, 13: 770-780. 10.1038/nrg3296.
Paeschke K, Simonsson T, Postberg J, Rhodes D, Lipps HJ: Telomere end-binding proteins control the formation of G-quadruplex DNA structures in vivo. Nat Struct Mol Biol. 2005, 12: 847-854. 10.1038/nsmb982.
Paeschke K, Juranek S, Simonsson T, Hempel A, Rhodes D, Lipps HJ: Telomerase recruitment by the telomere end binding protein-beta facilitates G-quadruplex DNA unfolding in ciliates. Nat Struct Mol Biol. 2008, 15: 598-604. 10.1038/nsmb.1422.
Huppert JL, Balasubramanian S: Prevalence of quadruplexes in the human genome. Nucleic Acids Res. 2005, 33: 2908-2916. 10.1093/nar/gki609.
Todd AK, Johnston M, Neidle S: Highly prevalent putative quadruplex sequence motifs in human DNA. Nucleic Acids Res. 2005, 33: 2901-2907. 10.1093/nar/gki553.
Biffi G, Tannahill D, McCafferty J, Balasubramanian S: Quantitative visualization of DNA G-quadruplex structures in human cells. Nat Chem. 2013, 5: 182-186. 10.1038/nchem.1548.
Henderson A, Wu Y, Huang YC, Chavez EA, Platt J, Johnson FB, Brosh RM, Sen D, Lansdorp PM: Detection of G-quadruplex DNA in mammalian cells. Nucleic Acids Res. 2013, 35: 406-413.
Capra JA, Paeschke K, Singh M, Zakian VA: G-quadruplex DNA sequences are evolutionarily conserved and associated with distinct genomic features in Saccharomyces cerevisiae. PLoS Comput Biol. 2010, 6: e1000861-10.1371/journal.pcbi.1000861.
Rawal P, Kummarasetti VB, Ravindran J, Kumar N, Halder K, Sharma R, Mukerji M, Das SK, Chowdhury S: Genome-wide prediction of G4 DNA as regulatory motifs: role in Escherichia coli global regulation. Genome Res. 2006, 16: 644-655. 10.1101/gr.4508806.
Huppert JL, Balasubramanian S: G-quadruplexes in promoters throughout the human genome. Nucleic Acids Res. 2007, 35: 406-413. 10.1093/nar/gkl1057.
Paeschke K, Capra JA, Zakian VA: DNA replication through G-quadruplex motifs is promoted by the Saccharomyces cerevisiae Pif1 DNA helicase. Cell. 2011, 145: 678-691. 10.1016/j.cell.2011.04.015.
Paeschke K, Bochman ML, Garcia PD, Cejka P, Friedman KL, Kowalczykowski SC, Zakian VA: Pif1 family helicases suppress genome instability at G-quadruplex motifs. Nature. 2013, 497: 458-462. 10.1038/nature12149.
Piazza A, Serero A, Boule JB, Legoix-Ne P, Lopes J, Nicolas A: Stimulation of gross chromosomal rearrangements by the human CEB1 and CEB25 minisatellites in Saccharomyces cerevisiae depends on G-quadruplexes or Cdc13. PLoS Genet. 2012, 8: e1003033-10.1371/journal.pgen.1003033.
Rodriguez R, Miller KM, Forment JV, Bradshaw CR, Nikan M, Britton S, Oelschlaegel T, Xhemalce B, Balasubramanian S, Jackson SP: Small-molecule-induced DNA damage identifies alternative DNA structures in human genes. Nat Chem Biol. 2012, 8: 301-310. 10.1038/nchembio.780.
Rhind N, Chen Z, Yassour M, Thompson DA, Haas BJ, Habib N, Wapinski I, Roy S, Lin MF, Heiman DI, Young SK, Furuya K, Guo Y, Pidoux A, Chen HM, Robbertse B, Goldberg JM, Aoki K, Bayne EH, Berlin AM, Desjardins CA, Dobbs E, Dukaj L, Fan L, FitzGerald MG, French C, Gujja S, Hansen K, Keifenheim D, Levin JZ, et al: Comparative functional genomics of the fission yeasts. Science 2011, 332:930–936. http://www.broadinstitute.org/annotation/genome/schizosaccharomyces_group/MultiDownloads.html.,
Sugawara NF: DNA sequences at the telomeres of the fission yeast S. pombe. PhD. thesis. Harvard University; 1989.,
Sanchez JA, Kim SM, Huberman JA: Ribosomal DNA replication in the fission yeast, Schizosaccharomyces pombe. Exp Cell Res. 1998, 238: 220-230. 10.1006/excr.1997.3835.
Naehring J, Kiefer S, Wolf K: Nucleotide sequence of the Schizosaccharomyces japonicus var. versatilis ribosomal RNA gene cluster and its phylogenetic implications. Curr Genet. 1995, 28: 353-359. 10.1007/BF00326433.
Hershman SG, Chen Q, Lee JY, Kozak ML, Yue P, Wang LS, Johnson FB: Genomic distribution and functional analyses of potential G-quadruplex-forming sequences in Saccharomyces cerevisiae. Nucleic Acids Res. 2008, 36: 144-156. 10.1093/nar/gkm986.
Hanakahi LA, Sun H, Maizels N: High affinity interactions of nucleolin with G-G-paired rDNA. J Biol Chem. 1999, 274: 15908-15912. 10.1074/jbc.274.22.15908.
Ben-Shem A, GarreaudeLoubresse N, Melnikov S, Jenner L, Yusupova G, Yusupov M: The structure of the eukaryotic ribosome at 3.0 A resolution. Science 2011, 334:1524–1529.,
Doniger SW, Fay JC: Frequent gain and loss of functional transcription factor binding sites. PLoS Comput Biol. 2007, 3: e99-10.1371/journal.pcbi.0030099.
Tuch BB, Li H, Johnson AD: Evolution of eukaryotic transcription circuits. Science. 2008, 319: 1797-1799. 10.1126/science.1152398.
Tuch BB, Galgoczy DJ, Hernday AD, Li H, Johnson AD: The evolution of combinatorial gene regulation in fungi. PLoS Biol. 2008, 6: e38-10.1371/journal.pbio.0060038.
Storey JD: A direct approach to false discovery rates. J Royal Stat Soc Ser B-Stat Method. 2002, 64: 479-498. 10.1111/1467-9868.00346.
McDonald KR, Sabouri N, Webb CJ, Zakian VA: The Pif1 family helicase Pfh1 facilitates telomere replication and has an RPA-dependent role during telomere lengthening. DNA Repair (Amst). 2014, 24: 80-86. 10.1016/j.dnarep.2014.09.008.
Azvolinsky A, Giresi PG, Lieb JD, Zakian VA: Highly transcribed RNA polymerase II genes are impediments to replication fork progression in Saccharomyces cerevisiae. Mol Cell. 2009, 34: 722-734. 10.1016/j.molcel.2009.05.022.
Rozenzhak S, Mejia-Ramirez E, Williams JS, Schaffer L, Hammond JA, Head SR, Russell P: Rad3 decorates critical chromosomal domains with gammaH2A to protect genome integrity during S-Phase in fission yeast. PLoS Genet. 2010, 6: e1001032-10.1371/journal.pgen.1001032.
Forsburg SL: The best yeast?. Trends Genet. 1999, 15: 340-344. 10.1016/S0168-9525(99)01798-9.
Eddy J, Maizels N: Gene function correlates with potential for G4 DNA formation in the human genome. Nucleic Acids Res. 2006, 34: 3887-3896. 10.1093/nar/gkl529.
Lantermann AB, Straub T, Stralfors A, Yuan GC, Ekwall K, Korber P: Schizosaccharomyces pombe genome-wide nucleosome mapping reveals positioning mechanisms distinct from those of Saccharomyces cerevisiae. Nat Struct Mol Biol. 2010, 17: 251-257. 10.1038/nsmb.1741.
Drygin D, Siddiqui-Jain A, O’Brien S, Schwaebe M, Lin A, Bliesath J, Ho CB, Proffitt C, Trent K, Whitten JP, Lim JK, Von Hoff D, Anderes K, Rice WG: Anticancer activity of CX-3543: a direct inhibitor of rRNA biogenesis. Cancer Res. 2009, 69: 7653-7661. 10.1158/0008-5472.CAN-09-1304.
Oganesian L, Moon IK, Bryan TM, Jarstfer MB: Extension of G-quadruplex DNA by ciliate telomerase. EMBO J. 2006, 25: 1148-1159. 10.1038/sj.emboj.7601006.
Zahler AM, Williamson JR, Cech TR, Prescott DM: Inhibition of telomerase by G-quartet DNA structures. Nature. 1991, 350: 718-720. 10.1038/350718a0.
Zaug AJ, Podell ER, Cech TR: Human POT1 disrupts telomeric G-quadruplexes allowing telomerase extension in vitro. Proc Natl Acad Sci U S A. 2005, 102: 10864-10869. 10.1073/pnas.0504744102.
Teytelman L, Thurtle DM, Rine J, van Oudenaarden A: Highly expressed loci are vulnerable to misleading ChIP localization of multiple unrelated proteins. Proc Natl Acad Sci U S A. 2013, 110: 18602-18607. 10.1073/pnas.1316064110.
Adelman K, Lis JT: Promoter-proximal pausing of RNA polymerase II: emerging roles in metazoans. Nat Rev Genet. 2012, 13: 720-731. 10.1038/nrg3293.
Boule JB, Zakian VA: The yeast Pif1p DNA helicase preferentially unwinds RNA DNA substrates. Nucleic Acids Res. 2007, 35: 5809-5818. 10.1093/nar/gkm613.
Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, Mangan M, Nekrutenko A, Taylor J: Galaxy: a web-based genome analysis tool for experimentalists. Curr Protoc Mol Biol 2010, Chapter 19:Unit 19 10 11–21.,
Wood V, Harris MA, McDowall MD, Rutherford K, Vaughan BW, Staines DM, Aslett M, Lock A, Bahler J, Kersey PJ, Oliver SG: PomBase: a comprehensive online resource for fission yeast. Nucleic Acids Res. 2012, 40: D695-D699. 10.1093/nar/gkr853.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.
Siow CC, Nieduszynska SR, Muller CA, Nieduszynski CA: OriDB, the DNA replication origin database updated and extended. Nucleic Acids Res. 2012, 40: D682-D686. 10.1093/nar/gkr1091.
Sadeghi L, Bonilla C, Stralfors A, Ekwall K, Svensson JP: Podbat: a novel genomic tool reveals Swr1-independent H2A.Z incorporation at gene coding sequences through epigenetic meta-analysis. PLoS Comput Biol. 2011, 7: e1002163-10.1371/journal.pcbi.1002163.
Fowler KR, Gutierrez-Velasco S, Martin-Castellanos C, Smith GR: Protein determinants of meiotic DNA break hot spots. Mol Cell. 2013, 49: 983-996. 10.1016/j.molcel.2013.01.008.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.
Work in the Zakian lab is supported by NIH grant GM26938. JAC was supported by development funds from Vanderbilt University, NS by the Wenner-Gren Foundations, the Kempe Foundations (SMK-1246 and SMK-1325), and the Swedish Society for Medical Research. We are thankful to C. Redon for sending us anti- γ-H2A antibodies, B. Diner for constructing the epitope tagged Cdc20-3HA strain, K. Petzold for expert comments on the ribosome structure, and P. Svensson for help with the genomic location list of NDRs. We thank K. Paeschke for critical comments on the manuscript and Princeton sequencing core facility for their expertise in sequencing. We thank P.D. Garcia for help with some experiments.
The authors declare that they have no competing interests.
All authors designed the experiments; NS created yeast strains, prepared samples for ChIP-seq experiments, and processed and mapped the reads; JAC performed genome wide and statistical analysis; NS, JAC, and VAZ discussed and contextualized the results and wrote the manuscript; VAZ supervised the study. All authors read and approved the final manuscript.