Interspecies systems biology uncovers metabolites affecting C. elegans gene expression and life history traits.

Diet greatly influences gene expression and physiology. In mammals, elucidating the effects and mechanisms of individual nutrients is challenging due to the complexity of both the animal and its diet. Here, we used an interspecies systems biology approach with Caenorhabditis elegans and two of its bacterial diets, Escherichia coli and Comamonas aquatica, to identify metabolites that affect the animal's gene expression and physiology. We identify vitamin B12 as the major dilutable metabolite provided by Comamonas aq. that regulates gene expression, accelerates development, and reduces fertility but does not affect lifespan. We find that vitamin B12 has a dual role in the animal: it affects development and fertility via the methionine/S-Adenosylmethionine (SAM) cycle and breaks down the short-chain fatty acid propionic acid, preventing its toxic buildup. Our interspecies systems biology approach provides a paradigm for understanding complex interactions between diet and physiology.

A cell cycle-dependent regulatory circuit composed of 53BP1-RIF1 and BRCA1-CtIP controls DNA repair pathway choice.

DNA double-strand break (DSB) repair pathway choice is governed by the opposing activities of 53BP1 and BRCA1. 53BP1 stimulates nonhomologous end joining (NHEJ), whereas BRCA1 promotes end resection and homologous recombination (HR). Here we show that 53BP1 is an inhibitor of BRCA1 accumulation at DSB sites, specifically in the G1 phase of the cell cycle. ATM-dependent phosphorylation of 53BP1 physically recruits RIF1 to DSB sites, and we identify RIF1 as the critical effector of 53BP1 during DSB repair. Remarkably, RIF1 accumulation at DSB sites is strongly antagonized by BRCA1 and its interacting partner CtIP. Lastly, we show that depletion of RIF1 is able to restore end resection and RAD51 loading in BRCA1-depleted cells. This work therefore identifies a cell cycle-regulated circuit, underpinned by RIF1 and BRCA1, that governs DSB repair pathway choice to ensure that NHEJ dominates in G1 and HR is favored from S phase onward.

Fitness and Productivity Increase with Ecotypic Diversity among Escherichia coli Strains That Coevolved in a Simple, Constant Environment.

The productivity of a biological community often correlates with its diversity. In the microbial world this phenomenon can sometimes be explained by positive, density-dependent interactions such as cross-feeding and syntrophy. These metabolic interactions help account for the astonishing variety of microbial life and drive many of the biogeochemical cycles without which life as we know it could not exist. While it is difficult to recapitulate experimentally how these interactions evolved among multiple taxa, we can explore in the laboratory how they arise within one. These experiments provide insight into how different bacterial ecotypes evolve and from these, possibly new "species." We have previously shown that in a simple, constant environment a single clone of Escherichia coli can give rise to a consortium of genetically and phenotypically differentiated strains, in effect, a set of ecotypes, that coexist by cross-feeding. We marked these different ecotypes and their shared ancestor by integrating fluorescent protein into their genomes and then used flow cytometry to show that each evolved strain is more fit than the shared ancestor, that pairs of evolved strains are fitter still, and that the entire consortium is the fittest of all. We further demonstrate that the rank order of fitness values agrees with estimates of yield, indicating that an experimentally evolved consortium more efficiently converts primary and secondary resources to offspring than its ancestor or any member acting in isolation.IMPORTANCE Polymicrobial consortia occur in both environmental and clinical settings. In many cases, diversity and productivity correlate in these consortia, especially when sustained by positive, density-dependent interactions. However, the evolutionary history of such entities is typically obscure, making it difficult to establish the relative fitness of consortium partners and to use those data to illuminate the diversity-productivity relationship. Here, we dissect an Escherichia coli consortium that evolved under continuous glucose limitation in the laboratory from a single common ancestor. We show that a partnership consisting of cross-feeding ecotypes is better able to secure primary and secondary resources and to convert those resources to offspring than the ancestral clone. Such interactions may be a prelude to a special form of syntrophy and are likely determinants of microbial community structure in nature, including those having clinical significance such as chronic infections.

Drosophila larvae synthesize the putative oncometabolite L-2-hydroxyglutarate during normal developmental growth.

