CHD4 and SMYD1 repress common transcriptional programs in the developing heart

Regulation of chromatin states is essential for proper temporal and spatial gene expression. Chromatin states are modulated by remodeling complexes composed of components that have enzymatic activities. CHD4 is the catalytic core of the Nucleosome Remodeling and Deacetylase (NuRD) complex that represses gene transcription. However, it remains to be determined how CHD4, a ubiquitous enzyme that remodels chromatin structure, functions in cardiomyocytes to maintain heart development. Particularly, there exists controversy as to whether other proteins besides the NuRD components interact with CHD4 in the heart. Using quantitative proteomics, we identified that CHD4 interacts with SMYD1, a striated muscle-restricted histone methyltransferase that is essential for cardiomyocyte differentiation and cardiac morphogenesis. Comprehensive transcriptomic and chromatin accessibility studies of Smyd1 and Chd4 null embryonic hearts revealed that SMYD1 and CHD4 repress a group of common genes and pathways that included glycolysis, response to hypoxia, and angiogenesis. Our study reveals a novel mechanism by which CHD4 functions during heart development, and we have revealed an uncharacterized mechanism regarding how SMYD1 represses cardiac transcription in the developing heart.

scifi-ATAC-seq: massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called scifi-ATAC-seq, single-cell combinatorial fluidic indexing ATAC-sequencing, which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using the 10X Genomics platform. With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples can be indexed in a single emulsion reaction, representing an approximately 20-fold increase in throughput compared to the standard 10X Genomics workflow.

ATAC-clock: An aging clock based on chromatin accessibility

The establishment of aging clocks highlighted the strong link between changes in DNA methylation and aging. Yet, it is not known if other epigenetic features could be used to predict age accurately. Furthermore, previous studies have observed a lack of effect of age-related changes in DNA methylation on gene expression, putting the interpretability of DNA methylation-based aging clocks into question. In this study, we explore the use of chromatin accessibility to construct aging clocks. We collected blood from 159 human donors and generated chromatin accessibility, transcriptomic, and cell composition data. We investigated how chromatin accessibility changes during aging and constructed a novel aging clock with a median absolute error of 5.27 years. The changes in chromatin accessibility used by the clock were strongly related to transcriptomic alterations, aiding clock interpretation. We additionally show that our chromatin accessibility clock performs significantly better than a transcriptomic clock trained on matched samples. In conclusion, we demonstrate that the clock relies on cell-intrinsic chromatin accessibility alterations rather than changes in cell composition. Further, we present a new approach to construct epigenetic aging clocks based on chromatin accessibility, which bear a direct link to age-related transcriptional alterations, but which allow for more accurate age predictions than transcriptomic clocks.

Transcriptome and open chromatin analysis reveals the process of myocardial cell development and key pathogenic target proteins in Long QT syndrome type 7

OBJECTIVE

Long QT syndrome type 7 (Andersen-Tawil syndrome, ATS), which is caused by KCNJ2 gene mutation, often leads to ventricular arrhythmia, periodic paralysis and skeletal malformations. The development, differentiation and electrophysiological maturation of cardiomyocytes (CMs) changes promote the pathophysiology of Long QT syndrome type 7(LQT7). We aimed to specifically reproduce the ATS disease phenotype and study the pathogenic mechanism.

METHODS AND RESULTS

We established a cardiac cell model derived from human induced pluripotent stem cells (hiPSCs) to the phenotypes and electrophysiological function, and the establishment of a human myocardial cell model that specifically reproduces the symptoms of ATS provides a reliable platform for exploring the mechanism of this disease or potential drugs. The spontaneous pulsation rate of myocardial cells in the mutation group was significantly lower than that in the repair CRISPR group, the action potential duration was prolonged, and the Kir2.1 current of the inward rectifier potassium ion channel was decreased, which is consistent with the clinical symptoms of ATS patients. Only ZNF528, a chromatin-accessible TF related to pathogenicity, was continuously regulated beginning from the cardiac mesodermal precursor cell stage (day 4), and continued to be expressed at low levels, which was identified by WGCNA method and verified with ATAC-seq data in the mutation group. Subsequently, it indicated that seven pathways were downregulated (all p < 0.05) by used single sample Gene Set Enrichment Analysis to evaluate the overall regulation of potassium-related pathways enriched in the transcriptome and proteome of late mature CMs. Among them, the three pathways (GO: 0008076, GO: 1990573 and GO: 0030007) containing the mutated gene KCNJ2 is involved that are related to the whole process by which a potassium ion enters the cell via the inward rectifier potassium channel to exert its effect were inhibited. The other four pathways are related to regulation of the potassium transmembrane pathway and sodium:potassium exchange ATPase (p < 0.05). ZNF528 small interfering (si)-RNA was applied to hiPSC-derived cardiomyocytes for CRISPR group to explore changes in potassium ion currents and growth and development related target protein levels that affect disease phenotype. Three consistently downregulated proteins (KCNJ2, CTTN and ATP1B1) associated with pathogenicity were verificated through correlation and intersection analysis.

CONCLUSION

This study uncovers TFs and target proteins related to electrophysiology and developmental pathogenicity in ATS myocardial cells, obtaining novel targets for potential therapeutic candidate development that does not rely on gene editing.

txci-ATAC-seq: a massive-scale single-cell technique to profile chromatin accessibility

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16-/-) model. Our study of CC16-/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model.

Chromatin accessibility and H3K9me3 landscapes reveal long-term epigenetic effects of fetal-neonatal iron deficiency in rat hippocampus

BACKGROUND

Iron deficiency (ID) during the fetal-neonatal period results in long-term neurodevelopmental impairments associated with pervasive hippocampal gene dysregulation. Prenatal choline supplementation partially normalizes these effects, suggesting an interaction between iron and choline in hippocampal transcriptome regulation. To understand the regulatory mechanisms, we investigated epigenetic marks of genes with altered chromatin accessibility (ATAC-seq) or poised to be repressed (H3K9me3 ChIP-seq) in iron-repleted adult rats having experienced fetal-neonatal ID exposure with or without prenatal choline supplementation.

RESULTS

Fetal-neonatal ID was induced by limiting maternal iron intake from gestational day (G) 2 through postnatal day (P) 7. Half of the pregnant dams were given supplemental choline (5.0 g/kg) from G11-18. This resulted in 4 groups at P65 (Iron-sufficient [IS], Formerly Iron-deficient [FID], IS with choline [ISch], and FID with choline [FIDch]). Hippocampi were collected from P65 iron-repleted male offspring and analyzed for chromatin accessibility and H3K9me3 enrichment. 22% and 24% of differentially transcribed genes in FID- and FIDch-groups, respectively, exhibited significant differences in chromatin accessibility, whereas 1.7% and 13% exhibited significant differences in H3K9me3 enrichment. These changes mapped onto gene networks regulating synaptic plasticity, neuroinflammation, and reward circuits. Motif analysis of differentially modified genomic sites revealed significantly stronger choline effects than early-life ID and identified multiple epigenetically modified transcription factor binding sites.

CONCLUSIONS

This study reveals genome-wide, stable epigenetic changes and epigenetically modifiable gene networks associated with specific chromatin marks in the hippocampus, and lays a foundation to further elucidate iron-dependent epigenetic mechanisms that underlie the long-term effects of fetal-neonatal ID, choline, and their interactions.

From degraded to deciphered: ATAC-seq's application potential in forensic diagnosis

