Chromosomal Mcm2-7 distribution and the genome replication program in species from yeast to humans.

The spatio-temporal program of genome replication across eukaryotes is thought to be driven both by the uneven loading of pre-replication complexes (pre-RCs) across the genome at the onset of S-phase, and by differences in the timing of activation of these complexes during S phase. To determine the degree to which distribution of pre-RC loading alone could account for chromosomal replication patterns, we mapped the binding sites of the Mcm2-7 helicase complex (MCM) in budding yeast, fission yeast, mouse and humans. We observed similar individual MCM double-hexamer (DH) footprints across the species, but notable differences in their distribution: Footprints in budding yeast were more sharply focused compared to the other three organisms, consistent with the relative sequence specificity of replication origins in S. cerevisiae. Nonetheless, with some clear exceptions, most notably the inactive X-chromosome, much of the fluctuation in replication timing along the chromosomes in all four organisms reflected uneven chromosomal distribution of pre-replication complexes.

Sir2 suppresses transcription-mediated displacement of Mcm2-7 replicative helicases at the ribosomal DNA repeats.

Repetitive DNA sequences within eukaryotic heterochromatin are poorly transcribed and replicate late in S-phase. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). Despite the widespread association between transcription and replication timing, it remains unclear how transcription might impinge on replication, or vice versa. Here we show that, when silencing of an RNA polymerase II (RNA Pol II)-transcribed non-coding RNA at the rDNA is disrupted by SIR2 deletion, RNA polymerase pushes and thereby relocalizes replicative Mcm2-7 helicases away from their loading sites to an adjacent region with low nucleosome occupancy, and this relocalization is associated with increased rDNA origin efficiency. Our results suggest a model in which two of the major defining features of heterochromatin, transcriptional silencing and late replication, are mechanistically linked through suppression of polymerase-mediated displacement of replication initiation complexes.

SIR2 suppresses replication gaps and genome instability by balancing replication between repetitive and unique sequences.

Replication gaps that persist into mitosis likely represent important threats to genome stability, but experimental identification of these gaps has proved challenging. We have developed a technique that allows us to explore the dynamics by which genome replication is completed before mitosis. Using this approach, we demonstrate that excessive allocation of replication resources to origins within repetitive regions, induced by SIR2 deletion, leads to persistent replication gaps and genome instability. Conversely, the weakening of replication origins in repetitive regions suppresses these gaps. Given known age- and cancer-associated changes in chromatin accessibility at repetitive sequences, we suggest that replication gaps resulting from misallocation of replication resources underlie age- and disease-associated genome instability.

Screen for reactivation of MeCP2 on the inactive X chromosome identifies the BMP/TGF-ß superfamily as a regulator of XIST expression.

Rett syndrome (RS) is a debilitating neurological disorder affecting mostly girls with heterozygous mutations in the gene encoding the methyl-CpG-binding protein MeCP2 on the X chromosome. Because restoration of MeCP2 expression in a mouse model reverses neurologic deficits in adult animals, reactivation of the wild-type copy of MeCP2 on the inactive X chromosome (Xi) presents a therapeutic opportunity in RS. To identify genes involved in MeCP2 silencing, we screened a library of 60,000 shRNAs using a cell line with a MeCP2 reporter on the Xi and found 30 genes clustered in seven functional groups. More than half encoded proteins with known enzymatic activity, and six were members of the bone morphogenetic protein (BMP)/TGF-ß pathway. shRNAs directed against each of these six genes down-regulated X-inactive specific transcript (XIST), a key player in X-chromosome inactivation that encodes an RNA that coats the silent X chromosome, and modulation of regulators of this pathway both in cell culture and in mice demonstrated robust regulation of XIST. Moreover, we show that Rnf12, an X-encoded ubiquitin ligase important for initiation of X-chromosome inactivation and XIST transcription in ES cells, also plays a role in maintenance of the inactive state through regulation of BMP/TGF-ß signaling. Our results identify pharmacologically suitable targets for reactivation of MeCP2 on the Xi and a genetic circuitry that maintains XIST expression and X-chromosome inactivation in differentiated cells.

A natural polymorphism in rDNA replication origins links origin activation with calorie restriction and lifespan.