L-2-hydroxyglutarate (L-2HG) has emerged as a putative oncometabolite that is capable of inhibiting enzymes involved in metabolism, chromatin modification, and cell differentiation. However, despite the ability of L-2HG to interfere with a broad range of cellular processes, this molecule is often characterized as a metabolic waste product. Here, we demonstrate that Drosophila larvae use the metabolic conditions established by aerobic glycolysis to both synthesize and accumulate high concentrations of L-2HG during normal developmental growth. A majority of the larval L-2HG pool is derived from glucose and dependent on the Drosophila estrogen-related receptor (dERR), which promotes L-2HG synthesis by up-regulating expression of the Drosophila homolog of lactate dehydrogenase (dLdh). We also show that dLDH is both necessary and sufficient for directly synthesizing L-2HG and the Drosophila homolog of L-2-hydroxyglutarate dehydrogenase (dL2HGDH), which encodes the enzyme that breaks down L-2HG, is required for stage-specific degradation of the L-2HG pool. In addition, dLDH also indirectly promotes L-2HG accumulation via synthesis of lactate, which activates a metabolic feed-forward mechanism that inhibits dL2HGDH activity and stabilizes L-2HG levels. Finally, we use a genetic approach to demonstrate that dLDH and L-2HG influence position effect variegation and DNA methylation, suggesting that this compound serves to coordinate glycolytic flux with epigenetic modifications. Overall, our studies demonstrate that growing animal tissues synthesize L-2HG in a controlled manner, reveal a mechanism that coordinates glucose catabolism with L-2HG synthesis, and establish the fly as a unique model system for studying the endogenous functions of L-2HG during cell growth and proliferation.

Functional genetic discovery of enzymes using full-scan mass spectrometry metabolomics 1.

Our understanding of metabolic networks is incomplete, and new enzymatic activities await discovery in well-studied organisms. Mass spectrometric measurement of cellular metabolites reveals compounds inside cells that are unexplained by current maps of metabolic reactions, and existing computational models are unable to account for all activities observed within cells. Additional large-scale genetic and biochemical approaches are required to elucidate metabolic gene function. We have used full-scan mass spectrometry metabolomics of polar small molecules to examine deletion mutants of candidate enzymes in the model yeast Saccharomyces cerevisiae. We report the identification of 25 genes whose deletion results in focal metabolic changes consistent with loss of enzymatic activity and describe the informatic approaches used to enrich for candidate enzymes from uncharacterized open reading frames. Triumphs and pitfalls of metabolic phenotyping screens are discussed, including estimates of the frequency of uncharacterized eukaryotic genes that affect metabolism and key issues to consider when searching for new enzymatic functions in other organisms.

Mutational landscape of EGFR-, MYC-, and Kras-driven genetically engineered mouse models of lung adenocarcinoma.

Genetically engineered mouse models (GEMMs) of cancer are increasingly being used to assess putative driver mutations identified by large-scale sequencing of human cancer genomes. To accurately interpret experiments that introduce additional mutations, an understanding of the somatic genetic profile and evolution of GEMM tumors is necessary. Here, we performed whole-exome sequencing of tumors from three GEMMs of lung adenocarcinoma driven by mutant epidermal growth factor receptor (EGFR), mutant Kirsten rat sarcoma viral oncogene homolog (Kras), or overexpression of MYC proto-oncogene. Tumors from EGFR- and Kras-driven models exhibited, respectively, 0.02 and 0.07 nonsynonymous mutations per megabase, a dramatically lower average mutational frequency than observed in human lung adenocarcinomas. Tumors from models driven by strong cancer drivers (mutant EGFR and Kras) harbored few mutations in known cancer genes, whereas tumors driven by MYC, a weaker initiating oncogene in the murine lung, acquired recurrent clonal oncogenic Kras mutations. In addition, although EGFR- and Kras-driven models both exhibited recurrent whole-chromosome DNA copy number alterations, the specific chromosomes altered by gain or loss were different in each model. These data demonstrate that GEMM tumors exhibit relatively simple somatic genotypes compared with human cancers of a similar type, making these autochthonous model systems useful for additive engineering approaches to assess the potential of novel mutations on tumorigenesis, cancer progression, and drug sensitivity.

Heritability and genetic basis of protein level variation in an outbred population.

