Rapid Irreversible Transcriptional Reprogramming in Human Stem Cells Accompanied by Discordance between Replication Timing and Chromatin Compartment

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

Three-dimensional nuclear organisation and the DNA replication timing program

In eukaryotic cells, genome duplication is temporally organised according to a program referred to as the replication-timing (RT) program. The RT of individual genomic domains strikingly parallels the three-dimensional architecture of their chromatin contacts and subnuclear distribution. However, it is unclear whether this correspondence is coincidental or whether it indicates a causal and regulatory relationship. In either case, the nature of the molecular mechanisms ensuring this spatio-temporal coordination is still unknown. Here, we review recent evidence that begins to uncover the existence of a shared molecular machinery at the core of the spatio-temporal co-regulation of DNA replication and genome architecture. Finally, we discuss the outstanding, key question of the biological role of their coordination.

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

Decades of work on the spatiotemporal organization of mammalian DNA replication timing (RT) continues to unveil novel correlations with aspects of transcription and chromatin organization but, until recently, mechanisms regulating RT and the biological significance of the RT program had been indistinct. We now know that the RT program is both influenced by and necessary to maintain chromatin structure, forming an epigenetic positive feedback loop. Moreover, the discovery of specific cis-acting elements regulating mammalian RT at both the domain and the whole-chromosome level has revealed multiple cell-type-specific and developmentally regulated mechanisms of RT control. We review recent evidence for diverse mechanisms employed by different cell types to regulate their RT programs and the biological significance of RT regulation during development.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Genome architecture plasticity underlies DNA replication timing dynamics in cell differentiation

During the S-phase of eukaryotic cell cycle, DNA is replicated in a dedicatedly regulated temporal order, with regions containing active and inactive genes replicated early and late, respectively. Recent advances in sequencing technology allow us to explore the connection between replication timing (RT), histone modifications, and three-dimensional (3D) chromatin structure in diverse cell types. To characterize the dynamics during cell differentiation, corresponding sequencing data for human embryonic stem cells and four differentiated cell types were collected. By comparing RT and its extent of conservation before and after germ layer specification, the human genome was partitioned into distinct categories. Each category is then subject to comparisons on genomic, epigenetic, and chromatin 3D structural features. As expected, while constitutive early and late replication regions showed active and inactive features, respectively, dynamic regions with switched RT showed intermediate features. Surprisingly, although early-to-late replication and late-to-early replication regions showed similar histone modification patterns in hESCs, their structural preferences were opposite. Specifically, in hESCs, early-to-late replication regions tended to appear in the B compartment and large topologically associated domains, while late-to-early replication regions showed the opposite. Our results uncover the coordinated regulation of RT and 3D genome structure that underlies the loss of pluripotency and lineage commitment and indicate the importance and potential roles of genome architecture in biological processes.

Quantifying Chromosome Structural Reorganizations during Differentiation, Reprogramming, and Transdifferentiation

We developed a nonequilibrium model to study chromosome structural reorganizations within a simplified cell developmental system. From the chromosome structural perspective, we predicted that the neural progenitor cell is on the neural developmental path and very close to the transdifferentiation path from the fibroblast to the neuron cell. We identified an early bifurcation of stem cell differentiation processes and the cell-of-origin-specific reprogramming pathways. Our theoretical results are in good agreement with available experimental evidence, promoting future applications of our approach.

HiCRes: a computational method to estimate and predict the genomic resolution of Hi-C libraries

Three-dimensional (3D) conformation of the chromatin is crucial to stringently regulate gene expression patterns and DNA replication in a cell-type specific manner. Hi-C is a key technique for measuring 3D chromatin interactions genome wide. Estimating and predicting the resolution of a library is an essential step in any Hi-C experimental design. Here, we present the mathematical concepts to estimate the resolution of a dataset and predict whether deeper sequencing would enhance the resolution. We have developed HiCRes, a docker pipeline, by applying these concepts to several Hi-C libraries.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Nuclear organisation and replication timing are coupled through RIF1-PP1 interaction

Three-dimensional genome organisation and replication timing are known to be correlated, however, it remains unknown whether nuclear architecture overall plays an instructive role in the replication-timing programme and, if so, how. Here we demonstrate that RIF1 is a molecular hub that co-regulates both processes. Both nuclear organisation and replication timing depend upon the interaction between RIF1 and PP1. However, whereas nuclear architecture requires the full complement of RIF1 and its interaction with PP1, replication timing is not sensitive to RIF1 dosage. The role of RIF1 in replication timing also extends beyond its interaction with PP1. Availing of this separation-of-function approach, we have therefore identified in RIF1 dual function the molecular bases of the co-dependency of the replication-timing programme and nuclear architecture.

Exons and introns exhibit transcriptional strand asymmetry of dinucleotide distribution, damage formation and DNA repair

Recent cancer sequencing efforts have uncovered asymmetry in DNA damage induced mutagenesis between the transcribed and non-transcribed strands of genes. Here, we investigate the major type of damage induced by ultraviolet (UV) radiation, the cyclobutane pyrimidine dimers (CPDs), which are formed primarily in TT dinucleotides. We reveal that a transcriptional asymmetry already exists at the level of TT dinucleotide frequency and therefore also in CPD damage formation. This asymmetry is conserved in vertebrates and invertebrates and is completely reversed between introns and exons. We show the asymmetry in introns is linked to the transcription process itself, and is also found in enhancer elements. In contrast, the asymmetry in exons is not correlated to transcription, and is associated with codon usage preferences. Reanalysis of nucleotide excision repair, normalizing repair to the underlying TT frequencies, we show repair of CPDs is more efficient in exons compared to introns, contributing to the maintenance and integrity of coding regions. Our results highlight the importance of considering the primary sequence of the DNA in determining DNA damage sensitivity and mutagenic potential.

