Chromatin Profiling Research Articles

Histone lactylation mediated by Fam172a in POMC neurons regulates energy balance

Glycolysis-derived lactate was identified as substrate for histone lactylation, which has been regarded as a significant role in transcriptional regulation in many tissues. However, the role of histone lactylation in the metabolic center, the hypothalamus, is still unknown. Here, we show that hypothalamic pro-opiomelanocortin (POMC) neuron-specific deletion of family with sequence similarity 172, member A (Fam172a) can increase histone lactylation and protect mice against diet-induced obesity (DIO) and related metabolic disorders. Conversely, overexpression of Fam172a in POMC neurons led to an obesity-like phenotype. Using RNA-seq and CUT&Tag chromatin profiling analyzes, we find that knockdown of Fam172a activates the glycolytic process and increases peptidylglycine α-amidating monooxygenase (PAM), which affects the synthesis of α-MSH, via H4K12la (histone lactylation). In addition, pharmacological inhibition of lactate production clearly abrogates the anti-obesity effect of PFKO (POMC-Cre, Fam172aloxP/loxP, POMC neurons Fam172a knockout). These findings highlight the importance of Fam172a and lactate in the development of obesity and our results mainly concern male mice.

Turnover of retroelements and satellite DNA drives centromere reorganization over short evolutionary timescales in Drosophila.

Centromeres reside in rapidly evolving, repeat-rich genomic regions, despite their essential function in chromosome segregation. Across organisms, centromeres are rich in selfish genetic elements such as transposable elements and satellite DNAs that can bias their transmission through meiosis. However, these elements still need to cooperate at some level and contribute to, or avoid interfering with, centromere function. To gain insight into the balance between conflict and cooperation at centromeric DNA, we take advantage of the close evolutionary relationships within the Drosophila simulans clade-D. simulans, D. sechellia, and D. mauritiana-and their relative, D. melanogaster. Using chromatin profiling combined with high-resolution fluorescence in situ hybridization on stretched chromatin fibers, we characterize all centromeres across these species. We discovered dramatic centromere reorganization involving recurrent shifts between retroelements and satellite DNAs over short evolutionary timescales. We also reveal the recent origin (<240 Kya) of telocentric chromosomes in D. sechellia, where the X and fourth centromeres now sit on telomere-specific retroelements. Finally, the Y chromosome centromeres, which are the only chromosomes that do not experience female meiosis, do not show dynamic cycling between satDNA and TEs. The patterns of rapid centromere turnover in these species are consistent with genetic conflicts in the female germline and have implications for centromeric DNA function and karyotype evolution. Regardless of the evolutionary forces driving this turnover, the rapid reorganization of centromeric sequences over short evolutionary timescales highlights their potential as hotspots for evolutionary innovation.

KAS-ATAC reveals the genome-wide single-stranded accessible chromatin landscape of the human genome.

Gene regulation in most eukaryotes involves two fundamental physical processes -- alterations in the packaging of the genome by nucleosomes, with active cis-regulatory elements (CREs) generally characterized by an open-chromatin configuration, and the activation of transcription. Mapping these physical properties and biochemical activities genome-wide -- through profiling chromatin accessibility and active transcription -- are key tools used to understand the logic and mechanisms of transcription and its regulation. However, the relationship between these two states has until now not been accessible to simultaneous measurement. To address this, we developed KAS-ATAC, a combination of the KAS-seq (Kethoxal-Assisted SsDNA sequencing and ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) methods for mapping single-stranded DNA (and thus active transcription) and chromatin accessibility, respectively, enabling the genome-wide identification of DNA fragments that are simultaneously accessible and contain ssDNA. We use KAS-ATAC to evaluate levels of active transcription over different classes of regulatory elements in the human genome, to estimate the absolute levels of transcribed accessible DNA over CREs, to map the nucleosomal configurations associated with RNA polymerase activities, and to assess transcription factor association with transcribed DNA through transcription factor binding site (TFBS) footprinting. We observe lower levels of transcription over distal enhancers compared to promoters and distinct nucleosomal configurations around transcription initiation sites associated with active transcription. Most TFs associate equally with transcribed and nontranscribed DNA but a few factors specifically do not exhibit footprints over ssDNA-containing fragments. We anticipate KAS-ATAC to continue to derive useful insights into chromatin organization and transcriptional regulation in other contexts in the future.

