Wild-type isogenic reference strains WT 3 and WT 5 were used to determine relative gene-expression changes in trt1^- samples. Before averaging the expression values from the two reference strains, the similarity of their expression profiles was assessed. The dye ratios measured by microarray for each strain were plotted against each other (Figure 1a). All genes had expression values that varied less than twofold between the two samples, indicating that the samples were highly similar. The wild-type values used in this paper are thus the average expression values of strains WT 3 and WT 5.

To learn whether changes in gene expression would result from subjecting cells to the continuous growth program for 15 days, gene-expression values from strain WT 5 on day 1 of the growth curve were compared with those of the same strain harvested on day 15 (Figure 1b). Only three genes (SPBC354.08c, atp8⁺ and cox1⁺) changed their expression values by more than twofold, and they were only slightly greater; thus, the vast majority of genes do not have altered expression as a result of long-term growth in culture, provided that expression is measured while the cells are in early log phase (see Materials and methods). These three genes also had expression changes of more than twofold in one or more conditions measured for trt1^- cells, but given their variable expression in wild-type cells, these changes were most probably unrelated to the absence of telomerase.

Diploid S. pombe cells that were heterozygous for trt1⁺ and able to maintain full-length telomeres were sporulated, and the resulting trt1⁺ and trt1^- cells propagated through a 15-day growth curve (Figure 2a). Cells lacking telomerase gave rise to survivors after day 8 concomitant with heterogeneous amplified telomeric repeats and telomere-associated sequence (TAS) (Figures 2b-d), indicative of linear chromosomes 14. By day 15, the culture was dominated by faster-growing cells with linear chromosomes. The linear structure of these chromosomes was confirmed by their ability to enter a pulsed-field gel (Figure 3b, lane g), and the existence of terminal chromosome fragments C, I, L and M after digestion of chromosomes with NotI (Figure 3a-d, lane e) 1415. Cells passing through crisis (days 7 and 9) also had weak hybridization signals for the C+M and I+L fragments (Figure 3d, lanes c-d), suggesting a mix of cells with either linear or circular chromosomes, or perhaps cells containing both linear and circular chromosomes. The inability to detect intact chromosomal DNA at day 7 (Figure 3b, lane e) may have resulted from the presence of cells with circularized chromosomes (Figure 3d, lane c) that do not enter pulsed-field gels.

Strains C1 and C5 had circular chromosomes as evidenced by lack of telomeric repeats (data not shown), lack of TAS2 sequence (data not shown), the inability of chromosomes to enter a pulsed-field gel (Figure 3b, lanes b-c), the lack of terminal chromosome fragments C, I, L and M (Figures 3c,d, lanes g-h) 1415, and hybridization signals to fragments C+M and I+L (Figure 3d, lanes g-h).

Two waves of altered gene expression were seen during the growth curve (Figure 4a), the first with a peak at day 7, consisting of around 110 genes with expression upregulated twofold or more, and the second with a peak at day 9, consisting of three microarray signals that appear to represent a single ORF (see below) (Figure 4a). The peak of the first wave (day 7) was nearly coincident with crisis in the cell population (day 8) (Figure 2a) and the time when telomeres were shortest (near day 7) (Figure 2c,d). The second peak of gene expression at day 9 was coincident with the emergence of survivors (Figure 2a-d).

The vast majority of expression changes involved upregulation, and only seven genes had downregulated expression of twofold or greater on two or more days of the growth curve. Notably, there were three cases of reduction in expression greater than tenfold: trt1⁺ (intentionally knocked out), SPAC2E1P3.04 (a predicted copper amine oxidase) and SPAC2E1P3.05c (unknown function). Hybridizations of genomic DNA to microarrays (data not shown) revealed that genes SPAC2E1P3.04 and SPAC2E1P3.05c were deleted from the genome in all strains except WT 3, WT 5 and C1. Interestingly, these two genes are within about 4 kb of transposable element SPAC167.08 (Tf2-2), suggesting a hotspot for DNA excision. In no case was gene amplification detected by genomic hybridization (data not shown), so the observed increases in expression were most probably due to transcriptional or post-transcriptional regulation, as opposed to changes in gene copy number.