In recent years, molecular biology-based diagnostic techniques have made remarkable strides and are now extensively utilized in clinical practice, providing invaluable insights for disease diagnosis and treatment. However, forensic medicine, especially forensic pathology, has witnessed relatively limited progress in the application of molecular biology technologies. A significant challenge in employing molecular techniques for forensic diagnoses lies in the quantitative and qualitative changes observed in diagnostic markers due to sample degradation-a recognized and formidable obstacle. Inspired by the success of DNA sequencing in forensic practices, which enables accurate individual identification even in cases involving degraded and deteriorated tissues and organs, we propose the application of the assay for transposase-accessible chromatin with sequencing (ATAC-seq) to identify targets at the transcriptional onset, exploring chromatin and DNA-level alterations for injury and disease inference in forensic samples. This study employs ATAC-seq to explore alterations in chromatin accessibility post-injury and their subsequent changes over a 2-h degradation period, employing traumatic brain injury (TBI) as a representative model. Our findings reveal high sensitivity of chromatin accessibility sites to injury, evidenced by shifts in thousands of peak positions post-TBI. Remarkably, these alterations remain largely unaffected by early degradation. Our results robustly endorse the notion that integrating and incorporating these specific loci for injury and disease diagnosis in forensic samples holds tremendous promise for practical application. We further validated the above results using human cortical tissue, which supported that early degradation did not significantly affect chromatin accessibility. This pioneering advancement in molecular diagnostic techniques may revolutionize the field of forensic science, especially forensic pathology.

Integrative analysis of chromatin accessibility and transcriptome landscapes in the induction of peritoneal fibrosis by high glucose

BACKGROUND

Peritoneal fibrosis is the prevailing complication induced by prolonged exposure to high glucose in patients undergoing peritoneal dialysis.

METHODS

To elucidate the molecular mechanisms underlying this process, we conducted an integrated analysis of the transcriptome and chromatin accessibility profiles of human peritoneal mesothelial cells (HMrSV5) during high-glucose treatment.

RESULTS

Our study identified 2775 differentially expressed genes (DEGs) related to high glucose-triggered pathological changes, including 1164 upregulated and 1611 downregulated genes. Genome-wide DEGs and network analysis revealed enrichment in the epithelial-mesenchymal transition (EMT), inflammatory response, hypoxia, and TGF-beta pathways. The enriched genes included VEGFA, HIF-1α, TGF-β1, EGF, TWIST2, and SNAI2. Using ATAC-seq, we identified 942 hyper (higher ATAC-seq signal in high glucose-treated HMrSV5 cells than in control cells) and 714 hypo (lower ATAC-seq signal in high glucose-treated HMrSV5 cells versus control cells) peaks with differential accessibility in high glucose-treated HMrSV5 cells versus controls. These differentially accessible regions were positively correlated (R = 0.934) with the nearest DEGs. These genes were associated with 566 up- and 398 downregulated genes, including SNAI2, TGF-β1, HIF-1α, FGF2, VEGFA, and VEGFC, which are involved in critical pathways identified by transcriptome analysis. Integrated ATAC-seq and RNA-seq analysis also revealed key transcription factors (TFs), such as HIF-1α, ARNTL, ELF1, SMAD3 and XBP1. Importantly, we demonstrated that HIF-1α is involved in the regulation of several key genes associated with EMT and the TGF-beta pathway. Notably, we predicted and experimentally validated that HIF-1α can exacerbate the expression of TGF-β1 in a high glucose-dependent manner, revealing a novel role of HIF-1α in high glucose-induced pathological changes in human peritoneal mesothelial cells (HPMCs).

CONCLUSIONS

In summary, our study provides a comprehensive view of the role of transcriptome deregulation and chromosome accessibility alterations in high glucose-induced pathological fibrotic changes in HPMCs. This analysis identified hub genes, signaling pathways, and key transcription factors involved in peritoneal fibrosis and highlighted the novel glucose-dependent regulation of TGF-β1 by HIF-1α. This integrated approach has offered a deeper understanding of the pathogenesis of peritoneal fibrosis and has indicated potential therapeutic targets for intervention.

Anp32e protects against accumulation of H2A.Z at Sox motif containing promoters during zebrafish gastrulation

Epigenetic regulation of chromatin states is crucial for proper gene expression programs and progression during development, but precise mechanisms by which epigenetic factors influence differentiation remain poorly understood. Here we find that the histone variant H2A.Z accumulates at Sox motif-containing promoters during zebrafish gastrulation while neighboring genes become transcriptionally active. These changes coincide with reduced expression of anp32e, the H2A.Z histone removal chaperone, suggesting that loss of Anp32e may lead to increases in H2A.Z binding during differentiation. Remarkably, genetic removal of Anp32e in embryos leads to H2A.Z accumulation prior to gastrulation and developmental genes become precociously active. Accordingly, H2A.Z accumulation occurs most extensively at Sox motif-associated genes, including many which are normally activated following gastrulation. Altogether, our results provide compelling evidence for a mechanism in which Anp32e preferentially restricts H2A.Z accumulation at Sox motifs to regulate the initial phases of developmental differentiation in zebrafish.

Identification of Regulatory Elements in Primary Sensory Neurons Involved in Trauma-Induced Neuropathic Pain

Chronic pain is a significant public health issue that is often refractory to existing therapies. Here we use a multiomic approach to identify cis-regulatory elements that show differential chromatin accessibility and reveal transcription factor (TF) binding motifs with functional regulation in the rat dorsal root ganglion (DRG), which contain cell bodies of primary sensory neurons, after nerve injury. We integrated RNA-seq to understand how differential chromatin accessibility after nerve injury may influence gene expression. Using TF protein arrays and chromatin immunoprecipitation-qPCR, we confirmed C/EBPγ binding to a differentially accessible sequence and used RNA-seq to identify processes in which C/EBPγ plays an important role. Our findings offer insights into TF motifs that are associated with chronic pain. These data show how interactions between chromatin landscapes and TF expression patterns may work together to determine gene expression programs in rat DRG neurons after nerve injury.

Sall4 regulates posterior trunk mesoderm development by promoting mesodermal gene expression and repressing neural genes in the mesoderm

The trunk axial skeleton develops from paraxial mesoderm cells. Our recent study demonstrated that conditional knockout of the stem cell factor Sall4 in mice by TCre caused tail truncation and a disorganized axial skeleton posterior to the lumbar level. Based on this phenotype, we hypothesized that, in addition to the previously reported role of Sall4 in neuromesodermal progenitors, Sall4 is involved in the development of the paraxial mesoderm tissue. Analysis of gene expression and SALL4 binding suggests that Sall4 directly or indirectly regulates genes involved in presomitic mesoderm differentiation, somite formation and somite differentiation. Furthermore, ATAC-seq in TCre; Sall4 mutant posterior trunk mesoderm shows that Sall4 knockout reduces chromatin accessibility. We found that Sall4-dependent open chromatin status drives activation and repression of WNT signaling activators and repressors, respectively, to promote WNT signaling. Moreover, footprinting analysis of ATAC-seq data suggests that Sall4-dependent chromatin accessibility facilitates CTCF binding, which contributes to the repression of neural genes within the mesoderm. This study unveils multiple mechanisms by which Sall4 regulates paraxial mesoderm development by directing activation of mesodermal genes and repression of neural genes.

scAuto as a comprehensive framework for single-cell chromatin accessibility data analysis

Interpreting single-cell chromatin accessibility data is crucial for understanding intercellular heterogeneity regulation. Despite the progress in computational methods for analyzing this data, there is still a lack of a comprehensive analytical framework and a user-friendly online analysis tool. To fill this gap, we developed a pre-trained deep learning-based framework, single-cell auto-correlation transformers (scAuto), to overcome the challenge. Following DNABERT's methodology of pre-training and fine-tuning, scAuto learns a general understanding of DNA sequence's grammar by being pre-trained on unlabeled human genome via self-supervision; it is then transferred to the single-cell chromatin accessibility analysis task of scATAC-seq data for supervised fine-tuning. We extensively validated scAuto on the Buenrostro2018 dataset, demonstrating its superior performance on chromatin accessibility prediction, single-cell clustering, and data denoising. Based on scAuto, we further developed an interactive web server for single-cell chromatin accessibility data analysis. It integrates tutorial-style interfaces for those with limited programming skills. The platform is accessible at http://zhanglab.icaup.cn. To our knowledge, this work is expected to help analyze single-cell chromatin accessibility data and facilitate the development of precision medicine.

