Chromatin Domains Research Articles

LAP2β transmits force to upregulate genes via chromatin domain stretching but not compression.

There is increasing evidence that force impacts almost every aspect of cells and tissues in physiology and disease including gene regulation. However, the molecular pathway of force transmission from the nuclear lamina to the chromatin remain largely elusive. Here we employ two different approaches of a local stress on cell apical surface via an RGD (Arg-Gly-Asp)-coated magnetic bead and whole cell deformation at cell basal surface via uniaxial or biaxial deformation of a fibronectin-coated flexible polydimethylsiloxane substrate. We find that nuclear protein LAP2β mediates force transmission from the nuclear lamina to the chromatin. Knocking down LAP2β increases spontaneous movements of the chromatin by reducing tethering of the chromatin and substantially inhibits the magnetic bead-stress or the substrate-deformation induced chromatin domain stretching and the ensuing dihydrofolate reductase (DHFR) gene upregulation. Analysis of DHFR gene-containing chromatin domain alignments along or perpendicular to the direction of the stretching/compressing reveals that the chromatin domain must be stretched and not compressed in order for the gene to be rapidly upregulated. Together these results suggest that external-load induced rapid transcription upregulation originates from chromatin domain stretching but not compressing and depends on the molecular force transmission pathway of LAP2β. STATEMENT OF SIGNIFICANCE: How force regulates gene expression has been elusive. Here we show that the orientation of the chromatin domain relative to the stress direction is crucial in determining if the chromatin domain will be stretched or compressed in response to a cell surface loading. We also show that nuclear protein Lap2b is a critical molecule that mediates force transmission from the nuclear laminar to the chromatin to regulate gene transcription. This study reveals the molecular force transmission pathway for force-induced gene regulation.

DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Modeling of ionizing radiation-induced chromosome aberration and tumor prevalence based on two classes of DNA double-strand breaks clustering in chromatin domains

There has been some controversy over the use of radiobiological models when modeling the dose-response curves of ionizing radiation (IR)-induced chromosome aberration and tumor prevalence, as those curves usually show obvious non-targeted effects (NTEs) at low doses of high linear energy transfer (LET) radiation. The lack of understanding the contribution of NTEs to IR-induced carcinogenesis can lead to distinct deviations of relative biological effectiveness (RBE) estimations of carcinogenic potential, which are widely used in radiation risk assessment and radiation protection. In this work, based on the initial pattern of two classes of IR-induced DNA double-strand breaks (DSBs) clustering in chromatin domains and the subsequent incorrect repair processes, we proposed a novel radiobiological model to describe the dose-response curves of two carcinogenic-related endpoints within the same theoretical framework. The representative experimental data was used to verify the consistency and validity of the present model. The fitting results indicated that, compared with targeted effect (TE) and NTE models, the current model has better fitting ability when dealing with the experimental data of chromosome aberration and tumor prevalence induced by multiple types of IR with different LETs. Notably, the present model without introducing an NTE term was adequate to describe the dose-response curves of IR-induced chromosome aberration and tumor prevalence with NTEs in low-dose regions. Based on the fitting parameters, the LET-dependent RBE values were calculated for three given low doses. Our results showed that the RBE values predicted by the current model gradually decrease with the increase of doses for the endpoints of chromosome aberration and tumor prevalence. In addition, the calculated RBE was also compared with those evaluated from other models. These analyses show that the proposed model can be used as an alternative tool to well describe dose-response curves of multiple carcinogenic-related endpoints and effectively estimate RBE in low-dose regions.

SUMO control of centromere homeostasis.

