Modeling the relationship of epigenetic modifications to transcription factor binding

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.

H2AK119ub dynamics controls hair follicle stem cell quiescence

The transition of stem cells from a quiescent state to an active state is a finely tuned process that requires the dismantling of the quiescence program and the establishment of a cell cycle-promoting transcriptional landscape. Whether epigenetic processes control stem cell states to promote the regeneration of adult tissues remains elusive. In this study, we show that a repressive histone modification, H2AK119ub, is dynamic between quiescent and active hair follicle stem cells (HFSCs) in the adult murine skin. Ablation of H2AK119ub in HFSCs leads to impaired quiescence leading to premature activation and an eventual exhaustion of HFSC pool. Transcriptional and chromatin studies revealed that H2AK119ub directly represses a proliferation promoting transcriptional program in the HFSCs to preserve quiescence. Lastly, we identify that the inhibitory FGF signaling produced by the hair follicle niche keratinocytes maintains H2AK119ub in quiescent HFSCs. Together, these findings reveal that a repressive histone mark, H2AK119ub, is under the dynamic regulation of inhibitory niche signaling to prevent the untimely establishment of an activated state to preserve SC function and longevity.

Schlank orchestrates insect developmental transition by switching H3K27 acetylation to trimethylation in the prothoracic gland

Insect developmental transitions are precisely coordinated by ecdysone and juvenile hormone (JH). We previously revealed that accumulated H3K27 trimethylation (H3K27me3) at the locus encoding JH signal transducer Hairy is involved in the larval-pupal transition in insects, but the underlying mechanism remains to be fully defined. Here, we show in Drosophila and Bombyx that Rpd3-mediated H3K27 deacetylation in the prothoracic gland during the last larval instar promotes ecdysone biosynthesis and the larval-pupal transition by enabling H3K27me3 accumulation at the Hairy locus to induce its transcriptional repression. Importantly, we find that the homeodomain transcription factor Schlank acts to switch active H3K27 acetylation (H3K27ac) to repressive H3K27me3 at the Hairy locus by directly binding to the Hairy promoter and then recruiting the histone deacetylase Rpd3 and the histone methyltransferase PRC2 component Su(z)12 through physical interactions. Moreover, Schlank inhibits Hairy transcription to facilitate the larval-pupal transition, and the Schlank signaling cascade is suppressed by JH but regulated in a positive feedback manner by ecdysone. Together, our data uncover that Schlank mediates epigenetic reprogramming of H3K27 modifications in hormone actions during insect developmental transition.

Deciphering the dynamic code: DNA recognition by transcription factors in the ever-changing genome

Transcription factors (TFs) intricately navigate the vast genomic landscape to locate and bind specific DNA sequences for the regulation of gene expression programs. These interactions occur within a dynamic cellular environment, where both DNA and TF proteins experience continual chemical and structural perturbations, including epigenetic modifications, DNA damage, mechanical stress, and post-translational modifications (PTMs). While many of these factors impact TF-DNA binding interactions, understanding their effects remains challenging and incomplete. This review explores the existing literature on these dynamic changes and their potential impact on TF-DNA interactions.

Dynamic DNA methylation turnover in gene bodies is associated with enhanced gene expression plasticity in plants

BACKGROUND

In several eukaryotes, DNA methylation occurs within the coding regions of many genes, termed gene body methylation (GbM). Whereas the role of DNA methylation on the silencing of transposons and repetitive DNA is well understood, gene body methylation is not associated with transcriptional repression, and its biological importance remains unclear.

RESULTS

We report a newly discovered type of GbM in plants, which is under constitutive addition and removal by dynamic methylation modifiers in all cells, including the germline. Methylation at Dynamic GbM genes is removed by the DRDD demethylation pathway and added by an unknown source of de novo methylation, most likely the maintenance methyltransferase MET1. We show that the Dynamic GbM state is present at homologous genes across divergent lineages spanning over 100 million years, indicating evolutionary conservation. We demonstrate that Dynamic GbM is tightly associated with the presence of a promoter or regulatory chromatin state within the gene body, in contrast to other gene body methylated genes. We find Dynamic GbM is associated with enhanced gene expression plasticity across development and diverse physiological conditions, whereas stably methylated GbM genes exhibit reduced plasticity. Dynamic GbM genes exhibit reduced dynamic range in drdd mutants, indicating a causal link between DNA demethylation and enhanced gene expression plasticity.

CONCLUSIONS

We propose a new model for GbM in regulating gene expression plasticity, including a novel type of GbM in which increased gene expression plasticity is associated with the activity of DNA methylation writers and erasers and the enrichment of a regulatory chromatin state.

Regulation of human ZNF687, a gene associated with Paget's disease of bone

Mutations in Zinc finger 687 (ZNF687) were associated with Paget's disease of bone (PDB), a disease characterized by increased bone resorption and excessive bone formation. It was suggested that ZNF687 plays a role in bone differentiation and development. However, the mechanisms involved in ZNF687 regulation remain unknown. This study aimed to obtain novel knowledge regarding ZNF687 transcriptional and epigenetic regulation. Through in silico analysis, we hypothesized three ZNF687 promoter regions located upstream exon 1 A, 1B, and 1 C and denominated promoter regions 1, 2, and 3, respectively. Their functionality was confirmed by luciferase activity assays and positive/negative regulatory regions were identified using promoter deletions constructs. In silico analysis revealed a high density of CpG islands in these promoter regions and in vitro methylation suppressed promoters' activity. Using bioinformatic approaches, bone-associated transcription factor binding sites containing CpG dinucleotides were identified, including those for NFκB, PU.1, DLX5, and SOX9. By co-transfection in HEK293 and hFOB cells, we found that DLX5 specifically activated ZNF687 promoter region 1, and its methylation impaired DLX5-driven promoter stimulation. NFκB repressed and activated promoter regions 1 and 2, respectively, and these activities were affected by methylation. PU.1 induced ZNF687 promoter region 1 which was affected by methylation. SOX9 differentially regulated ZNF687 promoters in HEK293 and hFOB cells that were impaired after methylation. In conclusion, this study provides novel insights into ZNF687 regulation by demonstrating that NFκB, PU.1, DLX5, and SOX9 are regulators of ZNF687 promoters, and DNA methylation influences their activity. The contribution of the dysregulation of these mechanisms in PDB should be further elucidated.

Integrated Analysis of Gene Expression and Methylation Data to Identify Potential Biomarkers Related to Atherosclerosis Onset

Atherosclerosis is a kind of chronic inflammatory cardiovascular disease. Epigenetic regulation plays a crucial role in atherosclerosis. Our study was aimed at finding potential biomarkers associated with the occurrence of atherosclerosis. Two datasets were downloaded from the Gene Expression Omnibus (GEO) database. The epigenome-wide association study (EWAS) analysis was performed on methylation data using CpGassoc package. The differential expression analysis was conducted on mRNA data using limma package. The GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) functional enrichment was done in clusterProfiler package. Finally, the logistic regression model was constructed using generalized linear model (glm) function. Between atherosclerotic vs. nonatherosclerotic samples, totally 4980 cytosine-phosphate-guanine (CpG) sites (annotated to 2860 genes) and 132 differentially expressed genes (DEGs) related to atherosclerosis were identified. The annotated 2860 genes and 132 DEGs were significantly enriched in 9 and 4 KEGG pathways and 289 and 132 GO terms, respectively. After cross-analysis, 6 crucial CpG sites were screened to build the model, including cg01187920, cg03422911, cg08018825, cg10967350, cg14473924, and cg25313204. The diagnostic model could reliably separate the atherosclerosis samples from nonatherosclerotic samples. In conclusion, the 6 CpG sites are probably potential diagnostic biomarkers for atherosclerosis, including cg01187920, cg03422911, cg08018825, cg10967350, cg14473924, and cg25313204.