The dynamic landscape of chromatin accessibility and active regulatory elements in the mediobasal hypothalamus influences the seasonal activation of the reproductive axis in the male quail under long light exposure

BACKGROUND

In cold and temperate zones, seasonal reproduction plays a crucial role in the survival and reproductive success of species. The photoperiod influences reproductive processes in seasonal breeders through the hypothalamic-pituitary-gonadal (HPG) axis, in which the mediobasal hypothalamus (MBH) serves as the central region responsible for transmitting light information to the endocrine system. However, the cis-regulatory elements and the transcriptional activation mechanisms related to seasonal activation of the reproductive axis in MBH remain largely unclear. In this study, an artificial photoperiod program was used to induce the HPG axis activation in male quails, and we compared changes in chromatin accessibility changes during the seasonal activation of the HPG axis.

RESULTS

Alterations in chromatin accessibility occurred in the mediobasal hypothalamus (MBH) and stabilized at LD7 during the activation of the HPG axis. Most open chromatin regions (OCRs) are enriched mainly in introns and distal intergenic regions. The differentially accessible regions (DARs) showed enrichment of binding motifs of the RFX, NKX, and MEF family of transcription factors that gained-loss accessibility under long-day conditions, while the binding motifs of the nuclear receptor (NR) superfamily and BZIP family gained-open accessibility. Retinoic acid signaling and GTPase-mediated signal transduction are involved in adaptation to long days and maintenance of the HPG axis activation. According to our footprint analysis, three clock-output genes (TEF, DBP, and HLF) and the THRA were the first responders to long days in LD3. THRB, NR3C2, AR, and NR3C1 are the key players associated with the initiation and maintenance of the activation of the HPG axis, which appeared at LD7 and tended to be stable under long-day conditions. By integrating chromatin and the transcriptome, three genes (DIO2, SLC16A2, and PDE6H) involved in thyroid hormone signaling showed differential chromatin accessibility and expression levels during the seasonal activation of the HPG axis. TRPA1, a target of THRB identified by DAP-seq, was sensitive to photoactivation and exhibited differential expression levels between short- and long-day conditions.

CONCLUSION

Our data suggest that trans effects were the main factors affecting gene expression during the seasonal activation of the HPG axis. This study could lead to further research on the seasonal reproductive behavior of birds, particularly the role of MBH in controlling seasonal reproductive behavior.

Massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called single-cell combinatorial fluidic indexing ATAC-sequencing ("scifi-ATAC-seq"), which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using a widely commercialized microfluidics platform (10X Genomics). With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples in a single emulsion reaction can be indexed, representing a ~20-fold increase in throughput compared to the standard 10X Genomics workflow.

Integration analysis of ATAC-seq and RNA-seq provides insight into fatty acid biosynthesis in Schizochytrium limacinum under nitrogen limitation stress

BACKGROUND

Schizochytrium limacinum holds significant value utilized in the industrial-scale synthesis of natural DHA. Nitrogen-limited treatment can effectively increase the content of fatty acids and DHA, but there is currently no research on chromatin accessibility during the process of transcript regulation. The objective of this research was to delve into the workings of fatty acid production in S. limacinum by examining the accessibility of promoters and profiling gene expressions.

RESULTS

Results showed that differentially accessible chromatin regions (DARs)-associated genes were enriched in fatty acid metabolism, signal transduction mechanisms, and energy production. By identifying and annotating DARs-associated motifs, the study obtained 54 target transcription factor classes, including BPC, RAMOSA1, SPI1, MYC, and MYB families. Transcriptomics results revealed that several differentially expressed genes (DEGs), including SlFAD2, SlALDH, SlCAS1, SlNSDHL, and SlDGKI, are directly related to the biosynthesis of fatty acids, meanwhile, SlRPS6KA, SlCAMK1, SlMYB3R1, and SlMYB3R5 serve as transcription factors that could potentially influence the regulation of fatty acid production. In the integration analysis of DARs and ATAC-seq, 13 genes were identified, which were shared by both DEGs and DARs-associated genes, including SlCAKM, SlRP2, SlSHOC2, SlTN, SlSGK2, SlHMP, SlOGT, SlclpB, and SlDNAAF3.

CONCLUSIONS

SlCAKM may act as a negative regulator of fatty acid and DHA synthesis, while SlSGK2 may act as a positive regulator, which requires further study in the future. These insights enhance our comprehension of the processes underlying fatty acid and DHA production in S. limacinum. They also supply a foundational theoretical framework and practical assistance for the development of strains rich in fatty acids and DHA.

Dynamic chromatin accessibility landscapes of osteoblast differentiation and mineralization

Bone acts as a self-healing organ, which undergoes continuous regeneration process that is tightly regulated by the cooperation of osteoclasts with the capability of bone resorption and osteoblasts with the capability of bone formation. Generally, bone marrow derived mesenchymal stem cells (BMSCs) differentiated to final osteoblasts have been considered as critical role in bone remodeling. In this regard, several transcription factors (TFs) whose binding sites are initially hidden deep within accessible chromatin that participate in modulating osteoblast differentiation and bone matrix mineralization. Then, it is necessary to explore further the dynamic changes about the epigenetic transcription machinery during osteoblastogenesis. Here, we performed the chromatin accessibility and transcriptomic landscape of osteoblast differentiation and mineralization by using transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-Seq). Our data found that global chromatin accessibility during osteoblastogenesis was extensively improved. Above this, it is shown that key target genes including Col6a3, Serpina3n, Ms4a4d, Lyz2, Phf11b and Grin3a were enriched in differential loci RNA-seq and ATAC-Seq peaks with continuous changed tendency during osteoblasts differentiation and mineralization. In addition, Analysis of Motif Enrichment (AME) was used to elucidate TFs which modulated these target genes. In this study, it was shown for the first time that these important TFs including MEF2A, PRRX1, Shox2 and HOXB13 could alter promoter accessibility of target genes during osteoblastogenesis. This helps us understand how TF binding motif accessibility influences osteoblast differentiation. In addition, it also suggests that modulating the chromatin accessibility of osteogenesis could be developed as the promising strategies to regulate bone regeneration.

MYH1G-AS is a chromatin-associated lncRNA that regulates skeletal muscle development in chicken

BACKGROUND

Skeletal muscle development is pivotal for animal growth and health. Recently, long noncoding RNAs (lncRNAs) were found to interact with chromatin through diverse roles. However, little is known about how lncRNAs act as chromatin-associated RNAs to regulate skeletal muscle development. Here, we aim to investigate the regulation of chromatin-associated RNA (MYH1G-AS) during skeletal muscle development.

METHODS

We provided comprehensive insight into the RNA profile and chromatin accessibility of different myofibers, combining RNA sequencing (RNA-seq) with an assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq). The dual-luciferase reporter assay and chromatin immunoprecipitation (ChIP) assay were used to analyze the transcriptional regulation mechanism of MYH1G-AS. ALKBH5-mediated MYH1G-AS N6-methyladenosine (m6A) demethylation was assessed by a single-base elongation and ligation-based qPCR amplification method (SELECT) assay. Functions of MYH1G-AS were investigated through a primary myoblast and lentivirus/cholesterol-modified antisense oligonucleotide (ASO)-mediated animal model. To validate the interaction of MYH1G-AS with fibroblast growth factor 18 (FGF18) protein, RNA pull down and an RNA immunoprecipitation (RIP) assay were performed. Specifically, the interaction between FGF18 and SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A member 5 (SMARCA5) protein was analyzed by coimmunoprecipitation (Co-IP) and a yeast two-hybrid assay.

