Chromatin state marks cell-type- and gender-specific replication of the Drosophila genome

Genomics tools for unraveling chromosome architecture

The spatial organization of chromosomes inside the cell nucleus is still poorly understood. This organization is guided by intra- and interchromosomal contacts and by interactions of specific chromosomal loci with relatively fixed nuclear 'landmarks' such as the nuclear envelope and the nucleolus. Researchers have begun to use new molecular genome-wide mapping techniques to uncover both types of molecular interactions, providing insights into the fundamental principles of interphase chromosome folding.

Systematic protein location mapping reveals five principal chromatin types in Drosophila cells

Chromatin is important for the regulation of transcription and other functions, yet the diversity of chromatin composition and the distribution along chromosomes are still poorly characterized. By integrative analysis of genome-wide binding maps of 53 broadly selected chromatin components in Drosophila cells, we show that the genome is segmented into five principal chromatin types that are defined by unique yet overlapping combinations of proteins and form domains that can extend over > 100 kb. We identify a repressive chromatin type that covers about half of the genome and lacks classic heterochromatin markers. Furthermore, transcriptionally active euchromatin consists of two types that differ in molecular organization and H3K36 methylation and regulate distinct classes of genes. Finally, we provide evidence that the different chromatin types help to target DNA-binding factors to specific genomic regions. These results provide a global view of chromatin diversity and domain organization in a metazoan cell.

Evaluating genome-scale approaches to eukaryotic DNA replication

Mechanisms regulating where and when eukaryotic DNA replication initiates remain a mystery. Recently, genome-scale methods have been brought to bear on this problem. The identification of replication origins and their associated proteins in yeasts is a well-integrated investigative tool, but corresponding data sets from multicellular organisms are scarce. By contrast, standardized protocols for evaluating replication timing have generated informative data sets for most eukaryotic systems. Here, I summarize the genome-scale methods that are most frequently used to analyse replication in eukaryotes, the kinds of questions each method can address and the technical hurdles that must be overcome to gain a complete understanding of the nature of eukaryotic replication origins.

E(nos)/CG4699 required for nanos function in the female germ line of Drosophila

The translational repressor Nanos is required in the germ line stem cells of the Drosophila ovary to maintain their capacity for self-renewal. Following division of the stem cells, Nanos is inhibited in the daughters that differentiate into cysts and ultimately become mature oocytes. The control of Nanos activity is thus an important aspect of the switch from self-renewal to differentiation. In this report, we describe a genetic interaction between nanos and Enhancer of nos, an allele of the previously uncharacterized locus CG4699. We find that E(nos) protein is required for normal accumulation of Nanos in the ovary and thus for maintenance of the germ line. The mechanism by which E(nos)/CG4699 protein acts is not clear, although it has been found in a complex with Mof acetylase. Consistent with the finding that E(nos) interacts with Mof, we observe that nanos and mof also interact genetically to maintain normal oogenesis.

Gene density profile reveals the marking of late replicated domains in the Drosophila melanogaster genome

Regulation of replication timing has been a focus of many studies. It has been shown that numerous chromosomal regions switch their replication timing on cell differentiation in Drosophila and mice. However, it is not clear which features of these regions are essential for such regulation. In this study, we examined the organization of late underreplicated regions (URs) of the Drosophila melanogaster genome. When compared with their flanks, these regions showed decreased gene density. A detailed view revealed that these regions originate from unusual combination of short genes and long intergenic spacers. Furthermore, gene expression study showed that this pattern is mostly contributed by short testis-specific genes abundant in the URs. Based on these observations, we developed a genome scanning algorithm and identified 110 regions possessing similar gene density and transcriptional profiles. According to the published data, replication of these regions has been significantly shifted towards late S-phase in two Drosophila cell lines and in polytene chromosomes. Our results suggest that genomic organization of the underreplicated areas of Drosophila polytene chromosomes may be associated with the regulation of their replication timing.

Eukaryotic chromosome DNA replication: where, when, and how?