GoPeaks: histone modification peak calling for CUT&Tag

Genome-wide mapping of histone modifications is critical to understanding transcriptional regulation. CUT&Tag is a new method for profiling histone modifications, offering improved sensitivity and decreased cost compared with ChIP-seq. Here, we present GoPeaks, a peak calling method specifically designed for histone modification CUT&Tag data. We compare the performance of GoPeaks against commonly used peak calling algorithms to detect histone modifications that display a range of peak profiles and are frequently used in epigenetic studies. We find that GoPeaks robustly detects genome-wide histone modifications and, notably, identifies a substantial number of H3K27ac peaks with improved sensitivity compared to other standard algorithms.

An Update in Epigenetics in Metabolic-Associated Fatty Liver Disease

Metabolic-associated fatty liver disease (MAFLD) is characterized by hepatic steatosis accompanied by one of three features: overweight or obesity, T2DM, or lean or normal weight with evidence of metabolic dysregulation. It is distinguished by excessive fat accumulation in hepatocytes, and a decrease in the liver's ability to oxidize fats, the accumulation of ectopic fat, and the activation of proinflammatory pathways. Chronic damage will keep this pathophysiologic cycle active causing progression from hepatic steatosis to cirrhosis and eventually, hepatocarcinoma. Epigenetics affecting gene expression without altering DNA sequence allows us to study MAFLD pathophysiology from a different perspective, in which DNA methylation processes, histone modifications, and miRNAs expression have been closely associated with MAFLD progression. However, these considerations also faced us with the circumstance that modifying those epigenetics patterns might lead to MAFLD regression. Currently, epigenetics is an area of great interest because it could provide new insights in therapeutic targets and non-invasive biomarkers. This review comprises an update on the role of epigenetic patterns, as well as innovative therapeutic targets and biomarkers in MAFLD.

CRISPR activation of endogenous genes reprograms fibroblasts into cardiovascular progenitor cells for myocardial infarction therapy

Fibroblasts can be reprogrammed into cardiovascular progenitor cells (CPCs) using transgenic approaches, although the underlying mechanism remains unclear. We determined whether activation of endogenous genes such as Gata4, Nkx2.5, and Tbx5 can rapidly establish autoregulatory loops and initiate CPC generation in adult extracardiac fibroblasts using a CRISPR activation system. The induced fibroblasts (>80%) showed phenotypic changes as indicated by an Nkx2.5 cardiac enhancer reporter. The progenitor characteristics were confirmed by colony formation and expression of cardiovascular genes. Cardiac sphere induction segregated the early and late reprogrammed cells that can generate functional cardiomyocytes and vascular cells in vitro. Therefore, they were termed CRISPR-induced CPCs (ciCPCs). Transcriptomic analysis showed that cell cycle and heart development pathways were important to accelerate CPC formation during the early reprogramming stage. The CRISPR system opened the silenced chromatin locus, thereby allowing transcriptional factors to access their own promoters and eventually forming a positive feedback loop. The regenerative potential of ciCPCs was assessed after implantation in mouse myocardial infarction models. The engrafted ciCPCs differentiated into cardiovascular cells in vivo but also significantly improved contractile function and scar formation. In conclusion, multiplex gene activation was sufficient to drive CPC reprogramming, providing a new cell source for regenerative therapeutics.

Epigenetic Regulation of Vascular Smooth Muscle Cell Phenotype Switching in Atherosclerotic Artery Remodeling: A Mini-Review

Atherosclerosis is a chronic inflammatory disease characterized by extensive remodeling of medium and large-sized arteries. Inward remodeling (=lumen shrinkage) of the vascular walls is the underlying cause for ischemia in target organs. Therefore, inward remodeling can be considered the predominant feature of atherosclerotic pathology. Outward remodeling (=lumen enlargement) is a physiological response compensating for lumen shrinkage caused by neointimal hyperplasia, but as a pathological response to changes in blood flow, outward remodeling leads to substantial arterial wall thinning. Thinned vascular walls are prone to rupture, and subsequent thrombus formation accounts for the majority of acute cardiovascular events. Pathological remodeling is driven by inflammatory cells which induce vascular smooth muscle cells to switch from quiescent to a proliferative and migratory phenotype. After decades of intensive research, the molecular mechanisms of arterial remodeling are starting to unfold. In this mini-review, we summarize the current knowledge of the epigenetic and transcriptional regulation of vascular smooth muscle cell phenotype switching from the contractile to the synthetic phenotype involved in arterial remodeling and discuss potential therapeutic options.

STAT1 epigenetically regulates LCP2 and TNFAIP2 by recruiting EP300 to contribute to the pathogenesis of inflammatory bowel disease

Background

The aetiology of inflammatory bowel disease (IBD) is related to genetics and epigenetics. Epigenetic regulation of the pathogenesis of IBD has not been well defined. Here, we investigated the role of H3K27ac events in the pathogenesis of IBD. Based on previous ChIP-seq and RNA-seq assays, we studied signal transducer and activator of transcription 1 (STAT1) as a transcription factor (TF) and investigated whether the STAT1–EP300–H3K27ac axis contributes to the development of IBD. We performed ChIP-PCR to investigate the interaction between STAT1 and H3K27ac, and co-IP assays were performed to investigate the crosstalk between STAT1 and EP300.

Results

Lymphocyte cytosolic protein 2 (LCP2) and TNF-α‐inducible protein 2 (TNFAIP2) are target genes of STAT1. p-STAT1 binds to the enhancer loci of the two genes where H3K27ac is enriched, and EP300 subsequently binds to regulate their expression. In mice with dextran sulfate sodium (DSS)-induced acute colitis, an EP300 inhibitor significantly inhibited colitis.

Conclusions

p-STAT1 and EP300 promote TNFAIP2 and LCP2 expression through an increase in H3K27ac enrichment on their enhancers and contribute to the pathogenesis of chronic inflammation.

Graphic abstract

PCSK9 genetic (rs11591147) and epigenetic (DNA methylation) modifications associated with PCSK9 expression and serum proteins in CAD patients

INTRODUCTION

Proprotein convertase subtilisin/kexin type 9 (PCSK9) genetic polymorphisms play a significant role in cholesterol homeostasis. Therefore, we aimed to investigate the association of PCSK9 genetic variations NM_174936.3:c.137G>T (R46L, rs11591147) and NM_174936.3:c.1120G>T (D374Y, rs137852912), as well as promoter DNA methylation status, with mRNA expression and circulating serum protein levels in coronary artery disease (CAD) patients.

METHODS

The present study includes 300 CAD cases and 300 controls from South India. Biochemical assays were performed using commercially available kits. PCSK9 rs11591147 and rs137852912 polymorphisms were analyzed by the polymerase chain reaction (PCR)-restriction fragment length polymorphism method, whereas promoter DNA methylation status and gene expression were determined using methylation specific PCR and quantitative PCR respectively.

RESULTS