RESULTS

A total of 45 differentially expressed (DE) lncRNAs, with DE ATAC-seq peaks in their promoter region, were classified as open chromatin-associated lncRNAs. A skeletal muscle-specific lncRNA (MSTRG.15576.9; MYH1G-AS), which is one of the open chromatin-associated lncRNA, was identified. MYH1G-AS transcription is coordinately regulated by transcription factors (TF) SMAD3 and SP2. Moreover, SP2 represses ALKBH5 transcription to weaken ALKBH5-mediated m6A demethylation of MYH1G-AS, thus destroying MYH1G-AS RNA stability. MYH1G-AS accelerates myoblast proliferation but restrains myoblast differentiation. Moreover, MYH1G-AS drives a switch from slow-twitch to fast-twitch fibers and causes muscle atrophy. Mechanistically, MYH1G-AS inhibits FGF18 protein stabilization to reduce the interaction of FGF18 to SMARCA5, thus repressing chromatin accessibility of the SMAD4 promoter to activate the SMAD4-dependent pathway.

CONCLUSIONS

Our results reveal a new pattern of the regulation of lncRNA expression at diverse levels and help expound the regulation of m6A methylation on chromatin status.

Epigenetic inhibition of CTCF by HN1 promotes dedifferentiation and stemness of anaplastic thyroid cancer

Anaplastic thyroid cancer (ATC) is one of the deadliest cancers, whose important malignant feature is dedifferentiation. Chromatin remodeling is critical for tumorigenesis and progression, while its roles and regulator in facilitating dedifferentiation of ATC had been poorly understood. In our study, an emerging function of hematological and neurological expressed 1 (HN1) in promoting dedifferentiation of ATC cells was uncovered. HN1 expression was negatively correlated with the thyroid differentiation markers both at mRNA and protein level. Knockdown of HN1 in ATC cells effectively upregulated the thyroid differentiation markers and impeded the sphere formation capacity, accompanying with the loss of cancer stemness. In contrast, overexpression of HN1 drove the gain of stemness and the loss of thyroid differentiation markers. Nude mouse and zebrafish xenograft models showed that inhibition of HN1 in ATC cells effectively hindered tumor growth due to the loss of cancer stemness. Further study showed that HN1 was negatively correlated with CTCF in an independent thyroid-cancer cohort, and inhibition of HN1 enhanced the expression of CTCF in ATC cells. Overexpression of CTCF significantly reversed the dedifferentiation phenotypes of ATC cells, whereas simultaneously inhibiting HN1 and CTCF was unable to recover the level of thyroid differentiation markers. The combination of ATAC-seq and ChIP-seq analysis confirmed that CTCF regulated genes relating with thyroid gland development through influencing their chromatin accessibility. HN1 inhibited the acetylation of H3K27 at the promoter of CTCF by recruiting HDAC2, thereby inhibiting the transcriptional activation of CTCF. These findings demonstrated an essential role of HN1 in regulating the chromatin accessibility of thyroid differentiation genes during ATC dedifferentiation.

Single-Nucleus ATAC-seq for Mapping Chromatin Accessibility in Individual Cells of Murine Hearts

During the last decade a wide range of single-cell and single-nucleus next-generation sequencing techniques have been developed, which revolutionized detection of rare cell populations, enabling creation of comprehensive cell atlases of complex organs and tissues. State-of-the-art methods do not only allow classical transcriptomics of individual cells but also comprise a number of epigenetic approaches, including assessment of chromatin accessibility by single-nucleus Assay for Transposase Accessible Chromatin ATAC-seq (snATAC-seq). The snATAC-seq assay detects "open chromatin," a term for low nucleosome occupancy of genomic regions, which is a prerequisite for effective transcription factor binding. Information about open chromatin at the single-nucleus level helps to recognize epigenetic changes, sometimes before transcription of respective genes occurs. snATAC-seq detects cellular heterogeneity in otherwise still transcriptionally and/or morphologically homogeneous cell populations. Chromatin accessibility assays may be used to detect epigenetic changes in cardiac lineages during heart development, chromatin landscape changes during aging, and epigenetic alterations in heart diseases. Here, we provide an optimized protocol for snATAC-seq of murine hearts. We describe isolation of single nuclei from snap-frozen hearts, provide hints for preparation of libraries suitable for snATAC-seq next-generation sequencing (NGS) using the Chromium 10× platform, and give general recommendations for downstream analysis using conventional bioinformatic pipelines and packages. The protocol should serve as a beginner's guide to generate high-quality snATAC-seq datasets and to perform chromatin accessibility analysis of individual heart-derived cell nuclei.

Combined Analysis of mRNA Expression and Open Chromatin in Microglia

The advance of single-cell RNA-sequencing technologies in the past years has enabled unprecedented insights into the complexity and heterogeneity of microglial cell states in the homeostatic and diseased brain. This includes rather complex proteomic, metabolomic, morphological, transcriptomic, and epigenetic adaptations to external stimuli and challenges resulting in a novel concept of core microglia properties and functions. To uncover the regulatory programs facilitating the rapid transcriptomic adaptation in response to changes in the local microenvironment, the accessibility of gene bodies and gene regulatory elements can be assessed. Here, we describe the application of a previously published method for simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq) on microglia nuclei isolated from frozen mouse brain tissue.

Epigenetic landscape analysis reveals the significance of early reduced chromatin accessibility in osteoclastogenesis

Recently improved techniques could provide snapshots of chromatin structure generated based on chromatin accessibility. Since chromatin accessibility determines transcriptional potential, it has been attempted in a variety of cell systems. However, there has been no genome-wide analysis of chromatin accessibility for the entire murine osteoclast (OC) differentiation process. We performed an Assay for Transposase-Accessible Chromatin (ATAC)-sequencing (seq) during RANKL-induced OC differentiation and found that global chromatin accessibility decreased, especially early in OC differentiation. The global histone H3K27Ac level, an active histone modification mark, was diminished during OC differentiation by western blot and histone extract experiments. Its genomic enrichment was also reduced based on publicly available H3K27Ac chromatin immunoprecipitation (ChIP)-seq data. ATAC-seq and H3K27Ac ChIP-seq data demonstrated that RANKL induced a less accessible chromatin state during OC differentiation. Restoration of reduced H3K27Ac, presumably representing accessible states upon acetate treatment, suppresses OC differentiation by provoking immune-related gene expression. Subsequential integrative analysis of ATAC-seq, RNA-seq after acetate treatment, and H3K27Ac ChIP-seq reveals that Irf8 and its downstream targets are the most vulnerable to chromatin accessibility changes and acetate supplementation. Taken together, our study generated chromatin accessibility maps during the whole OC differentiation and suggested perturbation of chromatin accessibility might be a potential therapeutic strategy for excessive OC diseases.

3-D chromatin conformation, accessibility, and gene expression profiling of triple-negative breast cancer

OBJECTIVES

Triple-negative breast cancer (TNBC) is a highly aggressive breast cancer subtype with limited treatment options. Unlike other breast cancer subtypes, the scarcity of specific therapies and greater frequencies of distant metastases contribute to its aggressiveness. We aimed to find epigenetic changes that aid in the understanding of the dissemination process of these cancers.

DATA DESCRIPTION