Centromeres are unique chromosomal loci that form the anchorage point for the mitotic spindle during mitosis and meiosis. Their position and function are specified by a unique chromatin domain featuring the histone H3 variant CENP-A. While typically formed on centromeric satellite arrays, CENP-A nucleosomes are maintained and assembled by a strong self-templated feedback mechanism that can propagate centromeres even at non-canonical sites. Central to the epigenetic chromatin-based transmission of centromeres is the stable inheritance of CENP-A nucleosomes. While long-lived at centromeres, CENP-A can turn over rapidly at non-centromeric sites and even erode from centromeres in non-dividing cells. Recently, SUMO modification of the centromere complex has come to the forefront as a mediator of centromere complex stability, including CENP-A chromatin. We review evidence from different models and discuss the emerging view that limited SUMOylation appears to play a constructive role in centromere complex formation, while polySUMOylation drives complex turnover. The deSUMOylase SENP6/Ulp2 and the proteins segregase p97/Cdc48 constitute the dominant opposing forces that balance CENP-A chromatin stability. This balance may be key to ensuring proper kinetochore strength at the centromere while preventing ectopic centromere formation.

HMGXB4 Targets Sleeping Beauty Transposition to Germinal Stem Cells.

Transposons are parasitic genetic elements that frequently hijack vital cellular processes of their host. HMGXB4 is a known Wnt signaling-regulating HMG-box protein, previously identified as a host-encoded factor of Sleeping Beauty (SB) transposition. Here, we show that HMGXB4 is predominantly maternally expressed, and marks both germinal progenitor and somatic stem cells. SB piggybacks HMGXB4 to activate transposase expression and target transposition to germinal stem cells, thereby potentiating heritable transposon insertions. The HMGXB4 promoter is located within an active chromatin domain, offering multiple looping possibilities with neighboring genomic regions. HMGXB4 is activated by ERK2/MAPK1, ELK1 transcription factors, coordinating pluripotency and self-renewal pathways, but suppressed by the KRAB-ZNF/TRIM28 epigenetic repression machinery, also known to regulate transposable elements. At the post-translational level, SUMOylation regulates HMGXB4, which modulates binding affinity to its protein interaction partners and controls its transcriptional activator function via nucleolar compartmentalization. When expressed, HMGXB4 can participate in nuclear-remodeling protein complexes and transactivate target gene expression in vertebrates. Our study highlights HMGXB4 as an evolutionarily conserved host-encoded factor that assists Tc1/Mariner transposons to target the germline, which was necessary for their fixation and may explain their abundance in vertebrate genomes.

Abstract 3558: Role of FLI portion of EWS::FLI in transcription regulation via modulation of chromatin 3D landscape in Ewing sarcoma

Abstract Ewing sarcoma is an aggressive bone-associated tumor currently treated with dose-intense chemotherapy, radiation, and surgery and it affects adolescents and young adults. The hallmark of Ewing sarcoma is a translocated fusion transcription factor named EWS::FLI that drives the oncogenic process. We hypothesize that FLI portion of EWS::FLI containing a crucial alpha-helix plays a novel role in transcription regulation of thousands of genes by modulating chromatin looping. The hypothesis is based on recent evidences from our lab: an alpha-helix immediately downstream of DNA binding domain as an important player in regulating transcriptional activity and that EWS::FLI has a substantial role in shaping the chromatin landscape of Ewing sarcoma cells. The objective of this study is to elucidate the mechanism underlying transcriptional regulation by FLI. We utilized knockdown/rescue experiments in which EWS::FLI was depleted with shRNA and replaced with constructs containing the ETS DNA Binding Domain (DBD) alone or DBD+ (DBD and 4th alpha-helix). The binding pattern and transcriptional regulation were assessed with CUT&Tag and RNA-Seq respectively. The extent of roles each of these construct play in organization, structure, and function of chromatin are being assessed using Micro-C technique, a variation of Hi-C with improved resolution, higher signal-to-noise ratio and more information on chromatin domain boundaries and chromatin looping. The DNA binding and genomic localization of EWS::FLI was unaltered by the deletion surrounding the DNA binding domain (which contains a 4th alpha-helix) in A673 and TTC466 cells. Despite this similarity in genomic localization and binding, the transcriptional output driven by EWS::FLI was significantly diminished by the deletion. With this current study, we hope to understand if the flanking region neighboring the DNA binding domain contributes to the chromatin architecture remodeling function of EWS::FLI and whether this function could be attributed to the transcriptional output differences in the DBD and DBD+ conditions. Citation Format: Ariunaa Bayanjargal, Cenny Taslim, Jesse Crow, Julia Selich-Anderson, Stephen L. Lessnick. Role of FLI portion of EWS::FLI in transcription regulation via modulation of chromatin 3D landscape in Ewing sarcoma. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3558.