Because a relatively large number of trt1^- strains were studied, the identification of genes with consistently altered expression was facilitated by selecting those genes with expression changes of twofold or more in two or more days of the growth curve or, alternatively, in both strains C1 and C5. This criterion was met by 123 genes, of which 54 (44%) overlapped between the growth curve and survivors with circularized chromosomes. In addition, of the 67 genes that had their expression changed twofold or more exclusively in the growth curve, many displayed altered expression just below the cutoff in survivors with circularized chromosomes. Two genes - SPBC1683.06c (a predicted uridine ribohydrolase) and SPBC1198.01 (a predicted formaldehyde dehydrogenase) - had expression changes of twofold or more in both strains C1 and C5, but no significant changes during the growth curve. As a measure of confidence, 84 of the 123 genes (approximately 68%) met a more stringent criterion requiring a gene to change its expression in three or more of the 17 conditions. Additional confidence that expression changes scored as significant were not false positives came from the remarkably continuous manner in which gene expression changed throughout the growth curve (Figure 4a).

The 123 genes with altered expression encompass a broad range of functions, but were especially enriched in genes associated with energy production and carbohydrate metabolism (Table 1). There were seven pseudogenes and 29 predicted genes that did not have assigned functions at the time of writing. For nearly all the gene-type categories, there was a larger number of genes with altered expression in the growth curve than in the survivors with circular chromosomes (Table 1). This difference may be attributable to the fact that cells in the growth curve were experiencing crisis whereas strains C1 and C5 were survivors, presumably with established mechanisms to cope with the absence of or the loss of telomeres.

The telomerase-deletion response had a large overlap with genes that changed expression in response to environmental stresses. Fission yeast stress-response genes can be separated into a CESR, in which genes changed expression in all or most of the stresses studied (oxidative stress, heavy metals, heat shock, osmotic stress and DNA damage), and into more specific stress responses 16. Of the 123 genes with altered expression in trt1^- cells, 48 (about 39%) also had upregulated expression among a conservative list of CESR genes (P ~ 10^-77) 16, and two genes had downregulated expression in the CESR and in this study. Of the 110 genes with expression upregulated twofold or more on day 7 of the growth curve, 44% overlapped with the CESR. Comparison with a less conservative list of CESR genes 16 suggested that 54% of the 123 genes with altered expression in trt1^- cells had overlap with the CESR (P ~ 10^-81). With respect to specific stress responses 16, there were 17/123 genes in common with the oxidative stress response (P ~ 10^-32), and 11/123 genes in common with the heat stress response (P ~ 10^-24). The stress response study found that the DNA damage response and the oxidative stress response have substantial overlap 16. Therefore, the genes with altered expression in this study that overlap with the oxidative-stress response may represent a DNA damage response to short telomeres.

Comparisons of all the gene-expression profiles in this study revealed striking differences between the profiles of survivors with linear chromosomes versus those with circular chromosomes. Survivors with linear chromosomes (days 12-15 of the growth curve) had gene-expression patterns similar to those of cells with native telomeres in the first two days of the growth curve. To illustrate, by day 12 of the growth curve, the gene-expression profiles of survivors became relatively constant and remained so through day 15. The profiles of days 12-15 appear most similar to days 1 and 2 of the growth curve, immediately after cells lost telomerase and were experiencing shortening telomeres (Figure 4b). This observation was confirmed by hierarchical clustering (Figure 4c). Conversely, survivors with circular chromosomes had gene-expression profiles that most resembled those of cells in crisis during days 5-8 of the growth curve (Figure 4b,c).

There were 54 genes with clearly altered expression (twofold or more) mainly during crisis in the growth curve that also had altered expression in the survivors with circular chromosomes (Table 2, Figure 5). The expression of all but three of these 54 genes was not altered in survivors with linear chromosomes (growth curve days 12-15) (Table 2). Of the 54 genes, 30 (56%) overlapped with the conservative list of CESR genes (P ~ 10^-46), and eight genes (15%) overlapped with the oxidative stress response (P ~ 10^-14). There were 8/54 genes (15%) that overlapped with the heat stress response (P ~ 10^-17). Because of the extensive overlap of the 54 genes with the CESR, we conclude that survivors with circular chromosomes had a sustained stress response.