The genotypic distribution of PCSK9 rs11591147 revealed that individuals with the TT-genotype and T-allele have a reduced risk for CAD. Furthermore, patients with the PCSK9 rs11591147 TT genotype have a significantly lower total cholesterol and low-density lipoprotein-cholesterol levels and also higher high-density lipoprotein-cholesterol levels than individuals with the GG genotype. Logistic regression analysis has shown that the GG and GT (p = 1.51 × 10^-8 , p = 1.47 × 10^-9 ) genotypes predicted the risk for CAD with an odds ratio of 5.8 and 7.3 respectively. In addition, individuals with the TT genotype were hypermethylated at promoter DNA of PCSK9, resulting in lower mRNA expression and circulating serum proteins than in individuals with the GG genotype. In silico analyses revealed that rs11591147 T-allele has protein destabilizing capacity.

CONCLUSIONS

In conclusion, the present study indicates that the PCSK9 gene expression and circulating serum protein levels are not only associated with rs11591147 genotype, but also with promoter DNA methylation. Furthermore, the findings with respect to both single nucleotide polymorphism and promoter DNA methylation may open avenues for novel treatment possibilities targeting PCSK9 for CAD management.

CpG content-dependent associations between transcription factors and histone modifications

Understanding the factors that underlie the epigenetic regulation of genes is crucial to understand the gene regulatory machinery as a whole. Several experimental and computational studies examined the relationship between different factors involved. Here we investigate the relationship between transcription factors (TFs) and histone modifications (HMs), based on ChIP-seq data in cell lines. As it was shown that gene regulation by TFs differs depending on the CpG class of a promoter, we study the impact of the CpG content in promoters on the associations between TFs and HMs. We suggest an approach based on sparse linear regression models to infer associations between TFs and HMs with respect to CpG content. A study of the partial correlation of HMs for the two classes of high and low CpG content reveals possible CpG dependence and potential candidates for confounding factors in our models. We show that the models are accurate, inferred associations reflect known biological relationships, and we give new insight into associations with respect to CpG content. Moreover, analysis of a ChIP-seq dataset in HepG2 cells of the HM H3K122ac, an HM about little is known, reveals novel TF associations and supports a previously established link to active transcription.

intePareto: an R package for integrative analyses of RNA-Seq and ChIP-Seq data

BACKGROUND

RNA-Seq, the high-throughput sequencing (HT-Seq) of mRNAs, has become an essential tool for characterizing gene expression differences between different cell types and conditions. Gene expression is regulated by several mechanisms, including epigenetically by post-translational histone modifications which can be assessed by ChIP-Seq (Chromatin Immuno-Precipitation Sequencing). As more and more biological samples are analyzed by the combination of ChIP-Seq and RNA-Seq, the integrated analysis of the corresponding data sets becomes, theoretically, a unique option to study gene regulation. However, technically such analyses are still in their infancy.

RESULTS

Here we introduce intePareto, a computational tool for the integrative analysis of RNA-Seq and ChIP-Seq data. With intePareto we match RNA-Seq and ChIP-Seq data at the level of genes, perform differential expression analysis between biological conditions, and prioritize genes with consistent changes in RNA-Seq and ChIP-Seq data using Pareto optimization.

CONCLUSION

intePareto facilitates comprehensive understanding of high dimensional transcriptomic and epigenomic data. Its superiority to a naive differential gene expression analysis with RNA-Seq and available integrative approach is demonstrated by analyzing a public dataset.

Intronic regions of the human coagulation factor VIII gene harboring transcription factor binding sites with a strong bias towards the short-interspersed elements

Increasing data show that intronic derived regulatory elements, such as transcription factor binding sites (TFBs), play key roles in gene regulation, and malfunction. Accordingly, characterizing the sequence context of the intronic regions of the human coagulation factor VIII (hFVIII) gene can be important. In this study, the intronic regions of the hFVIII gene were scrutinized based on in-silico methods. The results disclosed that these regions harbor a rich array of functional elements such as repetitive elements (REs), splicing sites, and transcription factor binding sites (TFBs). Among these elements, TFBs and REs showed a significant distribution and correlation to each other. This survey indicated that 31% of TFBs are localized in the intronic regions of the gene. Moreover, TFBs indicate a strong bias in the regions far from splice sites of introns with mapping to different REs. Accordingly, TFBs showed highly bias toward Short Interspersed Elements (SINEs), which in turn they covering about 12% of the total of REs. However, the distribution pattern of TFBs-REs showed different bias in the intronic regions, spatially into the Introns 13 and 25. The rich array of SINE-TFBs and CR1-TFBs were situated within 5′UTR of the gene that may be an important driving force for regulatory innovation of the hFVIII gene. Taken together, these data may lead to revealing intronic regions with the capacity to renewing gene regulatory networks of the hFVIII gene. On the other hand, these correlations might provide the novel idea for a new hypothesis of molecular evolution of the FVIII gene, and treatment of Hemophilia A which should be considered in future studies.

Hypermethylation of Auxin-Responsive Motifs in the Promoters of the Transcription Factor Genes Accompanies the Somatic Embryogenesis Induction in Arabidopsis

The auxin-induced embryogenic reprogramming of plant somatic cells is associated with extensive modulation of the gene expression in which epigenetic modifications, including DNA methylation, seem to play a crucial role. However, the function of DNA methylation, including the role of auxin in epigenetic regulation of the SE-controlling genes, remains poorly understood. Hence, in the present study, we analysed the expression and methylation of the TF genes that play a critical regulatory role during SE induction (LEC1, LEC2, BBM, WUS and AGL15) in auxin-treated explants of Arabidopsis. The results showed that auxin treatment substantially affected both the expression and methylation patterns of the SE-involved TF genes in a concentration-dependent manner. The auxin treatment differentially modulated the methylation of the promoter (P) and gene body (GB) sequences of the SE-involved genes. Relevantly, the SE-effective auxin treatment (5.0 µM of 2,4-D) was associated with the stable hypermethylation of the P regions of the SE-involved genes and a significantly higher methylation of the P than the GB fragments was a characteristic feature of the embryogenic culture. The presence of auxin-responsive (AuxRE) motifs in the hypermethylated P regions suggests that auxin might substantially contribute to the DNA methylation-mediated control of the SE-involved genes.

Prediction of transcription factors binding events based on epigenetic modifications in different human cells

Aim: We aim to predict transcription factor (TF) binding events from knowledge of gene expression and epigenetic modifications. Materials & methods: TF-binding events based on the Encode project and The Cancer Genome Atlas data were analyzed by the random forest method. Results: We showed the high performance of TF-binding predictive models in GM12878, HeLa, HepG2 and K562 cell lines and applied them to other cell lines and tissues. The genes bound by the top TFs (MAX and MAZ) were significantly associated with cancer-related processes such as cell proliferation and DNA repair. Conclusion: We successfully constructed TF-binding predictive models in cell lines and applied them in tissues.

Integrative Functional Annotation of 52 Genetic Loci Influencing Myocardial Mass Identifies Candidate Regulatory Variants and Target Genes

BACKGROUND

Regulatory elements may be involved in the mechanisms by which 52 loci influence myocardial mass, reflected by abnormal amplitude and duration of the QRS complex on the ECG. Functional annotation thus far did not take into account how these elements are affected in disease context.

METHODS

We generated maps of regulatory elements on hypertrophic cardiomyopathy patients (ChIP-seq N=14 and RNA-seq N=11) and nondiseased hearts (ChIP-seq N=4 and RNA-seq N=11). We tested enrichment of QRS-associated loci on elements differentially acetylated and directly regulating differentially expressed genes between hypertrophic cardiomyopathy patients and controls. We further performed functional annotation on QRS-associated loci using these maps of differentially active regulatory elements.

