Nonparametric methods for analyzing replication origins in genomewide data

Due to the advent of high-throughput genomic technology, it has become possible to monitor cellular activities on a genomewide basis. With these new methods, scientists can begin to address important biological questions. One such question involves the identification of replication origins, which are regions in the chromosomes where DNA replication is initiated. One hypothesis is that their locations are nonrandom throughout the genome. In this article, we analyze data from a recent yeast study in which candidate replication origins were profiled using cDNA microarrays to test this hypothesis. We find no evidence for such clustering.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/47937/1/10142_2004_Article_122.pd

Distinct Differences in Chromatin Structure at Subtelomeric X and Y' Elements in Budding Yeast

In Saccharomyces cerevisiae, all ends of telomeric DNA contain telomeric repeats of (TG1–3), but the number and position of subtelomeric X and Y' repeat elements vary. Using chromatin immunoprecipitation and genome-wide analyses, we here demonstrate that the subtelomeric X and Y' elements have distinct structural and functional properties. Y' elements are transcriptionally active and highly enriched in nucleosomes, whereas X elements are repressed and devoid of nucleosomes. In contrast to X elements, the Y' elements also lack the classical hallmarks of heterochromatin, such as high Sir3 and Rap1 occupancy as well as low levels of histone H4 lysine 16 acetylation. Our analyses suggest that the presence of X and Y' elements govern chromatin structure and transcription activity at individual chromosome ends

Missing value imputation improves clustering and interpretation of gene expression microarray data

<p>Abstract</p> <p>Background</p> <p>Missing values frequently pose problems in gene expression microarray experiments as they can hinder downstream analysis of the datasets. While several missing value imputation approaches are available to the microarray users and new ones are constantly being developed, there is no general consensus on how to choose between the different methods since their performance seems to vary drastically depending on the dataset being used.</p> <p>Results</p> <p>We show that this discrepancy can mostly be attributed to the way in which imputation methods have traditionally been developed and evaluated. By comparing a number of advanced imputation methods on recent microarray datasets, we show that even when there are marked differences in the measurement-level imputation accuracies across the datasets, these differences become negligible when the methods are evaluated in terms of how well they can reproduce the original gene clusters or their biological interpretations. Regardless of the evaluation approach, however, imputation always gave better results than ignoring missing data points or replacing them with zeros or average values, emphasizing the continued importance of using more advanced imputation methods.</p> <p>Conclusion</p> <p>The results demonstrate that, while missing values are still severely complicating microarray data analysis, their impact on the discovery of biologically meaningful gene groups can – up to a certain degree – be reduced by using readily available and relatively fast imputation methods, such as the Bayesian Principal Components Algorithm (BPCA).</p

Diversity of Eukaryotic DNA Replication Origins Revealed by Genome-Wide Analysis of Chromatin Structure

Eukaryotic DNA replication origins differ both in their efficiency and in the characteristic time during S phase when they become active. The biological basis for these differences remains unknown, but they could be a consequence of chromatin structure. The availability of genome-wide maps of nucleosome positions has led to an explosion of information about how nucleosomes are assembled at transcription start sites, but no similar maps exist for DNA replication origins. Here we combine high-resolution genome-wide nucleosome maps with comprehensive annotations of DNA replication origins to identify patterns of nucleosome occupancy at eukaryotic replication origins. On average, replication origins contain a nucleosome depleted region centered next to the ACS element, flanked on both sides by arrays of well-positioned nucleosomes. Our analysis identified DNA sequence properties that correlate with nucleosome occupancy at replication origins genome-wide and that are correlated with the nucleosome-depleted region. Clustering analysis of all annotated replication origins revealed a surprising diversity of nucleosome occupancy patterns. We provide evidence that the origin recognition complex, which binds to the origin, acts as a barrier element to position and phase nucleosomes on both sides of the origin. Finally, analysis of chromatin reconstituted in vitro reveals that origins are inherently nucleosome depleted. Together our data provide a comprehensive, genome-wide view of chromatin structure at replication origins and suggest a model of nucleosome positioning at replication origins in which the underlying sequence occludes nucleosomes to permit binding of the origin recognition complex, which then (likely in concert with nucleosome modifiers and remodelers) positions nucleosomes adjacent to the origin to promote replication origin function