Of the 54 genes, 51 represent a gene-expression signature that differentiates survivors with circular chromosomes from those with linear chromosomes. As an independent test of whether these 51 genes can serve as a signature for cells with circularized chromosomes, two additional cultures (strains H1 and H2, see Materials and methods) with circularized chromosomes were grown and analyzed by microarray. Both strains clearly displayed altered expression of the 51 genes whereas survivors with linear chromosomes did not (Figure 5), thus validating this gene signature.

One feature of microarray studies is that genes not previously recognized to be under the control of a common regulator can often be associated by similar expression patterns 17. On the basis of this hypothesis, a list of genes known to be involved in telomere maintenance and recombination was inspected. However, the expression patterns of all these genes were not substantially changed throughout the course of the study (data not shown). Genes investigated included pku70⁺ and lig4⁺, which encode components of the non-homologous end-joining pathway 18; taz1⁺ 19 and pot1⁺ 20 encoding telomere DNA-binding proteins; telomerase component est1⁺ 21; homologous recombination-related genes rad22⁺ 22, rhp54⁺ 23, rad32⁺ 24 and rhp51⁺ 25; RecQ helicase gene rqh1⁺ 26; silencing component clr4⁺ 27; and telomere maintenance components pof3⁺ 28 and rad3⁺ 15. Interestingly, even though pof3⁺ and clr4⁺ expression did not change, the genes with altered expression in this study had a statistically significant overlap with the lists of genes with induced expression in pof3 mutants (P < 10^-45) 28 and clr4 mutants (P < 10^-45) 29; a significant correlation was also observed with genes that changed expression in the RNA interference (RNAi)-machinery mutants dcr1⁺, ago1⁺ and rdp1⁺ (P ~ 10^-22) 29. These genes with altered expression may act in common pathways downstream of trt1⁺, clr4⁺, pof3⁺ and the RNAi machinery.

The second wave of gene-expression changes during the growth curve (Figure 4a) consisted of three microarray signals: SPAC212.11 (largest magnitude), SPAC212.06 (second largest magnitude) and the reverse transcript of centromeric dh repeats 30. Inspection of the sequences revealed that the microarray signals from SPAC212.06 and centromeric dh repeats most probably resulted from cross-hybridization with the SPAC212.11 transcript (see Materials and methods).

A BLAST search of the SPAC212.11 predicted protein sequence found that the ORF has the most similarity to RecQ DNA helicases of superfamily II (Figure 6) (reviewed in 31). We report a role for the helicase in cells passing through crisis in a separate study (J.G.M., K.J. Goodrich, J.B. and T.R.C., unpublished work) and investigate its transcriptional regulation here.

SPAC212.11 is the last sequenced ORF on the left arm of chromosome I. The sub-telomeric regions of chromosomes I and II have significant similarity 32. A BLAST search performed with the SPAC212.11 DNA sequence (5.6 kb) revealed a paralog, SPBCPT2R1.08c (6.3 kb), located on the right arm of chromosome II (the microarray had no probe for SPBCPT2R1.08c), and partial homology on the right arm of chromosome I. The annotated sequence of SPBCPT2R1.08c includes the entirety of the SPAC212.11 sequence with only a single base change. The SPAC212.11 sequence does not contain a stop codon because the ORF is located at the end of the sequencing contig, which ended before a stop codon was reached. Comparison with the annotated SPBCPT2R1.08c sequence suggests that SPAC212.11 has an additional 95 bp before the stop codon.

Both SPBCPT2R1.08c and SPAC212.11 are the last predicted genes on their respective sub-telomeric sequencing contigs. Analysis of contig pT2R1 revealed that the 3' end of SPBCPT2R1.08c is approximately 2.8 kb upstream from the start of TAS3 (Figure 2b). Since TAS3 is around 7 kb from the chromosome end, the 3' end of SPBCPT2R1.08c is approximately 10 kb from the telomeric repeats.

It is not known which of the paralogs contributed to the SPAC212.11 microarray signal. For the sake of simplicity, further references in the text to 'the putative helicase' are meant to include SPAC212.11, SPBCPT2R1.08c and any paralogs, collectively.