RESULTS

Regions differentially affected in disease showed a stronger enrichment ( P=8.6×10^-5) for QRS-associated variants than those not showing differential activity ( P=0.01). Promoters of genes differentially regulated between hypertrophic cardiomyopathy patients and controls showed more enrichment ( P=0.001) than differentially acetylated enhancers ( P=0.8) and super-enhancers ( P=0.025). We also identified 74 potential causal variants overlapping these differential regulatory elements. Eighteen of the genes mapped confirmed previous findings, now also pinpointing the potentially affected regulatory elements and candidate causal variants. Fourteen new genes were also mapped.

CONCLUSIONS

Our results suggest differentially active regulatory elements between hypertrophic cardiomyopathy patients and controls can offer more insights into the mechanisms of QRS-associated loci than elements not affected by disease.

Comprehensive functional annotation of susceptibility variants associated with asthma

Genome-wide association studies (GWAS) have identified hundreds of primarily non-coding disease-susceptibility variants that further need functional interpretation to prioritize and discriminate the disease-relevant variants. We present a comprehensive genome-wide non-coding variant prioritization scheme followed by validation using Pyrosequencing and TaqMan assays in asthma. We implemented a composite Functional Annotation Score (cFAS) to investigate over 32,000 variants consisting of 1525 GWAS-lead asthma-susceptibility variants and their LD proxies (r² ≥ 0.80). Functional annotation pipeline in cFAS revealed 274 variants with significant score at 1% false discovery rate. This study implicates a novel locus 4p16 (SLC26A1) with eQTL variant (rs11936407) and known loci in 17q12-21 and 5q22 which encode ORM1-like protein 3 (ORMDL3, rs406527, and rs12936231) and thymic stromal lymphopoietin (TSLP, rs3806932 and rs10073816) epithelial gene, respectively. Follow-up validation analysis through pyrosequencing of CpG sites in and nearby rs4065275 and rs11936407 showed genotype-dependent hypomethylation on asthma cases compared with healthy controls. Prioritized variants are enriched for asthma-specific histone modification associated with active chromatin (H3K4me1 and H3K27ac) in T cells, B cells, lung, and immune-related interferon gamma signaling pathways. Our findings, together with those from prior studies, suggest that SNPs can affect asthma by regulating enhancer activity, and our comprehensive bioinformatics and functional analysis could lead to biological insights into asthma pathogenesis.Graphic abstract.

Systematical identification of cell-specificity of CTCF-gene binding based on epigenetic modifications

The CCCTC-binding factor (CTCF) mediates transcriptional regulation and implicates epigenetic modifications in cancers. However, the systematically unveiling inverse regulatory relationship between CTCF and epigenetic modifications still remains unclear, especially the mechanism by which histone modification mediates CTCF binding. Here, we developed a systematic approach to investigate how epigenetic changes affect CTCF binding. Through integration analysis of CTCF binding in 30 cell lines, we concluded that CTCF generally binds with higher intensity in normal cell lines than that in cancers, and higher intensity in genome regions closed to transcription start sites. To facilitate the better understanding of their associations, we constructed linear mixed-effect models to analyze the effects of the epigenetic modifications on CTCF binding in four cancer cell lines and six normal cell lines, and identified seven epigenetic modifications as potential epigenetic patterns that influence CTCF binding intensity in promoter regions and six epigenetic modifications in enhancer regions. Further analysis of the effects in different locations revealed that the epigenetic regulation of CTCF binding was location-specific and cancer cell line-specific. Moreover, H3K4me2 and H3K9ac showed the potential association with immune regulation of disease. Taken together, our method can contribute to improve the understanding of the epigenetic regulation of CTCF binding and provide potential therapeutic targets for treating tumors associated with CTCF.

A computational method to predict topologically associating domain boundaries combining histone Marks and sequence information

Background

The three-dimensional (3D) structure of chromatins plays significant roles during cell differentiation and development. Hi-C and other 3C-based technologies allow us to look deep into the chromatin architectures. Many studies have suggested that topologically associating domains (TAD), as the structure and functional unit, are conserved across different organs. However, our understanding about the underlying mechanism of the TAD boundary formation is still limited.

Results

We developed a computational method, TAD–Lactuca, to infer this structure by taking the contextual information of the epigenetic modification signals and the primary DNA sequence information on the genome. TAD–Lactuca is found stable in the case of multi-resolutions and different datasets. It could achieve high accuracy and even outperforms the state-of-art methods when the sequence patterns were incorporated. Moreover, several transcript factor binding motifs, besides the well-known CCCTC-binding factor (CTCF) motif, were found significantly enriched on the boundaries.

Conclusions

We provided a low cost, effective method to predict TAD boundaries. Above results suggested the incorporation of sequence features could significantly improve the performance. The sequence motif enrichment analysis indicates several gene regulation motifs around the boundaries, which is consistent with TADs may serve as the functional units of gene regulation and implies the sequence patterns would be important in chromatin folding.

Individual DNA Methylation Profile is Correlated with Age and can be Targeted to Modulate Healthy Aging and Longevity

:

Patterns of DNA methylation, the best characterized epigenetic modification, are modulated by aging. In humans, different studies at both site-specific and genome-wide levels have reported that modifications of DNA methylation are associated with the chronological aging process but also with the quality of aging (or biological aging), providing new perspectives for establishing powerful biomarkers of aging.

:

In this article, the role of DNA methylation in aging and longevity has been reviewed by analysing literature data about DNA methylation variations occurring during the lifetime in response to environmental factors and genetic background, and their association with the aging process and, in particular, with the quality of aging. Special attention has been devoted to the relationship between nuclear DNA methylation patterns, mitochondrial DNA epigenetic modifications, and longevity. Mitochondrial DNA has recently been reported to modulate global DNA methylation levels of the nuclear genome during the lifetime, and, in spite of the previous belief, it has been found to be the target of methylation modifications.

:

Analysis of DNA methylation profiles across lifetime shows that a remodeling of the methylome occurs with age and/or with age-related decline. Thus, it can be an excellent biomarker of aging and of the individual decline and frailty status. The knowledge about the mechanisms underlying these modifications is crucial since it might allow the opportunity for targeted treatment to modulate the rate of aging and longevity.

Labelled regulatory elements are pervasive features of the macrophage genome and are dynamically utilized by classical and alternative polarization signals

The concept of tissue-specific gene expression posits that lineage-determining transcription factors (LDTFs) determine the open chromatin profile of a cell via collaborative binding, providing molecular beacons to signal-dependent transcription factors (SDTFs). However, the guiding principles of LDTF binding, chromatin accessibility and enhancer activity have not yet been systematically evaluated. We sought to study these features of the macrophage genome by the combination of experimental (ChIP-seq, ATAC-seq and GRO-seq) and computational approaches. We show that Random Forest and Support Vector Regression machine learning methods can accurately predict chromatin accessibility using the binding patterns of the LDTF PU.1 and four other key TFs of macrophages (IRF8, JUNB, CEBPA and RUNX1). Any of these TFs alone were not sufficient to predict open chromatin, indicating that TF binding is widespread at closed or weakly opened chromatin regions. Analysis of the PU.1 cistrome revealed that two-thirds of PU.1 binding occurs at low accessible chromatin. We termed these sites labelled regulatory elements (LREs), which may represent a dormant state of a future enhancer and contribute to macrophage cellular plasticity. Collectively, our work demonstrates the existence of LREs occupied by various key TFs, regulating specific gene expression programs triggered by divergent macrophage polarizing stimuli.

