Takayo Sasaki

An encyclopedia of mouse DNA elements (Mouse ENCODE)

To complement the human Encyclopedia of DNA Elements (ENCODE) project and to enable a broad range of mouse genomics efforts, the Mouse ENCODE Consortium is applying the same experimental pipelines developed for human ENCODE to annotate the mouse genome.

Replication Timing: A Fingerprint for Cell Identity and Pluripotency

Many types of epigenetic profiling have been used to classify stem cells, stages of cellular differentiation, and cancer subtypes. Existing methods focus on local chromatin features such as DNA methylation and histone modifications that require extensive analysis for genome-wide coverage. Replication timing has emerged as a highly stable cell type-specific epigenetic feature that is regulated at the megabase-level and is easily and comprehensively analyzed genome-wide. Here, we describe a cell classification method using 67 individual replication profiles from 34 mouse and human cell lines and stem cell-derived tissues, including new data for mesendoderm, definitive endoderm, mesoderm and smooth muscle. Using a Monte-Carlo approach for selecting features of replication profiles conserved in each cell type, we identify “replication timing fingerprints” unique to each cell type and apply a k nearest neighbor approach to predict known and unknown cell types. Our method correctly classifies 67/67 independent replication-timing profiles, including those derived from closely related intermediate stages. We also apply this method to derive fingerprints for pluripotency in human and mouse cells. Interestingly, the mouse pluripotency fingerprint overlaps almost completely with previously identified genomic segments that switch from early to late replication as pluripotency is lost. Thereafter, replication timing and transcription within these regions become difficult to reprogram back to pluripotency, suggesting these regions highlight an epigenetic barrier to reprogramming. In addition, the major histone cluster Hist1 consistently becomes later replicating in committed cell types, and several histone H1 genes in this cluster are downregulated during differentiation, suggesting a possible instrument for the chromatin compaction observed during differentiation. Finally, we demonstrate that unknown samples can be classified independently using site-specific PCR against fingerprint regions. In sum, replication fingerprints provide a comprehensive means for cell characterization and are a promising tool for identifying regions with cell type-specific organization.

Highly stable loading of Mcm proteins onto chromatin in living cells requires replication to unload

Components of the minichromosome maintenance complex (Mcm2-7) remain indefinitely bound to chromatin during G1 phase and replication arrest.

Pre-replication complex proteins assemble at regions of low nucleosome occupancy within the Chinese hamster dihydrofolate reductase initiation zone

Genome-scale mapping of pre-replication complex proteins has not been reported in mammalian cells. Poor enrichment of these proteins at specific sites may be due to dispersed binding, poor epitope availability or cell cycle stage-specific binding. Here, we have mapped sites of biotin-tagged ORC and MCM protein binding in G1-synchronized populations of Chinese hamster cells harboring amplified copies of the dihydrofolate reductase (DHFR) locus, using avidin-affinity purification of biotinylated chromatin followed by high-density microarray analysis across the DHFR locus. We have identified several sites of significant enrichment for both complexes distributed throughout the previously identified initiation zone. Analysis of the frequency of initiations across stretched DNA fibers from the DHFR locus confirmed a broad zone of de-localized initiation activity surrounding the sites of ORC and MCM enrichment. Mapping positions of mononucleosomal DNA empirically and computing nucleosome-positioning information in silico revealed that ORC and MCM map to regions of low measured and predicted nucleosome occupancy. Our results demonstrate that specific sites of ORC and MCM enrichment can be detected within a mammalian intitiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400-800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.

The Chinese Hamster Dihydrofolate Reductase Replication Origin Decision Point Follows Activation of Transcription and Suppresses Initiation of Replication within Transcription Units

ABSTRACT Chinese hamster ovary (CHO) cells select specific replication origin sites within the dihydrofolate reductase (DHFR) locus at a discrete point during G1 phase, the origin decision point (ODP). Origin selection is sensitive to transcription but not protein synthesis inhibitors, implicating a pretranslational role for transcription in origin specification. We have constructed a DNA array covering 121 kb surrounding the DHFR locus, to comprehensively investigate replication initiation and transcription in this region. When nuclei isolated within the first 3 h of G1 phase were stimulated to initiate replication in Xenopus egg extracts, replication initiated without any detectable preference for specific sites. At the ODP, initiation became suppressed from within the Msh3, DHFR, and 2BE2121 transcription units. Active transcription was mostly confined to these transcription units, and inhibition of transcription by alpha-amanitin resulted in the initiation of replication within transcription units, indicating that transcription is necessary to limit initiation events to the intergenic region. However, the resumption of DHFR transcription after mitosis took place prior to the ODP and so is not on its own sufficient to suppress initiation of replication. Together, these results demonstrate a remarkable flexibility in sequence selection for initiating replication and implicate transcription as one important component of origin specification at the ODP.

The Distribution of Genomic Variations in Human iPSCs Is Related to Replication-Timing Reorganization during Reprogramming

Cell-fate change involves significant genome reorganization, including changes in replication timing, but how these changes are related to genetic variation has not been examined. To study how a change in replication timing that occurs during reprogramming impacts the copy-number variation (CNV) landscape, we generated genome-wide replication-timing profiles of induced pluripotent stem cells (iPSCs) and their parental fibroblasts. A significant portion of the genome changes replication timing as a result of reprogramming, indicative of overall genome reorganization. We found that early- and late-replicating domains in iPSCs are differentially affected by copy-number gains and losses and that in particular, CNV gains accumulate in regions of the genome that change to earlier replication during the reprogramming process. This differential relationship was present irrespective of reprogramming method. Overall, our findings reveal a functional association between reorganization of replication timing and the CNV landscape that emerges during reprogramming.

Murine esBAF chromatin remodeling complex subunits BAF250a and Brg1 are necessary to maintain and reprogram pluripotency-specific replication timing of select replication domains

BackgroundCellular differentiation and reprogramming are accompanied by changes in replication timing and 3D organization of large-scale (400 to 800 Kb) chromosomal domains (‘replication domains’), but few gene products have been identified whose disruption affects these properties.ResultsHere we show that deletion of esBAF chromatin-remodeling complex components BAF250a and Brg1, but not BAF53a, disrupts replication timing at specific replication domains. Also, BAF250a-deficient fibroblasts reprogrammed to a pluripotency-like state failed to reprogram replication timing in many of these same domains. About half of the replication domains affected by Brg1 loss were also affected by BAF250a loss, but a much larger set of domains was affected by BAF250a loss. esBAF binding in the affected replication domains was dependent upon BAF250a but, most affected domains did not contain genes whose transcription was affected by loss of esBAF.ConclusionsLoss of specific esBAF complex subunits alters replication timing of select replication domains in pluripotent cells.

Discovery of the tail tube gene of bacteriophage Mu and sequence analysis of the sheath and tube genes

The nucleotide sequence was determined for 2.75 kbp of phage Mu DNA encoding the contractile tail sheath protein L. N-terminal sequence analysis of Mu tail tube and sheath proteins identified the open reading frame just downstream of gene L as the tube gene. This clustering and order of the sheath and tube genes appear to be common among the myoviridae. Database homology searches revealed high similarity between the Mu sheath and tube proteins and two proteins in a Haemophilus influenzae Mu-like prophage, suggesting that they are the sheath and tube proteins of that prophage.

Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq

This protocol is an extension to: Nat. Protoc. 6, 870–895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.