Serinc5 Regulates Sequential Chondrocyte Differentiation by Inhibiting Sox9 Function in Pre-Hypertrophic Chondrocytes.

The growth plate is the primary site of longitudinal bone growth with chondrocytes playing a pivotal role in endochondral bone development. Chondrocytes undergo a series of differentiation steps, resulting in the formation of a unique hierarchical columnar structure comprising round, proliferating, pre-hypertrophic, and hypertrophic chondrocytes. Pre-hypertrophic chondrocytes, which exist in the transitional stage between proliferating and hypertrophic stages, are a critical cell population in the growth plate. However, the molecular basis of pre-hypertrophic chondrocytes remains largely undefined. Here, we employed scRNA-seq analysis on fluorescently labeled growth plate chondrocytes for their molecular characterization. Serine incorporator 5 (Serinc5) was identified as a marker gene for pre-hypertrophic chondrocytes. Histological analysis revealed that Serinc5 is specifically expressed in pre-hypertrophic chondrocytes, overlapping with Indian hedgehog (Ihh). Serinc5 represses cell proliferation and Col2a1 and Acan expression by inhibiting the transcriptional activity of Sox9 in primary chondrocytes. Chromatin profiling using ChIP-seq and ATAC-seq revealed an active enhancer of Serinc5 located in intron 1, with its chromatin status progressively activated during chondrocyte differentiation. Collectively, our findings suggest that Serinc5 regulates sequential chondrocyte differentiation from proliferation to hypertrophy by inhibiting Sox9 function in pre-hypertrophic chondrocytes, providing novel insights into the mechanisms underlying chondrocyte differentiation in growth plates.

Expression and Purification of pAG-MNase for CUT&RUN.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a chromatin profiling strategy that releases specific DNA-protein complexes from in situ chromatin to the supernatant. CUT&RUN employs the hybrid protein pAG-MNase, composed of IgG binding proteins A and G, followed by micrococcal nuclease. Here, I describe the expression and purification of recombinant pAG-MNase as well as the assessment of its enzymatic activity. Purified pAG-MNase is hyperactive and enables large-scale analyses.

BPDCN MYB fusions regulate cell cycle genes, impair differentiation and induce myeloid-dendritic cell leukemia.

MYB fusions are recurrently found in select cancers, including blastic plasmacytoid dendritic cell neoplasm (BPDCN), an acute leukemia with poor prognosis. They are markedly enriched in BPDCN compared to other blood cancers, and in some patients are the only obvious somatic mutation detected. This suggests they may alone be sufficient to drive dendritic cell transformation. MYB fusions are hypothesized to alter the normal transcription factor activity of MYB, but mechanistically how they promote leukemogenesis is poorly understood. Using CUT&RUN chromatin profiling, we found that in BPDCN leukemogenesis, MYB switches from being a regulator of dendritic cell lineage genes to aberrantly regulating G2/M cell cycle control genes. MYB fusions found in BPDCN patients increased the magnitude of DNA binding at these locations, and this was linked to BPDCN-associated gene expression changes. Furthermore, expression of MYB fusions in vivo impaired dendritic cell differentiation and induced transformation to generate a mouse model of myeloid-dendritic acute leukemia. Therapeutically, we present evidence that all-trans retinoic acid (ATRA) may cause loss of MYB protein and cell death in BPDCN.

Origin of Ewing sarcoma by embryonic reprogramming of neural crest to mesoderm.

