A single-cell chromatin accessibility dataset of human primed and naïve pluripotent stem cell-derived teratoma.

Teratoma, due to its remarkable ability to differentiate into multiple cell lineages, is a valuable model for studying human embryonic development. The similarity of the gene expression and chromatin accessibility patterns in these cells to those observed in vivo further underscores its potential as a research tool. Notably, teratomas derived from human naïve (pre-implantation epiblast-like) pluripotent stem cells (PSCs) have larger embryonic cell diversity and contain extraembryonic lineages, making them more suitable to study developmental processes. However, the cell type-specific epigenetic profiles of naïve PSC teratomas have not been yet characterized. Using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we analyzed 66,384 cell profiles from five teratomas derived from human naïve PSCs and their post-implantation epiblast-like (primed) counterparts. We observed 17 distinct cell types from both embryonic and extraembryonic lineages, resembling the corresponding cell types in human fetal tissues. Additionally, we identified key transcription factors specific to different cell types. Our dataset provides a resource for investigating gene regulatory programs in a relevant model of human embryonic development.

nucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome.

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

SETDB1 regulates short interspersed nuclear elements and chromatin loop organization in mouse neural precursor cells.

BACKGROUND: Transposable elements play a critical role in maintaining genome architecture during neurodevelopment. Short Interspersed Nuclear Elements (SINEs), a major subtype of transposable elements, are known to harbor binding sites for the CCCTC-binding factor (CTCF) and pivotal in orchestrating chromatin organization. However, the regulatory mechanisms controlling the activity of SINEs in the developing brain remains elusive. RESULTS: In our study, we conduct a comprehensive genome-wide epigenetic analysis in mouse neural precursor cells using ATAC-seq, ChIP-seq, whole genome bisulfite sequencing, in situ Hi-C, and RNA-seq. Our findings reveal that the SET domain bifurcated histone lysine methyltransferase 1 (SETDB1)-mediated H3K9me3, in conjunction with DNA methylation, restricts chromatin accessibility on a selective subset of SINEs in neural precursor cells. Mechanistically, loss of Setdb1 increases CTCF access to these SINE elements and contributes to chromatin loop reorganization. Moreover, de novo loop formation contributes to differential gene expression, including the dysregulation of genes enriched in mitotic pathways. This leads to the disruptions of cell proliferation in the embryonic brain after genetic ablation of Setdb1 both in vitro and in vivo. CONCLUSIONS: In summary, our study sheds light on the epigenetic regulation of SINEs in mouse neural precursor cells, suggesting their role in maintaining chromatin organization and cell proliferation during neurodevelopment.

The chromatin accessibility and transcriptomic landscape of the aging mice cochlea and the identification of potential functional super-enhancers in age-related hearing loss.

BACKGROUND: Presbycusis, also referred to as age-related hearing loss (ARHL), is a condition that results from the cumulative effects of aging on an individual's auditory capabilities. Given the limited understanding of epigenetic mechanisms in ARHL, our research focuses on alterations in chromatin-accessible regions. METHODS: We employed assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) in conjunction with unique identifier (UID) mRNA-seq between young and aging cochleae, and conducted integrated analysis as well as motif/TF-gene prediction. Additionally, the essential role of super-enhancers (SEs) in the development of ARHL was identified by comparative analysis to previous research. Meanwhile, an ARHL mouse model and an aging mimic hair cell (HC) model were established with a comprehensive identification of senescence phenotypes to access the role of SEs in ARHL progression. RESULTS: The control cochlear tissue exhibited greater chromatin accessibility than cochlear tissue affected by ARHL. Furthermore, the levels of histone 3 lysine 27 acetylation were significantly depressed in both aging cochlea and aging mimic HEI-OC1 cells, highlighting the essential role of SEs in the development of ARHL. The potential senescence-associated super-enhancers (SASEs) of ARHL were identified, most of which exhibited decreased chromatin accessibility. The majority of genes related to the SASEs showed obvious decreases in mRNA expression level in aging HCs and was noticeably altered following treatment with JQ1 (a commonly used SE inhibitor). CONCLUSION: The chromatin accessibility in control cochlear tissue was higher than that in cochlear tissue affected by ARHL. Potential SEs involved in ARHL were identified, which might provide a basis for future therapeutics targeting SASEs related to ARHL.

Topological Spectra and Entropy of Chromatin Loop Networks.

