Gene regulation relies on the precise binding of transcription factors (TFs) at regulatory elements, but simultaneously detecting hundreds of TFs on chromatin is challenging. We developed cFOOT-seq, a cytosine deaminase-based TF footprinting assay, for high-resolution, quantitative genome-wide assessment of TF binding in both open and closed chromatin regions, even with small cell numbers. By utilizing the dsDNA deaminase SsdAtox, cFOOT-seq converts accessible cytosines to uracil while preserving genomic integrity, making it compatible with techniques like ATAC-seq for sensitive and cost-effective detection of TF occupancy at the single-molecule and single-cell level. Our approach enables the delineation of TF footprints, quantification of occupancy, and examination of chromatin influences on TF binding. Notably, cFOOT-seq, combined with FootTrack analysis, enables de novo prediction of TF binding sites and tracking of TF occupancy dynamics. We demonstrate its application in capturing cell type-specific TFs, analyzing TF dynamics during reprogramming, and revealing TF dependencies on chromatin remodelers. Overall, cFOOT-seq represents a robust approach for investigating the genome-wide dynamics of TF occupancy and elucidating the cis-regulatory architecture underlying gene regulation.
Transcription Factors/genetics* ; Humans ; Chromatin/genetics* ; Animals ; Binding Sites ; Mice ; DNA Footprinting/methods*

The p60 subunit of the chromatin assembly factor-1 complex, that is, chromatin assembly factor-1 subunit B (CHAF1B), is a histone H3/H4 chaperone crucial for the transcriptional regulation of cell differentiation and self-renewal. CHAF1B is overexpressed in several cancers and may represent a potential target for cancer therapy. However, its expression and clinical significance in lung squamous-cell carcinoma (LUSC) remain unclear. In this study, we performed weighted gene correlation network analysis to analyze the Gene Expression Omnibus GSE68793 LUSC dataset and identified CHAF1B as one of the most important driver gene candidates. Immunohistochemical analysis of 126 LUSC tumor samples and 80 adjacent normal lung tissues showed the marked upregulation of CHAF1B in tumor tissues and the negative association of its expression level with patient survival outcomes. Silencing of CHAF1B suppressed LUSC proliferation in vitro and LUSC tumor growth in vivo. Furthermore, bulk RNA sequencing of CHAF1B knockdown cells indicated SET domain containing 7 (SETD7) as a significant CHAF1B target gene. In addition, CHAF1B competitively binds to the SETD7 promoter region and represses its transcription. Altogether, these results imply that CHAF1B plays a vital role in LUSC tumorigenesis and may represent a potential molecular target for this deadly disease.
Humans ; Lung Neoplasms/metabolism* ; Histone-Lysine N-Methyltransferase/metabolism* ; Carcinoma, Squamous Cell/metabolism* ; Gene Expression Regulation, Neoplastic ; Disease Progression ; Cell Proliferation/genetics* ; Cell Line, Tumor ; Chromatin Assembly Factor-1/metabolism* ; Animals ; Mice ; Male ; Female

Objective To investigate the effects of folate deficiency on changes in histone H3 lysine 4 (H3K4) mono-methylation (me1)-marked enhancers and the molecular mechanism underpinning the folate deficiency-induced neural tube defects (NTD). Methods Mouse embryonic stem cells (mESCs) were cultured in the folate-free DMEM medium (folate-deficient group) and the DMEM medium containing 4 mg/L folate (normal control group),respectively.Chromatin immunoprecipitation sequencing (ChIP-seq) was performed for H3K4me1. The mouse model of folate-induced NTD was established,and transcriptome sequencing (RNA-seq) was performed for the brain tissue of fetal mice to reveal the differential expression profiles.The results were validated through real-time quantitative polymerase chain reaction (RT-qPCR).The activity of the differential peak regions of H3K4me1 was verified through the luciferase reporter assay. Results The folate content in the mESCs cultured in the folate-free medium reduced compared with that in the normal control group (P=0.008).The H3K4me1-maked enhancers in the mESCs cultured in the folate-free medium induced significant changes in intronic regions,and these changes were concentrated in metabolic and energy metabolism processes (q=9.56×10^-48,P=1.28×10^-47).The differentially expressed genes harboring H3K4me1-marked enhancers in mESCs were mainly enriched in the Wnt signaling pathway (q=0.004,P=0.004 7).ChIP-qPCR results confirmed that H3K4me1 binding decreased in the differential peak regions of the Ldlrap1 gene (P=0.008),Camta1 gene (P=0.002),and Apc2 gene (P=0.012).The H3K4 demethylase inhibitor T-448 effectively reversed the H3K4me1 binding in the differential peak regions of the aforementioned genes (P=0.01).The results of RNA-seq for the brain tissue of NTD fetal mice showed significant enrichment of the differentially expressed genes in the Wnt signaling pathway (P=1.52×10^-5).The enrichment of differential peak regions of H3K4me1-marked enhancers in Apc2,Ldlrap1,and Camta1 genes in the brain tissue also showed significant changes.The differential peak region in Apc2 exhibited transcription factor activity (P=0.020). Conclusion Folate deficiency may affect changes in H3K4me1-marked enhancers to participate in the regulation of neural tube closure genes,thereby inducing the occurrence of NTD.
Neural Tube Defects/genetics* ; Animals ; Mice ; Folic Acid Deficiency/complications* ; Histones/metabolism* ; Folic Acid/metabolism* ; Methylation ; Mouse Embryonic Stem Cells/metabolism* ; Wnt Signaling Pathway ; Lysine/metabolism* ; Chromatin Immunoprecipitation Sequencing

