Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency

Linear Interaction Between Replication and Transcription Shapes DNA Break Dynamics at Recurrent DNA Break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unraveled the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA revealed co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity revealed that head-to-head interactions between replication and transcription machinery resulted in 60% DSB contribution to the head-on compared to 40% for co-directional. Our findings revealed TTR as a novel fragile class and highlighted how the linear interaction between transcription and replication impacts genome stability.

Nuclear-localized, iron-bound superoxide dismutase-2 antagonizes epithelial lineage programs to promote stemness of breast cancer cells via a histone demethylase activity

The dichotomous behavior of superoxide dismutase-2 (SOD2) in cancer biology has long been acknowledged and more recently linked to different posttranslational forms of the enzyme. However, a distinctive activity underlying its tumor-promoting function is yet to be described. Here, we report that acetylation, one of such posttranslational modifications (PTMs), increases SOD2 affinity for iron, effectively changing the biochemical function of this enzyme from that of an antioxidant to a demethylase. Acetylated, iron-bound SOD2 localizes to the nucleus, promoting stem cell gene expression via removal of suppressive epigenetic marks such as H3K9me3 and H3K927me3. Particularly, H3K9me3 was specifically removed from regulatory regions upstream of Nanog and Oct-4, two pluripotency factors involved in cancer stem cell reprogramming. Phenotypically, cells expressing nucleus-targeted SOD2 (NLS-SOD2) have increased clonogenicity and metastatic potential. FeSOD2 operating as H3 demethylase requires H2O2 as substrate, which unlike cofactors of canonical demethylases (i.e., oxygen and 2-oxoglutarate), is more abundant in tumor cells than in normal tissue. Therefore, our results indicate that FeSOD2 is a demethylase with unique activities and functions in the promotion of cancer evolution toward metastatic phenotypes.

The Role of the Histone Variant H2A.Z in Metazoan Development

During the emergence and radiation of complex multicellular eukaryotes from unicellular ancestors, transcriptional systems evolved by becoming more complex to provide the basis for this morphological diversity. The way eukaryotic genomes are packaged into a highly complex structure, known as chromatin, underpins this evolution of transcriptional regulation. Chromatin structure is controlled by a variety of different epigenetic mechanisms, including the major mechanism for altering the biochemical makeup of the nucleosome by replacing core histones with their variant forms. The histone H2A variant H2A.Z is particularly important in early metazoan development because, without it, embryos cease to develop and die. However, H2A.Z is also required for many differentiation steps beyond the stage that H2A.Z-knockout embryos die. H2A.Z can facilitate the activation and repression of genes that are important for pluripotency and differentiation, and acts through a variety of different molecular mechanisms that depend upon its modification status, its interaction with histone and nonhistone partners, and where it is deposited within the genome. In this review, we discuss the current knowledge about the different mechanisms by which H2A.Z regulates chromatin function at various developmental stages and the chromatin remodeling complexes that determine when and where H2A.Z is deposited.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

Evolution of replication origins in vertebrate genomes: rapid turnover despite selective constraints

The replication program of vertebrate genomes is driven by the chromosomal distribution and timing of activation of tens of thousands of replication origins. Genome-wide studies have shown the association of origins with promoters and CpG islands, and their enrichment in G-quadruplex motifs (G4). However, the genetic determinants driving their activity remain poorly understood. To gain insight on the constraints operating on origins, we conducted the first evolutionary comparison of origins across vertebrates. We generated a genome-wide map of chicken origins (the first of a bird genome), and performed a comparison with human and mouse maps. The analysis of intra-species polymorphism revealed a strong depletion of genetic diversity at the core of replication initiation loci. This depletion is not linked to the presence of G4 motifs, promoters or CpG islands. In contrast, we show that origins experienced a rapid turnover during vertebrate evolution, since pairwise comparisons of origin maps revealed that <24% of them are conserved among vertebrates. This study unravels the existence of a novel determinant of origins, the precise functional role of which remains to be determined. Despite the importance of replication initiation for the fitness of organisms, the distribution of origins along vertebrate chromosomes is highly flexible.

Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions

The spatiotemporal program of metazoan DNA replication is regulated during development and altered in cancers. We have generated novel OK-seq, Repli-seq and RNA-seq data to compare the DNA replication and gene expression programs of twelve cancer and non-cancer human cell types. Changes in replication fork directionality (RFD) determined by OK-seq are widespread but more frequent within GC-poor isochores and largely disconnected from transcription changes. Cancer cell RFD profiles cluster with non-cancer cells of similar developmental origin but not with different cancer types. Importantly, recurrent RFD changes are detected in specific tumour progression pathways. Using a model for establishment and early progression of chronic myeloid leukemia (CML), we identify 1027 replication initiation zones (IZs) that progressively change efficiency during long-term expression of the BCR-ABL1 oncogene, being twice more often downregulated than upregulated. Prolonged expression of BCR-ABL1 results in targeting of new IZs and accentuation of previous efficiency changes. Targeted IZs are predominantly located in GC-poor, late replicating gene deserts and frequently silenced in late CML. Prolonged expression of BCR-ABL1 results in massive deletion of GC-poor, late replicating DNA sequences enriched in origin silencing events. We conclude that BCR-ABL1 expression progressively affects replication and stability of GC-poor, late-replicating regions during CML progression.

Evidence for DNA Sequence Encoding of an Accessible Nucleosomal Array across Vertebrates

Nucleosome-depleted regions around which nucleosomes order following the "statistical" positioning scenario were recently shown to be encoded in the DNA sequence in human. This intrinsic nucleosomal ordering strongly correlates with oscillations in the local GC content as well as with the interspecies and intraspecies mutation profiles, revealing the existence of both positive and negative selection. In this letter, we show that these predicted nucleosome inhibitory energy barriers (NIEBs) with compacted neighboring nucleosomes are indeed ubiquitous to all vertebrates tested. These 1 kb-sized chromatin patterns are widely distributed along vertebrate chromosomes, overall covering more than a third of the genome. We have previously observed in human deviations from neutral evolution at these genome-wide distributed regions, which we interpreted as a possible indication of the selection of an open, accessible, and dynamic nucleosomal array to constitutively facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. As a first, very appealing observation supporting this hypothesis, we report evidence of a strong association between NIEB borders and the poly(A) tails of Alu sequences in human. These results suggest that NIEBs provide adequate chromatin patterns favorable to the integration of Alu retrotransposons and, more generally to various transposable elements in the genomes of primates and other vertebrates.

Crosstalk between pluripotency factors and higher-order chromatin organization

Pluripotent cells are characterized by a globally open and accessible chromatin organization that is thought to contribute to cellular plasticity and developmental decision-making. We recently identified the pluripotency factor Nanog as a key regulator of this form of chromatin architecture in mouse embryonic stem cells. In particular, we demonstrated that the transcription factors Nanog and Sall1 co-dependently mediate the epigenetic state of pericentromeric heterochromatin to reinforce a more open and accessible organization in pluripotent cells. Here, we summarize our main findings and place the work into a broader context. We explore how heterochromatin domains could be targets of transcriptional networks in pluripotent cells and are coordinated with cell state. We propose this integration may be to balance the requirement for a dynamic and plastic chromatin organization in pluripotent cells, together with priming for a more restrictive nuclear compartmentalization that is triggered rapidly upon lineage commitment.

Multi-scale structural community organisation of the human genome

BACKGROUND

Structural interaction frequency matrices between all genome loci are now experimentally achievable thanks to high-throughput chromosome conformation capture technologies. This ensues a new methodological challenge for computational biology which consists in objectively extracting from these data the structural motifs characteristic of genome organisation.

RESULTS