The genetic basis of heritable traits has been studied for decades. Although recent mapping efforts have elucidated genetic determinants of transcript levels, mapping of protein abundance has lagged. Here, we analyze levels of 4084 GFP-tagged yeast proteins in the progeny of a cross between a laboratory and a wild strain using flow cytometry and high-content microscopy. The genotype of trans variants contributed little to protein level variation between individual cells but explained >50% of the variance in the population's average protein abundance for half of the GFP fusions tested. To map trans-acting factors responsible, we performed flow sorting and bulk segregant analysis of 25 proteins, finding a median of five protein quantitative trait loci (pQTLs) per GFP fusion. Further, we find that cis-acting variants predominate; the genotype of a gene and its surrounding region had a large effect on protein level six times more frequently than the rest of the genome combined. We present evidence for both shared and independent genetic control of transcript and protein abundance: More than half of the expression QTLs (eQTLs) contribute to changes in protein levels of regulated genes, but several pQTLs do not affect their cognate transcript levels. Allele replacements of genes known to underlie trans eQTL hotspots confirmed the correlation of effects on mRNA and protein levels. This study represents the first genome-scale measurement of genetic contribution to protein levels in single cells and populations, identifies more than a hundred trans pQTLs, and validates the propagation of effects associated with transcript variation to protein abundance.

TFIIH and P-TEFb coordinate transcription with capping enzyme recruitment at specific genes in fission yeast.

Cyclin-dependent kinases (CDKs) are subunits of transcription factor (TF) IIH and positive transcription elongation factor b (P-TEFb). To define their functions, we mutated the TFIIH-associated kinase Mcs6 and P-TEFb homologs Cdk9 and Lsk1 of fission yeast, making them sensitive to inhibition by bulky purine analogs. Selective inhibition of Mcs6 or Cdk9 blocks cell division, alters RNA polymerase (Pol) II carboxyl-terminal domain (CTD) phosphorylation, and represses specific, overlapping subsets of transcripts. At a common target gene, both CDKs must be active for normal Pol II occupancy, and Spt5-a CDK substrate and regulator of elongation-accumulates disproportionately to Pol II when either kinase is inhibited. In contrast, Mcs6 activity is sufficient-and necessary-to recruit the Cdk9/Pcm1 (mRNA cap methyltransferase) complex. In vitro, phosphorylation of the CTD by Mcs6 stimulates subsequent phosphorylation by Cdk9. We propose that TFIIH primes the CTD and promotes recruitment of P-TEFb/Pcm1, serving to couple elongation and capping of select pre-mRNAs.

Cell cycle-regulated oscillator coordinates core histone gene transcription through histone acetylation.

DNA replication occurs during the synthetic (S) phase of the eukaryotic cell cycle and features a dramatic induction of histone gene expression for concomitant chromatin assembly. Ectopic production of core histones outside of S phase is toxic, underscoring the critical importance of regulatory pathways that ensure proper expression of histone genes. Several regulators of histone gene expression in the budding yeast Saccharomyces cerevisiae are known, yet the key oscillator responsible for restricting gene expression to S phase has remained elusive. Here, we show that suppressor of Ty (Spt)10, a putative histone acetyltransferase, and its binding partner Spt21 are key determinants of S-phase-specific histone gene expression. We show that Spt21 abundance is restricted to S phase in part by anaphase promoting complex Cdc20-homologue 1 (APC(Cdh1)) and that it is recruited to histone gene promoters in S phase by Spt10. There, Spt21-Spt10 enables the recruitment of a cascade of regulators, including histone chaperones and the histone-acetyltransferase general control nonderepressible (Gcn) 5, which we hypothesize lead to histone acetylation and consequent transcription activation.

Mapping the pericentric heterochromatin by comparative genomic hybridization analysis and chromosome deletions in Drosophila melanogaster.