A change in biophysical properties accompanies heterochromatin formation in mouse embryos.

The majority of our genome is composed of repeated DNA sequences that assemble into heterochromatin, a highly compacted structure that constrains their mutational potential. How heterochromatin forms during development and how its structure is maintained are not fully understood. Here, we show that mouse heterochromatin phase-separates after fertilization, during the earliest stages of mammalian embryogenesis. Using high-resolution quantitative imaging and molecular biology approaches, we show that pericentromeric heterochromatin displays properties consistent with a liquid-like state at the two-cell stage, which change at the four-cell stage, when chromocenters mature and heterochromatin becomes silent. Disrupting the condensates results in altered transcript levels of pericentromeric heterochromatin, suggesting a functional role for phase separation in heterochromatin function. Thus, our work shows that mouse heterochromatin forms membrane-less compartments with biophysical properties that change during development and provides new insights into the self-organization of chromatin domains during mammalian embryogenesis.

P97/VCP drives turnover of SUMOylated centromeric CCAN proteins and CENP-A.

The centromere is a unique chromatin domain that links sister chromatids and forms the attachment site for spindle microtubules in mitosis. Centromere inheritance is largely DNA sequence-independent but strongly reliant on a self-propagating chromatin domain featuring nucleosomes containing the H3 variant CENP-A. Unlike other histones, CENP-A is maintained with unusually high stability in chromatin. Previously, we have shown that mitotic maintenance of CENP-A and other constitutive centromere-associated network (CCAN) proteins is controlled by a dynamic SUMO cycle and that the deSUMOylase SENP6 is necessary for stable maintenance of CENP-A at the centromere. Here, we discover that the removal of SENP6 leads to a rapid loss of the CCAN, followed by a delayed loss of centromeric CENP-A, indicating that the CCAN is the primary SUMO target. We found that the ATP-dependent segregase p97/VCP removes centromeric CENP-A in a SUMO-dependent manner and interacts physically with the CCAN and CENP-A chromatin. Our data suggest a direct role of p97 in removing centromeric CENP-A via SUMOylated CCAN proteins, thereby ensuring centromere homeostasis and potentially preventing ectopic CENP-A accumulation.

Whole-genome doubling drives oncogenic loss of chromatin segregation

Whole-genome doubling (WGD) is a recurrent event in human cancers and it promotes chromosomal instability and acquisition of aneuploidies1–8. However, the three-dimensional organization of chromatin in WGD cells and its contribution to oncogenic phenotypes are currently unknown. Here we show that in p53-deficient cells, WGD induces loss of chromatin segregation (LCS). This event is characterized by reduced segregation between short and long chromosomes, A and B subcompartments and adjacent chromatin domains. LCS is driven by the downregulation of CTCF and H3K9me3 in cells that bypassed activation of the tetraploid checkpoint. Longitudinal analyses revealed that LCS primes genomic regions for subcompartment repositioning in WGD cells. This results in chromatin and epigenetic changes associated with oncogene activation in tumours ensuing from WGD cells. Notably, subcompartment repositioning events were largely independent of chromosomal alterations, which indicates that these were complementary mechanisms contributing to tumour development and progression. Overall, LCS initiates chromatin conformation changes that ultimately result in oncogenic epigenetic and transcriptional modifications, which suggests that chromatin evolution is a hallmark of WGD-driven cancer.

Heterochromatin rewiring and domain disruption-mediated chromatin compaction during erythropoiesis.