Transcription Factors That Govern Development and Disease: An Achilles Heel in Cancer

Development requires the careful orchestration of several biological events in order to create any structure and, eventually, to build an entire organism. On the other hand, the fate transformation of terminally differentiated cells is a consequence of erroneous development, and ultimately leads to cancer. In this review, we elaborate how development and cancer share several biological processes, including molecular controls. Transcription factors (TF) are at the helm of both these processes, among many others, and are evolutionarily conserved, ranging from yeast to humans. Here, we discuss four families of TFs that play a pivotal role and have been studied extensively in both embryonic development and cancer-high mobility group box (HMG), GATA, paired box (PAX) and basic helix-loop-helix (bHLH) in the context of their role in development, cancer, and their conservation across several species. Finally, we review TFs as possible therapeutic targets for cancer and reflect on the importance of natural resistance against cancer in certain organisms, yielding knowledge regarding TF function and cancer biology.

Genetic Contexts Characterize Dynamic Histone Modification Patterns Among Cell Types

Histone modifications play critical roles in mammalian development, regulating chromatin structure and gene expression. Dynamic histone modifications among cell types have been shown to associate with changes in mammalian development. However, how to quantitatively measure the histone modification alterations and how histone modifications vary across cell types under different genetic contexts remain largely unexplored and whether these changes are related to the primary DNA sequence remains limited. Here, we employed an entropy-based method to measure histone modification alterations in six definite genomic regions across five cell types and identified lineage-specific histone modification genes. We observed that histone modification alterations prefer to enrich in 5'-UTR exons, and also in 3'-UTR exons and its downstream. Then we built a model to predict the histone modification patterns from the primary DNA sequence. We found that the frequencies of k-mer sequence compositions are predictive of histone modification patterns, suggesting that the primary DNA sequence correlated with the histone modification alterations among cell types. Additionally, the lineage-specific histone modification genes display a higher conservation and lower GC-content. Together, we performed a systematic analysis for histone modification alterations and demonstrated how to identify genomic region-specific elements of epigenetic and genetic regulation and histone modification patterns across different cell types.

Epigenetic Reprogramming of TGF-β Signaling in Breast Cancer

The Transforming Growth Factor-β (TGF-β) signaling pathway has a well-documented, context-dependent role in breast cancer development. In normal and premalignant cells, it acts as a tumor suppressor. By contrast, during the malignant phases of breast cancer progression, the TGF-β signaling pathway elicits tumor promoting effects particularly by driving the epithelial to mesenchymal transition (EMT), which enhances tumor cell migration, invasion and ultimately metastasis to distant organs. The molecular and cellular mechanisms that govern this dual capacity are being uncovered at multiple molecular levels. This review will focus on recent advances relating to how epigenetic changes such as acetylation and methylation control the outcome of TGF-β signaling and alter the fate of breast cancer cells. In addition, we will highlight how this knowledge can be further exploited to curb tumorigenesis by selective targeting of the TGF-β signaling pathway.

An Integrated Regulatory Network Based on Comprehensive Analysis of mRNA Expression, Gene Methylation and Expression of Long Non-coding RNAs (lncRNAs) in Myelodysplastic Syndromes

Myelodysplastic syndromes (MDS) are a heterogeneous group of disorders characterized by ineffective hematopoiesis, defective differentiation of hematopoietic precursors, and expansion of the abnormal clones. The prevalence of MDS has raised great concerns worldwide, but its pathogenetic mechanisms remain elusive. To provide insights on novel biomarkers for the diagnosis and therapy of MDS, we performed high-throughput genome-wide mRNA expression profiling, DNA methylation analysis, and long non-coding RNAs (lncRNA) analysis on bone marrows from four MDS patients and four age-matched healthy controls. We identified 1,937 differentially expressed genes (DEGs), 515 methylated genes, and 214 lncRNA that showed statistically significant differences. As the most significant module-related DEGs, TCL1A, PTGS2, and MME were revealed to be enriched in regulation of cell differentiation and cell death pathways. In addition, the GeneGo pathway maps identified by top DEGs were shown to converge on cancer, immunoregulation, apoptosis and regulation of actin cytoskeleton, most of which are known contributors in MDS etiology and pathogenesis. Notably, as potential biomarkers for diagnosis of MDS, four specific genes (ABAT, FADD, DAPP1, and SMPD3) were further subjected to detailed pathway analysis. Our integrative analysis on mRNA expression, gene methylation and lncRNAs profiling facilitates further understanding of the pathogenesis of MDS, and may promote the diagnosis and novel therapeutics for this disease.

Histone modifications underlie monocyte dysregulation in patients with systemic sclerosis, underlining the treatment potential of epigenetic targeting

BACKGROUND AND OBJECTIVE

Systemic sclerosis (SSc) is a severe autoimmune disease, in which the pathogenesis is dependent on both genetic and epigenetic factors. Altered gene expression in SSc monocytes, particularly of interferon (IFN)-responsive genes, suggests their involvement in SSc development. We investigated the correlation between epigenetic histone marks and gene expression in SSc monocytes.

METHODS

Chromatin immunoprecipitation followed by sequencing (ChIPseq) for histone marks H3K4me3 and H3K27ac was performed on monocytes of nine healthy controls and 14 patients with SSc. RNA sequencing was performed in parallel to identify aberrantly expressed genes and their correlation with the levels of H3K4me3 and H3K27ac located nearby their transcription start sites. ChIP-qPCR assays were used to verify the role of bromodomain proteins, H3K27ac and STATs on IFN-responsive gene expression.

RESULTS

1046 and 534 genomic loci showed aberrant H3K4me3 and H3K27ac marks, respectively, in SSc monocytes. The expression of 381 genes was directly and significantly proportional to the levels of such chromatin marks present near their transcription start site. Genes correlated to altered histone marks were enriched for immune, IFN and antiviral pathways and presented with recurrent binding sites for IRF and STAT transcription factors at their promoters. IFNα induced the binding of STAT1 and STAT2 at the promoter of two of these genes, while blocking acetylation readers using the bromodomain BET family inhibitor JQ1 suppressed their expression.

CONCLUSION

SSc monocytes have altered chromatin marks correlating with their IFN signature. Enzymes modulating these reversible marks may provide interesting therapeutic targets to restore monocyte homeostasis to treat or even prevent SSc.

Identifying Transcriptional Regulatory Modules Among Different Chromatin States in Mouse Neural Stem Cells

Gene expression regulation is a complex process involving the interplay between transcription factors and chromatin states. Significant progress has been made toward understanding the impact of chromatin states on gene expression. Nevertheless, the mechanism of transcription factors binding combinatorially in different chromatin states to enable selective regulation of gene expression remains an interesting research area. We introduce a nonparametric Bayesian clustering method for inhomogeneous Poisson processes to detect heterogeneous binding patterns of multiple proteins including transcription factors to form regulatory modules in different chromatin states. We applied this approach on ChIP-seq data for mouse neural stem cells containing 21 proteins and observed different groups or modules of proteins clustered within different chromatin states. These chromatin-state-specific regulatory modules were found to have significant influence on gene expression. We also observed different motif preferences for certain TFs between different chromatin states. Our results reveal a degree of interdependency between chromatin states and combinatorial binding of proteins in the complex transcriptional regulatory process. The software package is available on Github at - https://github.com/BSharmi/DPM-LGCP.

