Search

Features of higher-order chromatin organization-such as A/B compartments, topologically associating domains and chromatin loops-are temporarily disrupted during mitosis1,2. Because these structures are thought to influence gene regulation, it is important to understand how they are re-established after mitosis. Here we examine the dynamics of chromosome reorganization by Hi-C after mitosis in highly purified, synchronous mouse erythroid cell populations. We observed rapid establishment of A/B compartments, followed by their gradual intensification and expansion. Contact domains form from the 'bottom up'-smaller subTADs are formed initially, followed by convergence into multi-domain TAD structures. CTCF is partially retained on mitotic chromosomes and immediately resumes full binding in ana/telophase. By contrast, cohesin is completely evicted from mitotic chromosomes and regains focal binding at a slower rate. The formation of CTCF/cohesin co-anchored structural loops follows the kinetics of cohesin positioning. Stripe-shaped contact patterns-anchored by CTCF-grow in length, which is consistent with a loop-extrusion process after mitosis. Interactions between cis-regulatory elements can form rapidly, with rates exceeding those of CTCF/cohesin-anchored contacts. Notably, we identified a group of rapidly emerging transient contacts between cis-regulatory elements in ana/telophase that are dissolved upon G1 entry, co-incident with the establishment of inner boundaries or nearby interfering chromatin loops. We also describe the relationship between transcription reactivation and architectural features. Our findings indicate that distinct but mutually influential forces drive post-mitotic chromatin reconfiguration.

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability across dividing cells in metazoans. It is currently unknown how the location of replication origins and the timing of their activation is determined in the human genome. Here, we dissect the role for G1 phase topologically associating domains (TADs), subTADs, and loops in the activation of replication initiation zones (IZs). We identify twelve subtypes of self-interacting chromatin domains distinguished by their degree of nesting, the presence of corner dot structures indicative of loops, and their co-localization with A/B compartments. Early replicating IZs localize to boundaries of nested corner-dot TAD/subTADs anchored by high density arrays of co-occupied CTCF+cohesin binding sites with divergently oriented motifs. By contrast, late replicating IZs localize to weak TADs/subTAD boundaries devoid of corner dots and most often anchored by singlet CTCF+cohesin sites. Upon global knock-down of cohesin-mediated loops in G1, early wave focal IZs replicate later in S phase and convert to diffuse placement along the genome. Moreover, IZs in mid-late S phase are delayed to the final minutes before entry into G2 when cohesin-mediated dot-less boundaries are ablated. We also delete a specific loop anchor and observe a sharp local delay of an early wave IZ to replication in late S phase. Our data demonstrate that cohesin-mediated loops at genetically-encoded TAD/subTAD boundaries in G1 phase are an essential determinant of the precise genomic placement of human replication origins in S phase.

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3–6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

Mammalian genomes are folded in a hierarchy of compartments, topologically associating domains (TADs), subTADs and looping interactions. Here, we describe 3DNetMod, a graph theory-based method for sensitive and accurate detection of chromatin domains across length scales in Hi-C data. We identify nested, partially overlapping TADs and subTADs genome wide by optimizing network modularity and varying a single resolution parameter. 3DNetMod can be applied broadly to understand genome reconfiguration in development and disease.

Mammalian genomes are folded into a hierarchy of compartments, topologically associating domains (TADs), subTADs, and long-range looping interactions. The higher-order folding patterns of chromatin contacts within TADs and how they localize to disease-associated single nucleotide variants (daSNVs) remains an open area of investigation. Here, we analyze high-resolution Hi-C data with graph theory to understand possible mesoscale network architecture within chromatin domains. We identify a subset of TADs exhibiting strong core-periphery mesoscale structure in embryonic stem cells, neural progenitor cells, and cortical neurons. Hyper-connected core nodes co-localize with genomic segments engaged in multiple looping interactions and enriched for occupancy of the architectural protein CCCTC binding protein (CTCF). CTCF knockdown and in silico deletion of CTCF-bound core nodes disrupts core-periphery structure, whereas in silico mutation of cell type-specific enhancer or gene nodes has a negligible effect. Importantly, neuropsychiatric daSNVs are significantly more likely to localize with TADs folded into core-periphery networks compared to domains devoid of such structure. Together, our results reveal that a subset of TADs encompasses looping interactions connected into a core-periphery mesoscale network. We hypothesize that daSNVs in the periphery of genome folding networks might preserve global nuclear architecture but cause local topological and functional disruptions contributing to human disease. By contrast, daSNVs co-localized with hyper-connected core nodes might cause severe topological and functional disruptions. Overall, these findings shed new light into the mesoscale network structure of fine scale genome folding within chromatin domains and its link to common genetic variants in human disease.