Hypertrophic cardiomyopathy (HCM) is the most common inherited heart disease and is characterized by primary left ventricular hypertrophy usually caused by mutations in sarcomere genes. The mechanism underlying cardiac remodeling in HCM remains incompletely understood. An investigation of HCM through integrative analysis at multi-omics levels will be helpful for treating HCM. DNA methylation and chromatin accessibility, as well as gene expression, were assessed by nucleosome occupancy and methylome sequencing (NOMe-seq) and RNA-seq, respectively, using the cardiac tissues of HCM patients. Compared with those of the controls, the transcriptome, DNA methylome, and chromatin accessibility of the HCM myocardium showed multifaceted differences. At the transcriptome level, HCM hearts returned to the fetal gene program through decreased sarcomeric and metabolic gene expression and increased extracellular matrix gene expression. In the DNA methylome, hypermethylated and hypomethylated differentially methylated regions were identified in HCM. At the chromatin accessibility level, HCM hearts showed changes in different genome elements. Several transcription factors, including SP1 and EGR1, exhibited a fetal-like pattern of binding motifs in nucleosome-depleted regions in HCM. In particular, the inhibition of SP1 or EGR1 in an HCM mouse model harboring sarcomere mutations markedly alleviated the HCM phenotype of the mutant mice and reversed fetal gene reprogramming. Overall, this study not only provides a high-precision multi-omics map of HCM heart tissue but also sheds light on the therapeutic strategy by intervening in the fetal gene reprogramming in HCM.
Cardiomyopathy, Hypertrophic/metabolism* ; Humans ; Animals ; DNA Methylation ; Mice ; Transcriptome ; Chromatin/genetics* ; Early Growth Response Protein 1/metabolism* ; Male ; Epigenome ; Nucleosomes/genetics* ; Female ; Middle Aged ; Disease Models, Animal ; Adult

Previous studies have shown that long-term spermatogonial stem cells (SSCs) have the potential to spontaneously transform into pluripotent stem cells, which is speculated to be related to the tumorigenesis of testicular germ cells, especially when p53 is deficient in SSCs which shows a significant increase in the spontaneous transformation efficiency. Energy metabolism has been proved to be strongly associated with the maintenance and acquisition of pluripotency. Recently, we compared the difference in chromatin accessibility and gene expression profiles between wild-type (p53^+/+) and p53 deficient (p53^-/-) mouse SSCs using the Assay for Targeting Accessible-Chromatin with high-throughput sequencing (ATAC-seq) and transcriptome sequencing (RNA-seq) techniques, and revealed that SMAD3 is a key transcription factor in the transformation of SSCs into pluripotent cells. In addition, we also observed significant changes in the expression levels of many genes related to energy metabolism after p53 deletion. To further reveal the role of p53 in the regulation of pluripotency and energy metabolism, this paper explored the effects and mechanism of p53 deletion on energy metabolism during the pluripotent transformation of SSCs. The results of ATAC-seq and RNA-seq from p53^+/+ and p53^-/- SSCs revealed that gene chromatin accessibility related to positive regulation of glycolysis and electron transfer and ATP synthesis was increased, and the transcription levels of genes encoding key glycolytic enzymes and regulating electron transport-related enzymes were markedly increased. Furthermore, transcription factors SMAD3 and SMAD4 promoted glycolysis and energy homeostasis by binding to the chromatin of the Prkag2 gene which encodes the AMPK subunit. These results suggest that p53 deficiency activates the key enzyme genes of glycolysis in SSCs and enhances the chromatin accessibility of genes associated with glycolysis activation to improve glycolysis activity and promote transformation to pluripotency. Moreover, SMAD3/SMAD4-mediated transcription of the Prkag2 gene ensures the energy demand of cells in the process of pluripotency transformation and maintains cell energy homeostasis by promoting AMPK activity. These results shed light on the importance of the crosstalk between energy metabolism and stem cell pluripotency transformation, which might be helpful for clinical research of gonadal tumors.
Animals ; Mice ; AMP-Activated Protein Kinases ; Chromatin ; Energy Metabolism ; Gene Deletion ; Stem Cells ; Tumor Suppressor Protein p53/genetics* ; Spermatogonia/cytology* ; Male