The 3D folding of a mammalian gene can be studied by a polymer model, where the chromatin fiber is represented by a semiflexible polymer which interacts with multivalent proteins, representing complexes of DNA-binding transcription factors and RNA polymerases. This physical model leads to the natural emergence of clusters of proteins and binding sites, accompanied by the folding of chromatin into a set of topologies, each associated with a different network of loops. Here, we combine numerics and analytics to first classify these networks and then find their relative importance or statistical weight, when the properties of the underlying polymer are those relevant to chromatin. Unlike polymer networks previously studied, our chromatin networks have finite average distances between successive binding sites, and this leads to giant differences between the weights of topologies with the same number of edges and nodes but different wiring. These weights strongly favor rosettelike structures with a local cloud of loops with respect to more complicated nonlocal topologies. Our results suggest that genes should overwhelmingly fold into a small fraction of all possible 3D topologies, which can be robustly characterized by the framework we propose here.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia.

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Intrauterine growth in chromatinopathies: A long road for better understanding and for improving clinical management.

BACKGROUND: Chromatinopathies are a heterogeneous group of genetic disorders caused by pathogenic variants in genes coding for chromatin state balance proteins. Remarkably, many of these syndromes present unbalanced postnatal growth, both under- and over-, although little has been described in the literature. Fetal growth measurements are common practice in pregnancy management and values within normal ranges indicate proper intrauterine growth progression; on the contrary, abnormalities in intrauterine fetal growth open the discussion of possible pathogenesis affecting growth even in the postnatal period. METHODS: Among the numerous chromatinopathies, we have selected six of the most documented in the literature offering evidence about two fetal overgrowth (Sotos and Weaver syndrome) and four fetal undergrowth syndromes (Bohring Opitz, Cornelia de Lange, Floating-Harbor, and Meier Gorlin syndrome), describing their molecular characteristics, maternal biochemical results and early pregnancy findings, prenatal ultrasound findings, and postnatal characteristics. RESULTS/CONCLUSION: To date, the scarce data in the literature on prenatal findings are few and inconclusive, even though these parameters may contribute to a more rapid and accurate diagnosis, calling for a better and more detailed description of pregnancy findings.

Chromatin Immunoprecipitation Using Imbibed Seeds of Arabidopsis thaliana.

Chromatin immunoprecipitation (ChIP) is used to analyze the targeting of a protein to a specific region of chromatin in vivo. Here, we present an instructive ChIP protocol for Arabidopsis imbibed seeds. The protocol covers all steps, from the sampling of imbibed seeds to the reverse crosslinking of immunoprecipitated protein-DNA complexes, and includes experimental tips and notes. The targeting of the protein to DNA is determined by quantitative PCR (qPCR) using reverse crosslinked DNA. The protocol can be further scaled up for ChIP-sequencing (ChIP-seq) analysis. As an example of the protocol, we include a ChIP-quantitative PCR (ChIP-qPCR) analysis demonstrating the targeting of PIF1 to the ABI5 promoter.

Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish.

BACKGROUND: Cis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq). RESULTS: Using the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions. CONCLUSIONS: Thus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

Polymerization of ZBTB transcription factors regulates chromatin occupancy.

BCL6, an oncogenic transcription factor (TF), forms polymers in the presence of a small-molecule molecular glue that stabilizes a complementary interface between homodimers of BCL6's broad-complex, tramtrack, and bric-à-brac (BTB) domain. The BTB domains of other proteins, including a large class of TFs, have similar architectures and symmetries, raising the possibility that additional BTB proteins self-assemble into higher-order structures. Here, we surveyed 189 human BTB proteins with a cellular fluorescent reporter assay and identified 18 ZBTB TFs that show evidence of polymerization. Through biochemical and cryoelectron microscopy (cryo-EM) studies, we demonstrate that these ZBTB TFs polymerize into filaments. We found that BTB-domain-mediated polymerization of ZBTB TFs enhances chromatin occupancy within regions containing homotypic clusters of TF binding sites, leading to repression of target genes. Our results reveal a role of higher-order structures in regulating ZBTB TFs and suggest an underappreciated role for TF polymerization in modulating gene expression.

Roles of Tubulin Concentration during Prometaphase and Ran-GTP during Anaphase of Caenorhabditis elegans Meiosis.