DNA replication is central to cell proliferation. Studies in the past six decades since the proposal of a semiconservative mode of DNA replication have confirmed the high degree of conservation of the basic machinery of DNA replication from prokaryotes to eukaryotes. However, the need for replication of a substantially longer segment of DNA in coordination with various internal and external signals in eukaryotic cells has led to more complex and versatile regulatory strategies. The replication program in higher eukaryotes is under a dynamic and plastic regulation within a single cell, or within the cell population, or during development. We review here various regulatory mechanisms that control the replication program in eukaryotes and discuss future directions in this dynamic field.

Arabidopsis thaliana chromosome 4 replicates in two phases that correlate with chromatin state

DNA replication programs have been studied extensively in yeast and animal systems, where they have been shown to correlate with gene expression and certain epigenetic modifications. Despite the conservation of core DNA replication proteins, little is known about replication programs in plants. We used flow cytometry and tiling microarrays to profile DNA replication of Arabidopsis thaliana chromosome 4 (chr4) during early, mid, and late S phase. Replication profiles for early and mid S phase were similar and encompassed the majority of the euchromatin. Late S phase exhibited a distinctly different profile that includes the remaining euchromatin and essentially all of the heterochromatin. Termination zones were consistent between experiments, allowing us to define 163 putative replicons on chr4 that clustered into larger domains of predominately early or late replication. Early-replicating sequences, especially the initiation zones of early replicons, displayed a pattern of epigenetic modifications specifying an open chromatin conformation. Late replicons, and the termination zones of early replicons, showed an opposite pattern. Histone H3 acetylated on lysine 56 (H3K56ac) was enriched in early replicons, as well as the initiation zones of both early and late replicons. H3K56ac was also associated with expressed genes, but this effect was local whereas replication time correlated with H3K56ac over broad regions. The similarity of the replication profiles for early and mid S phase cells indicates that replication origin activation in euchromatin is stochastic. Replicon organization in Arabidopsis is strongly influenced by epigenetic modifications to histones and DNA. The domain organization of Arabidopsis is more similar to that in Drosophila than that in mammals, which may reflect genome size and complexity. The distinct patterns of association of H3K56ac with gene expression and early replication provide evidence that H3K56ac may be associated with initiation zones and replication origins.

During growth and development, all plants and animals must replicate their DNA. This process is regulated to ensure that all sequences are completely and accurately replicated and is limited to S phase of the cell cycle. In the cell, DNA is packaged with histone proteins into chromatin, and both DNA and histones are subject to epigenetic modifications that affect chromatin state. Euchromatin and heterochromatin are chromatin states marked by epigenetic modifications specifying open and closed conformations, respectively. Using the model plant Arabidopsis thaliana, we show that the time at which a DNA sequence replicates is influenced by the epigenetic modifications to the surrounding chromatin. DNA replication occurs in two phases, with euchromatin replicating in early and mid S phase and heterochromatin replicating late. DNA replication time has been linked to gene expression in other organisms, and this is also true in Arabidopsis because more genes are active in euchromatin when compared to heterochromatin. The earliest replicating DNA sequences are associated with acetylation of histone H3 on lysine 56 (H3K56ac). H3K56ac is also abundant in active genes, but the patterns of association of H3K56ac with gene expression and DNA replication are distinct, suggesting that H3K56ac is independently linked to both processes.

G2 phase chromatin lacks determinants of replication timing

Chromatin spatial organization helps establish the replication timing decision point at early G1. However, at G2, although retained, chromatin organization is no longer necessary or sufficient to maintain the replication timing program.

DNA replication in all eukaryotes follows a defined replication timing program, the molecular mechanism of which remains elusive. Using a Xenopus laevis egg extract replication system, we previously demonstrated that replication timing is established during early G1 phase of the cell cycle (timing decision point [TDP]), which is coincident with the repositioning and anchorage of chromatin in the newly formed nucleus. In this study, we use this same system to show that G2 phase chromatin lacks determinants of replication timing but maintains the overall spatial organization of chromatin domains, and we confirm this finding by genome-wide analysis of rereplication in vivo. In contrast, chromatin from quiescent cells retains replication timing but exhibits disrupted spatial organization. These data support a model in which events at the TDP, facilitated by chromatin spatial organization, establish determinants of replication timing that persist independent of spatial organization until the process of chromatin replication during S phase erases those determinants.