We deployed the fast multi-scale community mining algorithm based on spectral graph wavelets to characterise the networks of intra-chromosomal interactions in human cell lines. We observed that there exist structural domains of all sizes up to chromosome length and demonstrated that the set of structural communities forms a hierarchy of chromosome segments. Hence, at all scales, chromosome folding predominantly involves interactions between neighbouring sites rather than the formation of links between distant loci.

CONCLUSIONS

Multi-scale structural decomposition of human chromosomes provides an original framework to question structural organisation and its relationship to functional regulation across the scales. By construction the proposed methodology is independent of the precise assembly of the reference genome and is thus directly applicable to genomes whose assembly is not fully determined.

The Epigenetic Paradox of Pluripotent ES Cells

The propagation and maintenance of gene expression programs are at the foundation of the preservation of cell identity. A large and complex set of epigenetic mechanisms enables the long-term stability and inheritance of transcription states. A key property of authentic epigenetic regulation is being independent from the instructive signals used for its establishment. This makes epigenetic regulation, particularly epigenetic silencing, extremely robust and powerful to lock regulatory states and stabilise cell identity. In line with this, the establishment of epigenetic silencing during development restricts cell potency and maintains the cell fate choices made by transcription factors (TFs). However, how more immature cells that have not yet established their definitive fate maintain their transitory identity without compromising their responsiveness to signalling cues remains unclear. A paradigmatic example is provided by pluripotent embryonic stem (ES) cells derived from a transient population of cells of the blastocyst. Here, we argue that ES cells represent an interesting "epigenetic paradox": even though they are captured in a self-renewing state characterised by extremely efficient maintenance of their identity, which is a typical manifestation of robust epigenetic regulation, they seem not to heavily rely on classical epigenetic mechanisms. Indeed, self-renewal strictly depends on the TFs that previously instructed their undifferentiated identity and relies on a particular signalling-dependent chromatin state where repressive chromatin marks play minor roles. Although this "epigenetic paradox" may underlie their exquisite responsiveness to developmental cues, it suggests that alternative mechanisms to faithfully propagate gene regulatory states might be prevalent in ES cells.

Evidence of selection for an accessible nucleosomal array in human

BACKGROUND

Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.

RESULTS

We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.

CONCLUSION

The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.

Deciphering DNA replication dynamics in eukaryotic cell populations in relation with their averaged chromatin conformations

We propose a non-local model of DNA replication that takes into account the observed uncertainty on the position and time of replication initiation in eukaryote cell populations. By picturing replication initiation as a two-state system and considering all possible transition configurations, and by taking into account the chromatin’s fractal dimension, we derive an analytical expression for the rate of replication initiation. This model predicts with no free parameter the temporal profiles of initiation rate, replication fork density and fraction of replicated DNA, in quantitative agreement with corresponding experimental data from both S. cerevisiae and human cells and provides a quantitative estimate of initiation site redundancy. This study shows that, to a large extent, the program that regulates the dynamics of eukaryotic DNA replication is a collective phenomenon that emerges from the stochastic nature of replication origins initiation.

Metaboloepigenetic Regulation of Pluripotent Stem Cells

The differentiation of pluripotent stem cells is associated with extensive changes in metabolism, as well as widespread remodeling of the epigenetic landscape. Epigenetic regulation is essential for the modulation of differentiation, being responsible for cell type specific gene expression patterns through the modification of DNA and histones, thereby establishing cell identity. Each cell type has its own idiosyncratic pattern regarding the use of specific metabolic pathways. Rather than simply being perceived as a means of generating ATP and building blocks for cell growth and division, cellular metabolism can directly influence cellular regulation and the epigenome. Consequently, the significance of nutrients and metabolites as regulators of differentiation is central to understanding how cells interact with their immediate environment. This review serves to integrate studies on pluripotent stem cell metabolism, and the regulation of DNA methylation and acetylation and identifies areas in which current knowledge is limited.

Structural organization of human replication timing domains

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.