Impact of cytosine methylation on DNA binding specificities of human transcription factors

Using methylation data to improve transcription factor binding prediction

Modelling the regulatory mechanisms that determine cell fate, response to external perturbation, and disease state depends on measuring many factors, a task made more difficult by the plasticity of the epigenome. Scanning the genome for the sequence patterns defined by Position Weight Matrices (PWM) can be used to estimate transcription factor (TF) binding locations. However, this approach does not incorporate information regarding the epigenetic context necessary for TF binding. CpG methylation is an epigenetic mark influenced by environmental factors that is commonly assayed in human cohort studies. We developed a framework to score inferred TF binding locations using methylation data. We intersected motif locations identified using PWMs with methylation information captured in both whole-genome bisulfite sequencing and Illumina EPIC array data for six cell lines, scored motif locations based on these data, and compared with experimental data characterizing TF binding (ChIP-seq). We found that for most TFs, binding prediction improves using methylation-based scoring compared to standard PWM-scores. We also illustrate that our approach can be generalized to infer TF binding when methylation information is only proximally available, i.e. measured for nearby CpGs that do not directly overlap with a motif location. Overall, our approach provides a framework for inferring context-specific TF binding using methylation data. Importantly, the availability of DNA methylation data in existing patient populations provides an opportunity to use our approach to understand the impact of methylation on gene regulatory processes in the context of human disease.

Single-housing-induced islet epigenomic changes are related to polymorphisms in diabetic KK mice

A lack of social relationships is increasingly recognized as a type 2 diabetes (T2D) risk. To investigate the underlying mechanism, we used male KK mice, an inbred strain with spontaneous diabetes. Given the association between living alone and T2D risk in humans, we divided the non-diabetic mice into singly housed (KK-SH) and group-housed control mice. Around the onset of diabetes in KK-SH mice, we compared H3K27ac ChIP-Seq with RNA-Seq using pancreatic islets derived from each experimental group, revealing a positive correlation between single-housing-induced changes in H3K27ac and gene expression levels. In particular, single-housing-induced H3K27ac decreases revealed a significant association with islet cell functions and GWAS loci for T2D and related diseases, with significant enrichment of binding motifs for transcription factors representative of human diabetes. Although these H3K27ac regions were preferentially localized to a polymorphic genomic background, SNVs and indels did not cause sequence disruption of enriched transcription factor motifs in most of these elements. These results suggest alternative roles of genetic variants in environment-dependent epigenomic changes and provide insights into the complex mode of disease inheritance.

TET2-STAT3-CXCL5 nexus promotes neutrophil lipid transfer to fuel lung adeno-to-squamous transition

Phenotypic plasticity is a rising cancer hallmark, and lung adeno-to-squamous transition (AST) triggered by LKB1 inactivation is significantly associated with drug resistance. Mechanistic insights into AST are urgently needed to identify therapeutic vulnerability in LKB1-deficient lung cancer. Here, we find that ten-eleven translocation (TET)-mediated DNA demethylation is elevated during AST in KrasLSL-G12D/+; Lkb1L/L (KL) mice, and knockout of individual Tet genes reveals that Tet2 is required for squamous transition. TET2 promotes neutrophil infiltration through STAT3-mediated CXCL5 expression. Targeting the STAT3-CXCL5 nexus effectively inhibits squamous transition through reducing neutrophil infiltration. Interestingly, tumor-infiltrating neutrophils are laden with triglycerides and can transfer the lipid to tumor cells to promote cell proliferation and squamous transition. Pharmacological inhibition of macropinocytosis dramatically inhibits neutrophil-to-cancer cell lipid transfer and blocks squamous transition. These data uncover an epigenetic mechanism orchestrating phenotypic plasticity through regulating immune microenvironment and metabolic communication, and identify therapeutic strategies to inhibit AST.

Dynamic DNA methylation modifications in the cold stress response of cassava

Cassava, a crucial tropical crop, faces challenges from cold stress, necessitating an exploration of its molecular response. Here, we investigated the role of DNA methylation in moderating the response to moderate cold stress (10 °C) in cassava. Using whole-genome bisulfite sequencing, we examined DNA methylation patterns in leaf blades and petioles under control conditions, 5 h, and 48 h of cold stress. Tissue-specific responses were observed, with leaf blades exhibiting subtle changes, while petioles displayed a pronounced decrease in methylation levels under cold stress. We identified cold stress-induced differentially methylated regions (DMRs) that demonstrated both tissue and treatment specificity. Importantly, these DMRs were enriched in genes with altered expression, implying functional relevance. The cold-response transcription factor ERF105 associated with DMRs emerged as a significant and conserved regulator across tissues and treatments. Furthermore, we investigated DNA methylation dynamics in transposable elements, emphasizing the sensitivity of MITEs with bHLH binding motifs to cold stress. These findings provide insights into the epigenetic regulation of response to cold stress in cassava, contributing to an understanding of the molecular mechanisms underlying stress adaptation in this tropical plant.

Predicting regulatory mutations and their target genes by new computational integrative analysis: A study of follicular lymphoma

Mutations in DNA regulatory regions are increasingly being recognized as important drivers of cancer and other complex diseases. These mutations can regulate gene expression by affecting DNA-protein binding and epigenetic profiles, such as DNA methylation in genome regulatory elements. However, identifying mutation hotspots associated with expression regulation and disease progression in non-coding DNA remains a challenge. Unlike most existing approaches that assign a mutation score to individual single nucleotide polymorphisms (SNP), a mutation block (MB)-based approach was introduced in this study to assess the collective impact of a cluster of SNPs on transcription factor-DNA binding affinity, differential gene expression (DEG), and nearby DNA methylation. Moreover, the long-distance target genes of functional MBs were identified using a new permutation-based algorithm that assessed the significance of correlations between DNA methylation at regulatory regions and target gene expression. Two new Python packages were developed. The Differential Methylation Region (DMR-analysis) analysis tool was used to detect DMR and map them to regulatory elements. The second tool, an integrated DMR, DEG, and SNP analysis tool (DDS-analysis), was used to combine the omics data to identify functional MBs and long-distance target genes. Both tools were validated in follicular lymphoma (FL) cohorts, where not only known functional MBs and their target genes (BCL2 and BCL6) were recovered, but also novel genes were found, including CDCA4 and JAG2, which may be associated with FL development. These genes are linked to target gene expression and are significantly correlated with the methylation of nearby DNA sequences in FL. The proposed computational integrative analysis of multiomics data holds promise for identifying regulatory mutations in cancer and other complex diseases.

Uncovering the role of TET2-mediated ENPEP activation in trophoblast cell fate determination