The nucleotide BLAST search performed with the SPAC212.11 sequence also revealed that the ORF contains regions of homology to dh repeats (Figure 6), which are targeted for heterochromatin formation via an RNAi-mediated mechanism in S. pombe 3334. These repeats are typically located at centromeres and the K region of the mating-type locus 3033353637.

Centromeric repeats, previously thought to be transcriptionally silent, are transcribed in both the forward and reverse directions, leading to formation of double-stranded RNA (dsRNA). However, these transcripts do not accumulate in wild-type cells. Reverse-strand centromeric transcripts are synthesized and rapidly processed by the RNAi machinery, while forward-strand synthesis is silenced transcriptionally. RNA-dependent RNA polymerase (Rdp1) associates with centromeric repeat DNA and may use siRNAs corresponding to centromeric transcripts 38 to prime forward transcription from reverse-strand templates, thus resulting in dsRNA formation and maintenance of the heterochromatic state. In the RNAi mutants dcr1^-, ago1^- and rdp1^-, centromeric silencing is abolished and accumulation of both forward and reverse centromeric transcripts is observed 33.

Microarray, northern blot and reverse transcription (RT)-PCR analysis indicated that the putative helicase gene was robustly expressed in cells emerging from crisis, but was weakly (or not at all) expressed in wild-type cells, strains C1 and C5 and survivors with linear chromosomes (Figures 4a,b, 7a, and data not shown). As the putative helicase transcript was not detectable by northern blot in wild-type cells (data not shown), we hypothesized that this ORF could be silenced by its dh repeats, but that this silencing may have been disrupted in trt1^- cells as a result of genomic instability. Arguing against this hypothesis, however, Southern analysis with probe P_5' (Figure 6), which is specific for the helicase, did not reveal any DNA rearrangements during crisis close to the helicase that might have contributed to loss of silencing (data not shown). Nevertheless, the loss of silencing observed might lead to expression of both strands of the putative helicase, as was found for centromeric dh repeats in RNAi mutants.

To test for the presence of both strands, strand-specific RT-PCR was used with primers spanning the dh repeats of the putative helicase (region P_dh in Figure 6). The forward strand was expressed at levels higher than in wild type in cells from days 7, 9 and 15 of the growth curve. These results were consistent with microarray analysis that detected the 3' end of the forward transcript (Figure 7a). The reverse strand was weakly detectable in cells from days 7 and 9 of the growth curve (Figure 7a).

dsRNA arising from the repeats presumably could have formed on days 7 and 9 of the growth curve, but why such RNA was not all processed by the RNAi machinery is not clear. On days 7 and 9 of the growth curve, the RNAi machinery was not apparently affected by the mutation of telomerase as centromeric dh repeat transcripts were not detected by RT-PCR (Figure 7a).

We next hypothesized that if the RNAi machinery were involved in transcriptional silencing of the putative helicase in wild-type cells, transcript should accumulate in mutant RNAi strains. Strikingly, both ago1^- and dcr1^- strains displayed significant accumulation of the forward transcript of the putative helicase, and the rdp1^- strain showed slightly increased accumulation with respect to wild-type (Figure 7b). The reverse strand did not accumulate in these three strains. Thus, transcriptional silencing of the putative helicase appeared to be relieved in RNAi mutants, implicating RNAi in the control of expression of this ORF.

The genome-wide survey of expressed genes in this study provided an opportunity to investigate the cellular response to loss of the gene for the telomerase catalytic subunit Trt1. A major finding was the tight correlation between the structures of chromosomes in survivors and gene expression profiles. Survivors with linear chromosomes had expression profiles remarkably similar to cells with canonical - yet shortened - telomeres, whereas cells with circular chromosomes maintained the upregulated expression of a significant number of genes that also had upregulated expression during senescence.

The stress response in survivors with circular chromosomes had significant overlaps with the S. pombe CESR and with the heat and oxidative stress responses. The CESR consists of genes that had upregulated expression in all or most responses to oxidative stress, heavy metal stress, heat shock, osmotic stress and DNA damage 16. The stress response may persist in survivors with circularized chromosomes because of impaired DNA segregation and DNA breakage and rearrangement. Indeed, compared with wild-type cells, survivors with circular chromosomes are larger and have slower growth rates, indicating that functions related to cell division are impaired 14.