Revealing transcription factor and histone modification co-localization and dynamics across cell lines by integrating ChIP-seq and RNA-seq data

Background

Interactions among transcription factors (TFs) and histone modifications (HMs) play an important role in the precise regulation of gene expression. The context specificity of those interactions and further its dynamics in normal and disease remains largely unknown. Recent development in genomics technology enables transcription profiling by RNA-seq and protein’s binding profiling by ChIP-seq. Integrative analysis of the two types of data allows us to investigate TFs and HMs interactions both from the genome co-localization and downstream target gene expression.

Results

We propose a integrative pipeline to explore the co-localization of 55 TFs and 11 HMs and its dynamics in human GM12878 and K562 by matched ChIP-seq and RNA-seq data from ENCODE. We classify TFs and HMs into three types based on their binding enrichment around transcription start site (TSS). Then a set of statistical indexes are proposed to characterize the TF-TF and TF-HM co-localizations. We found that Rad21, SMC3, and CTCF co-localized across five cell lines. High resolution Hi-C data in GM12878 shows that they associate most of the Hi-C peak loci with a specific CTCF-motif “anchor” and supports that CTCF, SMC3, and RAD2 co-localization serves important role in 3D chromatin structure. Meanwhile, 17 TF-TF pairs are highly dynamic between GM12878 and K562. We then build SVM models to correlate high and low expression level of target genes with TF binding and HM strength. We found that H3k9ac, H3k27ac, and three TFs (ELF1, TAF1, and POL2) are predictive with the accuracy about 85~92%.

Conclusion

We propose a pipeline to analyze the co-localization of TF and HM and their dynamics across cell lines from ChIP-seq, and investigate their regulatory potency by RNA-seq. The integrative analysis of two level data reveals new insight for the cooperation of TFs and HMs and is helpful in understanding cell line specificity of TF/HM interactions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5278-5) contains supplementary material, which is available to authorized users.

Cancer epigenetics and the potential of epigenetic drugs for treating solid tumors

Introduction: Epigenetic modification without DNA sequence mutation plays an important role in cancer development. Some small molecular inhibitors targeting key epigenetic molecules have been approved by the Food and Drug Administration to treat hematological malignancies. However, the anticancer effects of these drugs on solid tumors are not satisfactory, and the mechanisms of action remain largely unknown. Areas covered: The review summarizes the latest research on cancer epigenetics and discusses the potentials and limitations of using epigenetic drugs to treat solid tumors. An analysis of possible reasons for epigenetic drug treatment failure in solid tumors in some clinical trials is discussed along with prospects for future development. Expert commentary: Next-generation small molecule inhibitors will target novel epigenetic regulators with high cancer specificity. Combined modalities exploiting epigenetic drugs with chemo-/radiotherapy, molecular-targeting drugs, and immunotherapy will be able to effectively treat solid tumors in the near future.

Alternative splicing links histone modifications to stem cell fate decision

BACKGROUND

Understanding the embryonic stem cell (ESC) fate decision between self-renewal and proper differentiation is important for developmental biology and regenerative medicine. Attention has focused on mechanisms involving histone modifications, alternative pre-messenger RNA splicing, and cell-cycle progression. However, their intricate interrelations and joint contributions to ESC fate decision remain unclear.

RESULTS

We analyze the transcriptomes and epigenomes of human ESC and five types of differentiated cells. We identify thousands of alternatively spliced exons and reveal their development and lineage-dependent characterizations. Several histone modifications show dynamic changes in alternatively spliced exons and three are strongly associated with 52.8% of alternative splicing events upon hESC differentiation. The histone modification-associated alternatively spliced genes predominantly function in G2/M phases and ATM/ATR-mediated DNA damage response pathway for cell differentiation, whereas other alternatively spliced genes are enriched in the G1 phase and pathways for self-renewal. These results imply a potential epigenetic mechanism by which some histone modifications contribute to ESC fate decision through the regulation of alternative splicing in specific pathways and cell-cycle genes. Supported by experimental validations and extended datasets from Roadmap/ENCODE projects, we exemplify this mechanism by a cell-cycle-related transcription factor, PBX1, which regulates the pluripotency regulatory network by binding to NANOG. We suggest that the isoform switch from PBX1a to PBX1b links H3K36me3 to hESC fate determination through the PSIP1/SRSF1 adaptor, which results in the exon skipping of PBX1.

CONCLUSION

We reveal the mechanism by which alternative splicing links histone modifications to stem cell fate decision.

Targeting Transcription Factors for Cancer Treatment

Transcription factors are involved in a large number of human diseases such as cancers for which they account for about 20% of all oncogenes identified so far. For long time, with the exception of ligand-inducible nuclear receptors, transcription factors were considered as “undruggable” targets. Advances knowledge of these transcription factors, in terms of structure, function (expression, degradation, interaction with co-factors and other proteins) and the dynamics of their mode of binding to DNA has changed this postulate and paved the way for new therapies targeted against transcription factors. Here, we discuss various ways to target transcription factors in cancer models: by modulating their expression or degradation, by blocking protein/protein interactions, by targeting the transcription factor itself to prevent its DNA binding either through a binding pocket or at the DNA-interacting site, some of these inhibitors being currently used or evaluated for cancer treatment. Such different targeting of transcription factors by small molecules is facilitated by modern chemistry developing a wide variety of original molecules designed to specifically abort transcription factor and by an increased knowledge of their pathological implication through the use of new technologies in order to make it possible to improve therapeutic control of transcription factor oncogenic functions.

Genome and epigenome analysis of monozygotic twins discordant for congenital heart disease

BACKGROUND

Congenital heart disease (CHD) is the leading non-infectious cause of death in infants. Monozygotic (MZ) twins share nearly all of their genetic variants before and after birth. Nevertheless, MZ twins are sometimes discordant for common complex diseases. The goal of this study is to identify genomic and epigenomic differences between a pair of twins discordant for a form of congenital heart disease, double outlet right ventricle (DORV).

RESULTS

A monoamniotic monozygotic (MZ) twin pair discordant for DORV were subjected to genome-wide sequencing and methylation analysis. We identified few genomic differences but 1566 differentially methylated regions (DMRs) between the MZ twins. Twenty percent (312/1566) of the DMRs are located within 2 kb upstream of transcription start sites (TSS), containing 121 binding sites of transcription factors. Particularly, ZIC3 and NR2F2 are found to have hypermethylated promoters in both the diseased twin and additional patients suffering from DORV.

CONCLUSIONS

The results showed a high correlation between hypermethylated promoters at ZIC3 and NR2F2 and down-regulated gene expression levels of these two genes in patients with DORV compared to normal controls, providing new insight into the potential mechanism of this rare form of CHD.

Annotating function to differentially expressed LincRNAs in myelodysplastic syndrome using a network-based method

Motivation

Long non-coding RNAs (lncRNAs) have been implicated in the regulation of diverse biological functions. The number of newly identified lncRNAs has increased dramatically in recent years but their expression and function have not yet been described from most diseases. To elucidate lncRNA function in human disease, we have developed a novel network based method (NLCFA) integrating correlations between lncRNA, protein coding genes and noncoding miRNAs. We have also integrated target gene associations and protein-protein interactions and designed our model to provide information on the combined influence of mRNAs, lncRNAs and miRNAs on cellular signal transduction networks.

