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.

SMARCB1/INI1-deficient undifferentiated tumour of the thorax: a case report and review of the literature.

The 5th WHO classification of thoracic tumours includes thoracic SMARCA4-deficient undifferentiated tumour (SMARCA4-UT) among the "other epithelial tumours of the lung" chapter. Herein, we present a case of undifferentiated thoracic neoplasm with retention of SMARCA4 expression, lack of NUT fusion protein and loss of SMARCB1/INI1 expression. After presenting the clinical and pathological features of the tumour, we carried out a review of the literature on the same topic. Albeit very rare, we believe this entity should be included in the heterogeneous group of undifferentiated neoplasms of the thorax.

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.

TFscope: systematic analysis of the sequence features involved in the binding preferences of transcription factors.

Characterizing the binding preferences of transcription factors (TFs) in different cell types and conditions is key to understand how they orchestrate gene expression. Here, we develop TFscope, a machine learning approach that identifies sequence features explaining the binding differences observed between two ChIP-seq experiments targeting either the same TF in two conditions or two TFs with similar motifs (paralogous TFs). TFscope systematically investigates differences in the core motif, nucleotide environment and co-factor motifs, and provides the contribution of each key feature in the two experiments. TFscope was applied to > 305 ChIP-seq pairs, and several examples are discussed.

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.

ARID1A safeguards the canalization of the cell fate decision during osteoclastogenesis.

Chromatin remodeler ARID1A regulates gene transcription by modulating nucleosome positioning and chromatin accessibility. While ARID1A-mediated stage and lineage-restricted gene regulation during cell fate canalization remains unresolved. Using osteoclastogenesis as a model, we show that ARID1A transcriptionally safeguards the osteoclast (OC) fate canalization during proliferation-differentiation switching at single-cell resolution. Notably, ARID1A is indispensable for the transcriptional apparatus condensates formation with coactivator BRD4/lineage-specifying transcription factor (TF) PU.1 at Nfatc1 super-enhancer during safeguarding the OC fate canalization. Besides, the antagonist function between ARID1A-cBAF and BRD9-ncBAF complex during osteoclastogenesis has been validated with in vitro assay and compound mutant mouse model. Furthermore, the antagonistic function of ARID1A-"accelerator" and BRD9-"brake" both depend on coactivator BRD4-"clutch" during osteoclastogenesis. Overall, these results uncover sophisticated cooperation between chromatin remodeler ARID1A, coactivator, and lineage-specifying TF at super-enhancer of lineage master TF in a condensate manner, and antagonist between distinct BAF complexes in the proper and balanced cell fate canalization.

Hepatitis B Virus Covalently Closed Circular DNA Chromatin Immunoprecipitation Assay.

Hepatitis B virus (HBV) is an obligate human hepatotropic DNA virus causing both transient and chronic infection. The livers of chronic hepatitis B patients have a high risk of developing liver fibrosis, cirrhosis, and hepatocellular carcinoma. The nuclear episomal viral DNA intermediate, covalently closed circular DNA (cccDNA), forms a highly stable complex with host and viral proteins to serve as a transcription template and support HBV infection chronicity. Thus, characterization of the composition and dynamics of cccDNA nucleoprotein complexes providing cccDNA stability and gene regulation is of high importance for both basic and medical research. The presented method for chromatin immunoprecipitation coupled with qPCR (ChIP-qPCR) allows to assess provisional physical interaction of the protein of interest (POI) with cccDNA using POI-specific antibody, the level of enrichment of a POI on cccDNA versus control/background is characterized quantitatively using qPCR.

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.

Single gene analysis in yeast suggests nonequilibrium regulatory dynamics for transcription.

Fluctuations in the initiation rate of transcription, the first step in gene expression, ensue from the stochastic behavior of the molecular process that controls transcription. In steady state, the regulatory process is often assumed to operate reversibly, i.e., in equilibrium. However, reversibility imposes fundamental limits to information processing. For instance, the assumption of equilibrium is difficult to square with the precision with which the regulatory process executes its task in eukaryotes. Here we provide evidence - from microscopic analyses of the transcription dynamics at a single gene copy of yeast - that the regulatory process for transcription is cyclic and irreversible (out of equilibrium). The necessary coupling to reservoirs of free energy occurs via sequence-specific transcriptional activators and the recruitment, in part, of ATP-dependent chromatin remodelers. Our findings may help explain how eukaryotic cells reconcile the dual but opposing requirements for fast regulatory kinetics and high regulatory specificity.