Modeling the evolution of a classic genetic switch

Abstract Background The regulatory network underlying the yeast galactose-use pathway has emerged as a model system for the study of regulatory network evolution. Evidence has recently been provided for adaptive evolution in this network following a whole genome duplication event. An ancestral gene encoding a bi-functional galactokinase and co-inducer protein molecule has become subfunctionalized as paralogous genes (GAL1 and GAL3) in Saccharomyces cerevisiae, with most fitness gains being attributable to changes in cis- regulatory elements. However, the quantitative functional implications of the evolutionary changes in this regulatory network remain unexplored. Results We develop a modeling framework to examine the evolution of the GAL regulatory network. This enables us to translate molecular changes in the regulatory network to changes in quantitative network function. We computationally reconstruct an inferred ancestral version of the network and trace the evolutionary paths in the lineage leading to S. cerevisiae. We explore the evolutionary landscape of possible regulatory networks and find that the operation of intermediate networks leading to S. cerevisiae differs substantially depending on the order in which evolutionary changes accumulate; in particular, we systematically explore evolutionary paths and find that some network features cannot be optimized simultaneously. Conclusions We find that a computational modeling approach can be used to analyze the evolution of a well-studied regulatory network. Our results are consistent with several experimental studies of the evolutionary of the GAL regulatory network, including increased fitness in Saccharomyces due to duplication and adaptive regulatory divergence. The conceptual and computational tools that we have developed may be applicable in further studies of regulatory network evolution

Silent but Not Static: Accelerated Base-Pair Substitution in Silenced Chromatin of Budding Yeasts

Subtelomeric DNA in budding yeasts, like metazoan heterochromatin, is gene poor, repetitive, transiently silenced, and highly dynamic. The rapid evolution of subtelomeric regions is commonly thought to arise from transposon activity and increased recombination between repetitive elements. However, we found evidence of an additional factor in this diversification. We observed a surprising level of nucleotide divergence in transcriptionally silenced regions in inter-species comparisons of Saccharomyces yeasts. Likewise, intra-species analysis of polymorphisms also revealed increased SNP frequencies in both intergenic and synonymous coding positions of silenced DNA. This analysis suggested that silenced DNA in Saccharomyces cerevisiae and closely related species had increased single base-pair substitution that was likely due to the effects of the silencing machinery on DNA replication or repair

Nuclear Mitochondrial DNA Activates Replication in Saccharomyces cerevisiae

The nuclear genome of eukaryotes is colonized by DNA fragments of mitochondrial origin, called NUMTs. These insertions have been associated with a variety of germ-line diseases in humans. The significance of this uptake of potentially dangerous sequences into the nuclear genome is unclear. Here we provide functional evidence that sequences of mitochondrial origin promote nuclear DNA replication in Saccharomyces cerevisiae. We show that NUMTs are rich in key autonomously replicating sequence (ARS) consensus motifs, whose mutation results in the reduction or loss of DNA replication activity. Furthermore, 2D-gel analysis of the mrc1 mutant exposed to hydroxyurea shows that several NUMTs function as late chromosomal origins. We also show that NUMTs located close to or within ARS provide key sequence elements for replication. Thus NUMTs can act as independent origins, when inserted in an appropriate genomic context or affect the efficiency of pre-existing origins. These findings show that migratory mitochondrial DNAs can impact on the replication of the nuclear region they are inserted in

Epigenetically-Inherited Centromere and Neocentromere DNA Replicates Earliest in S-Phase