Telomeric repeats contribute to recruiting the molecular components collectively involved in the protective capping of chromosome ends 203940. These repeats are maintained in the absence of telomerase in cells from diverse organisms that normally use telomerase (reviewed in 3). Interestingly, the survivors with linear chromosomes abated their stress response concomitant with the appearance of amplified telomeric and TAS repeats as rare survivors took over the population, suggesting that the repeats helped to ameliorate the stress response.

Neither cells in the growth curve that experienced shortened telomeres nor survivors with long telomeres displayed upregulation of telomeric gene expression, supporting the notion that telomeric length changes alone do not affect gene expression in S. pombe 19. In addition, in survivors with circular chromosomes, only eight microarray signals, corresponding to as few as two genes (due to cross-hybridization) near former telomeres had altered expression, although such changes might have been expected as a result of the large alterations in chromosome structure at these sites.

As in fission yeast, genes with changed expression in the budding yeast response to loss of telomerase had significant overlaps with genes whose expression was altered by environmental stresses such as heat shock, osmotic shock, dithiothreitol (DTT), nitrogen starvation and peroxide (1141 see also 42). A difference in the stress responses between the two yeasts was that in budding yeast a large but specific subset of the environmental stress-response genes persisted in survivors with linear chromosomes four days after crisis, whereas in fission yeast survivors with linear chromosomes, the stress response mostly abated by the fourth day after crisis (Figure 4b, day 12). The different yeast responses may be due to a fission yeast telomere structure that was not as strongly recognized as aberrant, perhaps mitigating a DNA-damage response. It is also possible that had budding yeast survivors been followed longer, providing a period for adaptation, the stress response would have subsided.

In fission yeast, the expression of a number of mitochondrial ATP synthase genes was upregulated (Table 1) with orthologs similarly induced in budding yeast. In both cases, the changes did not overlap with the DNA-damage responses of the yeasts, further supporting a link between short telomeres and alterations in the metabolic program suggested by Nautiyal et al. 11.

RecQ helicases have recently been implicated in telomerase-independent telomere maintenance in both S. cerevisiae and human ALT cells. BLM and WRN, human RecQ helicases associated with cancer and disease 31, have both been shown to associate with duplex telomere repeat binding protein TRF2 in vivo, and BLM co-localizes to telomeric foci exclusively in ALT cells 434445. The S. cerevisiae ortholog of human WRN and BLM, Sgs1, was also shown to be required for telomere elongation of type II survivors in the absence of telomerase 464748. The long, heterogeneous telomeres of S. pombe survivors with linear chromosomes are similar to those of S. cerevisiae survivors and human ALT cells, suggesting a role for RecQ helicases in fission yeast telomerase-independent telomere maintenance.

This is the first report to our knowledge of naturally occurring dh repeats outside of the centromeric and mating-type regions in fission yeast. We have presented several results that suggest that sub-telomeric dh repeats promote heterochromatin formation at the helicase locus. First, transcript from this ORF was only weakly expressed in wild-type cells as determined by RT-PCR (Figure 7a) (and was not detectable at all by northern hybridization, data not shown), consistent with transcriptional regulation of this ORF by heterochromatin. Second, expression of the putative helicase was robust in ago1^- and dcr1^- mutants, which would be expected if RNAi has a role in transcriptionally silencing this ORF. In trt1^- mutants experiencing genomic instability, we detected both forward and reverse transcripts of sub-telomeric dh repeats (Figure 7a). The presence of these complementary transcripts suggests the existence of dsRNA that had not been processed by the RNAi machinery, consistent with a lack of silencing at this locus. Intriguingly, after maximal expression of both strands on day 9 of the growth curve, subsequent downregulation was observed by day 15 (Figure 7a), consistent with restoration of silencing.