METTL3 and METTL14 are two components that form the core heterodimer of the main RNA m6A methyltransferase complex (MTC) that installs m6A. Surprisingly, depletion of METTL3 or METTL14 displayed distinct effects on stemness maintenance of mouse embryonic stem cell (mESC). While comparable global hypo-methylation in RNA m6A was observed in Mettl3 or Mettl14 knockout mESCs, respectively. Mettl14 knockout led to a globally decreased nascent RNA synthesis, whereas Mettl3 depletion resulted in transcription upregulation, suggesting that METTL14 might possess an m6A-independent role in gene regulation. We found that METTL14 colocalizes with the repressive H3K27me3 modification. Mechanistically, METTL14, but not METTL3, binds H3K27me3 and recruits KDM6B to induce H3K27me3 demethylation independent of METTL3. Depletion of METTL14 thus led to a global increase in H3K27me3 level along with a global gene suppression. The effects of METTL14 on regulation of H3K27me3 is essential for the transition from self-renewal to differentiation of mESCs. This work reveals a regulatory mechanism on heterochromatin by METTL14 in a manner distinct from METTL3 and independently of m6A, and critically impacts transcriptional regulation, stemness maintenance, and differentiation of mESCs.
Animals ; Mice ; Methylation ; Chromatin ; Histones/metabolism* ; RNA, Messenger/genetics* ; Methyltransferases/metabolism* ; RNA/metabolism*

The cranial neural crest plays a fundamental role in orofacial development and morphogenesis. Accordingly, mutations with impact on the cranial neural crest and its development lead to orofacial malformations such as cleft lip and palate. As a pluripotent and dynamic cell population, the cranial neural crest undergoes vast transcriptional and epigenomic alterations throughout the formation of facial structures pointing to an essential role of factors regulating chromatin state or transcription levels. Using CRISPR/Cas9-guided genome editing and conditional mutagenesis in the mouse, we here show that inactivation of Kat5 or Ep400 as the two essential enzymatic subunits of the Tip60/Ep400 chromatin remodeling complex severely affects carbohydrate and amino acid metabolism in cranial neural crest cells. The resulting decrease in protein synthesis, proliferation and survival leads to a drastic reduction of cranial neural crest cells early in fetal development and a loss of most facial structures in the absence of either protein. Following heterozygous loss of Kat5 in neural crest cells palatogenesis was impaired. These findings point to a decisive role of the Tip60/Ep400 chromatin remodeling complex in facial morphogenesis and lead us to conclude that the orofacial clefting observed in patients with heterozygous KAT5 missense mutations is at least in part due to disturbances in the cranial neural crest.
Animals ; Mice ; Chromatin Assembly and Disassembly ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; DNA Helicases/metabolism* ; DNA-Binding Proteins ; Neural Crest/metabolism* ; Skull ; Transcription Factors/metabolism*

Genome-wide association studies (GWASs) are the most widely used method to identify genetic risk loci associated with orofacial clefts (OFC). However, despite the increasing size of cohort, GWASs are still insufficient to detect all the heritability, suggesting there are more associations under the current stringent statistical threshold. In this study, we obtained an integrated epigenomic dataset based on the chromatin conformation of a human oral epithelial cell line (HIOEC) using RNA-seq, ATAC-seq, H3K27ac ChIP-seq, and DLO Hi-C. Presumably, this epigenomic dataset could reveal the missing functional variants located in the oral epithelial cell active enhancers/promoters along with their risk target genes, despite relatively less-stringent statistical association with OFC. Taken a non-syndromic cleft palate only (NSCPO) GWAS data of the Chinese Han population as an example, 3664 SNPs that cannot reach the strict significance threshold were subjected to this functional identification pipeline. In total, 254 potential risk SNPs residing in active cis-regulatory elements interacting with 1 718 promoters of oral epithelium-expressed genes were screened. Gapped k-mer machine learning based on enhancers interacting with epithelium-expressed genes along with in vivo and in vitro reporter assays were employed as functional validation. Among all the potential SNPs, we chose and confirmed that the risk alleles of rs560789 and rs174570 reduced the epithelial-specific enhancer activity by preventing the binding of transcription factors related to epithelial development. In summary, we established chromatin conformation datasets of human oral epithelial cells and provided a framework for testing and understanding how regulatory variants impart risk for clefts.
Chromatin ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; Epithelium ; Genome-Wide Association Study ; Humans