Results

We have generated lncRNA expression profiles from the CD34+ haematopoietic stem and progenitor cells (HSPCs) from patients with Myelodysplastic syndromes (MDS) and healthy donors. We report, for the first time, aberrantly expressed lncRNAs in MDS and further prioritize biologically relevant lncRNAs using the NLCFA. Taken together, our data suggests that aberrant levels of specific lncRNAs are intimately involved in network modules that control multiple cancer-associated signalling pathways and cellular processes. Importantly, our method can be applied to prioritize aberrantly expressed lncRNAs for functional validation in other diseases and biological contexts.

Availability and implementation

The method is implemented in R language and Matlab.

Contact

xizhou@wakehealth.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Homocysteine induces vascular inflammatory response via SMAD7 hypermethylation in human umbilical vein smooth muscle cells

Homocysteine (Hcy) can induce atherosclerosis through the inflammatory response and DNA methylation disorder. Our recent study has reported a novel epigenetic modified gene related to atherosclerosis -SMAD7. To further understand the pathogenesis of atherosclerosis, the current study was designed to investigate an inflammatory role of Hcy in human umbilical vein smooth muscle cells (HUVSMCs) through interfering with SMAD7 methylation. Using MALDI-TOF MS, we found that Hcy increased DNA methylation levels of SMAD7 promoter in a dose and time-dependent manner in HUVSMCs. Meanwhile, both SMAD7 mRNA and protein levels were decreased along with the increase of Hcy concentrations and treating time. Decreased SMAD7 levels led to up regulation of pro-inflammatory cytokines (TNF-α and IL-1β) expression in HUVSMCs. Furthermore, we found that activation of NF-κB pathway was the mechanism by which reduced Smad7 levels enhanced vascular inflammation. Thus, Hcy is able to activate NF-κB-mediated vascular inflammatory response via inducing hypermethylation of SMAD7 promoter in HUVSMCs. The in vitro findings supplement our recent clinical study that SMAD7 methylation as a novel marker in atherosclerosis and further elucidate the role of Hcy in atherogenesis.

Recognition of the long range enhancer-promoter interactions by further adding DNA structure properties and transcription factor binding motifs in human cell lines

The enhancer-promoter interactions (EPIs) with strong tissue-specificity play an important role in cis-regulatory mechanism of human cell lines. However, it still remains a challenging work to predict these interactions so far. Due to that these interactions are regulated by the cooperativeness of diverse functional genomic signatures, DNA spatial structure and DNA sequence elements. In this paper, by adding DNA structure properties and transcription factor binding motifs, we presented an improved computational method to predict EPIs in human cell lines. In comparison with the results of other group on the same datasets, our best accuracies by cross-validation test were about 15%-24% higher in the same cell lines, and the accuracies by independent test were about 11%-15% higher in new cell lines. Meanwhile, we found that transcription factor binding motifs and DNA structure properties have important information that would largely determine long range EPIs prediction. From the distribution comparisons, we also found their distinct differences between interacting and non-interacting sets in each cell line. Then, the correlation analysis and network models for relationships among top-ranked functional genomic signatures indicated that diverse genomic signatures would cooperatively establish a complex regulatory network to facilitate long range EPIs. The experimental results provided additional insights about the roles of DNA intrinsic properties and functional genomic signatures in EPIs prediction.

Transcriptional activator DOT1L putatively regulates human embryonic stem cell differentiation into the cardiac lineage

Background

Commitment of pluripotent stem cells into differentiated cells and associated gene expression necessitate specific epigenetic mechanisms that modify the DNA and corresponding histone proteins to render the chromatin in an open or closed state. This in turn dictates the associated genetic machinery, including transcription factors, acknowledging the cellular signals provided. Activating histone methyltransferases represent crucial enzymes in the epigenetic machinery that cause transcription initiation by delivering the methyl mark on histone proteins. A number of studies have evidenced the vital role of one such histone modifier, DOT1L, in transcriptional regulation. Involvement of DOT1L in differentiating pluripotent human embryonic stem (hES) cells into the cardiac lineage has not yet been investigated.

Methods

The study was conducted on in-house derived (KIND1) and commercially available (HES3) human embryonic stem cell lines. Chromatin immunoprecipitation (ChIP) was performed followed by sequencing to uncover the cardiac genes harboring the DOT1L specific mark H3K79me2. Following this, dual immunofluorescence was employed to show the DOT1L co-occupancy along with the cardiac progenitor specific marker. DOT1L was knocked down by siRNA to further confirm its role during cardiac differentiation.

Results

ChIP sequencing revealed a significant number of peaks characterizing H3K79me2 occupancy in the proximity of the transcription start site. This included genes like MYOF, NR2F2, NKX2.5, and HAND1 in cardiac progenitors and cardiomyocytes, and POU5F1 and NANOG in pluripotent hES cells. Consistent with this observation, we also show that DOT1L co-localizes with the master cardiac transcription factor NKX2.5, suggesting its direct involvement during gene activation. Knockdown of DOT1L did not alter the pluripotency of hES cells, but it led to the disruption of cardiac differentiation observed morphologically as well as at transcript and protein levels.

Conclusions

Collectively, our data suggests the crucial role of H3K79me2 methyltransferase DOT1L for activation of NKX2.5 during the cardiac differentiation of hES cells.

Electronic supplementary material

The online version of this article (10.1186/s13287-018-0810-8) contains supplementary material, which is available to authorized users.

Relationship between histone modifications and transcription factor binding is protein family specific

The very small fraction of putative binding sites (BSs) that are occupied by transcription factors (TFs) in vivo can be highly variable across different cell types. This observation has been partly attributed to changes in chromatin accessibility and histone modification (HM) patterns surrounding BSs. Previous studies focusing on BSs within DNA regulatory regions found correlations between HM patterns and TF binding specificities. However, a mechanistic understanding of TF-DNA binding specificity determinants is still not available. The ability to predict in vivo TF binding on a genome-wide scale requires the identification of features that determine TF binding based on evolutionary relationships of DNA binding proteins. To reveal protein family-dependent mechanisms of TF binding, we conducted comprehensive comparisons of HM patterns surrounding BSs and non-BSs with exactly matched core motifs for TFs in three cell lines: 33 TFs in GM12878, 37 TFs in K562, and 18 TFs in H1-hESC. These TFs displayed protein family-specific preferences for HM patterns surrounding BSs, with high agreement among cell lines. Moreover, compared to models based on DNA sequence and shape at flanking regions of BSs, HM-augmented quantitative machine-learning methods resulted in increased performance in a TF family-specific manner. Analysis of the relative importance of features in these models indicated that TFs, displaying larger HM pattern differences between BSs and non-BSs, bound DNA in an HM-specific manner on a protein family-specific basis. We propose that TF family-specific HM preferences reveal distinct mechanisms that assist in guiding TFs to their cognate BSs by altering chromatin structure and accessibility.

MEF2C transcription factor is associated with the genetic and epigenetic risk architecture of schizophrenia and improves cognition in mice