While the finding of homology with dh repeats at the sub-telomere was unexpected, dh repeats have been shown to function in silencing at sites outside of centromeres and the mating-type locus. Reporter genes fused to centromeric repeat fragments as short as 580 bp were silenced when integrated at ectopic locations in the genome 4950 and this silencing required the RNAi machinery 5152. The longest (nearly continuous) stretch of sequence with homology to dh repeats found in the helicase ORF was about 600 bp (Figure 6), presumably long enough to promote heterochromatin formation. In addition, RNAi-mediated silencing triggered by both a synthetic hairpin RNA and transposon long terminal repeats have been shown to induce heterochromatin formation away from centromeres and the mating-type locus 53.

In a separate study, telomeric silencing of a reporter gene and binding of Swi6 at the telomere were not affected in dcr1^-, ago1^- and rdp1^- mutants 54. The lack of an observed effect may have been due to the ability of telomeric repeats to recruit silencing factors. Indeed, telomeric heterochromatin is largely promoted by telomeric repeats. However, the study by Hall and co-workers 54 did report defective mitotic and meiotic telomere clustering in RNAi mutants, supporting a role for RNAi at telomeres.

Given the correlation between disruption of telomeric heterochromatin and expression of the helicase ORF, events other than telomere erosion that disrupt heterochromatin might also induce helicase expression.

The trt1⁺ and trt1^- cells used in this study were generated by sporulating S. pombe diploid strain G4 (h^-/h⁺ ade6-M210/ade6-M216 trt1⁺/trt1^-) on ME plates 18. The parent diploid strain was made heterozygous for trt1⁺ by using a standard two-step integration procedure 55 with a linearized plasmid containing about 1 kb each of the 5' and 3' flanking regions of the trt1⁺ ORF separated by HSV1-tk and KanMX4 56. The plasmid was linearized in the middle of the 3' flanking region with FseI and transformed using the lithium acetate method 57 into a diploid strain created by crossing PP68 (h^- ade6-M210) and PP69 (h⁺ ade6-M216). Cells were re-streaked twice on yeast extract low adenine (YEA) + geneticin plates 18 to select for stable genomic integrants, which were subsequently confirmed by Southern hybridization to a uniquely sized EcoRI restriction fragment 3' of trt1⁺ which was present only in integrants. Cells were then plated on YEA + 50 μM 5-fluorodeoxyuridine (5-FUdR) plates to select for those that had excised HSV1-tk, KanMX4 and the XbaI-XhoI fragment (around 5 kb) of trt1⁺ from their genomes. Random surviving colonies were screened for heterozygous diploids by Southern hybridization to the 3' region of the trt1⁺ KpnI restriction fragment. The heterozygous state was evidenced by hybridization signals to both full-length trt1⁺ and a shortened, nonfunctional version. Loss of markers was confirmed by lack of a Southern hybridization signal to HSV1-tk, and by lack of growth on YEA + geneticin plates.

After germination of G4 and growth of spores at 32°C for three days on YEA plates 18, plates were stored at 4°C while the genotypes of random colonies were determined. A portion of single colonies was used for crossing and visual inspection to identify those that had an h^- ade6-M210 genotype, which were further screened by Southern hybridization for the presence or absence of trt1⁺ (performed as described in 'Strain construction' below). Colonies were subsequently used as described in 'Growth curve', or alternatively used to create strains C1, C5, H1 and H2.

Strains C1, C5, H1 and H2 were created from four separate trt1^- colonies that were each successively re-streaked on YEA plates 15 times (with growth for 2 to 3 days at 32°C between re-streaks), to permit colonies to form without competition from faster-growing survivors with linear chromosomes. During this time cells were presumed to senesce and give rise to survivors. After the last re-streak, a single colony from each strain was randomly selected and used to prepare freeze stocks.