In many animal species, the oocyte meiotic spindle, which is required for chromosome segregation, forms without centrosomes. In some systems, Ran-GEF on chromatin initiates spindle assembly. We found that in Caenorhabditis elegans oocytes, endogenously-tagged Ran-GEF dissociates from chromatin during spindle assembly but re-associates during meiotic anaphase. Meiotic spindle assembly occurred after auxin-induced degradation of Ran-GEF, but anaphase I was faster than controls and extrusion of the first polar body frequently failed. In search of a possible alternative pathway for spindle assembly, we found that soluble tubulin concentrates in the nuclear volume during germinal vesicle breakdown. We found that the concentration of soluble tubulin in the metaphase spindle region is enclosed by ER sheets which exclude cytoplasmic organelles including mitochondria and yolk granules. Measurement of the volume occupied by yolk granules and mitochondria indicated that volume exclusion would be sufficient to explain the concentration of tubulin in the spindle volume. We suggest that this concentration of soluble tubulin may be a redundant mechanism promoting spindle assembly near chromosomes.

CWAS-Plus: estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data.

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Evolutionary adaptation of an HP1-protein chromodomain integrates chromatin and DNA sequence signals.

Members of the diverse heterochromatin protein 1 (HP1) family play crucial roles in heterochromatin formation and maintenance. Despite the similar affinities of their chromodomains for di- and tri-methylated histone H3 lysine 9 (H3K9me2/3), different HP1 proteins exhibit distinct chromatin-binding patterns, likely due to interactions with various specificity factors. Previously, we showed that the chromatin-binding pattern of the HP1 protein Rhino, a crucial factor of the Drosophila PIWI-interacting RNA (piRNA) pathway, is largely defined by a DNA sequence-specific C2H2 zinc finger protein named Kipferl (Baumgartner et al., 2022). Here, we elucidate the molecular basis of the interaction between Rhino and its guidance factor Kipferl. Through phylogenetic analyses, structure prediction, and in vivo genetics, we identify a single amino acid change within Rhino's chromodomain, G31D, that does not affect H3K9me2/3 binding but disrupts the interaction between Rhino and Kipferl. Flies carrying the rhinoG31D mutation phenocopy kipferl mutant flies, with Rhino redistributing from piRNA clusters to satellite repeats, causing pronounced changes in the ovarian piRNA profile of rhinoG31D flies. Thus, Rhino's chromodomain functions as a dual-specificity module, facilitating interactions with both a histone mark and a DNA-binding protein.

CloudATAC: a cloud-based framework for ATAC-Seq data analysis.

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) generates genome-wide chromatin accessibility profiles, providing valuable insights into epigenetic gene regulation at both pooled-cell and single-cell population levels. Comprehensive analysis of ATAC-seq data involves the use of various interdependent programs. Learning the correct sequence of steps needed to process the data can represent a major hurdle. Selecting appropriate parameters at each stage, including pre-analysis, core analysis, and advanced downstream analysis, is important to ensure accurate analysis and interpretation of ATAC-seq data. Additionally, obtaining and working within a limited computational environment presents a significant challenge to non-bioinformatic researchers. Therefore, we present Cloud ATAC, an open-source, cloud-based interactive framework with a scalable, flexible, and streamlined analysis framework based on the best practices approach for pooled-cell and single-cell ATAC-seq data. These frameworks use on-demand computational power and memory, scalability, and a secure and compliant environment provided by the Google Cloud. Additionally, we leverage Jupyter Notebook's interactive computing platform that combines live code, tutorials, narrative text, flashcards, quizzes, and custom visualizations to enhance learning and analysis. Further, leveraging GPU instances has significantly improved the run-time of the single-cell framework. The source codes and data are publicly available through NIH Cloud lab https://github.com/NIGMS/ATAC-Seq-and-Single-Cell-ATAC-Seq-Analysis. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

The Upstream Sequence Transcription Complex dictates nucleosome positioning and promoter accessibility at piRNA genes in the C. elegans germ line.

The piRNA pathway is a conserved germline-specific small RNA pathway that ensures genomic integrity and continued fertility. In C. elegans and other nematodes, Type-I piRNAs are expressed from >10,000 independently transcribed genes clustered within two discrete domains of 1.5 and 3.5 MB on Chromosome IV. Clustering of piRNA genes contributes to their germline-specific expression, but the underlying mechanisms are unclear. We analyze isolated germ nuclei to demonstrate that the piRNA genomic domains are located in a heterochromatin-like environment. USTC (Upstream Sequence Transcription Complex) promotes strong association of nucleosomes throughout piRNA clusters, yet organizes the local nucleosome environment to direct the exposure of individual piRNA genes. Localization of USTC to the piRNA domains depends upon the ATPase chromatin remodeler ISW-1, which maintains high nucleosome density across piRNA clusters and ongoing production of piRNA precursors. Overall, this work provides insight into how chromatin states coordinate transcriptional regulation over large genomic domains, with implications for global genome organization.