Understanding the regulatory networks for germ cell fate specification is necessary to developing strategies for improving the efficiency of germ cell production in vitro. In this study, we developed a coupled screening strategy that took advantage of an arrayed bi-molecular fluorescence complementation (BiFC) platform for protein-protein interaction screens and epiblast-like cell (EpiLC)-induction assays using reporter mouse embryonic stem cells (mESCs). Investigation of candidate interaction partners of core human pluripotent factors OCT4, NANOG, KLF4 and SOX2 in EpiLC differentiation assays identified novel primordial germ cell (PGC)-inducing factors including BEN-domain (BEND/Bend) family members. Through RNA-seq, ChIP-seq, and ATAC-seq analyses, we showed that Bend5 worked together with Bend4 and helped mark chromatin boundaries to promote EpiLC induction in vitro. Our findings suggest that BEND/Bend proteins represent a new family of transcriptional modulators and chromatin boundary factors that participate in gene expression regulation during early germline development.
Animals ; Cell Differentiation/genetics* ; Chromatin/metabolism* ; Embryonic Stem Cells ; Germ Cells/metabolism* ; Germ Layers/metabolism* ; Mice

OBJECTIVE: To observe the effect of Buyi Pishen acupuncture (acupuncture for invigorating spleen and kidney) on inflammatory factor and synovial cartilage matrix in adjuvant arthritis (AA) rats, and to explore the mechanism of acupuncture for rheumatoid arthritis (RA). METHODS: A total of 60 clean male Wistar rats were randomized into a normal group, a model group, a tripterygium wilfordii polyglycoside tablet (TWP) group and an acupuncture group, 15 rats in each group. Rats in the model group, the TWP group and the acupuncture group received intradermal injection of Freund's complete adjuvant (FCA) at right hind foot pad to induce the AA model. TWP suspension of 8 mg/kg was given by gavage in the TWP group. Acupuncture was applied at "Shenshu" (BL 23), "Pishu" (BL 20) and right "Housanli" (ST 36), "Sanyinjiao" (SP 6), "Yanglingquan" (GB 34) in the acupuncture group, 15 min a time, once a day. The intervention was given 15 days in both TWP group and acupuncture group. The foot-pad swelling degree before modeling, before and after intervention and the arthritis index (AI) score before and after intervention were calculated; the serum levels of interleukin (IL)-1β, IL-4, IL-10 and tumor necrosis factor-α (TNF-α) were detected by ELISA method; the ultrastructure and histomorphological changes of synovium issue were observed by transmission electron microscope and HE staining; the positive expression of matrix metalloproteinase (MMP)-3 and MMP-9 in synovium issue was detected by immunohistochemistry method. RESULTS: Before intervention, foot-pad swelling degree of the model group, the TWP group and the acupuncture group was increased compared with the normal group (P<0.01). After intervention, foot-pad swelling degree and AI score were increased compared with the normal group (P<0.01), foot-pad swelling degree and AI scores in the TWP group and the acupuncture group were lower than the model group (P<0.05), and those in the acupuncture group were decreased compared with the TWP group (P<0.05). The model group exhibited unclear nuclear membrane of synovial cells, chromatin pyknosis, massive inflammatory cell infiltration and hyperplasia in synovial tissue; the TWP group and the acupuncture group exhibited clear and smooth nuclear membrane of synovial cells, inapparent chromatin pyknosis, less inflammatory cell infiltration and hyperplasia in synovial tissue, the acupuncture group exhibited less matrix destruction as well. Compared with the normal group, serum levels of IL-1β and TNF-α and positive expression of MMP-3 and MMP-9 in synovium issue were increased (P<0.01), while serum levels of IL-4 and IL-10 were decreased (P<0.01) in the model group. Compared with the model group, serum levels of IL-1β and TNF-α and positive expression of MMP-3 and MMP-9 in synovium issue were decreased (P<0.05, P<0.01), while serum levels of IL-4 and IL-10 were increased (P<0.05) in the TWP group and the acupuncture group; compared with the TWP group, serum level of TNF-α and positive expression of MMP-3 and MMP-9 in synovium issue were decreased (P<0.05), while serum levels of IL-4 and IL-10 were increased (P<0.05) in the acupuncture group. CONCLUSION Buyi Pishen acupuncture can effectively improve the injury of articular cartilage in AA rats, its mechanism maybe related to reducing the inflammatory reaction in synovium and inhibiting the degradation of articular cartilage matrix.
Acupuncture Therapy ; Animals ; Arthritis, Experimental/therapy* ; Cartilage, Articular ; Chromatin ; Hyperplasia ; Interleukin-10 ; Interleukin-4 ; Male ; Matrix Metalloproteinase 3 ; Matrix Metalloproteinase 9 ; Rats ; Rats, Wistar ; Tumor Necrosis Factor-alpha/genetics*