Domain-wide regulation of DNA replication timing during mammalian development

Studies of replication timing provide a handle into previously impenetrable higher-order levels of chromosome organization and their plasticity during development. Although mechanisms regulating replication timing are not clear, novel genome-wide studies provide a thorough survey of the extent to which replication timing is regulated during most of the early cell fate transitions in mammals, revealing coordinated changes of a defined set of 400-800 kb chromosomal segments that involve at least half the genome. Furthermore, changes in replication time are linked to changes in sub-nuclear organization and domain-wide transcriptional potential, and tissue-specific replication timing profiles are conserved from mouse to human, suggesting that the program has developmental significance. Hence, these studies have provided a solid foundation for linking megabase level chromosome structure to function, and suggest a central role for replication in domain-level genome organization.

Genome-wide analysis of the replication program in mammals

Microarray technology has facilitated the research of eukaryotic DNA replication on a genome-wide scale. Recent studies have shed light on the association between time of replication and chromosome structure, on the organization principles of the replication program, and on the correlation between replication timing and transcription. In this review, we summarize various genomic measurement approaches and the biological insights achieved through applying them in the study of the mammalian replication program.

Paucity and preferential suppression of transgenes in late replication domains of the D. melanogaster genome

Background

Eukaryotic genomes are organized in extended domains with distinct features intimately linking genome structure, replication pattern and chromatin state. Recently we identified a set of long late replicating euchromatic regions that are underreplicated in salivary gland polytene chromosomes of D. melanogaster.

Results

Here we demonstrate that these underreplicated regions (URs) have a low density of P-element and piggyBac insertions compared to the genome average or neighboring regions. In contrast, Minos-based transposons show no paucity in URs but have a strong bias to testis-specific genes. We estimated the suppression level in 2,852 stocks carrying a single P-element by analysis of eye color determined by the mini-white marker gene and demonstrate that the proportion of suppressed transgenes in URs is more than three times higher than in the flanking regions or the genomic average. The suppressed transgenes reside in intergenic, genic or promoter regions of the annotated genes. We speculate that the low insertion frequency of P-elements and piggyBacs in URs partially results from suppression of transgenes that potentially could prevent identification of transgenes due to complete suppression of the marker gene. In a similar manner, the proportion of suppressed transgenes is higher in loci replicating late or very late in Kc cells and these loci have a lower density of P-elements and piggyBac insertions. In transgenes with two marker genes suppression of mini-white gene in eye coincides with suppression of yellow gene in bristles.

Conclusions

Our results suggest that the late replication domains have a high inactivation potential apparently linked to the silenced or closed chromatin state in these regions, and that such inactivation potential is largely maintained in different tissues.

Heterochromatin protein 1 (HP1) modulates replication timing of the Drosophila genome

The replication of a chromosomal region during S phase can be highly dynamic between cell types that differ in transcriptome and epigenome. Early replication timing has been positively correlated with several histone modifications that occur at active genes, while repressive histone modifications mark late replicating regions. This raises the question if chromatin modulates the initiating events of replication. To gain insights into this question, we have studied the function of heterochromatin protein 1 (HP1), which is a reader of repressive methylation at histone H3 lysine 9, in genome-wide organization of replication. Cells with reduced levels of HP1 show an advanced replication timing of centromeric repeats in agreement with the model that repressive chromatin mediates the very late replication of large clusters of constitutive heterochromatin. Surprisingly, however, regions with high levels of interspersed repeats on the chromosomal arms, in particular on chromosome 4 and in pericentromeric regions of chromosome 2, behave differently. Here, loss of HP1 results in delayed replication. The fact that these regions are bound by HP1 suggests a direct effect. Thus while HP1 mediates very late replication of centromeric DNA, it is also required for early replication of euchromatic regions with high levels of repeats. This observation of opposing functions of HP1 suggests a model where HP1-mediated repeat inactivation or replication complex loading on the chromosome arms is required for proper activation of origins of replication that fire early. At the same time, HP1-mediated repression at constitutive heterochromatin is required to ensure replication of centromeric repeats at the end of S phase.

Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types