Mammalian erythropoiesis involves progressive chromatin compaction and subsequent enucleation in terminal differentiation, but the mechanisms underlying the three-dimensional chromatin reorganization remain obscure. Here, we systematically analyze the higher-order chromatin in purified populations of primary human erythroblasts. Our results reveal that heterochromatin regions undergo substantial compression, with H3K9me3 markers relocalizing to the nuclear periphery and forming a significant number of long-range interactions, and that ~58% of the topologically associating domain (TAD) boundaries are disrupted, while certain TADs enriched for markers of the active transcription state and erythroid master regulators, GATA1 and KLF1, are selectively maintained during terminal erythropoiesis. Finally, we demonstrate that GATA1 is involved in safeguarding selected essential chromatin domains during terminal erythropoiesis. Our study therefore delineates the molecular characteristics of a development-driven chromatin compaction process, which reveals transcription competence as a key indicator of the selected domain maintenance to ensure appropriate gene expression during the extreme compaction of chromatin.

SLFN5-mediated chromatin dynamics sculpt higher-order DNA repair topology

Repair of DNA double-strand breaks (DSBs) elicits three-dimensional (3D) chromatin topological changes. A recent finding reveals that 53BP1 assembles into a 3D chromatin topology pattern around DSBs. How this formation of a higher-order structure is configured and regulated remains enigmatic. Here, we report that SLFN5 is a critical factor for 53BP1 topological arrangement at DSBs. Using super-resolution imaging, we find that SLFN5 binds to 53BP1 chromatin domains to assemble a higher-order microdomain architecture by driving damaged chromatin dynamics at both DSBs and deprotected telomeres. Mechanistically, we proposethat 53BP1 topology is shaped by two processes: (1) chromatin mobility driven by the SLFN5-LINC-microtubule axis and (2) the assembly of 53BP1 oligomers mediated by SLFN5. In mammals, SLFN5 deficiency disrupts the DSB repair topology and impairs non-homologous end joining, telomere fusions, class switch recombination, and sensitivity to poly (ADP-ribose) polymerase inhibitor. We establish a molecular mechanism that shapes higher-order chromatin topologies to safeguard genomic stability.

On the Advent of Super-Resolution Microscopy in the Realm of Polycomb Proteins.

The genomes of metazoans are organized at multiple spatial scales, ranging from the double helix of DNA to whole chromosomes. The intermediate genomic scale of kilobases to megabases, which corresponds to the 50-300 nm spatial scale, is particularly interesting, as the 3D arrangement of chromatin is implicated in multiple regulatory mechanisms. In this context, polycomb group (PcG) proteins stand as major epigenetic modulators of chromatin function, acting prevalently as repressors of gene transcription by combining chemical modifications of target histones with physical crosslinking of distal genomic regions and phase separation. The recent development of super-resolution microscopy (SRM) has strongly contributed to improving our comprehension of several aspects of nano-/mesoscale (10-200 nm) chromatin domains. Here, we review the current state-of-the-art SRM applied to PcG proteins, showing that the application of SRM to PcG activity and organization is still quite limited and mainly focused on the 3D assembly of PcG-controlled genomic loci. In this context, SRM approaches have mostly been applied to multilabel fluorescence in situ hybridization (FISH). However, SRM data have complemented the maps obtained from chromosome capture experiments and have opened a new window to observe how 3D chromatin topology is modulated by PcGs.

The Histone Variant H2A.Z C-Terminal Domain Has Locus-Specific Differential Effects on H2A.Z Occupancy and Nucleosome Localization.