Role of diffusion and reaction of the constituents in spreading of histone modification marks.

Cells switch genes ON or OFF by altering the state of chromatin via histone modifications at specific regulatory locations along the chromatin polymer. These gene regulation processes are carried out by a network of reactions in which the histone marks spread to neighboring regions with the help of enzymes. In the literature, this spreading has been studied as a purely kinetic, non-diffusive process considering the interactions between neighboring nucleosomes. In this work, we go beyond this framework and study the spreading of modifications using a reaction-diffusion (RD) model accounting for the diffusion of the constituents. We quantitatively segregate the modification profiles generated from kinetic and RD models. The diffusion and degradation of enzymes set a natural length scale for limiting the domain size of modification spreading, and the resulting enzyme limitation is inherent in our model. We also demonstrate the emergence of confined modification domains without the explicit requirement of a nucleation site. We explore polymer compaction effects on spreading and show that single-cell domains may differ from averaged profiles. We find that the modification profiles from our model are comparable with existing H3K9me3 data of S. pombe.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample.

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

CUT&Tag applied to zebrafish adult tail fins reveals a return of embryonic H3K4me3 patterns during regeneration.

Regenerative potential is governed by a complex process of transcriptional reprogramming, involving chromatin reorganization and dynamics in transcription factor binding patterns throughout the genome. The degree to which chromatin and epigenetic changes contribute to this process remains only partially understood. Here we provide a modified CUT&Tag protocol suitable for improved characterization and interrogation of changes in chromatin modifications during adult fin regeneration in zebrafish. Our protocol generates data that recapitulates results from previously published ChIP-Seq methods, requires far fewer cells as input, and significantly improves signal to noise ratios. We deliver high-resolution enrichment maps for H3K4me3 of uninjured and regenerating fin tissues. During regeneration, we find that H3K4me3 levels increase over gene promoters which become transcriptionally active and genes which lose H3K4me3 become silenced. Interestingly, these reprogramming events recapitulate the H3K4me3 patterns observed in developing fin folds of 24-h old zebrafish embryos. Our results indicate that changes in genomic H3K4me3 patterns during fin regeneration occur in a manner consistent with reactivation of developmental programs, demonstrating CUT&Tag to be an effective tool for profiling chromatin landscapes in regenerating tissues.

Determinants of p53 DNA binding, gene regulation, and cell fate decisions.

The extent to which transcription factors read and respond to specific information content within short DNA sequences remains an important question that the tumor suppressor p53 is helping us answer. We discuss recent insights into how local information content at p53 binding sites might control modes of p53 target gene activation and cell fate decisions. Significant prior work has yielded data supporting two potential models of how p53 determines cell fate through its target genes: a selective target gene binding and activation model and a p53 level threshold model. Both of these models largely revolve around an analogy of whether p53 is acting in a "smart" or "dumb" manner. Here, we synthesize recent and past studies on p53 decoding of DNA sequence, chromatin context, and cellular signaling cascades to elicit variable cell fates critical in human development, homeostasis, and disease.

The formation of aggregated chromatin/chromosomes in mouse oocytes treated with high concentration of IBMX as a model for a chromosome transfer in human.

The presence of cyclic adenosine monophosphate (cAMP) has been considered to be a fundamental factor in ensuring meiotic arrest prior to ovulation. cAMP is regarded as a key molecule in the regulation of oocyte maturation. However, it has been reported that increased levels of intracellular cAMP can result in abnormal cytokinesis, with some MI oocytes leading to symmetrically cleaved 2-cell MII oocytes. Consequently, we aimed to investigate the effects of elevated intracellular cAMP levels on abnormal cytokinesis and oocyte maturation during the meiosis of mouse oocytes. This study found that a high concentration of isobutylmethylxanthine (IBMX) also caused chromatin/chromosomes aggregation (AC) after the first meiosis. The rates of AC increased the greater the concentration of IBMX. In addition, AC formation was found to be reversible, showing that the re-formation of the spindle chromosome complex was possible after the IBMX was removed. In human oocytes, the chromosomes aggregate after the germinal vesicle breakdown and following the first and second polar body extrusions (the AC phase), while mouse oocytes do not have this AC phase. The results of our current study may indicate that the AC phase in human oocytes could be related to elevated levels of intracytoplasmic cAMP.