Heterochromatin represents a significant portion of eukaryotic genomes and has essential structural and regulatory functions. Its molecular organization is largely unknown due to difficulties in sequencing through and assembling repetitive sequences enriched in the heterochromatin. Here we developed a novel strategy using chromosomal rearrangements and embryonic phenotypes to position unmapped Drosophila melanogaster heterochromatic sequence to specific chromosomal regions. By excluding sequences that can be mapped to the assembled euchromatic arms, we identified sequences that are specific to heterochromatin and used them to design heterochromatin specific probes ("H-probes") for microarray. By comparative genomic hybridization (CGH) analyses of embryos deficient for each chromosome or chromosome arm, we were able to map most of our H-probes to specific chromosome arms. We also positioned sequences mapped to the second and X chromosomes to finer intervals by analyzing smaller deletions with breakpoints in heterochromatin. Using this approach, we were able to map >40% (13.9 Mb) of the previously unmapped heterochromatin sequences assembled by the whole-genome sequencing effort on arm U and arm Uextra to specific locations. We also identified and mapped 110 kb of novel heterochromatic sequences. Subsequent analyses revealed that sequences located within different heterochromatic regions have distinct properties, such as sequence composition, degree of repetitiveness, and level of underreplication in polytenized tissues. Surprisingly, although heterochromatin is generally considered to be transcriptionally silent, we detected region-specific temporal patterns of transcription in heterochromatin during oogenesis and early embryonic development. Our study provides a useful approach to elucidate the molecular organization and function of heterochromatin and reveals region-specific variation of heterochromatin.

Minotaur is critical for primary piRNA biogenesis.

Piwi proteins and their associated small RNAs are essential for fertility in animals. In part, this is due to their roles in guarding germ cell genomes against the activity of mobile genetic elements. piRNA populations direct Piwi proteins to silence transposon targets and, as such, form a molecular code that discriminates transposons from endogenous genes. Information ultimately carried by piRNAs is encoded within genomic loci, termed piRNA clusters. These give rise to long, single-stranded, primary transcripts that are processed into piRNAs. Despite the biological importance of this pathway, neither the characteristics that define a locus as a source of piRNAs nor the mechanisms that catalyze primary piRNA biogenesis are well understood. We searched an EMS-mutant collection annotated for fertility phenotypes for genes involved in the piRNA pathway. Twenty-seven homozygous sterile strains showed transposon-silencing defects. One of these, which strongly impacted primary piRNA biogenesis, harbored a causal mutation in CG5508, a member of the Drosophila glycerol-3-phosphate O-acetyltransferase (GPAT) family. These enzymes catalyze the first acylation step on the path to the production of phosphatidic acid (PA). Though this pointed strongly to a function for phospholipid signaling in the piRNA pathway, a mutant form of CG5508, which lacks the GPAT active site, still functions in piRNA biogenesis. We have named this new biogenesis factor Minotaur.

MOB rules: Antibiotic Exposure Reprograms Metabolism to Mobilize Bacillus subtilis in Competitive Interactions.

Antibiotics have dose-dependent effects on exposed bacteria. The medicinal use of antibiotics relies on their growth-inhibitory activities at sufficient concentrations. At subinhibitory concentrations, exposure effects vary widely among different antibiotics and bacteria. Bacillus subtilis responds to bacteriostatic translation inhibitors by mobilizing a population of cells (MOB-Mobilized Bacillus) to spread across a surface. How B. subtilis regulates the antibiotic-induced mobilization is not known. In this study, we used chloramphenicol to identify regulatory functions that B. subtilis requires to coordinate cell mobilization following subinhibitory exposure. We measured changes in gene expression and metabolism and mapped the results to a network of regulatory proteins that direct the mobile response. Our data reveal that several transcriptional regulators coordinately control the reprogramming of metabolism to support mobilization. The network regulates changes in glycolysis, nucleotide metabolism, and amino acid metabolism that are signature features of the mobilized population. Among the hundreds of genes with changing expression, we identified two, pdhA and pucA, where the magnitudes of their changes in expression, and in the abundance of associated metabolites, reveal hallmark metabolic features of the mobilized population. Using reporters of pdhA and pucA expression, we visualized the separation of major branches of metabolism in different regions of the mobilized population. Our results reveal a regulated response to chloramphenicol exposure that enables a population of bacteria in different metabolic states to mount a coordinated mobile response.

Daughter-specific transcription factors regulate cell size control in budding yeast.