Three strains were grown: two wild-type isolates (h^- ade6-M210 trt1⁺) designated WT 3 and WT 5, and a single mutant isolate (h^- ade6-M210 trt1^-) designated as 'GC Day X', where X represents the day of the growth curve that cells were collected. Single colonies were used to inoculate 5-ml starter cultures in yeast extract full supplements (YES) medium 18 and grown for 24 h with shaking at 32°C. Cells were counted and used to inoculate 200-ml YES cultures in 500-ml Erlenmeyer flasks at 2.5 × 10⁴ cells/ml, and were grown in an incubator (Innova 4430, New Brunswick Scientific) with continuous shaking at 200 rpm at 32°C. Cell density was monitored by periodic counting, and a portion of the cells was harvested for microarray analysis and Southern hybridization when the density reached 3-5 × 10⁶ cells/ml (early log phase). Cells harvested at this point were referred to as day 1 of the growth curve. The unharvested cells were permitted to continue growing until 24 h from the time of inoculation, at which time cells were counted and used to inoculate a fresh 200-ml YES culture at 2.5 × 10⁴ cells/ml, and the process repeated for 15 days. To harvest cells for microarray analysis, a volume of culture containing approximately 1.6 × 10⁸ cells was gently centrifuged at room temperature (2,000 rpm for 2 min), the supernatant removed, and the cell pellet snap-frozen in liquid N₂. For Southern hybridization, approximately 2 × 10⁸ cells were collected by centrifugation, washed twice in H₂O, and snap-frozen in liquid N₂. A portion of cells for pulsed-field gel analysis was also collected in the same manner as for Southern hybridization at the end of each 24 h period. trt1^- cells were collected daily, WT 3 and WT 5 on day 1, and WT 5 on day 15.

Cells were streaked onto YEA plates from freeze stocks, grown for 3 days at 32°C, and single colonies used to inoculate 5-ml starter cultures in YES medium. After 24 h, cells were counted and 200-ml YES cultures were inoculated at 2.5 × 10⁴ cells/ml, and cultures grown with constant shaking at 200 rpm at 32°C. When the cell density reached around 3 × 10⁶ cells/ml (early log phase), cells were collected as described in 'Growth curve' for microarray analysis, Southern hybridization and pulsed-field gel electrophoresis. Strains H1 and H2 are trt1^- isolates with circular chromosomes, as evidenced by pulsed-field gel electrophoresis (data not shown)

DNA from approximately 2 × 10⁸ S. pombe cells was prepared as described 18. After digestion with either EcoRI or HindIII, the DNA was subjected to electrophoresis on a 1% agarose gel in 1 × TBE (90 mM Tris, 90 mM borate, 2 mM EDTA pH 8.3). DNA was denatured by sodium hydroxide treatment and transferred to a nylon membrane (Hybond-N+ membrane, Amersham) by capillary transfer in 10 × SSC (1.5 M NaCl, 0.15 M sodium citrate). DNA was immobilized on the membrane by irradiation with 120 mJ/cm² at 254 nm in a UVStratalinker1800 (Stratagene). For molecular weight markers, the 1 kb DNA ladder (New England Biolabs) was labeled by filling in 5' overhangs with [α-³²P]dATP using DNA polymerase I Klenow fragment. Probes for pol1⁺, act1⁺, the putative helicase and the C, I, L and M chromosome fragments were generated by PCR amplification from a genomic DNA template and were gel purified. Probes were labeled by random-primed transcription of PCR products with the use of [α-³²P]dCTP and High Prime Mix (Boehringer Mannheim). Probes specific for the telomeric and telomere-associated sequences were created with the use of gel-purified fragments of pNSU70 58.

Cells (approximately 1 × 10⁸) were collected as described above. Plug preparation, chromosome digestion and electrophoresis were performed exactly as described 18. DNA was visualized by staining with ethidium bromide (1 μg/ml) for 30 min. The gel was then irradiated with 120 mJ/cm² at 254 nm in a UVStratalinker1800 to nick the DNA, treated with HCl, NaOH and neutralization buffer, and processed as described in 'Southern hybridization'.

RNA was prepared as for microarray analysis and used for RT-PCR (OneStep RT-PCR kit, Qiagen). First-strand cDNA synthesis was performed using primers complementary to either the forward or reverse strands. Both primers were present in subsequent cycles of PCR amplification after heat inactivation of reverse transcriptase at 95°C for 15 min. The control reaction lacking reverse transcriptase (act1⁺ sense, -RT) was not subjected to first-strand cDNA synthesis, but was otherwise treated identically.

The PCR primers used to generate probes C, I, L, and M have been published previously 14. The PCR primers spanning the regions described in Figure 6 were:

P_5': 5'-CTTCAAAAACTGCTAGAGATATCGCCGG-3' and