Ewing sarcoma is a malignant small round blue cell tumor of bones and soft tissues caused by chromosomal translocations that generate aberrant fusion oncogenes, most frequently EWSR1::FLI1. The cell of origin and mechanisms of EWSR1::FLI1-driven transformation have remained unresolved, largely due to lack of a representative animal model. By developing a zebrafish Ewing sarcoma model, we provide evidence for a neural crest origin of this cancer. Neural crest-derived cells uniquely tolerate expression of EWSR1::FLI1 and targeted expression of EWSR1::FLI1 in these cells generates Ewing sarcomas. Single-cell analysis of tumor initiation shows that EWSR1::FLI1 reprograms neural crest-derived cells to a mesoderm-like state, strikingly resulting in ectopic fins throughout the body. By profiling chromatin accessibility and genome-wide EWSR1::FLI1 binding, we find that the fusion oncogene hijacks developmental enhancers for neural crest to mesoderm reprogramming during cancer initiation. These findings show how a single mutation profoundly alters embryonic cell fate decisions to initiate a devastating childhood cancer.

Overloading And unpacKing (OAK) - droplet-based combinatorial indexing for ultra-high throughput single-cell multiomic profiling

Multiomic profiling of single cells by sequencing is a powerful technique for investigating cellular diversity. Existing droplet-based microfluidic methods produce many cell-free droplets, underutilizing bead barcodes and reagents. Combinatorial indexing on microplates is more efficient for barcoding but labor-intensive. Here we present Overloading And unpacKing (OAK), which uses a droplet-based barcoding system for initial compartmentalization followed by a second aliquoting round to achieve combinatorial indexing. We demonstrate OAK’s versatility with single-cell RNA sequencing as well as paired single-nucleus RNA sequencing and accessible chromatin profiling. We further showcase OAK’s performance on complex samples, including differentiated bronchial epithelial cells and primary retinal tissue. Finally, we examine transcriptomic responses of over 400,000 melanoma cells to a RAF inhibitor, belvarafenib, discovering a rare resistant cell population (0.12%). OAK’s ultra-high throughput, broad compatibility, high sensitivity, and simplified procedures make it a powerful tool for large-scale molecular analysis, even for rare cells.

Single-cell multimodal chromatin profiles revealing epigenetic regulations of cells in hepatocellular carcinoma.

AbstractHepatocellular carcinoma (HCC), with its increasing prevalence globally, is emerging as a major health challenge. Hepatocellular carcinoma cells exhibit both significant homogeneity and heterogeneity, governed by complex epigenomic regulations. Recently, numerous single‐cell unimodal techniques have been used to study HCC at various levelswhich have already revealed a complex interplay at genomic, transcriptomic, spatial and epigenomic levels. However, the use of single‐cell multimodal techniques combining different unimodal layers in HCC remains quite limited, necessitating further studies focusing on these methods to uncover novel markers, and mechanisms at epigenetic levels. In this commentary, we highlight how integrating multimodal approaches with epigenetic modifications can provide new insights into HCC and foster future therapeutic advancements.

Single-cell transcriptomes and chromatin accessibility of endothelial cells unravel transcription factors associated with dysregulated angiogenesis in systemic sclerosis

ObjectivesVasculopathy emerges early in systemic sclerosis (SSc) and links to endothelial cell (EC) injury and angiogenesis. Understanding EC transcriptomes and epigenomes is crucial for unravelling the mechanisms involved.MethodsTranscriptomes and chromatin...

Transcription factor networks disproportionately enrich for heritability of blood cell phenotypes.