To identify evolutionarily conserved features of replication timing and their relationship to epigenetic properties, we profiled replication timing genome-wide in four human embryonic stem cell (hESC) lines, hESC-derived neural precursor cells (NPCs), lymphoblastoid cells, and two human induced pluripotent stem cell lines (hiPSCs), and compared them with related mouse cell types. Results confirm the conservation of coordinately replicated megabase-sized "replication domains" punctuated by origin-suppressed regions. Differentiation-induced replication timing changes in both species occur in 400- to 800-kb units and are similarly coordinated with transcription changes. A surprising degree of cell-type-specific conservation in replication timing was observed across regions of conserved synteny, despite considerable species variation in the alignment of replication timing to isochore GC/LINE-1 content. Notably, hESC replication timing profiles were significantly more aligned to mouse epiblast-derived stem cells (mEpiSCs) than to mouse ESCs. Comparison with epigenetic marks revealed a signature of chromatin modifications at the boundaries of early replicating domains and a remarkably strong link between replication timing and spatial proximity of chromatin as measured by Hi-C analysis. Thus, early and late initiation of replication occurs in spatially separate nuclear compartments, but rarely within the intervening chromatin. Moreover, cell-type-specific conservation of the replication program implies conserved developmental changes in spatial organization of chromatin. Together, our results reveal evolutionarily conserved aspects of developmentally regulated replication programs in mammals, demonstrate the power of replication profiling to distinguish closely related cell types, and strongly support the hypothesis that replication timing domains are spatially compartmentalized structural and functional units of three-dimensional chromosomal architecture.

Mutational bias shaping fly copy number variation: implications for genome evolution

Copy number variants (CNVs) underlie several genomic disorders and are a major source of genetic innovation. Consequently, any bias affecting their placement in the genome will impact our understanding of human disease and genome evolution. Here we report a mutational bias affecting CNVs that generates different probabilities of duplication and deletion across the genome in association with DNA replication time. We show that this mutational bias has important consequences for genome evolution by leading to different probabilities of gene duplication for different classes of genes and by linking the probability of gene duplication with the transcriptional activity of genes.

Interplay between DNA replication and gene expression: a harmonious coexistence

Multicellular organisms have evolved highly sophisticated machinery to that their genomes are accurately duplicated and that the various gene expression programs are established correctly. Recent large-scale studies have shed light on how these fundamental processes interact. Although the machinery mediating these processes share similar cis-regulatory elements, they are not strictly coregulated. Furthermore, studies of the replisome show that highly transcribed genes present a major obstacle to its operation. Further studies will be needed to identify key regulators of the spatio-temporal program of DNA replication, for the elucidation of the complex interplay between replication and transcription.

MRC1-dependent scaling of the budding yeast DNA replication timing program

We describe the DNA replication timing programs of 14 yeast mutants with an extended S phase identified by a novel genome-wide screen. These mutants are associated with the DNA replication machinery, cell-cycle control, and dNTP synthesis and affect different parts of S phase. In 13 of the mutants, origin activation time scales with the duration of S phase. A limited number of origins become inactive in these strains, with inactive origins characterized by small replicons and distributed throughout S phase. In sharp contrast, cells deleted of MRC1, a gene implicated in replication fork stabilization and in the replication checkpoint pathway, maintained wild-type firing times despite over twofold lengthening of S phase. Numerous dormant origins were activated in this mutant. Our data suggest that most perturbations that lengthen S phase affect the entire program of replication timing, rather than a specific subset of origins, maintaining the relative order of origin firing time and delaying firing with relative proportions. Mrc1 emerges as a regulator of this robustness of the replication program.

The MSL complex: X chromosome and beyond