The incorporation of histone variant H2A.Z into nucleosomes creates specialized chromatin domains that regulate DNA-templated processes, such as gene transcription. In Saccharomyces cerevisiae, the diverging H2A.Z C terminus is thought to provide the H2A.Z exclusive functions. To elucidate the roles of this H2A.Z C terminus genome-wide, we used derivatives in which the C terminus was replaced with the corresponding region of H2A (ZA protein), or the H2A region plus a transcriptional activating peptide (ZA-rII'), with the intent of regenerating the H2A.Z-dependent regulation globally. The distribution of these H2A.Z derivatives indicates that the H2A.Z C-terminal region is crucial for both maintaining the occupation level of H2A.Z and the proper positioning of targeted nucleosomes. Interestingly, the specific contribution on incorporation efficiency versus nucleosome positioning varies enormously depending on the locus analyzed. Specifically, the role of H2A.Z in global transcription regulation relies on its C-terminal region. Remarkably, however, this mostly involves genes without a H2A.Z nucleosome in the promoter. Lastly, we demonstrate that the main chaperone complex which deposits H2A.Z to gene regulatory region (SWR1-C) is necessary to localize all H2A.Z derivatives at their specific loci, indicating that the differential association of these derivatives is not due to impaired interaction with SWR1-C. IMPORTANCE We provide evidence that the Saccharomyces cerevisiae C-terminal region of histone variant H2A.Z can mediate its special function in performing gene regulation by interacting with effector proteins and chaperones. These functional interactions allow H2A.Z not only to incorporate to very specific gene regulatory regions, but also to facilitate the gene expression process. To achieve this, we used a chimeric protein which lacks the native H2A.Z C-terminal region but contains an acidic activating region, a module that is known to interact with components of chromatin-remodeling entities and/or transcription modulators. We reasoned that because this activating region can fulfill the role of the H2A.Z C-terminal region, at least in part, the role of the latter would be to interact with these activating region targets.

Assembly of the 81.6 Mb centromere of pea chromosome 6 elucidates the structure and evolution of metapolycentric chromosomes.

Centromeres in the legume genera Pisum and Lathyrus exhibit unique morphological characteristics, including extended primary constrictions and multiple separate domains of centromeric chromatin. These so-called metapolycentromeres resemble an intermediate form between monocentric and holocentric types, and therefore provide a great opportunity for studying the transitions between different types of centromere organizations. However, because of the exceedingly large and highly repetitive nature of metapolycentromeres, highly contiguous assemblies needed for these studies are lacking. Here, we report on the assembly and analysis of a 177.6 Mb region of pea (Pisum sativum) chromosome 6, including the 81.6 Mb centromere region (CEN6) and adjacent chromosome arms. Genes, DNA methylation profiles, and most of the repeats were uniformly distributed within the centromere, and their densities in CEN6 and chromosome arms were similar. The exception was an accumulation of satellite DNA in CEN6, where it formed multiple arrays up to 2 Mb in length. Centromeric chromatin, characterized by the presence of the CENH3 protein, was predominantly associated with arrays of three different satellite repeats; however, five other satellites present in CEN6 lacked CENH3. The presence of CENH3 chromatin was found to determine the spatial distribution of the respective satellites during the cell cycle. Finally, oligo-FISH painting experiments, performed using probes specifically designed to label the genomic regions corresponding to CEN6 in Pisum, Lathyrus, and Vicia species, revealed that metapolycentromeres evolved via the expansion of centromeric chromatin into neighboring chromosomal regions and the accumulation of novel satellite repeats. However, in some of these species, centromere evolution also involved chromosomal translocations and centromere repositioning.

Super-resolved heterochromatin structure during DNA damage response revealed by CRISPR imaging.

Heterochromatin plays an important role in gene expression regulation and accurate chromosome segregation. Heterochromatin undergoes a DNA damage repair process distinct from that of euchromatin, involving large-scale rearrangement coupled to nuclear membrane and motor proteins. Revealing the nanoscale structure and dynamics of heterochromatin during its repair is essential for understanding the detailed mechanism of the wholeness DNA repair process, and the role of the distinct compaction and reorganization of heterochromatin. Here, we observed the novel structure and dynamics of heterochromatin in super resolution using CRISPR labeling techniques. Dense labeling of a pericentromeric region with CRISPR prevented cell cycle progression and induced DNA damage response accompanied by domain expansion and fiber-like extension beyond its chromosome territory. This structural change correlates with the DNA damage repair proteins. Due to persistent binding of an array of CRISPR complexes, the cell failed to repair the perceived damage, resulting in senescence-like phenotype. Use of advanced CRISPR labeling techniques allowed super-resolution imaging of the chromatin domain under repair in live cells. Our work highlights the DNA damage response in heterochromatin induced by molecular blockades and the resulting structural changes, providing a platform to study DNA repair process with high spatiotemporal resolution.

Optical nanoscopy approach to unveil chromatin alterations implicated in neuroblastoma transformation.

Chromatin organization is a well-orchestrated mechanism regulating cellular physiology and its changes are associated with the onset of many diseases, including cancer. Neuroblastoma (NB) is the most common extracranial tumor in childhood, and it shows remarkable heterogeneity in the phenotypes and evolutionary patterns. Although the molecular mechanisms of this heterogeneity are still unclear, the impairment of the “epigenetic machinery” and the consequent chromatin reorganization might play an important role. We aim to assess the possible alterations in chromatin nanoscale architecture correlated with NB transformation using NB cells genetically engineered to overexpress the non-coding RNA NDM29, which promotes differentiation of tumor cells into a neuronal lineage. In both malignant NB cells (Mock) and neuron-like cells (S1.1), we explored chromatin organization and nuclear architecture by optical nanoscale microscopy using both confocal and stimulation emission depletion (STED) imaging methods, in parallel with molecular biology analyses. We characterized the nuclear 3D morphology of Mock vs. S1.1 cells by measuring volume, shape, and circularity of the nuclei as well as chromatin compaction and distribution. We imaged the immunolabeled histone H3 modifications (H3K9 trimethylation and H3K9 acetylation) to characterize chromatin regions typically associated with highly and poorly condensed chromatin, respectively. In parallel, we quantified the global DNA methylation status by measuring the levels of 5-methylcytosine (5-mC) using the enzyme-linked immunosorbent assay (ELISA) MethylFlash Global DNA. Comparing Mock vs. S1.1 cells, we observed (i) significant remodeling of nuclear morphometry; (ii) marked rearrangement of chromatin domains rich in highly-acetylated and -methylated H3K9; (iii) changes in average chromatin compaction and DNA methylation. Altogether, the results point to an involvement of epigenetic mechanisms which may sustain the chromatin remodeling during NB malignant transformation, paving the way towards new prognostic and therapeutic approaches and understanding chromatin remodeling in physiological/pathological processes.

Illuminating dynamic structural states in hiPSC through endogenous gene tagging.

The Allen Institute for Cell Science aims to understand the principles by which human induced pluripotent stem cells (hiPSCs) establish and maintain robust dynamic localization of cellular structures, and how they transition between states during differentiation and disease. As a first step [JG1], we have created a collection of endogenously tagged hiPSC lines to illuminate cell organization. The tagged proteins include notable membrane-bound and membrane-less cellular organelles, signaling complexes, phase transition markers, transcription factors, and structural markers specific to differentiated cells. To date, the Allen Cell Collection consists of 54 single- or dual-edited lines that have undergone extensive quality control testing to ensure genomic, cell biological, and stem cell integrity. Here, we highlight our gene-editing and quality control workflow for mono- and biallelic editing of expressed or silent genes that are expressed specifically during differentiation and function as reporters of cellular state. Furthermore, we underscore our efforts to increase efficiency of the current workflow by utilizing Adeno-Associated Virus (AAV) to deliver donor DNAs and multiplexing transfection strategies for gene tagging at multiple loci. Our most recently released lines are CTCF-mEGFP and PCNA-mEGFP, for dynamic visualization of chromatin domains and replication foci, G3BP1-mEGFP to label stress granules, mono- and bi-allelic mEGFP-E-cadherin to visualize adherens junctions and DCP1a-mEGFP for P-bodies. We now have gene [RG2] [JG3]-edited cell populations for N- and VE-cadherin ready for single cell sorting. Our overarching goal is to make these tools openly available to the cell biology community to help accelerate biomedical research. In addition to our cell lines, the donor plasmids, segmented 3D images of about 200,000 cells from our lines, image analysis, visualization tools, integrated cell models and biological findings are all openly available to the research community (www.allencell.org).

Genomic glucocorticoid action in embryonic mouse neural stem cells

Prenatal exposure to synthetic glucocorticoids (sGCs) reprograms brain development and predisposes the developing fetus towards potential adverse neurodevelopmental outcomes. Using a mouse model of sGC administration, previous studies show that these changes are accompanied by sexually dimorphic alterations in the transcriptome of neural stem and progenitor cells (NSPCs) derived from the embryonic telencephalon. Because cell type-specific gene expression profiles tightly regulate cell fate decisions and are controlled by a flexible landscape of chromatin domains upon which transcription factors and enhancer elements act, we multiplexed data from four genome-wide assays: RNA-seq, ATAC-seq (assay for transposase accessible chromatin followed by genome wide sequencing), dual cross-linking ChIP-seq (chromatin immunoprecipitation followed by genome wide sequencing), and microarray gene expression to identify novel relationships between gene regulation, chromatin structure, and genomic glucocorticoid receptor (GR) action in NSPCs. These data reveal that GR binds preferentially to predetermined regions of accessible chromatin to influence gene programming and cell fate decisions. In addition, we identify SOX2 as a transcription factor that impacts the genomic response of select GR target genes to sGCs (i.e., dexamethasone) in NSPCs.

Deciphering Hierarchical Chromatin Domains and Preference of Genomic Position Forming Boundaries in Single Mouse Embryonic Stem Cells.

The exploration of single-cell 3D genome maps reveals that chromatin domains are indeed physical structures presenting in single cells, and domain boundaries vary from cell to cell. However, systematic analysis of the association between regulatory factor binding and elements and the formation of chromatin domains in single cells has not yet emerged. To this end, a hierarchical chromatin domain structure identification algorithm (named as HiCS) is first developed from individual single-cell Hi-C maps, with superior performance in both accuracy and efficiency. The results suggest that in addition to the known CTCF-cohesin complex, Polycomb, TrxG, pluripotent protein families, and other multiple factors also contribute to shaping chromatin domain boundaries in single embryonic stem cells. Different cooperation patterns of these regulatory factors drive genomic position categories with differential preferences forming boundaries, and the most extensive six types of retrotransposons are differentially distributed in these genomic position categories with preferential localization. The above results suggest that these different retrotransposons within genomic regions interplay with regulatory factors navigating the preference of genomic positions forming boundaries, driving the formation of higher-order chromatin structures, and thus regulating cell functions in single mouse embryonic stem cells.

Active enhancers strengthen insulation by RNA-mediated CTCF binding at chromatin domain boundaries.

Vertebrate genomes are partitioned into chromatin domains or topologically associating domains (TADs), which are typically bound by head-to-head pairs of CTCF binding sites. Transcription at domain boundaries correlates with better insulation; however, it is not known whether the boundary transcripts themselves contribute to boundary function. Here we characterize boundary-associated RNAs genome-wide, focusing on the disease-relevant INK4a/ARF and MYC TAD. Using CTCF site deletions and boundary-associated RNA knockdowns, we observe that boundary-associated RNAs facilitate recruitment and clustering of CTCF at TAD borders. The resulting CTCF enrichment enhances TAD insulation, enhancer-promoter interactions, and TAD gene expression. Importantly, knockdown of boundary-associated RNAs results in loss of boundary insulation function. Using enhancer deletions and CRISPRi of promoters, we show that active TAD enhancers, but not promoters, induce boundary-associated RNA transcription, thus defining a novel class of regulatory enhancer RNAs.