Most phenotype-associated genetic variants map to non-coding regulatory regions of the human genome. Moreover, variants associated with blood cell phenotypes are enriched in regulatory regions active during hematopoiesis. To systematically explore the nature of these regions, we developed a highly efficient strategy, Perturb-multiome, that makes it possible to simultaneously profile both chromatin accessibility and gene expression in single cells with CRISPR-mediated perturbation of a range of master transcription factors (TFs). This approach allowed us to examine the connection between TFs, accessible regions, and gene expression across the genome throughout hematopoietic differentiation. We discovered that variants within the TF-sensitive accessible chromatin regions, while representing less than 0.3% of the genome, show a ~100-fold enrichment in heritability across certain blood cell phenotypes; this enrichment is strikingly higher than for other accessible chromatin regions. Our approach facilitates large-scale mechanistic understanding of phenotype-associated genetic variants by connecting key cis-regulatory elements and their target genes within gene regulatory networks.

Genotype inference from aggregated chromatin accessibility data reveals genetic regulatory mechanisms.

Understanding the genetic causes for variability in chromatin accessibility can shed light on the molecular mechanisms through which genetic variants may affect complex traits. Thousands of ATAC-seq samples have been collected that hold information about chromatin accessibility across diverse cell types and contexts, but most of these are not paired with genetic information and come from diverse distinct projects and laboratories. We report here joint genotyping, chromatin accessibility peak calling, and discovery of quantitative trait loci which influence chromatin accessibility (caQTLs), demonstrating the capability of performing caQTL analysis on a large scale in a diverse sample set without pre-existing genotype information. Using 10,293 profiling samples representing 1,454 unique donor individuals across 653 studies from public databases, we catalog 23,381 caQTLs in total. After joint discovery analysis, we cluster samples based on accessible chromatin profiles to identify context-specific caQTLs. We find that caQTLs are strongly enriched for annotations of gene regulatory elements across diverse cell types and tissues and are often strongly linked with genetic variation associated with changes in expression (eQTLs), indicating that caQTLs can mediate genetic effects on gene expression. We demonstrate sharing of causal variants for chromatin accessibility and diverse complex human traits, enabling a more complete picture of the genetic mechanisms underlying complex human phenotypes. Our work provides a proof of principle for caQTL calling from previously ungenotyped samples, and represents one of the largest, most diverse caQTL resources currently available, informing mechanisms of genetic regulation of gene expression and contribution to disease.

Single-cell multi-modal chromatin profiles revealing epigenetic regulations of cells in hepatocellular carcinoma.

Various epigenetic regulations systematically govern gene expression in cells involving various biological processes. Dysregulation of the epigenome leads to aberrant transcriptional programs and subsequently results in diseases, such as cancer. Therefore, comprehensive profiling epigenomics is essential for exploring the mechanisms underlying gene expression regulation during development and disease. In this study, we developed single-cell chromatin proteins and accessibility tagmentation (scCPA-Tag), a multi-modal single-cell epigenetic profile capturing technique based on barcoded Tn5 transposases and a droplet microfluidics platform. scCPA-Tag enables the simultaneous capture of DNA profiles of histone modification and chromatin accessibility in the same cell. By applying scCPA-Tag to K562 cells and a hepatocellular carcinoma (HCC) sample, we found that the silence of several chromatin-accessible genes can be attributed to lysine-27-trimethylation of the histone H3 tail (H3K27me3) modification. We characterized the epigenetic features of the tumour cells and different immune cell types in the HCC tumour tissue by scCPA-Tag. Besides, a tumour cell subtype (C2) with more aggressive features was identified and characterized by high chromatin accessibility and a lower abundance of H3K27me3 on tumour-promoting genes. Our multi-modal scCPA-Tag provides a comprehensive approach for exploring the epigenetic landscapes of heterogeneous cell types and revealing the mechanisms of gene expression regulation during developmental and pathological processes at the single-cell level. scCPA-Tag offers a highly efficient and high throughput technique to simultaneously profile histone modification and chromatin accessibility within a single cell. scCPA-Tag enables to uncover multiple epigenetic modification features of cellular compositions within tumor tissues. scCPA-Tag facilitates the exploration of the epigenetic landscapes of heterogeneous cell types and provides the mechanisms governing gene expression regulation.