Aging and longevity are complex traits influenced by genetic and environmental factors. To identify quantitative trait loci (QTLs) that control replicative lifespan, we employed an outbred Saccharomyces cerevisiae model, generated by crossing a vineyard and a laboratory strain. The predominant QTL mapped to the rDNA, with the vineyard rDNA conferring a lifespan increase of 41%. The lifespan extension was independent of Sir2 and Fob1, but depended on a polymorphism in the rDNA origin of replication from the vineyard strain that reduced origin activation relative to the laboratory origin. Strains carrying vineyard rDNA origins have increased capacity for replication initiation at weak plasmid and genomic origins, suggesting that inability to complete genome replication presents a major impediment to replicative lifespan. Calorie restriction, a conserved mediator of lifespan extension that is also independent of Sir2 and Fob1, reduces rDNA origin firing in both laboratory and vineyard rDNA. Our results are consistent with the possibility that calorie restriction, similarly to the vineyard rDNA polymorphism, modulates replicative lifespan through control of rDNA origin activation, which in turn affects genome replication dynamics.

Genetic variation shapes protein networks mainly through non-transcriptional mechanisms.

Networks of co-regulated transcripts in genetically diverse populations have been studied extensively, but little is known about the degree to which these networks cause similar co-variation at the protein level. We quantified 354 proteins in a genetically diverse population of yeast segregants, which allowed for the first time construction of a coherent protein co-variation matrix. We identified tightly co-regulated groups of 36 and 93 proteins that were made up predominantly of genes involved in ribosome biogenesis and amino acid metabolism, respectively. Even though the ribosomal genes were tightly co-regulated at both the protein and transcript levels, genetic regulation of proteins was entirely distinct from that of transcripts, and almost no genes in this network showed a significant correlation between protein and transcript levels. This result calls into question the widely held belief that in yeast, as opposed to higher eukaryotes, ribosomal protein levels are regulated primarily by regulating transcript levels. Furthermore, although genetic regulation of the amino acid network was more similar for proteins and transcripts, regression analysis demonstrated that even here, proteins vary predominantly as a result of non-transcriptional variation. We also found that cis regulation, which is common in the transcriptome, is rare at the level of the proteome. We conclude that most inter-individual variation in levels of these particular high abundance proteins in this genetically diverse population is not caused by variation of their underlying transcripts.

Identification of 1600 replication origins in S. cerevisiae.

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

Sir2 and Fun30 regulate ribosomal DNA replication timing via Mcm helicase positioning and nucleosome occupancy.

The association between late replication timing and low transcription rates in eukaryotic heterochromatin is well-known, yet the specific mechanisms underlying this link remain uncertain. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). We have previously reported that in the absence of SIR2, a derepressed RNA PolII repositions MCM replicative helicases from their loading site at the ribosomal origin, where they abut well-positioned, high-occupancy nucleosomes, to an adjacent region with lower nucleosome occupancy. By developing a method that can distinguish activation of closely spaced MCM complexes, here we show that the displaced MCMs at rDNA origins have increased firing propensity compared to the non-displaced MCMs. Furthermore, we found that both, activation of the repositioned MCMs and low occupancy of the adjacent nucleosomes critically depend on the chromatin remodeling activity of FUN30. Our study elucidates the mechanism by which Sir2 delays replication timing, and it demonstrates, for the first time, that activation of a specific replication origin in vivo relies on the nucleosome context shaped by a single chromatin remodeler.

Identification of 1600 replication origins in S. cerevisiae.

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, "a site to which MCM is bound in G1" might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, a technique referred to as "Chromatin Endogenous Cleavage", we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least 3 orders of magnitude. Published data quantifying the production of ssDNA during S phase showed clear evidence of replication initiation among the most abundant 1600 of these sites, with replication activity decreasing in concert with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5,500 sites. Specifically, these sites (1) appeared in intergenic nucleosome-free regions that were flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. Furthermore, the high resolution of this technique allowed us to demonstrate a strong bias for detecting Mcm double-hexamers downstream rather than upstream of the ACS, which is consistent with the directionality of Mcm loading by Orc that has been observed in vitro. We conclude that, if sites at which Mcm double-hexamers are loaded can function as replication origins, then DNA replication origins are at least 3-fold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than is widely assumed.

Proteomic classification of acute leukemias by alignment-based quantitation of LC-MS/MS data sets.

Despite immense interest in the proteome as a source of biomarkers in cancer, mass spectrometry has yet to yield a clinically useful protein biomarker for tumor classification. To explore the potential of a particular class of mass spectrometry-based quantitation approaches, label-free alignment of liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) data sets, for the identification of biomarkers for acute leukemias, we asked whether a label-free alignment algorithm could distinguish known classes of leukemias on the basis of their proteomes. This approach to quantitation involves (1) computational alignment of MS1 peptide peaks across large numbers of samples; (2) measurement of the relative abundance of peptides across samples by integrating the area under the curve of the MS1 peaks; and (3) assignment of peptide IDs to those quantified peptide peaks on the basis of the corresponding MS2 spectra. We extracted proteins from blasts derived from four patients with acute myeloid leukemia (AML, acute leukemia of myeloid lineage) and five patients with acute lymphoid leukemia (ALL, acute leukemia of lymphoid lineage). Mobilized CD34+ cells purified from peripheral blood of six healthy donors and mononuclear cells (MNC) from the peripheral blood of two healthy donors were used as healthy controls. Proteins were analyzed by LC-MS/MS and quantified with a label-free alignment-based algorithm developed in our laboratory. Unsupervised hierarchical clustering of blinded samples separated the samples according to their known biological characteristics, with each sample group forming a discrete cluster. The four proteins best able to distinguish CD34+, AML, and ALL were all either known biomarkers or proteins whose biological functions are consistent with their ability to distinguish these classes. We conclude that alignment-based label-free quantitation of LC-MS/MS data sets can, at least in some cases, robustly distinguish known classes of leukemias, thus opening the possibility that large scale studies using such algorithms can lead to the identification of clinically useful biomarkers.