Using CRISPR/Cas9, our experimental approach led us to identify and disrupt an insulator element, IE8, whose activity seemed relevant for cell invasion. The experiments were performed in two well-established TNBC cellular models, the MDA-MB-231 and the MDA-MB-436. To gain insights into the underlying molecular mechanisms of TNBC invasion ability, we generated and characterized high-resolution chromatin interaction (Hi-C) and chromatin accessibility (ATAC-seq) maps in both cell models and complemented these datasets with gene expression profiling (RNA-seq) in MDA-MB-231, the cell line that showed more significant changes in chromatin accessibility. Altogether, our data provide a comprehensive resource for understanding the spatial organization of the genome in TNBC cells, which may contribute to accelerating the discovery of TNBC-specific alterations triggering advances for this devastating disease.

TEAD4 antagonizes cellular senescence by remodeling chromatin accessibility at enhancer regions

Dramatic alterations in epigenetic landscapes are known to impact genome accessibility and transcription. Extensive evidence demonstrates that senescent cells undergo significant changes in chromatin structure; however, the mechanisms underlying the crosstalk between epigenetic parameters and gene expression profiles have not been fully elucidated. In the present study, we delineate the genome-wide redistribution of accessible chromatin regions that lead to broad transcriptome effects during senescence. We report that distinct senescence-activated accessibility regions (SAAs) are always distributed in H3K27ac-occupied enhancer regions, where they are responsible for elevated flanking senescence-associated secretory phenotype (SASP) expression and aberrant cellular signaling relevant to SASP secretion. Mechanistically, a single transcription factor, TEAD4, moves away from H3K27ac-labled SAAs to allow for prominent chromatin accessibility reconstruction during senescence. The enhanced SAAs signal driven by TEAD4 suppression subsequently induces a robust increase in the expression of adjacent SASP genes and the secretion of downstream factors, which contribute to the progression of senescence. Our findings illustrate a dynamic landscape of chromatin accessibility following senescence entry, and further reveal an insightful function for TEAD4 in regulating the broad chromatin state that modulates the overall transcriptional program of SASP genes.

The evolution and mutational robustness of chromatin accessibility in Drosophila

BACKGROUND

The evolution of genomic regulatory regions plays a critical role in shaping the diversity of life. While this process is primarily sequence-dependent, the enormous complexity of biological systems complicates the understanding of the factors underlying regulation and its evolution. Here, we apply deep neural networks as a tool to investigate the sequence determinants underlying chromatin accessibility in different species and tissues of Drosophila.

RESULTS

We train hybrid convolution-attention neural networks to accurately predict ATAC-seq peaks using only local DNA sequences as input. We show that our models generalize well across substantially evolutionarily diverged species of insects, implying that the sequence determinants of accessibility are highly conserved. Using our model to examine species-specific gains in accessibility, we find evidence suggesting that these regions may be ancestrally poised for evolution. Using in silico mutagenesis, we show that accessibility can be accurately predicted from short subsequences in each example. However, in silico knock-out of these sequences does not qualitatively impair classification, implying that accessibility is mutationally robust. Subsequently, we show that accessibility is predicted to be robust to large-scale random mutation even in the absence of selection. Conversely, simulations under strong selection demonstrate that accessibility can be extremely malleable despite its robustness. Finally, we identify motifs predictive of accessibility, recovering both novel and previously known motifs.

CONCLUSIONS

These results demonstrate the conservation of the sequence determinants of accessibility and the general robustness of chromatin accessibility, as well as the power of deep neural networks to explore fundamental questions in regulatory genomics and evolution.

scIBD: a self-supervised iterative-optimizing model for boosting the detection of heterotypic doublets in single-cell chromatin accessibility data

Application of the widely used droplet-based microfluidic technologies in single-cell sequencing often yields doublets, introducing bias to downstream analyses. Especially, doublet-detection methods for single-cell chromatin accessibility sequencing (scCAS) data have multiple assay-specific challenges. Therefore, we propose scIBD, a self-supervised iterative-optimizing model for boosting heterotypic doublet detection in scCAS data. scIBD introduces an adaptive strategy to simulate high-confident heterotypic doublets and self-supervise for doublet-detection in an iteratively optimizing manner. Comprehensive benchmarking on various simulated and real datasets demonstrates the outperformance and robustness of scIBD. Moreover, the downstream biological analyses suggest the efficacy of doublet-removal by scIBD.

Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

BACKGROUND

Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states.

METHODS

We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing.

RESULTS

We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease.

CONCLUSION

This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

IL-33 priming and antigenic stimulation synergistically promote the transcription of proinflammatory cytokine and chemokine genes in human skin mast cells

BACKGROUND

Antigenic stimulation through cross-linking the IgE receptor and epithelial cell-derived cytokine IL-33 are potent stimuli of mast cell (MC) activation. Moreover, IL-33 primes a variety of cell types, including MCs to respond more vigorously to external stimuli. However, target genes induced by the combined IL-33 priming and antigenic stimulation have not been investigated in human skin mast cells (HSMCs) in a genome-wide manner. Furthermore, epigenetic changes induced by the combined IL-33 priming and antigenic stimulation have not been evaluated.

RESULTS

We found that IL-33 priming of HSMCs enhanced their capacity to promote transcriptional synergy of the IL1B and CXCL8 genes by 16- and 3-fold, respectively, in response to combined IL-33 and antigen stimulation compared to without IL-33 priming. We identified the target genes in IL-33-primed HSMCs in response to the combined IL-33 and antigenic stimulation using RNA sequencing (RNA-seq). We found that the majority of genes synergistically upregulated in the IL-33-primed HSMCs in response to the combined IL-33 and antigenic stimulation were predominantly proinflammatory cytokine and chemokine genes. Moreover, the combined IL-33 priming and antigenic stimulation increase chromatin accessibility in the synergy target genes but not synergistically. Transcription factor binding motif analysis revealed more binding sites for NF-κB, AP-1, GABPA, and RAP1 in the induced or increased chromatin accessible regions of the synergy target genes.

CONCLUSIONS

Our study demonstrates that IL-33 priming greatly potentiates MCs' ability to transcribe proinflammatory cytokine and chemokine genes in response to antigenic stimulation, shining light on how epithelial cell-derived cytokine IL-33 can cause exacerbation of skin MC-mediated allergic inflammation.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer's disease brains

BACKGROUND

The genetic underpinnings of late-onset Alzheimer's disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD.

METHODS

Here, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs).

RESULTS

Integrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes.

CONCLUSIONS

To our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis-trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

Inferring chromatin accessibility during murine hematopoiesis through phylogenetic analysis

OBJECTIVE

Diversification of cell types and changes in epigenetic states during cell differentiation processes are important for understanding development. Recently, phylogenetic analysis using DNA methylation and histone modification information has been shown useful for inferring these processes. The purpose of this study was to examine whether chromatin accessibility data can help infer these processes in murine hematopoiesis.

RESULTS

Chromatin accessibility data could partially infer the hematopoietic differentiation hierarchy. Furthermore, based on the ancestral state estimation of internal nodes, the open/closed chromatin states of differentiating progenitor cells could be predicted with a specificity of 0.86-0.99 and sensitivity of 0.29-0.72. These results suggest that the phylogenetic analysis of chromatin accessibility could offer important information on cell differentiation, particularly for organisms from which progenitor cells are difficult to obtain.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Ribosome biogenesis disruption mediated chromatin structure changes revealed by SRAtac, a customizable end to end analysis pipeline for ATAC-seq