Early trophoblast differentiation is crucial for embryo implantation, placentation and fetal development. Dynamic changes in DNA methylation occur during preimplantation development and are critical for cell fate determination. However, the underlying regulatory mechanism remains unclear. Recently, we derived morula-like expanded potential stem cells from human preimplantation embryos (hEPSC-em), providing a valuable tool for studying early trophoblast differentiation. Data analysis on published datasets showed differential expressions of DNA methylation enzymes during early trophoblast differentiation in human embryos and hEPSC-em derived trophoblastic spheroids. We demonstrated downregulation of DNA methyltransferase 3 members (DNMT3s) and upregulation of ten-eleven translocation methylcytosine dioxygenases (TETs) during trophoblast differentiation. While DNMT inhibitor promoted trophoblast differentiation, TET inhibitor hindered the process and reduced implantation potential of trophoblastic spheroids. Further integrative analysis identified that glutamyl aminopeptidase (ENPEP), a trophectoderm progenitor marker, was hypomethylated and highly expressed in trophoblast lineages. Concordantly, progressive loss of DNA methylation in ENPEP promoter and increased ENPEP expression were detected in trophoblast differentiation. Knockout of ENPEP in hEPSC-em compromised trophoblast differentiation potency, reduced adhesion and invasion of trophoblastic spheroids, and impeded trophoblastic stem cell (TSC) derivation. Importantly, TET2 was involved in the loss of DNA methylation and activation of ENPEP expression during trophoblast differentiation. TET2-null hEPSC-em failed to produce TSC properly. Collectively, our results illustrated the crucial roles of ENPEP and TET2 in trophoblast fate commitments and the unprecedented TET2-mediated loss of DNA methylation in ENPEP promoter.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Developmentally dynamic changes in DNA methylation in the human pancreas

Development of the human pancreas requires the precise temporal control of gene expression via epigenetic mechanisms and the binding of key transcription factors. We quantified genome-wide patterns of DNA methylation in human fetal pancreatic samples from donors aged 6 to 21 post-conception weeks. We found dramatic changes in DNA methylation across pancreas development, with > 21% of sites characterized as developmental differentially methylated positions (dDMPs) including many annotated to genes associated with monogenic diabetes. An analysis of DNA methylation in postnatal pancreas tissue showed that the dramatic temporal changes in DNA methylation occurring in the developing pancreas are largely limited to the prenatal period. Significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small proportion of sites showing sex-specific DNA methylation trajectories across pancreas development. Pancreas dDMPs were not distributed equally across the genome and were depleted in regulatory domains characterized by open chromatin and the binding of known pancreatic development transcription factors. Finally, we compared our pancreas dDMPs to previous findings from the human brain, identifying evidence for tissue-specific developmental changes in DNA methylation. This study represents the first systematic exploration of DNA methylation patterns during human fetal pancreas development and confirms the prenatal period as a time of major epigenomic plasticity.

Accurate structure prediction of biomolecular interactions with AlphaFold 3

The introduction of AlphaFold 21 has spurred a revolution in modelling the structure of proteins and their interactions, enabling a huge range of applications in protein modelling and design2-6. Here we describe our AlphaFold 3 model with a substantially updated diffusion-based architecture that is capable of predicting the joint structure of complexes including proteins, nucleic acids, small molecules, ions and modified residues. The new AlphaFold model demonstrates substantially improved accuracy over many previous specialized tools: far greater accuracy for protein-ligand interactions compared with state-of-the-art docking tools, much higher accuracy for protein-nucleic acid interactions compared with nucleic-acid-specific predictors and substantially higher antibody-antigen prediction accuracy compared with AlphaFold-Multimer v.2.37,8. Together, these results show that high-accuracy modelling across biomolecular space is possible within a single unified deep-learning framework.

Aberrant DNA methylation distorts developmental trajectories in atypical teratoid/rhabdoid tumors

Atypical teratoid/rhabdoid tumors (AT/RTs) are pediatric brain tumors known for their aggressiveness and aberrant but still unresolved epigenetic regulation. To better understand their malignancy, we investigated how AT/RT-specific DNA hypermethylation was associated with gene expression and altered transcription factor binding and how it is linked to upstream regulation. Medulloblastomas, choroid plexus tumors, pluripotent stem cells, and fetal brain were used as references. A part of the genomic regions, which were hypermethylated in AT/RTs similarly as in pluripotent stem cells and demethylated in the fetal brain, were targeted by neural transcriptional regulators. AT/RT-unique DNA hypermethylation was associated with polycomb repressive complex 2 and linked to suppressed genes with a role in neural development and tumorigenesis. Activity of the several NEUROG/NEUROD pioneer factors, which are unable to bind to methylated DNA, was compromised via the suppressed expression or DNA hypermethylation of their target sites, which was also experimentally validated for NEUROD1 in medulloblastomas and AT/RT samples. These results highlight and characterize the role of DNA hypermethylation in AT/RT malignancy and halted neural cell differentiation.

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis

BACKGROUND

The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants.

RESULTS

This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus.

CONCLUSIONS

NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.

Loci cg06256735 and cg15815843 in the MFAP5 gene regulatory regions are hypomethylated in varicose veins apparently due to active demethylation

Varicose vein disease (VVD) is a common health problem worldwide. Microfibril-associated protein 5 (MFAP5) is one of the potential key players in its pathogenesis. Our previous microarray analysis revealed the cg06256735 and cg15815843 loci in the regulatory regions of the MFAP5 gene as hypomethylated in varicose veins which correlated with its up-regulation. The aim of this work was to validate preliminary microarray data, estimate the level of 5-hydroxymethylcytosine (5hmC) at these loci, and determine the methylation status of one of them in different layers of the venous wall. For this, methyl- and hydroxymethyl-sensitive restriction techniques were used followed by real-time PCR and droplet digital PCR, correspondingly, as well as bisulfite pyrosequencing of +/- oxidized DNA. Our microarray data on hypomethylation at the cg06256735 and cg15815843 loci in whole varicose vein segments were confirmed and it was also demonstrated that the level of 5hmC at these loci is increased in VVD. Specifically, among other layers of the venous wall, tunica (t.) intima is the main contributor to hypomethylation at the cg06256735 locus in varicose veins. Thus, it was shown that hypomethylation at the cg06256735 and cg15815843 loci takes place in VVD, with evidence to suggest that it happens through their active demethylation leading to up-regulation of the MFAP5 gene, and t. intima is most involved in this biochemical process.

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Novel genome-wide DNA methylation profiling reveals distinct epigenetic landscape, prognostic model and cellular composition of early-stage lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) has been a leading cause of cancer-related mortality worldwide. Early intervention can significantly improve prognosis. DNA methylation could occur in the early stage of tumor. Comprehensive understanding the epigenetic landscape of early-stage LUAD is crucial in understanding tumorigenesis.

METHODS

Enzymatic methyl sequencing (EM-seq) was performed on 23 tumors and paired normal tissue to reveal distinct epigenetic landscape, for compared with The Cancer Genome Atlas (TCGA) 450K methylation microarray data. Then, an integrative analysis was performed combined with TCGA LUAD RNA-seq data to identify significant differential methylated and expressed genes. Subsequently, the prognostic risk model was constructed and cellular composition was analyzed.

RESULTS