X chromosomal regulation is a process that presents systematic problems of chromosome recognition and coordinated gene regulation. In Drosophila males, the ribonucleoprotein Male-Specific Lethal (MSL) complex plays an important role in hyperactivation of the X-linked genes to equalize gene dosage differences between the sexes. It appears that X chromosome recognition by the MSL complex may be mediated through a combination of sequence-specificity and transcriptional activities. The resulting transcriptional up-regulation also seems to involve several mechanisms, encompassing both gene-specific and chromosome-wide approaches. Interestingly the histone H4 lysine 16 specific MOF histone acetyl transferase, a key MSL member that hyper-acetylates the male X chromosome, is also involved in gene regulation beyond dosage compensation. A comparison of Drosophila and mammalian systems reveals intriguing parallels in MOF behavior, and highlights the multidisciplinary nature of this enzyme.

Drosophila ORC localizes to open chromatin and marks sites of cohesin complex loading

The origin recognition complex (ORC) is an essential DNA replication initiation factor conserved in all eukaryotes. In Saccharomyces cerevisiae, ORC binds to specific DNA elements; however, in higher eukaryotes, ORC exhibits little sequence specificity in vitro or in vivo. We investigated the genome-wide distribution of ORC in Drosophila and found that ORC localizes to specific chromosomal locations in the absence of any discernible simple motif. Although no clear sequence motif emerged, we were able to use machine learning approaches to accurately discriminate between ORC-associated sequences and ORC-free sequences based solely on primary sequence. The complex sequence features that define ORC binding sites are highly correlated with nucleosome positioning signals and likely represent a preferred nucleosomal landscape for ORC association. Open chromatin appears to be the underlying feature that is deterministic for ORC binding. ORC-associated sequences are enriched for the histone variant, H3.3, often at transcription start sites, and depleted for bulk nucleosomes. The density of ORC binding along the chromosome is reflected in the time at which a sequence replicates, with early replicating sequences having a high density of ORC binding. Finally, we found a high concordance between sites of ORC binding and cohesin loading, suggesting that, in addition to DNA replication, ORC may be required for the loading of cohesin on DNA in Drosophila.

Sequencing newly replicated DNA reveals widespread plasticity in human replication timing

Faithful transmission of genetic material to daughter cells involves a characteristic temporal order of DNA replication, which may play a significant role in the inheritance of epigenetic states. We developed a genome-scale approach--Repli Seq--to map temporally ordered replicating DNA using massively parallel sequencing and applied it to study regional variation in human DNA replication time across multiple human cell types. The method requires as few as 8,000 cytometry-fractionated cells for a single analysis, and provides high-resolution DNA replication patterns with respect to both cell-cycle time and genomic position. We find that different cell types exhibit characteristic replication signatures that reveal striking plasticity in regional replication time patterns covering at least 50% of the human genome. We also identified autosomal regions with marked biphasic replication timing that include known regions of monoallelic expression as well as many previously uncharacterized domains. Comparison with high-resolution genome-wide profiles of DNaseI sensitivity revealed that DNA replication typically initiates within foci of accessible chromatin comprising clustered DNaseI hypersensitive sites, and that replication time is better correlated with chromatin accessibility than with gene expression. The data collectively provide a unique, genome-wide picture of the epigenetic compartmentalization of the human genome and suggest that cell-lineage specification involves extensive reprogramming of replication timing patterns.

Genome-wide dynamics of replication timing revealed by in vitro models of mouse embryogenesis

Differentiation of mouse embryonic stem cells (mESCs) is accompanied by changes in replication timing. To explore the relationship between replication timing and cell fate transitions, we constructed genome-wide replication-timing profiles of 22 independent mouse cell lines representing 10 stages of early mouse development, and transcription profiles for seven of these stages. Replication profiles were cell-type specific, with 45% of the genome exhibiting significant changes at some point during development that were generally coordinated with changes in transcription. Comparison of early and late epiblast cell culture models revealed a set of early-to-late replication switches completed at a stage equivalent to the post-implantation epiblast, prior to germ layer specification and down-regulation of key pluripotency transcription factors [POU5F1 (also known as OCT4)/NANOG/SOX2] and coinciding with the emergence of compact chromatin near the nuclear periphery. These changes were maintained in all subsequent lineages (lineage-independent) and involved a group of irreversibly down-regulated genes, at least some of which were repositioned closer to the nuclear periphery. Importantly, many genomic regions of partially reprogrammed induced pluripotent stem cells (piPSCs) failed to re-establish ESC-specific replication-timing and transcription programs. These regions were enriched for lineage-independent early-to-late changes, which in female cells included the inactive X chromosome. Together, these results constitute a comprehensive "fate map" of replication-timing changes during early mouse development. Moreover, they support a model in which a distinct set of replication domains undergoes a form of "autosomal Lyonization" in the epiblast that is difficult to reprogram and coincides with an epigenetic commitment to differentiation prior to germ layer specification.