Progress in multifactorial single-cell chromatin profiling methods.

Chromatin states play a key role in shaping overall cellular states and fates. Building a complete picture of the functional state of chromatin in cells requires the co-detection of several distinct biochemical aspects. These span DNA methylation, chromatin accessibility, chromosomal conformation, histone posttranslational modifications, and more. While this certainly presents a challenging task, over the past few years many new and creative methods have been developed that now enable co-assay of these different aspects of chromatin at single cell resolution. This field is entering an exciting phase, where a confluence of technological improvements, decreased sequencing costs, and computational innovation are presenting new opportunities to dissect the diversity of chromatin states present in tissues, and how these states may influence gene regulation. In this review, I discuss the spectrum of current experimental approaches for multifactorial chromatin profiling, highlight some of the experimental and analytical challenges, as well as some areas for further innovation.

Single cell expression and chromatin accessibility of the Toxoplasma gondii lytic cycle identifies AP2XII-8 as an essential ribosome regulon driver

Sequential lytic cycles driven by cascading transcriptional waves underlie pathogenesis in the apicomplexan parasite Toxoplasma gondii. This parasite’s unique division by internal budding, short cell cycle, and jumbled up classically defined cell cycle stages have restrained in-depth transcriptional program analysis. Here, unbiased transcriptome and chromatin accessibility maps throughout the lytic cell cycle are established at the single-cell level. Correlated pseudo-timeline assemblies of expression and chromatin profiles maps transcriptional versus chromatin level transition points promoting the cell division cycle. Sequential clustering analysis identifies functionally related gene groups promoting cell cycle progression. Promoter DNA motif mapping reveals patterns of combinatorial regulation. Pseudo-time trajectory analysis reveals transcriptional bursts at different cell cycle points. The dominant burst in G1 is driven largely by transcription factor AP2XII-8, which engages a conserved DNA motif, and promotes the expression of 44 ribosomal proteins encoding regulon. Overall, the study provides integrated, multi-level insights into apicomplexan transcriptional regulation.

Epigenetic insights into GABAergic development in Dravet Syndrome iPSC and therapeutic implications

Dravet syndrome (DS) is a devastating early-onset refractory epilepsy syndrome caused by variants in the SCN1A gene. A disturbed GABAergic interneuron function is implicated in the progression to DS but the underlying developmental and pathophysiological mechanisms remain elusive, in particularly at the chromatin level. Induced pluripotent stem cells (iPSCs) derived from DS cases and healthy donors were used to model disease-associated epigenetic abnormalities of GABAergic development. Chromatin accessibility was assessed at multiple time points (Day 0, Day 19, Day 35, and Day 65) of GABAergic differentiation. Additionally, the effects of the commonly used anti-seizure drug valproic acid (VPA) on chromatin accessibility were elucidated in GABAergic cells. The distinct dynamics in the chromatin profile of DS iPSC predicted accelerated early GABAergic development, evident at D19, and diverged further from the pattern in control iPSC with continued differentiation, indicating a disrupted GABAergic maturation. Exposure to VPA at D65 reshaped the chromatin landscape at a variable extent in different iPSC-lines and rescued the observed dysfunctional development of some DS iPSC-GABA. The comprehensive investigation on the chromatin landscape of GABAergic differentiation in DS-patient iPSC offers valuable insights into the epigenetic dysregulations associated with interneuronal dysfunction in DS. Moreover, the detailed analysis of the chromatin changes induced by VPA in iPSC-GABA holds the potential to improve the development of personalized and targeted anti-epileptic therapies.

Epigenetic insights into GABAergic development in Dravet Syndrome iPSC and therapeutic implications.