Methylome analysis of EM-seq comparing tumor and normal tissues identified 25 million cytosine-phosphate-guanine (CpG) sites and 30,187 differentially methylated regions (DMR) with a greater number of untraditional types. EM-seq identified a significantly higher number of CpG sites and DMRs compared to the 450K microarray. By integrating the differentially methylated genes (DMGs) with LUAD-related differentially expressed genes (DEGs) from the TCGA database, we constructed prognostic model based on six differentially methylated-expressed genes (MEGs) and verified our prognostic model in GSE13213 and GSE42127 dataset. Finally, cell deconvolution based on the in-house EM-seq methylation profile was used to estimate cellular composition of early-stage LUAD.

CONCLUSIONS

This study firstly delves into novel pattern of epigenomic DNA methylation and provides a multidimensional analysis of the role of DNA methylation revealed by EM-seq in early-stage LUAD, providing distinctive insights into its potential epigenetic mechanisms.

Induced pluripotent stem cells (iPSCs): molecular mechanisms of induction and applications

The induced pluripotent stem cell (iPSC) technology has transformed in vitro research and holds great promise to advance regenerative medicine. iPSCs have the capacity for an almost unlimited expansion, are amenable to genetic engineering, and can be differentiated into most somatic cell types. iPSCs have been widely applied to model human development and diseases, perform drug screening, and develop cell therapies. In this review, we outline key developments in the iPSC field and highlight the immense versatility of the iPSC technology for in vitro modeling and therapeutic applications. We begin by discussing the pivotal discoveries that revealed the potential of a somatic cell nucleus for reprogramming and led to successful generation of iPSCs. We consider the molecular mechanisms and dynamics of somatic cell reprogramming as well as the numerous methods available to induce pluripotency. Subsequently, we discuss various iPSC-based cellular models, from mono-cultures of a single cell type to complex three-dimensional organoids, and how these models can be applied to elucidate the mechanisms of human development and diseases. We use examples of neurological disorders, coronavirus disease 2019 (COVID-19), and cancer to highlight the diversity of disease-specific phenotypes that can be modeled using iPSC-derived cells. We also consider how iPSC-derived cellular models can be used in high-throughput drug screening and drug toxicity studies. Finally, we discuss the process of developing autologous and allogeneic iPSC-based cell therapies and their potential to alleviate human diseases.

Integrative analysis in head and neck cancer reveals distinct role of miRNome and methylome as tumour epigenetic drivers

Head and neck cancer is the sixth most common malignancy worldwide, with the relatively low 5-year survival rate, mainly because it is diagnosed at a late stage. Infection with HPV is a well known aetiology, which affects the nature of these cancers and patients' survival. Besides, it is considered that the main driving force for this type of cancer could be epigenetics. In this study we aimed to find potential epigenetic biomarkers, by integrating miRNome, methylome, and transcriptome analyses. From the fresh head and neck cancer tissue samples, we chose a group for miRNome, methylome and transcriptome profiling, in comparison to adequate control samples. Bioinformatics analyses are performed in R v4.2.2. Count normalisation and group differential expression for mRNA and the previously obtained miRNA count data was performed with DESeq2 v1.36. Gene set enrichment analysis was performed and visualised using gProfiler2 v0.2.1 Identification of miRNA targets was performed by querying in miRTarBase using multiMiR v1.18.0. Annotation of CpG sites merging into islands was obtained from RnBeads.hg19 v1.28.0. package. For the integrative analysis we performed kmeans clustering using stats v4.2.2 package, using 8-12 clusters and nstart 100. We found that transcriptome analysis divides samples into cancers and controls clusters, with no relation to HPV status or cancer anatomical location. Differentially expressed genes (n = 2781) were predominantly associated with signalling pathways of tumour progression. We identified a cluster of genes under the control of the transcription factor E2F that are significantly underexpressed in cancer tissue, as well as T cell immunity genes and genes related to regulation of transcription. Among overexpressed genes in tumours we found those that belong to cell cycle regulation and vasculature. A small number of genes were found significantly differentially expressed in HPV-positive versus HPV-negative tumours (for example NEFH, ZFR2, TAF7L, ZNF541, and TYMS). In this comprehensive study on an overlapping set of samples where the integration of miRNome, methylome and transcriptome analysis were performed for head and neck cancer, we demonstrated that the majority of genes were associated exclusively with miRNome or methylome and, to a lesser extent, under the control of both epigenetic mechanisms.

Analysis of REST binding sites with canonical and non-canonical motifs in human cell lines

BACKGROUND

Repressor element 1 (RE1) silencing transcription factor (REST) is a transcriptional repressor abundantly expressed in aging human brains. It is known to regulate genes associated with oxidative stress, inflammation, and neurological disorders by binding to a canonical form of sequence motif and its non-canonical variations. Although analysis of genomic sequence motifs is crucial to understand transcriptional regulation by transcription factors (TFs), a comprehensive characterization of various forms of RE1 motifs in human cell lines has not been performed.

RESULTS

Here, we analyzed 23 ENCODE REST ChIP-seq datasets from diverse human cell lines and identified a non-redundant set of 68,975 loci with ChIP-seq peaks. Our systematic characterization of these binding sites revealed that the canonical form of REST binding motif was found primarily in ChIP-seq peaks shared across multiple cell lines, while non-canonical forms of motifs were identified in both cell-line-specific binding sites and those shared across cell lines. Remarkably, we observed a notable prevalence of non-canonical motifs that corresponded to half segments of the canonical motif. Furthermore, our analysis unveiled the presence of cell-line-specific REST binding patterns, as evidenced by the clustering of ChIP-seq experiments according to their respective cell lines. This observation underscores the cell-line specificity of REST binding at certain genomic loci, implying intricate cell-line-specific regulatory mechanisms.

CONCLUSIONS

Overall, our study provides a comprehensive characterization of REST binding motifs in human cell lines and genome-wide RE1 motif profiles. These findings contribute to a deeper understanding of REST-mediated transcriptional regulation and highlight the importance of considering cell-line-specific effects in future investigations.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

Integrative Single Cell Atlas Revealed Intratumoral Heterogeneity Generation from an Adaptive Epigenetic Cell State in Human Bladder Urothelial Carcinoma

Intratumor heterogeneity (ITH) of bladder cancer (BLCA) contributes to therapy resistance and immune evasion affecting clinical prognosis. The molecular and cellular mechanisms contributing to BLCA ITH generation remain elusive. It is found that a TM4SF1-positive cancer subpopulation (TPCS) can generate ITH in BLCA, evidenced by integrative single cell atlas analysis. Extensive profiling of the epigenome and transcriptome of all stages of BLCA revealed their evolutionary trajectories. Distinct ancestor cells gave rise to low-grade noninvasive and high-grade invasive BLCA. Epigenome reprograming led to transcriptional heterogeneity in BLCA. During early oncogenesis, epithelial-to-mesenchymal transition generated TPCS. TPCS has stem-cell-like properties and exhibited transcriptional plasticity, priming the development of transcriptionally heterogeneous descendent cell lineages. Moreover, TPCS prevalence in tumor is associated with advanced stage cancer and poor prognosis. The results of this study suggested that bladder cancer interacts with its environment by acquiring a stem cell-like epigenomic landscape, which might generate ITH without additional genetic diversification.