G2-seq: A High Throughput Sequencing-based Technique for Identifying Late Replicating Regions of the Genome.

Numerous techniques have been developed to follow the progress of DNA replication through the S phase of the cell cycle. Most of these techniques have been directed toward elucidation of the location and timing of initiation of genome duplication rather than its completion. However, it is critical that we understand regions of the genome that are last to complete replication, because these regions suffer elevated levels of chromosomal breakage and mutation, and they have been associated with both disease and aging. Here we describe how we have extended a technique that has been used to monitor replication initiation to instead identify those regions of the genome last to complete replication. This approach is based on a combination of flow cytometry and high throughput sequencing. Although this report focuses on the application of this technique to yeast, the approach can be used with any cells that can be sorted in a flow cytometer according to DNA content.

Perturbed maintenance of transcriptional repression on the inactive X-chromosome in the mouse brain after Xist deletion.

BACKGROUND: The long noncoding RNA Xist is critical for initiation and establishment of X-chromosome inactivation during embryogenesis in mammals, but it is unclear whether its continued expression is required for maintaining X-inactivation in vivo. RESULTS: By using an inactive X-chromosome-linked MeCP2-GFP reporter, which allowed us to enumerate reactivation events in the mouse brain even when they occur in very few cells, we found that deletion of Xist in the brain after establishment of X-chromosome inactivation leads to reactivation in 2-5% of neurons and in a smaller fraction of astrocytes. In contrast to global loss of both H3 lysine 27 trimethylation (H3K27m3) and histone H2A lysine 119 monoubiquitylation (H2AK119ub1) we observed upon Xist deletion, alterations in CpG methylation were subtle, and this was mirrored by only minor alterations in X-chromosome-wide gene expression levels, with highly expressed genes more prone to both derepression and demethylation compared to genes with low expression level. CONCLUSION: Our results demonstrate that Xist plays a role in the maintenance of histone repressive marks, DNA methylation and transcriptional repression on the inactive X-chromosome, but that partial loss of X-dosage compensation in the absence of Xist in the brain is well tolerated.

Pooled shRNA Screen for Reactivation of MeCP2 on the Inactive X Chromosome.

Forward genetic screens using reporter genes inserted into the heterochromatin have been extensively used to investigate mechanisms of epigenetic control in model organisms. Technologies including short hairpin RNAs (shRNAs) and clustered regularly interspaced short palindromic repeats (CRISPR) have enabled such screens in diploid mammalian cells. Here we describe a large-scale shRNA screen for regulators of X-chromosome inactivation (XCI), using a murine cell line with firefly luciferase and hygromycin resistance genes knocked in at the C-terminus of the methyl CpG binding protein 2 (MeCP2) gene on the inactive X-chromosome (Xi). Reactivation of the construct in the reporter cell line conferred survival advantage under hygromycin B selection, enabling us to screen a large shRNA library and identify hairpins that reactivated the reporter by measuring their post-selection enrichment using next-generation sequencing. The enriched hairpins were then individually validated by testing their ability to activate the luciferase reporter on Xi.

Genetic basis of proteome variation in yeast.

Proper regulation of protein levels is essential for health, and abnormal levels of proteins are hallmarks of many diseases. A number of studies have recently shown that messenger RNA levels vary among individuals of a species and that genetic linkage analysis can be used to identify quantitative trait loci that influence these levels. By contrast, little is known about the genetic basis of variation in protein levels in genetically diverse populations, in large part because techniques for large-scale measurements of protein abundance lag far behind those for measuring transcript abundance. Here we describe a label-free, mass spectrometry-based approach to measuring protein levels in total unfractionated cellular proteins, and we apply this approach to elucidate the genetic basis of variation in protein abundance in a cross between two diverse strains of yeast. Loci that influenced protein abundance differed from those that influenced transcript levels, emphasizing the importance of direct analysis of the proteome.