Dravet syndrome (DS) is a devastating early-onset refractory epilepsy syndrome caused by variants in the SCN1A gene. A disturbed GABAergic interneuron function is implicated in the progression to DS but the underlying developmental and pathophysiological mechanisms remain elusive, in particularly at the chromatin level. Induced pluripotent stem cells (iPSCs) derived from DS cases and healthy donors were used to model disease-associated epigenetic abnormalities of GABAergic development. Chromatin accessibility was assessed at multiple time points (Day 0, Day 19, Day 35, and Day 65) of GABAergic differentiation. Additionally, the effects of the commonly used anti-seizure drug valproic acid (VPA) on chromatin accessibility were elucidated in GABAergic cells. The distinct dynamics in the chromatin profile of DS iPSC predicted accelerated early GABAergic development, evident at D19, and diverged further from the pattern in control iPSC with continued differentiation, indicating a disrupted GABAergic maturation. Exposure to VPA at D65 reshaped the chromatin landscape at a variable extent in different iPSC-lines and rescued the observed dysfunctional development of some DS iPSC-GABA. The comprehensive investigation on the chromatin landscape of GABAergic differentiation in DS-patient iPSC offers valuable insights into the epigenetic dysregulations associated with interneuronal dysfunction in DS. Moreover, the detailed analysis of the chromatin changes induced by VPA in iPSC-GABA holds the potential to improve the development of personalized and targeted anti-epileptic therapies.

AChIP is an efficient and sensitive ChIP-seq technique for economically important plant organs.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is crucial for profiling histone modifications and transcription factor binding throughout the genome. However, its application in economically important plant organs (EIPOs) such as seeds, fruits and flowers is challenging due to their sturdy cell walls and complex constituents. Here we present advanced ChIP (aChIP), an optimized method that efficiently isolates chromatin from plant tissues while simultaneously removing cell walls and cellular constituents. aChIP precisely profiles histone modifications in all 14 tested EIPOs and identifies transcription factor and chromatin-modifying enzyme binding sites. In addition, aChIP enhances ChIP efficiency, revealing numerous novel modified sites compared with previous methods in vegetative tissues. aChIP reveals the histone modification landscape for rapeseed dry seeds, highlighting the intricate roles of chromatin dynamics during seed dormancy and germination. Altogether, aChIP is a powerful, efficient and sensitive approach for comprehensive chromatin profiling in virtually all plant tissues, especially in EIPOs.

AChIP for comprehensive chromatin profiling in economically important plant organs.

aChIP for comprehensive chromatin profiling in economically important plant organs.

Epigenetic heterogeneity shapes the transcriptional landscape of regional microglia.

Microglia, the innate immune cells in the central nervous system, exhibit distinct transcriptional profiles across brain regions that are important for facilitating their specialized function. There has been recent interest in identifying the epigenetic modifications associated with these distinct transcriptional profiles, as these may improve our understanding of the underlying mechanisms governing the functional specialization of microglia. One obstacle to achieving this goal is the large number of microglia required to obtain a genome-wide profile for a single histone modification. Given the cellular and regional heterogeneity of the brain, this would require pooling many samples which would impede biological applications that are limited by numbers of available animals. To overcome this obstacle, we have adapted a method of chromatin profiling known as Cleavage Under Targets and Tagmentation (CUT&Tag-Direct) to profile histone modifications associated with regional differences in gene expression throughout the brain reward system. Consistent with previous studies, we find that transcriptional profiles of microglia vary by brain region. However, here we report that these regional differences also exhibit transcriptional network signatures specific to each region. Additionally, we find that these region-dependent network signatures are associated with differential deposition of H3K27ac and H3K7me3, and while the H3K27me3 landscape is remarkably stable across brain regions, the H3K27ac landscape is most consistent with the anatomical location of microglia which explain their distinct transcriptional profiles. Altogether, these findings underscore the established role of H3K27me3 in cell fate determination and support the active role of H3K27ac in the dynamic regulation of microglial gene expression. In this study, we report a molecular and computational framework that can be applied to improve our understanding of the role of epigenetic regulation in microglia in both health and disease, using as few as 2,500 cells per histone mark.