Race-specific coregulatory and transcriptomic profiles associated with DNA methylation and androgen receptor in prostate cancer

BACKGROUND

Prostate cancer is a significant health concern, particularly among African American (AA) men who exhibit higher incidence and mortality compared to European American (EA) men. Understanding the molecular mechanisms underlying these disparities is imperative for enhancing clinical management and achieving better outcomes.

METHODS

Employing a multi-omics approach, we analyzed prostate cancer in both AA and EA men. Using Illumina methylation arrays and RNA sequencing, we investigated DNA methylation and gene expression in tumor and non-tumor prostate tissues. Additionally, Boolean analysis was utilized to unravel complex networks contributing to racial disparities in prostate cancer.

RESULTS

When comparing tumor and adjacent non-tumor prostate tissues, we found that DNA hypermethylated regions are enriched for PRC2/H3K27me3 pathways and EZH2/SUZ12 cofactors. Olfactory/ribosomal pathways and distinct cofactors, including CTCF and KMT2A, were enriched in DNA hypomethylated regions in prostate tumors from AA men. We identified race-specific inverse associations of DNA methylation with expression of several androgen receptor (AR) associated genes, including the GATA family of transcription factors and TRIM63. This suggests that race-specific dysregulation of the AR signaling pathway exists in prostate cancer. To investigate the effect of AR inhibition on race-specific gene expression changes, we generated in-silico patient-specific prostate cancer Boolean networks. Our simulations revealed prolonged AR inhibition causes significant dysregulation of TGF-β, IDH1, and cell cycle pathways specifically in AA prostate cancer. We further quantified global gene expression changes, which revealed differential expression of genes related to microtubules, immune function, and TMPRSS2-fusion pathways, specifically in prostate tumors of AA men. Enrichment of these pathways significantly correlated with an altered risk of disease progression in a race-specific manner.

CONCLUSIONS

Our study reveals unique signaling networks underlying prostate cancer biology in AA and EA men, offering potential insights for clinical management strategies tailored to specific racial groups. Targeting AR and associated pathways could be particularly beneficial in addressing the disparities observed in prostate cancer outcomes in the context of AA and EA men. Further investigation into these identified pathways may lead to the development of personalized therapeutic approaches to improve outcomes for prostate cancer patients across different racial backgrounds.

DNA Methylation in the Adaptive Response to Exercise

Emerging evidence published over the past decade has highlighted the role of DNA methylation in skeletal muscle function and health, including as an epigenetic transducer of the adaptive response to exercise. In this review, we aim to synthesize the latest findings in this field to highlight: (1) the shifting understanding of the genomic localization of altered DNA methylation in response to acute and chronic aerobic and resistance exercise in skeletal muscle (e.g., promoter, gene bodies, enhancers, intergenic regions, un-annotated regions, and genome-wide methylation); (2) how these global/regional methylation changes relate to transcriptional activity following exercise; and (3) the factors (e.g., individual demographic or genetic features, dietary, training history, exercise parameters, local epigenetic characteristics, circulating hormones) demonstrated to alter both the pattern of DNA methylation after exercise, and the relationship between DNA methylation and gene expression. Finally, we discuss the changes in non-CpG methylation and 5-hydroxymethylation after exercise, as well as the importance of emerging single-cell analyses to future studies-areas of increasing focus in the field of epigenetics. We anticipate that this review will help generate a framework for clinicians and researchers to begin developing and testing exercise interventions designed to generate targeted changes in DNA methylation as part of a personalized exercise regimen.

Transcriptional control of CCAAT/enhancer binding protein zeta gene in chicken adipose tissue

CCAAT/enhancer binding protein zeta (C/EBPZ) was differentially expressed in abdominal adipose tissues of fat and lean broilers and regulated adipogenesis in chicken. The objective of this study was to elucidate the transcriptional regulation of C/EBPZ gene in chicken adipose tissue. A 2,031-base pair (bp) chicken C/EBPZ sequence (2,025 nucleotides upstream to 6 nucleotides downstream from the initiator codon, -2,025/+6) was studied. The sequence exhibited a significant promoter activity (P < 0.05) and had some cis-acting elements, notably, a core promoter was identified in nucleotides -94 to +6. Additionally, DNA pull-down assay showed that proteins interacted with chicken C/EBPZ promoter (-173/+6) in preadipocytes were implicated in transcription, post-transcriptional regulation and translation. In addition, KLF2 facilitated the activities of chicken C/EBPZ promoter (-2,025/+6, -1,409/+6, -793/+6, -485/+6, -173/+6, and -94/+6) in preadipocytes (P < 0.05). The expression levels of KLF2 and C/EBPZ in chicken abdominal adipose tissue were substantially associated (r = 0.5978278, P < 0.0001), and KLF2 increased C/EBPZ expression in vitro (P < 0.05). Additionally, chromatin immunoprecipitation (ChIP)-PCR analysis revealed that KLF2 has the ability to interact with the chicken C/EBPZ promoter regions at least at the positions -1,245/-1,048 and -571/-397. Mutation analysis showed that the CGCAGCGCCCG motif located in the chicken C/EBPZ promoter at positions -45 to -35 is involved in regulating transcription and facilitates trans activation by KLF2. These results provided some information of transcription control of C/EBPZ in chicken adipose tissue.

DNA methylation signatures of early-life adversity are exposure-dependent in wild baboons

The early-life environment can profoundly shape the trajectory of an animal's life, even years or decades later. One mechanism proposed to contribute to these early-life effects is DNA methylation. However, the frequency and functional importance of DNA methylation in shaping early-life effects on adult outcomes is poorly understood, especially in natural populations. Here, we integrate prospectively collected data on fitness-associated variation in the early environment with DNA methylation estimates at 477,270 CpG sites in 256 wild baboons. We find highly heterogeneous relationships between the early-life environment and DNA methylation in adulthood: aspects of the environment linked to resource limitation (e.g., low-quality habitat, early-life drought) are associated with many more CpG sites than other types of environmental stressors (e.g., low maternal social status). Sites associated with early resource limitation are enriched in gene bodies and putative enhancers, suggesting they are functionally relevant. Indeed, by deploying a baboon-specific, massively parallel reporter assay, we show that a subset of windows containing these sites are capable of regulatory activity, and that, for 88% of early drought-associated sites in these regulatory windows, enhancer activity is DNA methylation-dependent. Together, our results support the idea that DNA methylation patterns contain a persistent signature of the early-life environment. However, they also indicate that not all environmental exposures leave an equivalent mark and suggest that socioenvironmental variation at the time of sampling is more likely to be functionally important. Thus, multiple mechanisms must converge to explain early-life effects on fitness-related traits.

Structural insights into the recognition of the A/T-rich motif in target gene promoters by the LMX1a homeobox domain

