Juan Carlos Rivera-Mulia

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.

Copy Number Variation Is a Fundamental Aspect of the Placental Genome

Discovery of lineage-specific somatic copy number variation (CNV) in mammals has led to debate over whether CNVs are mutations that propagate disease or whether they are a normal, and even essential, aspect of cell biology. We show that 1,000N polyploid trophoblast giant cells (TGCs) of the mouse placenta contain 47 regions, totaling 138 Megabases, where genomic copies are underrepresented (UR). UR domains originate from a subset of late-replicating heterochromatic regions containing gene deserts and genes involved in cell adhesion and neurogenesis. While lineage-specific CNVs have been identified in mammalian cells, classically in the immune system where V(D)J recombination occurs, we demonstrate that CNVs form during gestation in the placenta by an underreplication mechanism, not by recombination nor deletion. Our results reveal that large scale CNVs are a normal feature of the mammalian placental genome, which are regulated systematically during embryogenesis and are propagated by a mechanism of underreplication.

Replication timing and transcriptional control: beyond cause and effect-part III.

DNA replication is essential for faithful transmission of genetic information and is intimately tied to chromosome structure and function. Genome duplication occurs in a defined temporal order known as the replication-timing (RT) program, which is regulated during the cell cycle and development in discrete units referred to as replication domains (RDs). RDs correspond to topologically-associating domains (TADs) and are spatio-temporally compartmentalized in the nucleus. While improvements in experimental tools have begun to reveal glimpses of causality, they have also unveiled complex context-dependent relationships that challenge long recognized correlations of RT to chromatin organization and gene regulation. In particular, RDs/TADs that switch RT during development march to the beat of a different drummer.

Replicating Large Genomes: Divide and Conquer

Complete duplication of large metazoan chromosomes requires thousands of potential initiation sites, only a small fraction of which are selected in each cell cycle. Assembly of the replication machinery is highly conserved and tightly regulated during the cell cycle, but the sites of initiation are highly flexible, and their temporal order of firing is regulated at the level of large-scale multi-replicon domains. Importantly, the number of replication forks must be quickly adjusted in response to replication stress to prevent genome instability. Here we argue that large genomes are divided into domains for exactly this reason. Once established, domain structure abrogates the need for precise initiation sites and creates a scaffold for the evolution of other chromosome functions.

Large-Scale Chromatin Structure–Function Relationships during the Cell Cycle and Development: Insights from Replication Timing

Chromosome architecture has received a lot of attention since the recent development of genome-scale methods to measure chromatin interactions (Hi-C), enabling the first sequence-based models of chromosome tertiary structure. A view has emerged of chromosomes as a string of structural units (topologically associating domains; TADs) whose boundaries persist through the cell cycle and development. TADs with similar chromatin states tend to aggregate, forming spatially segregated chromatin compartments. However, high-resolution Hi-C has revealed substructure within TADs (subTADs) that poses a challenge for models that attribute significance to structural units at any given scale. More than 20 years ago, the DNA replication field independently identified stable structural (and functional) units of chromosomes (replication foci) as well as spatially segregated chromatin compartments (early and late foci), but lacked the means to link these units to genomic map units. Genome-wide studies of replication timing (RT) have now merged these two disciplines by identifying individual units of replication regulation (replication domains; RDs) that correspond to TADs and are arranged in 3D to form spatiotemporally segregated subnuclear compartments. Furthermore, classifying RDs/TADs by their constitutive versus developmentally regulated RT has revealed distinct classes of chromatin organization, providing unexpected insight into the relationship between large-scale chromosome structure and function.

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.