The nucleolus is a large nuclear body that serves as the primary site for ribosome biogenesis. Recent studies have suggested that it also plays an important role in organizing chromatin architecture. However, to establish a causal relationship between nucleolar ribosome assembly and chromatin architecture, genetic tools are required to disrupt nucleolar ribosome biogenesis. In this study, we used ATAC-seq to investigate changes in chromatin accessibility upon specific depletion of two ribosome biogenesis components, RPOA-2 and GRWD-1, in the model organism Caenorhabditis elegans. To facilitate the analysis of ATAC-seq data, we introduced two tools: SRAlign, an extensible NGS data processing workflow, and SRAtac, a customizable end-to-end ATAC-seq analysis pipeline. Our results revealed highly comparable changes in chromatin accessibility following both RPOA-2 and GRWD-1 perturbations. However, we observed a weak correlation between changes in chromatin accessibility and gene expression. While our findings corroborate the idea of a feedback mechanism between ribosomal RNA synthesis, nucleolar ribosome large subunit biogenesis, and chromatin structure during the L1 stage of C. elegans development, they also prompt questions regarding the functional impact of these alterations on gene expression.

Subcutaneous adipose tissue alteration in aging process associated with thyroid hormone signaling

BACKGROUND

Functional changes in subcutaneous adipose tissue (SAT) occur earlier in the aging process and play an important role in the occurrence and development of age-related metabolic diseases. The mechanism of this phenomenon is still unclear, and the change in adipose tissue with age is poorly understood.

METHODS

We used transcriptome sequencing (RNA seq) to screen differentially expressed genes at the mRNA level, and analyzed the functional characteristics of the differential genes through GO and KEGG analysis in human SAT of all ages. In order to clarify the specific mechanism of the functional change, we analyzed the chromatin accessibility in the promoter region in the same SAT used in the RNA seq by the assay for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) and obtained the functional genes in SAT changed with age. To verify these changes, we enlarged our sample content of human SAT. The primary mice adipocytes were extracted and stimulated by thyroid hormone of different concentration to construct an animal model, and the expression of the genes were determined through real-time Polymerase Chain Reaction(RT-PCR). The oxygen consumption test and immunofluorescence staining were used to determine the mitochondrial function of SAT.

RESULTS

RNA-seq showed characteristic gene expression of young and old human SAT, in which 331 genes were up-regulated and 349 genes were down-regulated. ATAC-seq, RNA-seq, combined with the mouse prediction model, determined the functional changed characteristics of seven genes. All these genes expressed differently in SAT of different ages, in which, NCF1, NLRP3, DUOX1 showed positive correlation with age; The expression of IFI30, P2RX1, P2RX6, PRODH, however, decreased with age. And all these genes showed dose dependent alternations under treatment of triiodothyroxine in mice SAT. The oxygen consumption rate revealed significant changes of the mitochondrial function and ROS accumulation in human SAT of different ages.

CONCLUSION

In elderly individuals, the function, in addition to distribution, of SAT undergoes significant changes, primarily in mitochondria, which may be due to insensitivity to thyroid hormone signaling. These results identified seven novel genes regulated by thyroid hormone, exhibiting significant changes in SAT of different age, and are probably related to the dysfunction of the aged SAT due to the mitochondrial damage and ROS accumulation.

High-fat diet in early life triggers both reversible and persistent epigenetic changes in the medaka fish (Oryzias latipes)

BACKGROUND

The nutritional status during early life can have enduring effects on an animal's metabolism, although the mechanisms underlying these long-term effects are still unclear. Epigenetic modifications are considered a prime candidate mechanism for encoding early-life nutritional memories during this critical developmental period. However, the extent to which these epigenetic changes occur and persist over time remains uncertain, in part due to challenges associated with directly stimulating the fetus with specific nutrients in viviparous mammalian systems.

RESULTS

In this study, we used medaka as an oviparous vertebrate model to establish an early-life high-fat diet (HFD) model. Larvae were fed with HFD from the hatching stages (one week after fertilization) for six weeks, followed by normal chow (NC) for eight weeks until the adult stage. We examined the changes in the transcriptomic and epigenetic state of the liver over this period. We found that HFD induces simple liver steatosis, accompanied by drastic changes in the hepatic transcriptome, chromatin accessibility, and histone modifications, especially in metabolic genes. These changes were largely reversed after the long-term NC, demonstrating the high plasticity of the epigenetic state in hepatocytes. However, we found a certain number of genomic loci showing non-reversible epigenetic changes, especially around genes related to cell signaling, liver fibrosis, and hepatocellular carcinoma, implying persistent changes in the cellular state of the liver triggered by early-life HFD feeding.

CONCLUSION

In summary, our data show that early-life HFD feeding triggers both reversible and persistent epigenetic changes in medaka hepatocytes. Our data provide novel insights into the epigenetic mechanism of nutritional programming and a comprehensive atlas of the long-term epigenetic state in an early-life HFD model of non-mammalian vertebrates.

Co-effects of m6A and chromatin accessibility dynamics in the regulation of cardiomyocyte differentiation

BACKGROUND

Cardiomyocyte growth and differentiation rely on precise gene expression regulation, with epigenetic modifications emerging as key players in this intricate process. Among these modifications, N6-methyladenosine (m6A) stands out as one of the most prevalent modifications on mRNA, exerting influence over mRNA metabolism and gene expression. However, the specific function of m6A in cardiomyocyte differentiation remains poorly understood.

RESULTS

We investigated the relationship between m6A modification and cardiomyocyte differentiation by conducting a comprehensive profiling of m6A dynamics during the transition from pluripotent stem cells to cardiomyocytes. Our findings reveal that while the overall m6A modification level remains relatively stable, the m6A levels of individual genes undergo significant changes throughout cardiomyocyte differentiation. We discovered the correlation between alterations in chromatin accessibility and the binding capabilities of m6A writers, erasers, and readers. The changes in chromatin accessibility influence the recruitment and activity of m6A regulatory proteins, thereby impacting the levels of m6A modification on specific mRNA transcripts.

CONCLUSION

Our data demonstrate that the coordinated dynamics of m6A modification and chromatin accessibility are prominent during the cardiomyocyte differentiation.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo, mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Y-box binding protein 1 promotes chromatin accessibility to aggravate liver fibrosis

Y-box binding protein 1 (YBX1) has been reported to be involved in the transcriptional regulation of various pathophysiological processes, such as inflammation, oxidative stress, and epithelial-mesenchymal transformation. However, its precise role and mechanism in regulating hepatic fibrosis remain unclear. In this study, we aimed to investigate the effects of YBX1 on liver fibrosis and its potential mechanism. The expression of YBX1 in human liver microarray, mice tissues and primary mouse hepatic stellate cells (HSCs) was validated to be upregulated in several hepatic fibrosis models (CCl₄ injection, TAA injection, and BDL). Hepatic-specific Ybx1 overexpression exacerbated the liver fibrosis phenotypes in vivo and in vitro. Moreover, the knockdown of YBX1 significantly improved TGF-β-induced fibrosis in the LX2 cell (a hepatic stellate cell line). Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) of hepatic-specific Ybx1 overexpression (Ybx1-OE) mice with CCl₄ injection showed increasing chromatin accessibility than CCl₄ only group. Functional enrichments of open regions in the Ybx1-OE group indicated that extracellular matrix (ECM) accumulation, lipid purine metabolism, and oxytocin-related pathways were more accessible in the Ybx1-OE group. Accessible regions of the Ybx1-OE group in the promoter also suggested significant activation of genes related to liver fibrogenesis, such as response to oxidative stress and ROS, lipid localization, angiogenesis and vascular development, and inflammatory regulation. Moreover, we screened and validated the expression of candidate genes (Fyn, Axl, Acsl1, Plin2, Angptl3, Pdgfb, Ccl24, and Arg2), which might be potential targets of Ybx1 in the pathogenesis of liver fibrosis.

Dynamic chromatin conformation and accessibility changes mediate the spatial-specific gene regulatory network in Bombyx mori