LIM homeodomain transcription factor 1-alpha (LMX1a) is a neuronal lineage-specific transcription activator that plays an essential role during the development of midbrain dopaminergic (mDA) neurons. LMX1a induces the expression of multiple key genes, which ultimately determine the morphology, physiology, and functional identity of mDA neurons. This function of LMX1a is dependent on its homeobox domain. Here, we determined the structures of the LMX1a homeobox domain in complex with the promoter sequences of the Wnt family member 1 (WNT1) or paired like homeodomain 3 (Pitx3) gene, respectively. The complex structures revealed that the LMX1a homeobox domain employed its α3 helix and an N-terminal loop to achieve specific target recognition. The N-terminal loop (loop1) interacted with the minor groove of the double-stranded DNA (dsDNA), whereas the third α-helix (α3) was tightly packed into the major groove of the dsDNA. Structure-based mutations in the α3 helix of the homeobox domain significantly reduced the binding affinity of LMX1a to dsDNA. Moreover, we identified a nonsyndromic hearing loss (NSHL)-related mutation, R199, which yielded a more flexible loop and disturbed the recognition in the minor groove of dsDNA, consistent with the molecular dynamics (MD) simulations. Furthermore, overexpression of Lmx1a promoted the differentiation of SH-SY5Y cells and upregulated the transcription of WNT1 and PITX3 genes. Hence, our work provides a detailed elucidation of the specific recognition between the LMX1a homeobox domain and its specific dsDNA targets, which represents valuable information for future investigations of the functional pathways that are controlled by LMX1a during mDA neuron development.

BNIP3 and DAPK1 methylation in peripheral blood leucocytes are noninvasive biomarkers for gastric cancer

OBJECTIVE

The objective of this study is to comprehensively investigate the potential value of BNIP3 and DAPK1 methylation in peripheral blood leukocytes as a non-invasive biomarker for the detection of gastric cancer (GC), prediction of chemotherapy efficacy, and prognosis assessment.

PATIENTS AND METHODS

Initially, multiple bioinformatic analyses were employed to explore the genetic landscape and biological effects of BNIP3 and DAPK1 in GC tissues. Subsequently, case-control and prospective follow-up studies were conducted to compare the differences in BNIP3 and DAPK1 methylation levels in peripheral blood leukocytes among GC patients and healthy controls, as well as between patients exhibiting sensitivity and resistance to platinum plus fluorouracil treatment, and between patients with varying survival outcomes of GC. Additionally, several predictive nomograms were constructed based on the identified CpG sites and relevant clinical parameters to forecast the occurrence of GC, chemotherapy efficacy, and prognosis.

RESULTS

The upregulation of BNIP3 and DAPK1 was found to be associated with the development and poorer survival outcomes of GC. Furthermore, the expression of BNIP3/DAPK1 exhibited an inverse relationship with their DNA methylation levels and demonstrated a positive correlation with immune cell infiltration, as well as the IC50 values of 5-Fluorouracil and Cisplatin in GC tissues. Increased infiltration of macrophages in the high-expression groups was observed to be linked to unfavorable GC survival. In the case-control and follow-up studies, lower methylation levels of BNIP3 and DAPK1 were identified in the peripheral leukocytes of GC patients compared to healthy controls. Hypomethylation was also associated with more aggressive subtypes, diminished chemotherapy efficacy, and poorer survival outcomes in GC.

CONCLUSION

The DNA methylation of BNIP3 and DAPK1 in peripheral blood leukocytes holds promise as a novel non-invasive biomarker for predicting the occurrence of GC, chemotherapy efficacy, and prognosis assessment.

DNA methylation-dependent and -independent binding of CDX2 directs activation of distinct developmental and homeostatic genes

Precise spatiotemporal and cell type-specific gene expression is essential for proper tissue development and function. Transcription factors (TFs) guide this process by binding to developmental stage-specific targets and establishing an appropriate enhancer landscape. In turn, DNA and chromatin modifications direct the genomic binding of TFs. However, how TFs navigate various chromatin features and selectively bind a small portion of the millions of possible genomic target loci is still not well understood. Here we show that Cdx2 - a pioneer TF that binds distinct targets in developing versus adult intestinal epithelial cells - has a preferential affinity for a non-canonical CpG-containing motif in vivo. A higher frequency of this motif at embryonic and fetal Cdx2 target loci and the specifically methylated state of the CpG during development allows selective Cdx2 binding and activation of developmental enhancers and linked genes. Conversely, demethylation at these enhancers prohibits ectopic Cdx2 binding in adult cells, where Cdx2 binds its canonical motif without a CpG. This differential Cdx2 binding allows for corecruitment of Ctcf and Hnf4, facilitating the establishment of intestinal superenhancers during development and enhancers mediating adult homeostatic functions, respectively. Induced gain of DNA methylation in the adult mouse epithelium or cultured cells causes ectopic recruitment of Cdx2 to the developmental target loci and facilitates cobinding of the partner TFs. Together, our results demonstrate that the differential CpG motif requirements for Cdx2 binding to developmental versus adult target sites allow it to navigate different DNA methylation profiles and activate cell type-specific genes at appropriate times.

Identification of dihydromyricetin as a natural DNA methylation inhibitor with rejuvenating activity in human skin

Changes in DNA methylation patterning have been reported to be a key hallmark of aged human skin. The altered DNA methylation patterns are correlated with deregulated gene expression and impaired tissue functionality, leading to the well-known skin aging phenotype. Searching for small molecules, which correct the aged methylation pattern therefore represents a novel and attractive strategy for the identification of anti-aging compounds. DNMT1 maintains epigenetic information by copying methylation patterns from the parental (methylated) strand to the newly synthesized strand after DNA replication. We hypothesized that a modest inhibition of this process promotes the restoration of the ground-state epigenetic pattern, thereby inducing rejuvenating effects. In this study, we screened a library of 1800 natural substances and 640 FDA-approved drugs and identified the well-known antioxidant and anti-inflammatory molecule dihydromyricetin (DHM) as an inhibitor of the DNA methyltransferase DNMT1. DHM is the active ingredient of several plants with medicinal use and showed robust inhibition of DNMT1 in biochemical assays. We also analyzed the effect of DHM in cultivated keratinocytes by array-based methylation profiling and observed a moderate, but significant global hypomethylation effect upon treatment. To further characterize DHM-induced methylation changes, we used published DNA methylation clocks and newly established age predictors to demonstrate that the DHM-induced methylation change is associated with a reduction in the biological age of the cells. Further studies also revealed re-activation of age-dependently hypermethylated and silenced genes in vivo and a reduction in age-dependent epidermal thinning in a 3-dimensional skin model. Our findings thus establish DHM as an epigenetic inhibitor with rejuvenating effects for aged human skin.

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

Early developmental carry-over effects on exploratory behaviour and DNA methylation in wild great tits (Parus major)