Higher-order chromatin organization such as A/B compartments, TADs and chromatin loops are temporarily disrupted during mitosis. These structures are thought to organize aspects of gene regulation, and thus it is important to understand how they are re-established after mitosis. We examined the dynamics of chromosome reorganization by Hi-C at defined time points following exit from mitosis in highly purified, synchronous cell populations. We observed that A/B compartments are rapidly established and progressively gain in strength following mitotic exit. Contact domain formation occurs from the “bottom-up” with smaller sub-TADs forming initially, followed by convergence into multi-domain TAD structures. CTCF is strongly retained at a significant fraction of sites on mitotic chromosomes and immediately resumes full binding at ana/telophase, the earliest tested time point. In contrast, cohesin is completely evicted from mitotic chromosomes and resumes focal binding with delayed kinetics. The formation of CTCF/cohesin co-anchored structural loops follows the kinetics of cohesin positioning. Stripe-shaped contacts anchored by CTCF grow in length, consistent with a loop extrusion process after mitosis. Interactions between cis-regulatory elements can form rapidly, preceding CTCF/cohesin anchored structural loops. Strikingly, we identified a group of rapidly emerging transient contacts between cis-regulatory elements in ana/telophase, that are dissolved upon G1 entry, co-incident with the establishment of inner boundaries or nearby interfering loops. Our findings indicate that distinct but mutually influential forces drive post-mitotic chromatin re-configuration to shape compartments, contact domains, cis-element contacts, and CTCF/cohesin dependent loops.

More than 25 inherited human disorders are causedby the unstable expansion of repetitive DNA sequences termed short tandem repeats (STRs). Afundamental unresolved question is why some STRs are susceptible to pathologic expansion, whereas thousands of repeat tracts across the human genome are relatively stable. Here, we discover that nearly all disease-associated STRs (daSTRs) are located at boundaries demarcating 3D chromatin domains. We identify a subset of boundaries with markedly higher CpG island density compared to the rest of the genome. daSTRs specifically localize to ultra-high-density CpG island boundaries, suggesting they might be hotspots for epigenetic misregulation or topological disruption linked to STR expansion. Fragile X syndrome patients exhibit severe boundary disruption in a manner that correlates with local loss of CTCF occupancy and the degree of FMR1 silencing. Our data uncover higher-order chromatin architecture as a new dimension in understanding repeat expansion disorders.

Methods for using exogenous fluorophore and general anesthetic 1-aminanthracene (1-AMA) and its photoactive derivative 1-azidoanthracene (1-AZA) are provided. 1-AMA potentiates GABA(A) chloride currents and immobilizes Xenopus laevis tadpoles. Cellular and tissue anesthetic distribution can be imaged for quantifying “on-pathway” and “off-pathway” targets. 1-AZA shares targets with 1-AMA and offers further optoanesthetic spatial and temporal control upon near-UV laser irradiation. Furthermore, 1-AZA adduction provides screening of possible relevant anesthetic protein targets and binding site characterization. We highlight several useful imaging and binding assays to demonstrate utility of 1-AMA and its derivative 1-AZA.

More than 25 inherited neurological disorders are caused by the unstable expansion of repetitive DNA sequences termed short tandem repeats (STRs). A fundamental unresolved question is why specific STRs are susceptible to unstable expansion leading to severe pathology, whereas tens of thousands of normal-length repeat tracts across the human genome are relatively stable. Here, we unexpectedly discover that nearly all STRs associated with repeat expansion diseases are located at boundaries demarcating 3-D chromatin domains. We find that boundaries exhibit markedly higher CpG island density compared to loci internal to domains. Importantly, disease-associated STRs are specifically localized to ultra-dense CpG island-rich boundaries, suggesting that these loci might be hotspots for epigenetic instability and topological disruption upon unstable expansion. In Fragile X Syndrome, mutation-length expansion at the Fmr1 gene results in severe disruption of the boundary between TADs. Our data uncover higher-order chromatin architecture as a new dimension in understanding the mechanistic basis of repeat expansion disorders.

Bromodomain and extraterminal motif (BET) proteins are pharmacologic targets for the treatment of diverse diseases, yet the roles of individual BET family members remain unclear. We find that BRD2, but not BRD4, co-localizes with the architectural/insulator protein CCCTC-binding factor (CTCF) genome-wide. CTCF recruits BRD2 to co-bound sites whereas BRD2 is dispensable for CTCF occupancy. Disruption of a CTCF/BRD2-occupied element positioned between two unrelated genes enables regulatory influence to spread from one gene to another, suggesting that CTCF and BRD2 form a transcriptional boundary. Accordingly, single-molecule mRNA fluorescence in situ hybridization (FISH) reveals that, upon site-specific CTCF disruption or BRD2 depletion, expression of the two genes becomes increasingly correlated. HiC shows that BRD2 depletion weakens boundaries co-occupied by CTCF and BRD2, but not those that lack BRD2. These findings indicate that BRD2 supports boundary activity, and they raise the possibility that pharmacologic BET inhibitors can influence gene expression in part by perturbing domain boundary function.