Large-scale consortia mapping the genomic risk architectures of schizophrenia provide vast amounts of molecular information, with largely unexplored therapeutic potential. We harnessed publically available information from the Psychiatric Genomics Consortium, and report myocyte enhancer factor 2C (MEF2C) motif enrichment in sequences surrounding the top scoring single-nucleotide polymorphisms within risk loci contributing by individual small effect to disease heritability. Chromatin profiling at base-pair resolution in neuronal nucleosomes extracted from prefrontal cortex of 34 subjects, including 17 cases diagnosed with schizophrenia, revealed MEF2C motif enrichment within cis-regulatory sequences, including neuron-specific promoters and superenhancers, affected by histone H3K4 hypermethylation in disease cases. Vector-induced short- and long-term Mef2c upregulation in mouse prefrontal projection neurons consistently resulted in enhanced cognitive performance in working memory and object recognition paradigms at baseline and after psychotogenic drug challenge, in conjunction with remodeling of local connectivity. Neuronal genome tagging in vivo by Mef2c-Dam adenine methyltransferase fusion protein confirmed the link between cognitive enhancement and MEF2C occupancy at promoters harboring canonical and variant MEF2C motifs. The multilayered integrative approaches presented here provide a roadmap to uncover the therapeutic potential of transcriptional regulators for schizophrenia and related disorders.

Carbon Catabolite Repression in Filamentous Fungi

Carbon Catabolite Repression (CCR) has fascinated scientists and researchers around the globe for the past few decades. This important mechanism allows preferential utilization of an energy-efficient and readily available carbon source over relatively less easily accessible carbon sources. This mechanism helps microorganisms to obtain maximum amount of glucose in order to keep pace with their metabolism. Microorganisms assimilate glucose and highly favorable sugars before switching to less-favored sources of carbon such as organic acids and alcohols. In CCR of filamentous fungi, CreA acts as a transcription factor, which is regulated to some extent by ubiquitination. CreD-HulA ubiquitination ligase complex helps in CreA ubiquitination, while CreB-CreC deubiquitination (DUB) complex removes ubiquitin from CreA, which causes its activation. CCR of fungi also involves some very crucial elements such as Hexokinases, cAMP, Protein Kinase (PKA), Ras proteins, G protein-coupled receptor (GPCR), Adenylate cyclase, RcoA and SnfA. Thorough study of molecular mechanism of CCR is important for understanding growth, conidiation, virulence and survival of filamentous fungi. This review is a comprehensive revision of the regulation of CCR in filamentous fungi as well as an updated summary of key regulators, regulation of different CCR-dependent mechanisms and its impact on various physical characteristics of filamentous fungi.

Deep learning of the splicing (epi)genetic code reveals a novel candidate mechanism linking histone modifications to ESC fate decision

Alternative splicing (AS) is a genetically and epigenetically regulated pre-mRNA processing to increase transcriptome and proteome diversity. Comprehensively decoding these regulatory mechanisms holds promise in getting deeper insights into a variety of biological contexts involving in AS, such as development and diseases. We assembled splicing (epi)genetic code, DeepCode, for human embryonic stem cell (hESC) differentiation by integrating heterogeneous features of genomic sequences, 16 histone modifications with a multi-label deep neural network. With the advantages of epigenetic features, DeepCode significantly improves the performance in predicting the splicing patterns and their changes during hESC differentiation. Meanwhile, DeepCode reveals the superiority of epigenomic features and their dominant roles in decoding AS patterns, highlighting the necessity of including the epigenetic properties when assembling a more comprehensive splicing code. Moreover, DeepCode allows the robust predictions across cell lineages and datasets. Especially, we identified a putative H3K36me3-regulated AS event leading to a nonsense-mediated mRNA decay of BARD1. Reduced BARD1 expression results in the attenuation of ATM/ATR signalling activities and further the hESC differentiation. These results suggest a novel candidate mechanism linking histone modifications to hESC fate decision. In addition, when trained in different contexts, DeepCode can be expanded to a variety of biological and biomedical fields.

Oxidative Stress Alters the Profile of Transcription Factors Related to Early Development on In Vitro Produced Embryos

High oxygen levels during in vitro culture (IVC) can induce oxidative stress through accumulation of reactive oxygen species (ROS), negatively affecting embryo development. This study evaluated the effect of different O₂ tensions during IVC on bovine blastocyst development and transcriptional status, considering transcription factors that play an essential role during early embryo development. For this purpose, embryos were produced in vitro by conventional protocols and cultured in two different oxygen tensions, physiological (5%) and atmospheric (20%). Expanded blastocysts were subjected to transcript quantitation analysis by RT-qPCR with Biomark™ HD System (Fluidigm, US), using 67 TaqMan assays specific for Bos taurus. Differences were observed in genes related to oxidation-reduction processes, DNA-dependent transcription factors, and factors related to important functional pathways for embryo development. Blastocyst rate was higher in the 5% O₂ group and the number of cells was assessed, with the 5% O₂ group having a higher number of cells. ROS concentration was evaluated, with a higher ROS presence in the 20% O₂ group. Taken together, these results allow us to conclude that IVC of embryos at atmospheric O₂ tension affects the expression of important transcription factors involved in multiple cell biology pathways that can affect embryo development, quality, and viability.

Mendelian Randomization Analysis Identifies CpG Sites as Putative Mediators for Genetic Influences on Cardiovascular Disease Risk

The extent to which genetic influences on cardiovascular disease risk are mediated by changes in DNA methylation levels has not been systematically explored. We developed an analytical framework that integrates genetic fine mapping and Mendelian randomization with epigenome-wide association studies to evaluate the causal relationships between methylation levels and 14 cardiovascular disease traits. We identified ten genetic loci known to influence proximal DNA methylation which were also associated with cardiovascular traits after multiple-testing correction. Bivariate fine mapping provided evidence that the individual variants responsible for the observed effects on cardiovascular traits at the ADCY3 and ADIPOQ loci were potentially mediated through changes in DNA methylation, although we highlight that we are unable to reliably separate causality from horizontal pleiotropy. Estimates of causal effects were replicated with results from large-scale consortia. Genetic variants and CpG sites identified in this study were enriched for histone mark peaks in relevant tissue types and gene promoter regions. Integrating our results with expression quantitative trait loci data, we provide evidence that variation at these regulatory regions is likely to also influence gene expression levels at these loci.

Plant-Specific Histone Deacetylases HDT1/2 Regulate GIBBERELLIN 2-OXIDASE2 Expression to Control Arabidopsis Root Meristem Cell Number

Root growth is modulated by environmental factors and depends on cell production in the root meristem (RM). New cells in the meristem are generated by stem cells and transit-amplifying cells, which together determine RM cell number. Transcription factors and chromatin-remodeling factors have been implicated in regulating the switch from stem cells to transit-amplifying cells. Here, we show that two Arabidopsis thaliana paralogs encoding plant-specific histone deacetylases, HDT1 and HDT2, regulate a second switch from transit-amplifying cells to expanding cells. Knockdown of HDT1/2 (hdt1,2i) results in an earlier switch and causes a reduced RM cell number. Our data show that HDT1/2 negatively regulate the acetylation level of the C₁₉-GIBBERELLIN 2-OXIDASE2 (GA2ox2) locus and repress the expression of GA2ox2 in the RM and elongation zone. Overexpression of GA2ox2 in the RM phenocopies the hdt1,2i phenotype. Conversely, knockout of GA2ox2 partially rescues the root growth defect of hdt1,2i These results suggest that by repressing the expression of GA2ox2, HDT1/2 likely fine-tune gibberellin metabolism and they are crucial for regulating the switch from cell division to expansion to determine RM cell number. We propose that HDT1/2 function as part of a mechanism that modulates root growth in response to environmental factors.