Adverse, postnatal conditions experienced during development are known to induce lingering effects on morphology, behaviour, reproduction and survival. Despite the importance of early developmental stress for shaping the adult phenotype, it is largely unknown which molecular mechanisms allow for the induction and maintenance of such phenotypic effects once the early environmental conditions are released. Here we aimed to investigate whether lasting early developmental phenotypic changes are associated with post-developmental DNA methylation changes. We used a cross-foster and brood size experiment in great tit (Parus major) nestlings, which induced post-fledging effects on biometric measures and exploratory behaviour, a validated personality trait. We investigated whether these post-fledging effects are associated with DNA methylation levels of CpG sites in erythrocyte DNA. Individuals raised in enlarged broods caught up on their developmental delay after reaching independence and became more explorative as days since fledging passed, while the exploratory scores of individuals that were raised in reduced broods remained stable. Although we previously found that brood enlargement hardly affected the pre-fledging methylation levels, we found 420 CpG sites that were differentially methylated between fledged individuals that were raised in small versus large sized broods. A considerable number of the affected CpG sites were located in or near genes involved in metabolism, growth, behaviour and cognition. Since the biological functions of these genes line up with the observed post-fledging phenotypic effects of brood size, our results suggest that DNA methylation provides organisms the opportunity to modulate their condition once the environmental conditions allow it. In conclusion, this study shows that nutritional stress imposed by enlarged brood size during early development associates with variation in DNA methylation later in life. We propose that treatment-associated DNA methylation differences may arise in relation to pre- or post-fledging phenotypic changes, rather than that they are directly induced by the environment during early development.

Epigenetic Reprogramming Drives Epithelial Disruption in Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) remains a major public health challenge that contributes greatly to mortality and morbidity worldwide. Although it has long been recognized that the epithelium is altered in COPD, there has been little focus on targeting it to modify the disease course. Therefore, mechanisms that disrupt epithelial cell function in patients with COPD are poorly understood. In this study, we sought to determine whether epigenetic reprogramming of the cell-cell adhesion molecule E-cadherin, encoded by the CDH1 gene, disrupts epithelial integrity. By reducing these epigenetic marks, we can restore epithelial integrity and rescue alveolar airspace destruction. We used differentiated normal and COPD-derived primary human airway epithelial cells, genetically manipulated mouse tracheal epithelial cells, and mouse and human precision-cut lung slices to assess the effects of epigenetic reprogramming. We show that the loss of CDH1 in COPD is due to increased DNA methylation site at the CDH1 enhancer D through the downregulation of the ten-eleven translocase methylcytosine dioxygenase (TET) enzyme TET1. Increased DNA methylation at the enhancer D region decreases the enrichment of RNA polymerase II binding. Remarkably, treatment of human precision-cut slices derived from patients with COPD with the DNA demethylation agent 5-aza-2'-deoxycytidine decreased cell damage and reduced air space enlargement in the diseased tissue. Here, we present a novel mechanism that targets epigenetic modifications to reverse the tissue remodeling in human COPD lungs and serves as a proof of concept for developing a disease-modifying target.

The potential of epigenetic therapy to target the 3D epigenome in endocrine-resistant breast cancer

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth.

Aspirin Challenge-Induced Genome-Wide DNA Methylation Profile of Peripheral Blood Lymphocytes in Aspirin-Exacerbated Respiratory Disease

Genetic variation and epigenetic factors are thought to contribute to the development of hypersensitivity to aspirin. DNA methylation fluctuates dynamically throughout the day. To discover new CpG methylation in lymphocytes associated with aspirin-exacerbated respiratory disease (AERD), we evaluated changes in global CpG methylation profiles from before to after an oral aspirin challenge in patients with AERD and aspirin-tolerant asthma (ATA). Whole-genome CpG methylation levels of peripheral blood mononuclear cells were quantified with an Illumina 860K Infinium Methylation EPIC BeadChip array and then adjusted for inferred lymphocyte fraction (ILF) with GLINT and Tensor Composition Analysis. Among the 866,091 CpGs in the array, differentially methylated CpGs (DMCs) were found in 6 CpGs in samples from all 12 patients with asthma included in the study (AERD, n = 6; ATA, n = 6). DMCs were found in 3 CpGs in the 6 ATA samples and in 615 CpGs in the 6 AERD samples. A total of 663 DMCs in 415 genes and 214 intergenic regions differed significantly in the AERD compared with the ATA. In promoters, 126 CpG loci were predicted to bind to 38 transcription factors (TFs), many of which were factors already known to be involved in the pathogenesis of asthma and immune responses. In conclusion, we identified 615 new CpGs methylated in peripheral blood lymphocytes by oral aspirin challenge in AERD but not in ATA. These findings indicate that oral aspirin challenge induces epigenetic changes in ILFs, specifically in AERD patients, possibly via changes in TF binding, which may have epigenetic effects on the development of AERD.

DNA methylation dynamics during yak adipocyte differentiation

In mammals, epigenetic modifications involving DNA methylation are necessary for the completion of the cell differentiation process. However, the global DNA methylation landscape and its dynamics during yak adipocyte differentiation remain unexplored. Here, we performed whole-genome bisulfite sequencing (WGBS) to asses DNA methylation in yak preadipocytes and adipocytes, combining these results with those of our previous studies on changes in chromatin accessibility and gene expression during yak adipogenesis. The results showed that CG methylation levels were lower in promoter than in exon and intron, and gradually decreasing from the distal regions to transcription start site (TSS). There was a significant negative correlation between CG methylation levels located in promoter and gene expression levels. The 46 genes shared by CG differentially methylated regions (DMRs) and differential chromatin accessibility were significantly enriched in Hedgehog and PI3K-Akt signaling pathways. ATAC-seq peaks with high chromatin accessibility located in both promoter (≤ 2 kb from TSS) and distal (> 2 kb from TSS) regions corresponded to low methylation levels. Taken together, these findings demonstrated that DNA methylation and its interactions with chromatin accessibility regulate gene expression during yak adipocyte differentiation, contributing to the understanding of mechanisms of various epigenetic modifications and their interactions in adipogenesis.

Structural basis for specific DNA sequence recognition by the transcription factor NFIL3

The CCAAT/enhancer-binding proteins (C/EBPs) constitute a family of pivotal transcription factors involved in tissue development, cellular function, proliferation, and differentiation. NFIL3, as one of them, plays an important role in regulating immune cell differentiation, circadian clock system, and neural regeneration, yet its specific DNA recognition mechanism remains enigmatic. In this study, we showed by the ITC binding experiments that NFIL3 prefers to bind to the TTACGTAA DNA motif. Our structural studies revealed that the α-helical NFIL3 bZIP domain dimerizes through its leucine zipper region, and binds to DNA via its basic region. The two basic regions of the NFIL3 bZIP dimer were pushed apart upon binding to DNA, facilitating the snug accommodation of the two basic regions within the major grooves of the DNA. Remarkably, our binding and structural data also revealed that both NFIL3 and C/EBPα/β demonstrate a shared preference for the TTACGTAA sequence. Furthermore, our study revealed that disease-associated mutations within the NFIL3 bZIP domain result in either reduction or complete disruption of its DNA binding ability. These discoveries not only provide valuable insights into the DNA binding mechanisms of NFIL3 but also elucidate the causal role of NFIL3 mutations in disease pathogenesis.

Transcription factors BZR2/MYC2 modulate brassinosteroid and jasmonic acid crosstalk during pear dormancy