Allele-specific control of replication timing and genome organization during development

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program. RT is regulated during development in discrete chromosomal units, coordinated with transcriptional activity and 3D genome organization. Here, we derived distinct cell types from F1 hybrid musculus × castaneus mouse crosses and exploited the high single-nucleotide polymorphism (SNP) density to characterize allelic differences in RT (Repli-seq), genome organization (Hi-C and promoter-capture Hi-C), gene expression (total nuclear RNA-seq), and chromatin accessibility (ATAC-seq). We also present HARP, a new computational tool for sorting SNPs in phased genomes to efficiently measure allele-specific genome-wide data. Analysis of six different hybrid mESC clones with different genomes (C57BL/6, 129/sv, and CAST/Ei), parental configurations, and gender revealed significant RT asynchrony between alleles across ∼12% of the autosomal genome linked to subspecies genomes but not to parental origin, growth conditions, or gender. RT asynchrony in mESCs strongly correlated with changes in Hi-C compartments between alleles but not as strongly with SNP density, gene expression, imprinting, or chromatin accessibility. We then tracked mESC RT asynchronous regions during development by analyzing differentiated cell types, including extraembryonic endoderm stem (XEN) cells, four male and female primary mouse embryonic fibroblasts (MEFs), and neural precursor cells (NPCs) differentiated in vitro from mESCs with opposite parental configurations. We found that RT asynchrony and allelic discordance in Hi-C compartments seen in mESCs were largely lost in all differentiated cell types, accompanied by novel sites of allelic asynchrony at a considerably smaller proportion of the genome, suggesting that genome organization of homologs converges to similar folding patterns during cell fate commitment.

DNA replication timing alterations identify common markers between distinct progeroid diseases

Significance We show that the temporal order of replication (replication timing, RT), normally an extremely stable cell type-specific chromosomal property, is altered in cells from two different premature aging (progeroid) diseases. By converting patient cells to stem cells and redifferentiating them as a model of disease progression, we identified the TP63 gene as one of the earliest RT alterations and altered RT was associated with abnormal TP63 gene expression. TP63 mutations have been linked to other diseases that share clinical features of progeroid syndromes. These findings introduce an approach for disease marker discovery, identify molecular abnormalities distinguishing progeroid diseases from natural aging, and point to TP63 as a molecular link to the pathophysiological manifestations of progeroid diseases. Progeroid syndromes are rare genetic disorders that phenotypically resemble natural aging. Different causal mutations have been identified, but no molecular alterations have been identified that are in common to these diseases. DNA replication timing (RT) is a robust cell type-specific epigenetic feature highly conserved in the same cell types from different individuals but altered in disease. Here, we characterized DNA RT program alterations in Hutchinson–Gilford progeria syndrome (HGPS) and Rothmund–Thomson syndrome (RTS) patients compared with natural aging and cellular senescence. Our results identified a progeroid-specific RT signature that is common to cells from three HGPS and three RTS patients and distinguishes them from healthy individuals across a wide range of ages. Among the RT abnormalities, we identified the tumor protein p63 gene (TP63) as a gene marker for progeroid syndromes. By using the redifferentiation of four patient-derived induced pluripotent stem cells as a model for the onset of progeroid syndromes, we tracked the progression of RT abnormalities during development, revealing altered RT of the TP63 gene as an early event in disease progression of both HGPS and RTS. Moreover, the RT abnormalities in progeroid patients were associated with altered isoform expression of TP63. Our findings demonstrate the value of RT studies to identify biomarkers not detected by other methods, reveal abnormal TP63 RT as an early event in progeroid disease progression, and suggest TP63 gene regulation as a potential therapeutic target.

Identification of cis elements for spatio-temporal control of DNA replication

The temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, but cis-elements regulating spatio-temporal control of replication have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discrete cis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function. Highlights cis-elements (ERCEs) regulate large scale chromosome structure and function Multiple ERCEs cooperatively control domain-wide replication ERCEs harbor prominent active chromatin features and form CTCF-independent loops ERCEs enable genetic dissection of large-scale chromosome structure-function.

Integrative detection and analysis of structural variation in cancer genomes

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.The authors present an integrative framework for identifying structural variants (SVs) in cancer that applies optical mapping, Hi-C, and whole-genome sequencing. They find SVs affecting distal regulatory sequences, DNA replication, and three-dimensional chromatin structure.