5'-GTACTGGTAGTCCTCTGATGTATGGG-3'

P_3': 5'-ATGCCCCGTACGCTTATCTA-3' and 5'-TTTGCCTTTCTAGCCCATGA-3'

P_dh: 5'-CAACACCAATACTGACGATGATG-3' and 5'-GCAATAGAACCAGCGGTTTG-3'

Primers for centromeric dh repeats have been published previously 33.

Whole-cell RNA was isolated from S. pombe cell pellets (~1.6 × 10⁸ cells) by hot-phenol extraction and purification with RNeasy columns (Qiagen) following a published protocol 59. Aliquots (10 μg) were made (henceforth referred to as 'sample RNA') and RNA quality was assessed by UV absorbance, by agarose gel electrophoresis to confirm intact rRNA bands, and by northern hybridization to act1⁺. A reference pool consisting of RNA from each sample was made, comprising 76% trt1^- cells and 24% trt1⁺ cells. This pool was divided into 10 μg aliquots (henceforth referred to as 'reference RNA') and used as the reference RNA in all hybridization experiments reported here.

A single large batch of YES medium was made at the start of the study and used to culture all cells analyzed by microarrays to prevent batch-to-batch medium variations that might yield artifactual microarray results.

The procedures performed and the S. pombe microarrays used have been described previously 59. Whole-cell RNA (10 μg) was labeled by directly incorporating either Cy3-dCTP (reference RNA) or Cy5-dCTP (sample RNA) through reverse transcription. The resulting cDNA was hybridized onto DNA microarrays containing spotted PCR products for over 5,269 different genes and genomic elements printed in duplicate on glass slides representing 99.9% of all known and predicted fission yeast genes. Microarrays were scanned using a GenePix 4000B laser scanner (Axon Instruments) and analyzed with GenePix Pro software. Low-quality signals were filtered out, and data were normalized using a customized Perl script (local adjustment of median of ratios to one within running windows of 1,000 spots).

Normalized data (Cy5/Cy3 ratios) were evaluated using GeneSpring (Silicon Genetics). All gene-expression values were normalized to the average of two trt1⁺ biological replicates (strains WT 3 and WT 5) collected on day 1 of the growth curve. Experiments and genes were clustered in GeneSpring using the Pearson correlation around zero (termed the Standard correlation in GeneSpring) with a minimum distance of 0.001 and a separation ratio of 1. Gene annotations were taken from GeneDB at the Wellcome Trust Sanger Institute 60. Lists of genes whose expression changed in the fission yeast stress response 16 were taken from the authors' website 61. BLAST searches were performed using the NCBI BLAST server 62.

The density of genes with changed regulation along the chromosome was determined by using a running window of 20 consecutive genes along each chromosome 63. For each window, the probability of obtaining the observed results by chance was calculated using the hypergeometric distribution.

There were two microarray signals - SPAC212.06 (a pseudogene) and the reverse transcript of centromeric dh repeats - that we believe were due to cross-hybridization with the SPAC212.11 transcript (or transcripts from identical ORFs, see text). Cross-hybridization becomes apparent with array element sequence identities higher than about 70% 59. The SPAC212.11 transcript is capable of hybridizing to the entire SPAC212.06 microarray probe (99% sequence identity), but the SPAC212.06 transcript does not contain the sequence required to hybridize to the SPAC212.11 microarray probe. The SPAC212.11 transcript also has a high degree of homology with the dh repeat microarray probe sequence (83% sequence identity), and both the SPAC212.11 and dh repeat transcripts are expected to be capable of hybridizing to each other's microarray probes. Significant levels of forward and reverse centromeric dh repeat transcripts could not be detected using RT-PCR with RNA from days 1, 7, 9 and 15 of the growth curve (Figure 7a) (indicating they could not hybridize to the microarray), although helicase RNA was detected by both RT-PCR and northern hybridization (Figure 7a and data not shown).

Twenty-one microarrays were used in this study, representing two wild-type biological repeats, 15 days of the growth curve, and four strains with circularized chromosomes. The complete raw and normalized data sets are available from ArrayExpress 64 (Accession number: E-MEXP-201).