Bud dormancy is an important physiological process during winter. Its release requires a certain period of chilling. In pear (Pyrus pyrifolia), the abscisic acid (ABA)-induced expression of DORMANCY-ASSOCIATED MADS-box (DAM) genes represses bud break, whereas exogenous gibberellin (GA) promotes dormancy release. However, with the exception of ABA and GA, the regulatory effects of phytohormones on dormancy remain largely uncharacterized. In this study, we confirmed brassinosteroids (BRs) and jasmonic acid (JA) contribute to pear bud dormancy release. If chilling accumulation is insufficient, both 24-epibrassinolide (EBR) and methyl jasmonic acid (MeJA) can promote pear bud break, implying that they positively regulate dormancy release. BRASSINAZOLE RESISTANT 2 (BZR2), which is a BR-responsive transcription factor, inhibited PpyDAM3 expression and accelerated pear bud break. The transient overexpression of PpyBZR2 increased endogenous GA, JA, and JA-Ile levels. In addition, the direct interaction between PpyBZR2 and MYELOCYTOMATOSIS 2 (PpyMYC2) enhanced the PpyMYC2-mediated activation of Gibberellin 20-oxidase genes PpyGA20OX1L1 and PpyGA20OX2L2 transcription, thereby increasing GA3 contents and accelerating pear bud dormancy release. Interestingly, treatment with 5 μm MeJA increased the bud break rate, while also enhancing PpyMYC2-activated PpyGA20OX expression and increasing GA3,4 contents. The 100 μm MeJA treatment decreased the PpyMYC2-mediated activation of the PpyGA20OX1L1 and PpyGA20OX2L2 promoters and suppressed the inhibitory effect of PpyBZR2 on PpyDAM3 transcription, ultimately inhibiting pear bud break. In summary, our data provide insights into the crosstalk between the BR and JA signaling pathways that regulate the BZR2/MYC2-mediated pathway in the pear dormancy release process.

Genome-wide methylation profiling of Peripheral T-cell lymphomas identifies TRIP13 as a critical driver of tumor proliferation and survival

Cytosine methylation of genomic DNA contributes to the regulation of gene expression and is involved in normal development including hematopoiesis in mammals. It is catalyzed by the family of DNA methyltransferases (DNMTs) that include DNMT1, DNMT3A, and DNMT3B. Peripheral T-cell lymphomas (PTCLs) represent a diverse group of aggressive mature T-cell malignancies accounting for approximately 10-15% of non-Hodgkin lymphoma cases in the US. PTCLs exhibit a broad spectrum of clinical, histological, and immunophenotypic features with poor prognosis and inadequately understood molecular pathobiology. To better understand the molecular landscape and identify candidate genes involved in disease maintenance, we used high-resolution Whole Genome Bisulfite Sequencing (WGBS) and RNA-seq to profile DNA methylation and gene expression of PTCLs and normal T-cells. We found that the methylation patterns in PTCLs are deregulated and heterogeneous but share 767 hypo- and 567 hypermethylated differentially methylated regions (DMRs) along with 231 genes up- and 91 genes downregulated in all samples suggesting a potential association with tumor development. We further identified 39 hypomethylated promoters associated with increased gene expression in the majority of PTCLs. This putative oncogenic signature included the TRIP13 (thyroid hormone receptor interactor 13) gene whose both genetic and pharmacologic inactivation, inhibited cellular growth of PTCL cell lines by inducing G2-M arrest accompanied by apoptosis suggesting that such an approach might be beneficial in human lymphoma treatment. Altogether we show that human PTCLs are characterized by a large number of recurrent methylation alterations, and demonstrated that TRIP13 is critical for PTCL maintenance in vitro .

Dnmt3b ablation affects fracture repair process by regulating apoptosis

PURPOSE

Previous studies have shown that DNA methyltransferase 3b (Dnmt3b) is the only Dnmt responsive to fracture repair and Dnmt3b ablation in Prx1-positive stem cells and chondrocyte cells both delayed fracture repair. Our study aims to explore the influence of Dnmt3b ablation in Gli1-positive stem cells in fracture healing mice and the underlying mechanism.

METHODS

We generated Gli1-CreERT2; Dnmt3bflox/flox (Dnmt3bGli1ER) mice to operated tibia fracture. Fracture callus tissues of Dnmt3bGli1ER mice and control mice were collected and analyzed by X-ray, micro-CT, biomechanical testing, histopathology and TUNEL assay.

RESULTS

The cartilaginous callus significantly decrease in ablation of Dnmt3b in Gli1-positive stem cells during fracture repair. The chondrogenic and osteogenic indicators (Sox9 and Runx2) in the fracture healing tissues in Dnmt3bGli1ER mice much less than control mice. Dnmt3bGli1ER mice led to delayed bone callus remodeling and decreased biomechanical properties of the newly formed bone during fracture repair. Both the expressions of Caspase-3 and Caspase-8 were upregulated in Dnmt3bGli1ER mice as well as the expressions of BCL-2.

CONCLUSIONS

Our study provides an evidence that Dnmt3b ablation Gli1-positive stem cells can affect fracture healing and lead to poor fracture healing by regulating apoptosis to decrease chondrocyte hypertrophic maturation.

DNA methylation-environment interactions in the human genome

Previously, we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here, we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

Cholestenoic acid as endogenous epigenetic regulator decreases hepatocyte lipid accumulation in vitro and in vivo

Cholestenoic acid (CA) has been reported as an important biomarker of many severe diseases, but its physiological and pathological roles remain unclear. This study aimed to investigate the potential role of CA in hepatic lipid homeostasis. Enzyme kinetic studies revealed that CA specifically activates DNA methyltransferases 1 (DNMT1) at low concentration with EC50 = 1.99 × 10-6 M and inhibits the activity at higher concentration with IC50 = 9.13 × 10-6 M, and specifically inhibits DNMT3a, and DNMT3b activities with IC50= 8.41 × 10-6 M and IC50= 4.89 × 10-6 M, respectively. In a human hepatocyte in vitro model of high glucose (HG)-induced lipid accumulation, CA significantly increased demethylation of 5mCpG in the promoter regions of over 7,000 genes, particularly those involved in master signaling pathways such as calcium-AMPK and 0.0027 at 6 h. RNA sequencing analysis showed that the downregulated genes are affected by CA encoding key enzymes, such as PCSK9, MVK, and HMGCR, which are involved in cholesterol metabolism and steroid biosynthesis pathways. In addition, untargeted lipidomic analysis showed that CA significantly reduced neutral lipid levels by 60% in the cells cultured in high-glucose media. Administration of CA in mouse metabolic dysfunction-associated steatotic liver disease (MASLD) models significantly decreases lipid accumulation, suppresses the gene expression involved in lipid biosynthesis in liver tissues, and alleviates liver function. This study shows that CA as an endogenous epigenetic regulator decreases lipid accumulation via epigenetic regulation. The results indicate that CA can be considered a potential therapeutic target for the treatment of metabolic disorders.NEW & NOTEWORTHY To our knowledge, this study is the first to identify the mitochondrial monohydroxy bile acid cholestenoic acid (CA) as an endogenous epigenetic regulator that regulates lipid metabolism through epigenome modification in human hepatocytes. The methods used in this study are all big data analysis, and the results of each part show the global regulation of CA on human hepatocytes rather than narrow point effects.