Silk gland genes of Bombyx mori can have strict spatial expression patterns, which impact their functions and silk quality; however, our understanding of their regulation mechanisms is currently insufficient. To address this, the middle silk gland (MSG) and posterior silk gland (PSG) of the silkworm were investigated. Gene ontology annotation showed that spatially specific expressed genes were involved in the formation of H3k9me and chromatin topology. Chromatin conformation data generated by Hi-C showed that the topologically associated domain boundaries around FibL and Sericin1 genes were significantly different between MSG and PSG. Changes in chromatin conformation led to changes in chromatin activity, which significantly affected the expression of nearby genes in silkworm. Chromatin accessibility regions of MSG and PSG were analyzed using FAIRE-seq, and 1006 transcription factor motifs were identified in open chromatin regions. Furthermore, the spatial-specific expression patterns of silk gland genes were mainly associated with homeobox-contained transcription factors, such as POU-M2, which was specifically bound and relatively highly expressed in the MSG. The regulatory network mediated by POU-M2 regulated most of the spatial-specific expressed genes in MSG, such as ADH1. These results can aid in improving silk performance, optimizing silkworm breeding, and improving the gene spatial regulatory model research for insects.

Co-profiling reveals distinct patterns of genomic chromatin accessibility and gene expression in pulmonary hypertension caused by chronic hypoxia

INTRODUCTION

Aberrant gene expression is a key mechanism underlying pulmonary hypertension (PH) development. The alterations of genomic chromatin accessibility and their relationship with the aberrant gene expressions in PH are poorly understood. We used bulk Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing (RNA-seq) in pulmonary artery smooth muscle cells (PASMCs) of chronic hypoxia-exposed rats mimicking group 3 human PH.

METHODS

Adult Sprague Dawley rats were commercially obtained from Hunan SJA (Hunan SJA Laboratory Animal Co., Changsha, China) and randomizedly allocated into four groups exposing to nomobaric hypoxia or normoxia for 1 or 28 days respectively. After the assessment of pulmonary hemodynamics, smooth muscle cells were isolated from intralobular arteries and simultaneously subjected to bulk Assay of ATAC-seq and RNA-seq.

RESULTS

Hypoxic exposure for continuous 28-days, but not for 1-day, induced established PH phenotypes in rats. ATAC-seq revealed a major distribution of differential accessibility regions (DARs) annotated to the genome in out-of-promoter regions, following 1-day or 28-days hypoxia. 1188 DAR-associated genes and 378 differentially expressed genes (DEGs) were identified in rats after exposure to 1-day hypoxia, while 238 DAR-associated genes and 452 DEGs for 28-days hypoxia. Most of the DAR-associated genes or DEGs in 1-day did not overlap with that of 28-days hypoxia. A Pearson correlation analysis indicated no significant correlation between ATAC-seq and RNA-seq.

CONCLUSIONS

The alterations in genomic chromatin accessibility and genes expression of PASMCs in the initial stage of hypoxia are distinct from the established stage of hypoxia-induced PH. The genomic differential accessibility regions may not be the main mechanisms directly underlying the differentially expressed genes observed either in the initial or established stages of PH. Thus the time-course alterations of gene expression and their possible indirect link with genomic chromatin accessibility warrant more attention in mechanistic study of pulmonary hypertension.

PlantCADB: A Comprehensive Plant Chromatin Accessibility Database

Chromatin accessibility landscapes are essential for detecting regulatory elements, illustrating the corresponding regulatory networks, and, ultimately, understanding the molecular basis underlying key biological processes. With the advancement of sequencing technologies, a large volume of chromatin accessibility data has been accumulated and integrated for humans and other mammals. These data have greatly advanced the study of disease pathogenesis, cancer survival prognosis, and tissue development. To advance the understanding of molecular mechanisms regulating plant key traits and biological processes, we developed a comprehensive plant chromatin accessibility database (PlantCADB) from 649 samples of 37 species. These samples are abiotic stress-related (such as heat, cold, drought, and salt; 159 samples), development-related (232 samples), and/or tissue-specific (376 samples). Overall, 18,339,426 accessible chromatin regions (ACRs) were compiled. These ACRs were annotated with genomic information, associated genes, transcription factor footprint, motif, and single-nucleotide polymorphisms (SNPs). Additionally, PlantCADB provides various tools to visualize ACRs and corresponding annotations. It thus forms an integrated, annotated, and analyzed plant-related chromatin accessibility resource, which can aid in better understanding genetic regulatory networks underlying development, important traits, stress adaptations, and evolution.PlantCADB is freely available at https://bioinfor.nefu.edu.cn/PlantCADB/.

Integrative epigenetic analysis reveals AP-1 promotes activation of tumor-infiltrating regulatory T cells in HCC

Regulatory T (Treg) cells that infiltrate human tumors exhibit stronger immunosuppressive activity compared to peripheral blood Treg cells (PBTRs), thus hindering the induction of effective antitumor immunity. Previous transcriptome studies have identified a set of genes that are conserved in tumor-infiltrating Treg cells (TITRs). However, epigenetic profiles of TITRs have not yet been completely deciphered. Here, we employed ATAC-seq and CUT&Tag assays to integrate transcriptome profiles and identify functional regulatory elements in TITRs. We observed a global difference in chromatin accessibility and enhancer landscapes between TITRs and PBTRs. We identified two types of active enhancer formation in TITRs. The H3K4me1-predetermined enhancers are poised to be activated in response to tumor microenvironmental stimuli. We found that AP-1 family motifs are enriched at the enhancer regions of TITRs. Finally, we validated that c-Jun binds at regulatory regions to regulate signature genes of TITRs and AP-1 is required for Treg cells activation in vitro. High c-Jun expression is correlated with poor survival in human HCC. Overall, our results provide insights into the mechanism of AP-1-mediated activation of TITRs and can hopefully be used to develop new therapeutic strategies targeting TITRs in liver cancer treatment.

NOMe-HiC: joint profiling of genetic variant, DNA methylation, chromatin accessibility, and 3D genome in the same DNA molecule

Cis-regulatory elements are coordinated to regulate the expression of their targeted genes. However, the joint measurement of cis-regulatory elements' activities and their interactions in spatial proximity is limited by the current sequencing approaches. We describe a method, NOMe-HiC, which simultaneously captures single-nucleotide polymorphisms, DNA methylation, chromatin accessibility (GpC methyltransferase footprints), and chromosome conformation changes from the same DNA molecule, together with the transcriptome, in a single assay. NOMe-HiC shows high concordance with state-of-the-art mono-omic assays across different molecular measurements and reveals coordinated chromatin accessibility at distal genomic segments in spatial proximity and novel types of long-range allele-specific chromatin accessibility.

Chromatin accessibility landscapes revealed the subgenome-divergent regulation networks during wheat grain development

Development of wheat (Triticum aestivum L.) grain mainly depends on the processes of starch synthesis and storage protein accumulation, which are critical for grain yield and quality. However, the regulatory network underlying the transcriptional and physiological changes of grain development is still not clear. Here, we combined ATAC-seq and RNA-seq to discover the chromatin accessibility and gene expression dynamics during these processes. We found that the chromatin accessibility changes are tightly associated with differential transcriptomic expressions, and the proportion of distal ACRs was increased gradually during grain development. Specific transcription factor (TF) binding sites were enriched at different stages and were diversified among the 3 subgenomes. We further predicted the potential interactions between key TFs and genes related with starch and storage protein biosynthesis and found different copies of some key TFs played diversified roles. Overall, our findings have provided numerous resources and illustrated the regulatory network during wheat grain development, which would shed light on the improvement of wheat yields and qualities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-023-00095-8.

Regulatory networks driving expression of genes critical for glioblastoma are controlled by the transcription factor c-Jun and the pre-existing epigenetic modifications