In budding yeast, asymmetric cell division yields a larger mother and a smaller daughter cell, which transcribe different genes due to the daughter-specific transcription factors Ace2 and Ash1. Cell size control at the Start checkpoint has long been considered to be a main regulator of the length of the G1 phase of the cell cycle, resulting in longer G1 in the smaller daughter cells. Our recent data confirmed this concept using quantitative time-lapse microscopy. However, it has been proposed that daughter-specific, Ace2-dependent repression of expression of the G1 cyclin CLN3 had a dominant role in delaying daughters in G1. We wanted to reconcile these two divergent perspectives on the origin of long daughter G1 times. We quantified size control using single-cell time-lapse imaging of fluorescently labeled budding yeast, in the presence or absence of the daughter-specific transcriptional regulators Ace2 and Ash1. Ace2 and Ash1 are not required for efficient size control, but they shift the domain of efficient size control to larger cell size, thus increasing cell size requirement for Start in daughters. Microarray and chromatin immunoprecipitation experiments show that Ace2 and Ash1 are direct transcriptional regulators of the G1 cyclin gene CLN3. Quantification of cell size control in cells expressing titrated levels of Cln3 from ectopic promoters, and from cells with mutated Ace2 and Ash1 sites in the CLN3 promoter, showed that regulation of CLN3 expression by Ace2 and Ash1 can account for the differential regulation of Start in response to cell size in mothers and daughters. We show how daughter-specific transcriptional programs can interact with intrinsic cell size control to differentially regulate Start in mother and daughter cells. This work demonstrates mechanistically how asymmetric localization of cell fate determinants results in cell-type-specific regulation of the cell cycle.

Identification of GIG1, a GlcNAc-induced gene in Candida albicans needed for normal sensitivity to the chitin synthase inhibitor nikkomycin Z.

The amino sugar N-acetylglucosamine (GlcNAc) is known to be an important structural component of cells from bacteria to humans, but its roles in cell signaling are less well understood. GlcNAc induces two pathways in the human fungal pathogen Candida albicans. One activates cyclic AMP (cAMP) signaling, which stimulates the formation of hyphal cells and the expression of virulence genes, and the other pathway induces genes needed to catabolize GlcNAc. Microarray analysis of gene expression was carried out under four different conditions in order to characterize the transcriptional changes induced by GlcNAc. The most highly induced genes include those that encode a GlcNAc transporter (NGT1) and the GlcNAc catabolic enzymes (HXK1, DAC1, and NAG1). GlcNAc also activated most of the genes whose expression is increased when cells are triggered with other stimuli to form hyphae. Surprisingly, GlcNAc also induced a subset of genes that are regulated by galactose (GAL1, GAL7, and GAL10), which may be due to cross talk between signaling pathways. A novel GlcNAc-induced gene, GIG1, which is not essential for GlcNAc catabolism or the induction of hyphae, was identified. However, a Gig1-green fluorescent protein (GFP) fusion protein was specifically induced by GlcNAc, and not by other sugars. Gig1-GFP localized to the cytoplasm, where GlcNAc metabolism occurs. Significantly, a gig1Δ mutant displayed increased resistance to nikkomycin Z, which inhibits chitin synthase from converting UDP-GlcNAc into cell wall chitin. Gig1 is highly conserved in fungi, especially those that contain GlcNAc catabolic genes. These results implicate Gig1 in GlcNAc metabolism.

Diverse conditions support near-zero growth in yeast: Implications for the study of cell lifespan.

Baker's yeast has a finite lifespan and ages in two ways: a mother cell can only divide so many times (its replicative lifespan), and a non-dividing cell can only live so long (its chronological lifespan). Wild and laboratory yeast strains exhibit natural variation for each type of lifespan, and the genetic basis for this variation has been generalized to other eukaryotes, including metazoans. To date, yeast chronological lifespan has chiefly been studied in relation to the rate and mode of functional decline among non-dividing cells in nutrient-depleted batch culture. However, this culture method does not accurately capture two major classes of long-lived metazoan cells: cells that are terminally differentiated and metabolically active for periods that approximate animal lifespan (e.g. cardiac myocytes), and cells that are pluripotent and metabolically quiescent (e.g. stem cells). Here, we consider alternative ways of cultivating Saccharomyces cerevisiae so that these different metabolic states can be explored in non-dividing cells: (i) yeast cultured as giant colonies on semi-solid agar, (ii) yeast cultured in retentostats and provided sufficient nutrients to meet minimal energy requirements, and (iii) yeast encapsulated in a semisolid matrix and fed ad libitum in bioreactors. We review the physiology of yeast cultured under each of these conditions, and explore their potential to provide unique insights into determinants of chronological lifespan in the cells of higher eukaryotes.

The cell cycle-regulated genes of Schizosaccharomyces pombe.