Predicting the DNA binding specificity of mutated transcription factors using family-level biophysically interpretable machine learning

Sequence-specific interactions of transcription factors (TFs) with genomic DNA underlie many cellular processes. High-throughput in vitro binding assays coupled with computational analysis have made it possible to accurately define such sequence recognition in a biophysically interpretable yet mechanism-agonistic way for individual TFs. The fact that such sequence-to-affinity models are now available for hundreds of TFs provides new avenues for predicting how the DNA binding specificity of a TF changes when its protein sequence is mutated. To this end, we developed an analytical framework based on a tetrahedron embedding that can be applied at the level of a given structural TF family. Using bHLH as a test case, we demonstrate that we can systematically map dependencies between the protein sequence of a TF and base preference within the DNA binding site. We also develop a regression approach to predict the quantitative energetic impact of mutations in the DNA binding domain of a TF on its DNA binding specificity, and perform SELEX-seq assays on mutated TFs to experimentally validate our results. Our results point to the feasibility of predicting the functional impact of disease mutations and allelic variation in the cell-wide TF repertoire by leveraging high-quality functional information across sets of homologous wild-type proteins.

Comparative analysis of models in predicting the effects of SNPs on TF-DNA binding using large-scale in vitro and in vivo data

Non-coding variants associated with complex traits can alter the motifs of transcription factor (TF)-deoxyribonucleic acid binding. Although many computational models have been developed to predict the effects of non-coding variants on TF binding, their predictive power lacks systematic evaluation. Here we have evaluated 14 different models built on position weight matrices (PWMs), support vector machines, ordinary least squares and deep neural networks (DNNs), using large-scale in vitro (i.e. SNP-SELEX) and in vivo (i.e. allele-specific binding, ASB) TF binding data. Our results show that the accuracy of each model in predicting SNP effects in vitro significantly exceeds that achieved in vivo. For in vitro variant impact prediction, kmer/gkm-based machine learning methods (deltaSVM_HT-SELEX, QBiC-Pred) trained on in vitro datasets exhibit the best performance. For in vivo ASB variant prediction, DNN-based multitask models (DeepSEA, Sei, Enformer) trained on the ChIP-seq dataset exhibit relatively superior performance. Among the PWM-based methods, tRap demonstrates better performance in both in vitro and in vivo evaluations. In addition, we find that TF classes such as basic leucine zipper factors could be predicted more accurately, whereas those such as C2H2 zinc finger factors are predicted less accurately, aligning with the evolutionary conservation of these TF classes. We also underscore the significance of non-sequence factors such as cis-regulatory element type, TF expression, interactions and post-translational modifications in influencing the in vivo predictive performance of TFs. Our research provides valuable insights into selecting prioritization methods for non-coding variants and further optimizing such models.

Pancreatic cancer acquires resistance to MAPK pathway inhibition by clonal expansion and adaptive DNA hypermethylation

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with poor prognosis. It is marked by extraordinary resistance to conventional therapies including chemotherapy and radiation, as well as to essentially all targeted therapies evaluated so far. More than 90% of PDAC cases harbor an activating KRAS mutation. As the most common KRAS variants in PDAC remain undruggable so far, it seemed promising to inhibit a downstream target in the MAPK pathway such as MEK1/2, but up to now preclinical and clinical evaluation of MEK inhibitors (MEKi) failed due to inherent and acquired resistance mechanisms. To gain insights into molecular changes during the formation of resistance to oncogenic MAPK pathway inhibition, we utilized short-term passaged primary tumor cells from ten PDACs of genetically engineered mice. We followed gain and loss of resistance upon MEKi exposure and withdrawal by longitudinal integrative analysis of whole genome sequencing, whole genome bisulfite sequencing, RNA-sequencing and mass spectrometry data.

RESULTS

We found that resistant cell populations under increasing MEKi treatment evolved by the expansion of a single clone but were not a direct consequence of known resistance-conferring mutations. Rather, resistant cells showed adaptive DNA hypermethylation of 209 and hypomethylation of 8 genomic sites, most of which overlap with regulatory elements known to be active in murine PDAC cells. Both DNA methylation changes and MEKi resistance were transient and reversible upon drug withdrawal. Furthermore, MEKi resistance could be reversed by DNA methyltransferase inhibition with remarkable sensitivity exclusively in the resistant cells.

CONCLUSION

Overall, the concept of acquired therapy resistance as a result of the expansion of a single cell clone with epigenetic plasticity sheds light on genetic, epigenetic and phenotypic patterns during evolvement of treatment resistance in a tumor with high adaptive capabilities and provides potential for reversion through epigenetic targeting.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

HOCOMOCO in 2024: a rebuild of the curated collection of binding models for human and mouse transcription factors

We present a major update of the HOCOMOCO collection that provides DNA binding specificity patterns of 949 human transcription factors and 720 mouse orthologs. To make this release, we performed motif discovery in peak sets that originated from 14 183 ChIP-Seq experiments and reads from 2554 HT-SELEX experiments yielding more than 400 thousand candidate motifs. The candidate motifs were annotated according to their similarity to known motifs and the hierarchy of DNA-binding domains of the respective transcription factors. Next, the motifs underwent human expert curation to stratify distinct motif subtypes and remove non-informative patterns and common artifacts. Finally, the curated subset of 100 thousand motifs was supplied to the automated benchmarking to select the best-performing motifs for each transcription factor. The resulting HOCOMOCO v12 core collection contains 1443 verified position weight matrices, including distinct subtypes of DNA binding motifs for particular transcription factors. In addition to the core collection, HOCOMOCO v12 provides motif sets optimized for the recognition of binding sites in vivo and in vitro, and for annotation of regulatory sequence variants. HOCOMOCO is available at https://hocomoco12.autosome.org and https://hocomoco.autosome.org.

MethMotif.Org 2024: a database integrating context-specific transcription factor-binding motifs with DNA methylation patterns

MethMotif (https://methmotif.org) is a publicly available database that provides a comprehensive repository of transcription factor (TF)-binding profiles, enriched with DNA methylation patterns. In this release, we have enhanced the platform, expanding our initial collection to over 700 position weight matrices (PWM), all of which include DNA methylation profiles. One of the key advancements in this release is the segregation of TF-binding motifs based on their cofactors and DNA methylation status. We have previously demonstrated that gene ontology (GO) enriched terms associated with TF target genes may differ based on their association with alternative cofactors and DNA methylation status. MethMotif provides precomputed GO annotations for each human TF of interest, as well as for TF-co-TF complexes, enabling a comprehensive analysis of TF functions in the context of their co-factors. Additionally, MethMotif has been updated to encompass data for two new species, Mus musculus and Arabidopsis thaliana, widening its applicability to a broader community. MethMotif stands out as the first and only TF-binding motifs database to incorporate context-specific PWM coupled with epigenetic information, thereby enlightening context-specific TF functions. This enhancement allows the community to explore and gain deeper insights into the regulatory mechanisms governing transcriptional processes.