Mammalian genomes are folded in a hierarchy of topologically associating domains (TADs), subTADs and looping interactions. The nested nature of chromatin domains has rendered it challenging to identify a sensitive and specific metric for detecting subTADs and quantifying their dynamic reconfiguration across cellular states. Here, we apply graph theoretic principles to quantify hierarchical folding patterns in high-resolution chromatin topology maps. We discover that TADs can be accurately detected using a Louvain-like locally greedy algorithm to maximize network modularity. By varying a resolution parameter in the modularity quality function, we accurately partition the mouse genome across length scales into a hierarchical nested structure of network communities exhibiting a wide range of sizes. To distinguish high probability subTADs from the full detected set, we developed and applied a new ‘hierarchical spatial variance minimization’ method. Moreover, we identified a large number of dynamically altered communities between pluripotent embryonic stem cells and multipotent neural progenitor cells. Cell type specific boundaries correlate with trends in dynamic occupancy of the architectural protein CTCF, thereby validating their biological relevance. Together, these data demonstrate the utility of metrics from network science in quantifying a nested hierarchy of dynamic 3D chromatin communities across length scales. Our findings are significant toward unraveling the link between higher-order genome folding and gene expression during healthy development and the deregulation of molecular pathways linked to disease.

Previous work by us showed that in replica exchange molecular dynamics, exchanges should be attempted extremely often, providing gains in efficiency and no undesired effects. Since that time some questions have been raised about the extendability of these claims to the general case. In this work, we answer this question in two ways. First, we perform a study measuring the effect of exchange attempt frequency in explicit solvent simulations including thousands of atoms. This shows, consistent with the previous assertion, that high exchange attempt frequency allows an optimal rate of exploration of configurational space. Second, we present an explanation of many theoretical and technical pitfalls when implementing replica exchange that cause "improper" exchanges resulting in erroneous data, exacerbated by high exchange attempt frequency.

Pharmacologic inhibitors of the bromodomain and extraterminal motif (BET) family of proteins have shown promise in the treatment of hematologic and other malignancies and are being developed for clinical use. However, BET inhibitors do not discriminate between the family members BRD2, BRD3, and BRD4, and thus the mechanistic basis for their therapeutic efficacy is not well understood. In addition, BRD2 and BRD4 are individually required for the activation of genes driven by the erythroid transcription factor GATA1 (Stonestrom et al., Blood 2015), yet how BRD2 in particular contributes to this process has not been studied. We examined BRD2 occupancy genome-wide in erythroid cells and find that BRD2 colocalizes extensively with the architectural/insulator protein CCCTC-binding factor (CTCF). We define a functional hierarchy whereby CTCF is required for BRD2 to occupy co-bound sites, while CTCF binding is BRD2-independent. Using CRISPR/Cas9-based genome editing we identify a boundary element occupied by CTCF and BRD2 that is adjacent to a GATA1-driven enhancer and ensures appropriate transcriptional regulation at the locus. Employing single-molecule RNA FISH we show that either site-specific CTCF disruption or BRD2 depletion leads to increased correlation in mature mRNA levels between the genes flanking this boundary, suggesting that they become inappropriately coregulated. Taken together these findings indicate that BRD2 collaborates with CTCF to constrain the activity of an erythroid enhancer and reveal a potential new role for BET proteins in chromatin domain boundary function. Disclosures No relevant conflicts of interest to declare.

Recently, we identified 1-aminoanthracene as a fluorescent general anesthetic. To investigate the mechanism of action, a photoactive analogue, 1-azidoanthracene, was synthesized. Administration of 1-azidoanthracene to albino stage 40-47 tadpoles was found to immobilize animals upon near-UV irradiation of the forebrain region. The immobilization was often reversible, but it was characterized by a longer duration consistent with covalent attachment of the ligand to functionally important targets. IEF/SDS-PAGE examination of irradiated tadpole brain homogenate revealed labeled protein, identified by mass spectrometry as β-tubulin. In vitro assays with aminoanthracene-cross-linked tubulin indicated inhibition of microtubule polymerization, similar to colchicine. Tandem mass spectrometry confirmed anthracene binding near the colchicine site. Stage 40-47 tadpoles were also incubated 1 h with microtubule stabilizing agents, epothilone D or discodermolide, followed by dosing with 1-aminoanthracene. The effective concentration of 1-aminoanthracene required to immobilize the tadpoles was significantly increased in the presence of either microtubule stabilizing agent. Epothilone D similarly mitigated the effects of a clinical neurosteroid general anesthetic, allopregnanolone, believed to occupy the colchicine site in tubulin. We conclude that neuronal microtubules are "on-pathway" targets for anthracene general anesthetics and may also represent functional targets for some neurosteroid general anesthetics.