BACKGROUND

Glioblastoma (GBM, WHO grade IV) is an aggressive, primary brain tumor. Despite extensive tumor resection followed by radio- and chemotherapy, life expectancy of GBM patients did not improve over decades. Several studies reported transcription deregulation in GBMs, but regulatory mechanisms driving overexpression of GBM-specific genes remain largely unknown. Transcription in open chromatin regions is directed by transcription factors (TFs) that bind to specific motifs, recruit co-activators/repressors and the transcriptional machinery. Identification of GBM-related TFs-gene regulatory networks may reveal new and targetable mechanisms of gliomagenesis.

RESULTS

We predicted TFs-regulated networks in GBMs in silico and intersected them with putative TF binding sites identified in the accessible chromatin in human glioma cells and GBM patient samples. The Cancer Genome Atlas and Glioma Atlas datasets (DNA methylation, H3K27 acetylation, transcriptomic profiles) were explored to elucidate TFs-gene regulatory networks and effects of the epigenetic background. In contrast to the majority of tumors, c-Jun expression was higher in GBMs than in normal brain and c-Jun binding sites were found in multiple genes overexpressed in GBMs, including VIM, FOSL2 or UPP1. Binding of c-Jun to the VIM gene promoter was stronger in GBM-derived cells than in cells derived from benign glioma as evidenced by gel shift and supershift assays. Regulatory regions of the majority of c-Jun targets have distinct DNA methylation patterns in GBMs as compared to benign gliomas, suggesting the contribution of DNA methylation to the c-Jun-dependent gene expression.

CONCLUSIONS

GBM-specific TFs-gene networks identified in GBMs differ from regulatory pathways attributed to benign brain tumors and imply a decisive role of c-Jun in controlling genes that drive glioma growth and invasion as well as a modulatory role of DNA methylation.

Simultaneous Measurement of DNA Methylation and Nucleosome Occupancy in Single Cells Using scNOMe-Seq

Single-cell Nucleosome Occupancy and Methylome sequencing (scNOMe-seq) is a multimodal assay that simultaneously measures endogenous DNA methylation and nucleosome occupancy (i.e., chromatin accessibility) in single cells. scNOMe-seq combines the activity of a GpC Methyltransferase, an enzyme which methylates cytosines in GpC dinucleotides, with bisulfite conversion, whereby unmethylated cytosines are converted into thymines. Because GpC Methyltransferase acts only on cytosines present in non-nucleosomal regions of the genome, the subsequent bisulfite conversion step not only detects the endogenous DNA methylation, but also reveals the genome-wide pattern of chromatin accessibility. Implementing this technology at the single-cell level helps to capture the dynamics governing methylation and accessibility vary across individual cells and cell types. Here, we provide a scalable plate-based protocol for preparing scNOMe-seq libraries from single nucleus suspensions.

Genome-Wide Mapping of Active Regulatory Elements Using ATAC-seq

Active cis-regulatory elements (cREs) in eukaryotes are characterized by nucleosomal depletion and, accordingly, higher accessibility. This property has turned out to be immensely useful for identifying cREs genome-wide and tracking their dynamics across different cellular states and is the basis of numerous methods taking advantage of the preferential enzymatic cleavage/labeling of accessible DNA. ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) has emerged as the most versatile and widely adaptable method and has been widely adopted as the standard tool for mapping open chromatin regions. Here, we discuss the current optimal practices and important considerations for carrying out ATAC-seq experiments, primarily in the context of mammalian systems.

Single-Molecule Mapping of Chromatin Accessibility Using NOMe-seq/dSMF

The bulk of gene expression regulation in most organisms is accomplished through the action of transcription factors (TFs) on cis-regulatory elements (CREs). In eukaryotes, these CREs are generally characterized by nucleosomal depletion and thus higher physical accessibility of DNA. Many methods exploit this property to map regions of high average accessibility, and thus putative active CREs, in bulk. However, these techniques do not provide information about coordinated patterns of accessibility along the same DNA molecule, nor do they map the absolute levels of occupancy/accessibility. SMF (Single-Molecule Footprinting) fills these gaps by leveraging recombinant DNA cytosine methyltransferases (MTase) to mark accessible locations on individual DNA molecules. In this chapter, we discuss current methods and important considerations for performing SMF experiments.

Simultaneous Single-Cell Profiling of the Transcriptome and Accessible Chromatin Using SHARE-seq

The ability to analyze the transcriptomic and epigenomic states of individual single cells has in recent years transformed our ability to measure and understand biological processes. Recent advancements have focused on increasing sensitivity and throughput to provide richer and deeper biological insights at the cellular level. The next frontier is the development of multiomic methods capable of analyzing multiple features from the same cell, such as the simultaneous measurement of the transcriptome and the chromatin accessibility of candidate regulatory elements. In this chapter, we discuss and describe SHARE-seq (Simultaneous high-throughput ATAC, and RNA expression with sequencing) for carrying out simultaneous chromatin accessibility and transcriptome measurements in single cells, together with the experimental and analytical considerations for achieving optimal results.

ATAC-See: A Tn5 Transposase-Mediated Assay for Detection of Chromatin Accessibility with Imaging

Assay of transposase-accessible chromatin with visualization (ATAC-see), a transposase-mediated imaging technology that enables direct imaging of the accessible genome in situ and deep sequencing to reveal the identity of the imaged elements. Here we image spatial organization of the accessible genome in HT1080 cells with this method.

Analysis of Chromatin Interaction and Accessibility by Trac-Looping

Spatial organization of the genome modulates pivotal biological processes. The emerging new technologies have provided novel insights into genome structure and its role in regulating cell activities. To examine the genome-wide chromatin interactions at accessible chromatin regions, we developed a DNA transposase-mediated analysis of chromatin looping (Trac-looping) method for simultaneously detecting chromatin interactions and chromatin accessibility. Here, we describe a detailed protocol of generating Trac-looping libraries.

The SAGA histone acetyltransferase module targets SMC5/6 to specific genes

BACKGROUND

Structural Maintenance of Chromosomes (SMC) complexes are molecular machines driving chromatin organization at higher levels. In eukaryotes, three SMC complexes (cohesin, condensin and SMC5/6) play key roles in cohesion, condensation, replication, transcription and DNA repair. Their physical binding to DNA requires accessible chromatin.

RESULTS

We performed a genetic screen in fission yeast to identify novel factors required for SMC5/6 binding to DNA. We identified 79 genes of which histone acetyltransferases (HATs) were the most represented. Genetic and phenotypic analyses suggested a particularly strong functional relationship between the SMC5/6 and SAGA complexes. Furthermore, several SMC5/6 subunits physically interacted with SAGA HAT module components Gcn5 and Ada2. As Gcn5-dependent acetylation facilitates the accessibility of chromatin to DNA-repair proteins, we first analysed the formation of DNA-damage-induced SMC5/6 foci in the Δgcn5 mutant. The SMC5/6 foci formed normally in Δgcn5, suggesting SAGA-independent SMC5/6 localization to DNA-damaged sites. Next, we used Nse4-FLAG chromatin-immunoprecipitation (ChIP-seq) analysis in unchallenged cells to assess SMC5/6 distribution. A significant portion of SMC5/6 accumulated within gene regions in wild-type cells, which was reduced in Δgcn5 and Δada2 mutants. The drop in SMC5/6 levels was also observed in gcn5-E191Q acetyltransferase-dead mutant.

CONCLUSION

Our data show genetic and physical interactions between SMC5/6 and SAGA complexes. The ChIP-seq analysis suggests that SAGA HAT module targets SMC5/6 to specific gene regions and facilitates their accessibility for SMC5/6 loading.