Many genes are regulated as an innate part of the eukaryotic cell cycle, and a complex transcriptional network helps enable the cyclic behavior of dividing cells. This transcriptional network has been studied in Saccharomyces cerevisiae (budding yeast) and elsewhere. To provide more perspective on these regulatory mechanisms, we have used microarrays to measure gene expression through the cell cycle of Schizosaccharomyces pombe (fission yeast). The 750 genes with the most significant oscillations were identified and analyzed. There were two broad waves of cell cycle transcription, one in early/mid G2 phase, and the other near the G2/M transition. The early/mid G2 wave included many genes involved in ribosome biogenesis, possibly explaining the cell cycle oscillation in protein synthesis in S. pombe. The G2/M wave included at least three distinctly regulated clusters of genes: one large cluster including mitosis, mitotic exit, and cell separation functions, one small cluster dedicated to DNA replication, and another small cluster dedicated to cytokinesis and division. S. pombe cell cycle genes have relatively long, complex promoters containing groups of multiple DNA sequence motifs, often of two, three, or more different kinds. Many of the genes, transcription factors, and regulatory mechanisms are conserved between S. pombe and S. cerevisiae. Finally, we found preliminary evidence for a nearly genome-wide oscillation in gene expression: 2,000 or more genes undergo slight oscillations in expression as a function of the cell cycle, although whether this is adaptive, or incidental to other events in the cell, such as chromatin condensation, we do not know.

Reporter-Based Synthetic Genetic Array Analysis: A Functional Genomics Approach for Investigating Transcript or Protein Abundance Using Fluorescent Proteins in Saccharomyces cerevisiae.

Fluorescent reporter genes have long been used to quantify various cell features such as transcript and protein abundance. Here, we describe a method, reporter synthetic genetic array (R-SGA) analysis, which allows for the simultaneous quantification of any fluorescent protein readout in thousands of yeast strains using an automated pipeline. R-SGA combines a fluorescent reporter system with standard SGA analysis and can be used to examine any array-based strain collection available to the yeast community. This protocol describes the R-SGA methodology for screening different arrays of yeast mutants including the deletion collection, a collection of temperature-sensitive strains for the assessment of essential yeast genes and a collection of inducible overexpression strains. We also present an alternative pipeline for the analysis of R-SGA output strains using flow cytometry of cells in liquid culture. Data normalization for both pipelines is discussed.

Checkpoint independence of most DNA replication origins in fission yeast.

BACKGROUND: In budding yeast, the replication checkpoint slows progress through S phase by inhibiting replication origin firing. In mammals, the replication checkpoint inhibits both origin firing and replication fork movement. To find out which strategy is employed in the fission yeast, Schizosaccharomyces pombe, we used microarrays to investigate the use of origins by wild-type and checkpoint-mutant strains in the presence of hydroxyurea (HU), which limits the pool of deoxyribonucleoside triphosphates (dNTPs) and activates the replication checkpoint. The checkpoint-mutant cells carried deletions either of rad3 (which encodes the fission yeast homologue of ATR) or cds1 (which encodes the fission yeast homologue of Chk2). RESULTS: Our microarray results proved to be largely consistent with those independently obtained and recently published by three other laboratories. However, we were able to reconcile differences between the previous studies regarding the extent to which fission yeast replication origins are affected by the replication checkpoint. We found (consistent with the three previous studies after appropriate interpretation) that, in surprising contrast to budding yeast, most fission yeast origins, including both early- and late-firing origins, are not significantly affected by checkpoint mutations during replication in the presence of HU. A few origins (approximately 3%) behaved like those in budding yeast: they replicated earlier in the checkpoint mutants than in wild type. These were located primarily in the heterochromatic subtelomeric regions of chromosomes 1 and 2. Indeed, the subtelomeric regions defined by the strongest checkpoint restraint correspond precisely to previously mapped subtelomeric heterochromatin. This observation implies that subtelomeric heterochromatin in fission yeast differs from heterochromatin at centromeres, in the mating type region, and in ribosomal DNA, since these regions replicated at least as efficiently in wild-type cells as in checkpoint-mutant cells. CONCLUSION: The fact that approximately 97% of fission yeast replication origins - both early and late - are not significantly affected by replication checkpoint mutations in HU-treated cells suggests that (i) most late-firing origins are restrained from firing in HU-treated cells by at least one checkpoint-independent mechanism, and (ii) checkpoint-dependent slowing of S phase in fission yeast when DNA is damaged may be accomplished primarily by the slowing of replication forks.