Long-range spreading of dosage compensation in Drosophila captures transcribed autosomal genes inserted on X

Dosage compensation in Drosophila melanogaster males is achieved via targeting of male-specific lethal (MSL) complex to X-linked genes. This is proposed to involve sequence-specific recognition of the X at approximately 150-300 chromatin entry sites, and subsequent spreading to active genes. Here we ask whether the spreading step requires transcription and is sequence-independent. We find that MSL complex binds, acetylates, and up-regulates autosomal genes inserted on X, but only if transcriptionally active. We conclude that a long-sought specific DNA sequence within X-linked genes is not obligatory for MSL binding. Instead, linkage and transcription play the pivotal roles in MSL targeting irrespective of gene origin and DNA sequence.

Localization and characteristics of DNA underreplication zone in the 75C region of intercalary heterochromatin in Drosophila melanogaster polytene chromosomes

In Drosophila polytene chromosomes, regions of intercalary heterochromatin are scattered throughout the euchromatic arms. Here, we present data on the first fine analysis of the individual intercalary heterochromatin region, 75C1-2, located in the 3L chromosome. By using electron microscopy, we demonstrated that this region appears as three closely adjacent condensed bands. Mapping of the region on the physical map by means of the chromosomal rearrangements with known breakpoints showed that the length of the region is about 445 kb. Although it seems that the SUUR protein binds to the whole 75C1-2 region, the proximal part of the region is fully polytenized, so the DNA underreplication zone is asymmetric and located in the distal half of the region. Finally, we speculate that intercalary heterochromatin regions of Drosophila polytene chromosomes are organized into three different types with respect to the localization of the underreplication zone.

Drosophila MSL complex globally acetylates H4K16 on the male X chromosome for dosage compensation

The Drosophila melanogaster male-specific lethal (MSL) complex binds the single male X chromosome to upregulate gene expression to equal that from the two female X chromosomes. However, it has been puzzling that approximately 25% of transcribed genes on the X chromosome do not stably recruit MSL complex. Here we find that almost all active genes on the X chromosome are associated with robust H4 Lys16 acetylation (H4K16ac), the histone modification catalyzed by the MSL complex. The distribution of H4K16ac is much broader than that of the MSL complex, and our results favor the idea that chromosome-wide H4K16ac reflects transient association of the MSL complex, occurring through spreading or chromosomal looping. Our results parallel those of localized Polycomb repressive complex and its more broadly distributed chromatin mark, trimethylated histone H3 Lys27 (H3K27me3), suggesting a common principle for the establishment of active and silenced chromatin domains.

Replication timing and transcriptional control: beyond cause and effect--part II

Replication timing is frequently discussed superficially in terms of its relationship to transcriptional activity via chromatin structure. However, so little is known about what regulates where and when replication initiates that it has been impossible to identify mechanistic and causal relationships. Moreover, much of our knowledge base has been anecdotal, derived from analyses of a few genes in unrelated cell lines. Recent studies have revisited long-standing hypotheses using genome-wide approaches. In particular, the foundation of this field was recently shored up with incontrovertible evidence that cellular differentiation is accompanied by coordinated changes in replication timing and transcription. These changes accompany subnuclear repositioning, and take place at the level of megabase-sized domains that transcend localized changes in chromatin structure or transcription. Inferring from these results, we propose that there exists a key transition during the middle of S-phase and that changes in replication timing traversing this period are associated with subnuclear repositioning and changes in the activity of certain classes of promoters.