Eukaryotic centromeres are maintained at specific chromosomal sites over many generations. In the budding yeast Saccharomyces cerevisiae, centromeres are genetic elements defined by a DNA sequence that is both necessary and sufficient for function; whereas, in most other eukaryotes, centromeres are maintained by poorly characterized epigenetic mechanisms in which DNA has a less definitive role. Here we use the pathogenic yeast Candida albicans as a model organism to study the DNA replication properties of centromeric DNA. By determining the genome-wide replication timing program of the C. albicans genome, we discovered that each centromere is associated with a replication origin that is the first to fire on its respective chromosome. Importantly, epigenetic formation of new ectopic centromeres (neocentromeres) was accompanied by shifts in replication timing, such that a neocentromere became the first to replicate and became associated with origin recognition complex (ORC) components. Furthermore, changing the level of the centromere-specific histone H3 isoform led to a concomitant change in levels of ORC association with centromere regions, further supporting the idea that centromere proteins determine origin activity. Finally, analysis of centromere-associated DNA revealed a replication-dependent sequence pattern characteristic of constitutively active replication origins. This strand-biased pattern is conserved, together with centromere position, among related strains and species, in a manner independent of primary DNA sequence. Thus, inheritance of centromere position is correlated with a constitutively active origin of replication that fires at a distinct early time. We suggest a model in which the distinct timing of DNA replication serves as an epigenetic mechanism for the inheritance of centromere position

Cdc45 Limits Replicon Usage from a Low Density of preRCs in Mammalian Cells

Little is known about mammalian preRC stoichiometry, the number of preRCs on chromosomes, and how this relates to replicon size and usage. We show here that, on average, each 100-kb of the mammalian genome contains a preRC composed of approximately one ORC hexamer, 4–5 MCM hexamers, and 2 Cdc6. Relative to these subunits, ∼0.35 total molecules of the pre-Initiation Complex factor Cdc45 are present. Thus, based on ORC availability, somatic cells contain ∼70,000 preRCs of this average total stoichiometry, although subunits may not be juxtaposed with each other. Except for ORC, the chromatin-bound complement of preRC subunits is even lower. Cdc45 is present at very low levels relative to the preRC subunits, but is highly stable, and the same limited number of stable Cdc45 molecules are present from the beginning of S-phase to its completion. Efforts to artificially increase Cdc45 levels through ectopic expression block cell growth. However, microinjection of excess purified Cdc45 into S-phase nuclei activates additional replication foci by three-fold, indicating that Cdc45 functions to activate dormant preRCs and is rate-limiting for somatic replicon usage. Paradoxically, although Cdc45 colocalizes in vivo with some MCM sites and is rate-limiting for DNA replication to occur, neither Cdc45 nor MCMs colocalize with active replication sites. Embryonic metazoan chromatin consists of small replicons that are used efficiently via an excess of preRC subunits. In contrast, somatic mammalian cells contain a low density of preRCs, each containing only a few MCMs that compete for limiting amounts of Cdc45. This provides a molecular explanation why, relative to embryonic replicon dynamics, somatic replicons are, on average, larger and origin efficiency tends to be lower. The stable, continuous, and rate-limiting nature of Cdc45 suggests that Cdc45 contributes to the staggering of replicon usage throughout S-phase, and that replicon activation requires reutilization of existing Cdc45 during S-phase

ATP-dependent chromatin remodeling shapes the DNA replication landscape.

The eukaryotic DNA replication machinery must traverse every nucleosome in the genome during S phase. As nucleosomes are generally inhibitory to DNA-dependent processes, chromatin structure must undergo extensive reorganization to facilitate DNA synthesis. However, the identity of chromatin-remodeling factors involved in replication and how they affect DNA synthesis is largely unknown. Here we show that two highly conserved ATP-dependent chromatin-remodeling complexes in Saccharomyces cerevisiae, Isw2 and Ino80, function in parallel to promote replication fork progression. As a result, Isw2 and Ino80 have especially important roles for replication of late-replicating regions during periods of replication stress. Both Isw2 and Ino80 complexes are enriched at sites of replication, suggesting that these complexes act directly to promote fork progression. These findings identify ATP-dependent chromatin-remodeling complexes that promote DNA replication and define a specific stage of replication that requires remodeling for normal function