Formation of a multi-layered 3-dimensional structure of the heterochromatin compartment during early mammalian development

During embryogenesis in mammals, the 3-dimensional (3D) genome organization changes globally in parallel with transcription changes in a cell-type specific manner. This involves the progressive formation of heterochromatin, the best example of which is the inactive X chromosome (Xi) in females, originally discovered as a compact 3D structure at the nuclear periphery known as the Barr body. The heterochromatin formation on the autosomes and the Xi is tightly associated with the differentiation state and the developmental potential of cells, making it an ideal readout of the cellular epigenetic state. At a glance, the heterochromatin appears to be uniform. However, recent studies are beginning to reveal a more complex picture, with multiple hierarchical levels co-existing within the heterochromatin compartment. Such hierarchical levels appear to exist in the heterochromatin compartment on autosomes as well as on the Xi. Here, we review recent progress in our understanding of the 3D genome organization changes during the period of differentiation surrounding pluripotency in vivo and in vitro, with a focus on the heterochromatin compartment. We first look at the whole genome, then focus on the Xi, and discuss their differences and similarities. Finally, we present a unified view of how the heterochromatin compartment is formed and regulated during early development. In particular, we emphasize that there are multiple layers within the heterochromatic compartment on both the autosomes and the Xi, with regulatory mechanisms common and specific to each layer.

Distinctive epigenomic alterations in NF1-deficient cutaneous and plexiform neurofibromas drive differential MKK/p38 signaling

Benign peripheral nerve sheath tumors are the clinical hallmark of Neurofibromatosis Type 1. They account for substantial morbidity and mortality in NF1. Cutaneous (CNF) and plexiform neurofibromas (PNF) share nearly identical histology, but maintain different growth rates and risk of malignant conversion. The reasons for this disparate clinical behavior are not well explained by recent genome or transcriptome profiling studies. We hypothesized that CNFs and PNFs are epigenetically distinct tumor types that exhibit differential signaling due to genome-wide and site-specific methylation events. We interrogated the methylation profiles of 45 CNFs and 17 PNFs from NF1 subjects with the Illumina EPIC 850K methylation array. Based on these profiles, we confirm that CNFs and PNFs are epigenetically distinct tumors with broad differences in higher-order chromatin states and specific methylation events altering genes involved in key biological and cellular processes, such as inflammation, RAS/MAPK signaling, actin cytoskeleton rearrangement, and oxytocin signaling. Based on our identification of two separate DMRs associated with alternative leading exons in MAP2K3, we demonstrate differential RAS/MKK3/p38 signaling between CNFs and PNFs. Epigenetic reinforcement of RAS/MKK/p38 was a defining characteristic of CNFs leading to pro-inflammatory signaling and chromatin conformational changes, whereas PNFs signaled predominantly through RAS/MEK. Tumor size also correlated with specific CpG methylation events. Taken together, these findings confirm that NF1 deficiency influences the epigenetic regulation of RAS signaling fates, accounting for observed differences in CNF and PNF clinical behavior. The extension of these findings is that CNFs may respond differently than PNFs to RAS-targeted therapeutics raising the possibility of targeting p38-mediated inflammation for CNF treatment.

Computational Stem Cell Biology: Open Questions and Guiding Principles

Computational biology is enabling an explosive growth in our understanding of stem cells and our ability to use them for disease modeling, regenerative medicine, and drug discovery. We discuss four topics that exemplify applications of computation to stem cell biology: cell typing, lineage tracing, trajectory inference, and regulatory networks. We use these examples to articulate principles that have guided computational biology broadly and call for renewed attention to these principles as computation becomes increasingly important in stem cell biology. We also discuss important challenges for this field with the hope that it will inspire more to join this exciting area.

Distinct molecular profile and restricted stem cell potential defines the prospective human cranial neural crest from embryonic stem cell state

Neural crest cells are an embryonic multipotent stem cell population. Recent studies in model organisms have suggested that neural crest cells are specified earlier than previously thought, at blastula stages. However, the molecular dynamics of early neural crest specification, and functional changes from pluripotent precursors to early specified NC, remain to be elucidated. In this report, we utilized a robust human model of cranial neural crest formation to address the distinct molecular character of the earliest stages of neural crest specification and assess the functional differences from its embryonic stem cell precursor. Our human neural crest model reveals a rapid change in the epigenetic state of neural crest and pluripotency genes, accompanied by changes in gene expression upon Wnt-based induction from embryonic stem cells. These changes in gene expression are directly regulated by the transcriptional activity of β-catenin. Furthermore, prospective cranial neural crest cells are characterized by restricted stem cell potential compared to embryonic stem cells. Our results suggest that human neural crest induced by Wnt/β-catenin signaling from human embryonic stem cells rapidly acquire a prospective neural crest cell state defined by a unique molecular signature and endowed with limited potential compared to pluripotent stem cells.

Chromosome Conformation Capture and Beyond: Toward an Integrative View of Chromosome Structure and Function

Rapidly developing technologies have recently fueled an exciting era of discovery in the field of chromosome structure and nuclear organization. In addition to chromosome conformation capture (3C) methods, new alternative techniques have emerged to study genome architecture and biological processes in the nucleus, often in single or living cells. This sets an unprecedented stage for exploring the mechanisms that link chromosome structure and biological function. Here we review popular as well as emerging approaches to study chromosome organization, focusing on the contribution of complementary methodologies to our understanding of structures revealed by 3C methods and their biological implications, and discuss the next technical and conceptual frontiers.