Core promoter structure and genomic context reflect histone 3 lysine 9 acetylation patterns

Global research trends on epigenetics and neuropathic pain: A bibliometric analysis

Objective

Neuropathic pain (NP) is a common disease that manifests with pathological changes in the somatosensory system. In recent years, the interactions of NP with the epigenetic mechanism have been increasingly elucidated. However, only a few studies have used bibliometric tools to systematically analyze knowledge in this field. The objective of this study is to visually analyze the trends, hotspots, and frontiers in epigenetics and NP research by using a bibliometric method.

Methods

Studies related to epigenetics and NP were searched from the Science Citation Index-Expanded of the Web of Science Core Collection database. Search time is from inception to November 30, 2022. No restrictions were placed on language. Only articles and reviews were included as document types. Data on institutions, countries, authors, journal distribution, and keywords were imported into CiteSpace software for visual analysis.

Results

A total of 867 publications met the inclusion criteria, which spanned the period from 2000 to 2022. Over the years, the number of publications and the frequency of citations exhibited a clear upward trend in general, reaching a peak in 2021. The major contributing countries in terms of the number of publications were China, the United States, and Japan. The top three institutions were Rutgers State University, Xuzhou Medical University, and Nanjing Medical University. Molecular Pain, Pain, and Journal of Neuroinflammation contributed significantly to the volume of issues. Among the top 10 authors in terms of the number of publications, Tao Yuan-Xiang contributed 30 entries, followed by Zhang Yi with 24 and Wu Shao-Gen with 20. On the basis of the burst and clusters of keywords, "DNA methylation," "Circular RNA," "acetylation," "long non-coding RNA," and "microglia" are global hotspots in the field.

Conclusion

The bibliometric analysis indicates that the number of publications related to epigenetics and NP is exhibiting a rapid increase. Keyword analysis shows that "DNA methylation," "Circular RNA," "acetylation," "long non-coding RNA" and "microglia" are the most interesting terms for researchers in the field. More rigorous clinical trials and additional studies that explore relevant mechanisms are required in the future.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

DNA methylation dynamics of long noncoding RNA during human fetal development

Aim: To determine whether the promoters of long noncoding RNAs (lncRNAs) undergo dynamic changes in DNA methylation during fetal development. Methods: ANOVA and the tissue specificity index were used to identify and validate tissue-specific methylation sites. Age-associated DNA methylation signatures were identified by applying the elastic net method. Results: The lncRNA methylome landscape was characterized in four types of fetal tissue and at three gestational time points, and specific characteristics relative to the tissue of origin and developmental age were identified. Higher levels of lncRNA methylation might be involved in tissue differentiation. LncRNAs harboring age-associated methylation signatures may participate in the fetal developmental process. Conclusion: This study provides novel insights into the role of lncRNA methylomes in fetal tissue specification and development.

Calcitonin gene-related peptide induces the histone H3 lysine 9 acetylation in astrocytes associated with neuroinflammation in rats with neuropathic pain

Aims

Calcitonin gene‐related peptide (CGRP) as a regulator of astrocyte activation may facilitate spinal nociceptive processing. Histone H3 lysine 9 acetylation (H3K9ac) is considered an important regulator of cytokine and chemokine gene expression after peripheral nerve injury. In this study, we explored the relationship between CGRP and H3K9ac in the activation of astrocytes, and elucidated the underlying mechanisms in the pathogenesis of chronic neuropathic pain.

Methods

Astroglial cells (C6) were treated with CGRP and differentially enrichments of H3K9ac on gene promoters were examined using ChIP‐seq. A chronic constriction injury (CCI) rat model was used to evaluate the role of CGRP on astrocyte activation and H3K9ac signaling in CCI‐induced neuropathic pain. Specific inhibitors were employed to delineate the involved signaling.

Results

Intrathecal injection of CGRP and CCI increased the number of astrocytes displaying H3K9ac in the spinal dorsal horn of rats. Treatment of CGRP was able to up‐regulate H3K9ac and glial fibrillary acidic protein (GFAP) expression in astroglial cells. ChIP‐seq data indicated that CGRP significantly altered H3K9ac enrichments on gene promoters in astroglial cells following CGRP treatment, including 151 gaining H3K9ac and 111 losing this mark, which mostly enriched in proliferation, autophagy, and macrophage chemotaxis processes. qRT‐PCR verified expressions of representative candidate genes (ATG12, ATG4C, CX3CR1, MMP28, MTMR14, HMOX1, RET) and RTCA verified astrocyte proliferation. Additionally, CGRP treatment increased the expression of H3K9ac, CX3CR1, and IL‐1β in the spinal dorsal horn. CGRP antagonist and HAT inhibitor attenuated mechanical and thermal hyperalgesia in CCI rats. Such analgesic effects were concurrently associated with the reduced levels of H3K9ac, CX3CR1, and IL‐1β in the spinal dorsal horn of CCI rats.

Conclusion

Our findings highly indicate that CGRP is associated with the development of neuropathic pain through astrocytes‐mediated neuroinflammatory responses via H3K9ac in spinal dorsa horn following nerve injury.

This study found that CGRP act on their astrocytic receptors and lead to H3K9 acetylation (H3K9ac), which are mainly associated with proliferation‐, autophagy‐, and inflammation‐related gene expression. The number of astrocytes with H3K9ac expression is increased after nerve injury. Inhibition of CGRP attenuates the development of neuropathic pain, which was accompanied by the suppression of H3K9ac, CX3CR1, and IL‐1β expression in CCI rats.

Aberrant gene expression induced by a high fat diet is linked to H3K9 acetylation in the promoter-proximal region

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disease, with an estimated global prevalence of 1 in 4 individuals. Aberrant transcriptional control of gene expression is central to the pathophysiology of metabolic diseases. However, the molecular mechanisms leading to gene dysregulation are not well understood. Histone modifications play important roles in the control of transcription. Acetylation of histone 3 at lysine 9 (H3K9ac) is associated with transcriptional activity and is implicated in transcript elongation by controlling RNA polymerase II (RNAPII) pause-release. Hence, changes in this histone modification may shed information on novel pathways linking transcription control and metabolic dysfunction. Here, we carried out genome-wide analysis of H3K9ac in the liver of mice fed a control or a high-fat diet (an animal model of NAFLD), and asked whether this histone mark associates with changes in gene expression. We found that over 70% of RNAPII peaks in promoter-proximal regions overlapped with H3K9ac, consistent with a role of H3K9ac in the regulation of transcription. When comparing high-fat with control diet, approximately 17% of the differentially expressed genes were associated with changes in H3K9ac in their promoters, showing a strong correlation between changes in H3K9ac signal and gene expression. Overall, our data indicate that in response to a high-fat diet, dysregulated gene expression of a subset of genes may be attributable to changes in transcription elongation driven by H3K9ac. Our results point at an added mechanism of gene regulation that may be important in the development of metabolic diseases.

Expression of clock gene Dbp in omental and mesenteric adipose tissue in patients with type 2 diabetes

INTRODUCTION

We previously reported in ob/ob mice, one of animal models of human type 2 diabetes mellitus (DM2), that (i) acetylation of histone H3 lysine 9 (H3K9) at the promoter region of clock gene Dbp and DBP mRNA expression are reduced in epididymal adipose tissue, (ii) binding of DBP to the promoter region of peroxisome proliferator-activated receptor (Ppar)-γ and mRNA expression of PPAR-γ1sv were decreased in preadipocytes and (iii) adiponectin secretion was decreased, leading to the impaired insulin sensitivity.

RESEARCH DESIGN AND METHODS

The present study was undertaken to evaluate whether such the changes in visceral adipose tissue were detected in patients with DM2. We obtained omental and mesenteric adipose tissue during surgery of lymph node dissection for gastric and colorectal cancers, and investigated these variables in adipose tissue (omental from gastric cancer; 13 non-DM, 12 DM2: mesenteric from colorectal cancer; 12 non-DM, 11 DM2).

RESULTS

Acetylation of histone H3K9 at the promoter region of Dbp and DBP mRNA expression in omental, but not in mesenteric adipose tissue were significantly lower in DM2 than in patients without DM. PPAR-γ mRNA expression in omental adipose tissue was also lower in patients with DM2, but not in mesenteric adipose tissue.

CONCLUSIONS

The changes in DBP-PPAR-γ axis observed in mice with diabetes were also detected in patients with DM2. Because adiponectin secretion is reported to be enhanced through the PPAR-γ-related mechanism, this study supports the hypothesis that omental adipose tissue is involved in the mechanism of DM2.

Bayesian integrative analysis of epigenomic and transcriptomic data identifies Alzheimer's disease candidate genes and networks

Biomedical research studies have generated large multi-omic datasets to study complex diseases like Alzheimer’s disease (AD). An important aim of these studies is the identification of candidate genes that demonstrate congruent disease-related alterations across the different data types measured by the study. We developed a new method to detect such candidate genes in large multi-omic case-control studies that measure multiple data types in the same set of samples. The method is based on a gene-centric integrative coefficient quantifying to what degree consistent differences are observed in the different data types. For statistical inference, a Bayesian hierarchical model is used to study the distribution of the integrative coefficient. The model employs a conditional autoregressive prior to integrate a functional gene network and to share information between genes known to be functionally related. We applied the method to an AD dataset consisting of histone acetylation, DNA methylation, and RNA transcription data from human cortical tissue samples of 233 subjects, and we detected 816 genes with consistent differences between persons with AD and controls. The findings were validated in protein data and in RNA transcription data from two independent AD studies. Finally, we found three subnetworks of jointly dysregulated genes within the functional gene network which capture three distinct biological processes: myeloid cell differentiation, protein phosphorylation and synaptic signaling. Further investigation of the myeloid network indicated an upregulation of this network in early stages of AD prior to accumulation of hyperphosphorylated tau and suggested that increased CSF1 transcription in astrocytes may contribute to microglial activation in AD. Thus, we developed a method that integrates multiple data types and external knowledge of gene function to detect candidate genes, applied the method to an AD dataset, and identified several disease-related genes and processes demonstrating the usefulness of the integrative approach.

Recent technological advances have led to a new generation of studies that interrogate multiple molecular levels in the same target tissue of a set of subjects, generating complex multi-omic datasets with which to study disease mechanism. These datasets of genetic, epigenomic, transcriptomic, and other data have the potential to reveal novel biological insights; however, integrative analyses remain challenging and require new computational methods. We developed an integrative Bayesian approach to detect genes with consistent differences between case and control samples across multiple data types. The method further integrates prior knowledge about gene function in the form of a gene functional similarity network to improve statistical inference by sharing information between related genes. We applied our method to an Alzheimer’s disease dataset of epigenomic and transcriptomic data and detected and then validated several novel and known candidate genes as well as three major disease-related biological processes. One of these processes reflected microglial activation and included the cytokine CSF1. Single-nucleus data revealed that CSF1 was primarily upregulated in astrocytes, implicating the involvement of this cell type in microglial activation. Hence, we demonstrated that integrative analysis approaches to multi-omic datasets can improve candidate gene detection and thereby generate new insights into complex diseases.

High-Mobility Group Box 1 Protein Signaling in Painful Diabetic Neuropathy

Diabetes is a global epidemic and more than 50% diabetic patients are also diagnosed with neuropathy, which greatly affects the quality of life of the patients. Available treatments are not always successful due to the limited efficacy and complications, such as addiction and dependency. Studies have implicated that high mobility group box1 (HMGB1) protein plays a crucial role in neuroinflammation and the development of neuropathic conditions. HMGB1 is a proinflammatory cytokine that can be released from necrotic cells in passive form or in response to inflammatory signals as an active form. HMGB1 is the ligand for the receptor for advanced glycation end products (RAGE), and toll-like receptors, (TLR)-2 and TLR4, which also indirectly activates C-X-C chemokine receptor type 4 (CXCR4). We investigated whether blocking of HMGB1 can reduce pain and inflammation in diabetic neuropathic animals to further understand the role of HMGB1 in diabetic neuropathy. Type 2 diabetic rats and mice were treated with natural inhibitor of HMGB1, glycyrrhizin (GLC) for five days/week for four weeks at a dose of 50 mg/kg per day by intraperitoneal injection. The animals were divided into three categories: naïve control, diabetic alone, diabetic with GLC treatment. All of the behavioral analyses were conducted before and after the treatment. The expression of inflammatory markers and changes in histone acetylation in the peripheral nervous system were measured by immunohistochemistry and Western blot analysis after the completion of the treatment. Our study revealed that TLR4, HMGB1, CXCR4, and Nod-like receptor protein 3 (NLRP3) levels were increased in the spinal and dorsal root ganglia (DRG) neurons of Type 2 diabetic mice and rats with painful neuropathy. GLC treatment inhibited the increases in TLR4, NLRP3, and CXCR4 expressions and improved the mechanical and thermal pain threshold in these animals. Immunohistochemical studies revealed that hyperglycemia mediated inflammation influenced HMGB1 acetylation and its release from the neurons. It also altered histone 3 acetylation in the microglial cells. The inhibition of HMGB1 by GLC prevented the release of HMGB1 as well as H3K9 acetylation. These findings indicate that the interruption of HMGB1 mediated inflammation could ameliorate diabetic neuropathy and might exhibit a unique target for the treatment.

Identification of Cholecystokinin by Genome-Wide Profiling as Potential Mediator of Serotonin-Dependent Behavioral Effects of Maternal Separation in the Amygdala

Converging evidence suggests a role of serotonin (5-hydroxytryptamine, 5-HT) and tryptophan hydroxylase 2 (TPH2), the rate-limiting enzyme of 5-HT synthesis in the brain, in modulating long-term, neurobiological effects of early-life adversity. Here, we aimed at further elucidating the molecular mechanisms underlying this interaction, and its consequences for socio-emotional behaviors, with a focus on anxiety and social interaction. In this study, adult, male Tph2 null mutant (Tph2 -/-) and heterozygous (Tph2 +/-) mice, and their wildtype littermates (Tph2 +/+) were exposed to neonatal, maternal separation (MS) and screened for behavioral changes, followed by genome-wide RNA expression and DNA methylation profiling. In Tph2 -/- mice, brain 5-HT deficiency profoundly affected socio-emotional behaviors, i.e., decreased avoidance of the aversive open arms in the elevated plus-maze (EPM) as well as decreased prosocial and increased rule breaking behavior in the resident-intruder test when compared to their wildtype littermates. Tph2 +/- mice showed an ambiguous profile with context-dependent, behavioral responses. In the EPM they showed similar avoidance of the open arm but decreased prosocial and increased rule breaking behavior in the resident-intruder test when compared to their wildtype littermates. Notably, MS effects on behavior were subtle and depended on the Tph2 genotype, in particular increasing the observed avoidance of EPM open arms in wildtype and Tph2 +/- mice when compared to their Tph2 -/- littermates. On the genomic level, the interaction of Tph2 genotype with MS differentially affected the expression of numerous genes, of which a subset showed an overlap with DNA methylation profiles at corresponding loci. Remarkably, changes in methylation nearby and expression of the gene encoding cholecystokinin, which were inversely correlated to each other, were associated with variations in anxiety-related phenotypes. In conclusion, next to various behavioral alterations, we identified gene expression and DNA methylation profiles to be associated with TPH2 inactivation and its interaction with MS, suggesting a gene-by-environment interaction-dependent, modulatory function of brain 5-HT availability.

DNA methylation regulates Microtubule-associated tumor suppressor 1 in human non-small cell lung carcinoma

Microtubule associated tumor suppressor 1 (MTUS1) has been recognized as a tumor suppressor gene in multiple cancers. However, the molecular mechanisms underlying the regulation of MTUS1 are yet to be investigated. This study aimed to clarify the significance of DNA methylation in silencing MTUS1 expression. We report that MTUS1 acts as tumor suppressor in non-small cell lung carcinoma (NSCLC). Analysis of in silico database and subsequent knockdown of DNMT1 suggested an inverse correlation between DNMT1 and MTUS1 function. Interestingly, increased methylation at MTUS1 promoter is associated with low expression of MTUS1. Treatment with DNA methyltransferases (DNMTs) inhibitor, 5-aza-2'-deoxycytidine (AZA) leads to both reduced promoter methylation accompanied with enrichment of H3K9Ac and enhanced MTUS1 expression. Remarkably, knockdown of MTUS1 showed increased proliferation and migration of NSCLC cells in contrast to diminished proliferation and migration, upon treatment with AZA. We concluded that low expression of MTUS1 correlates to DNA methylation and histone deacetylation in human NSCLC.

UHRF2 decreases H3K9ac expression by interacting with it through the PHD and SRA/YDG domain in HepG2 hepatocellular carcinoma cells

Ubiquitin-like with PHD and ring finger domains 2 (UHRF2) is a multi-domain E3 ubiquitin ligase which is involved in epigenetic regulation and plays an essential role in tumorigenesis. However, the role of UHRF2 in histone H3 acetylation has not yet been fully elucidated and few studies have reported its role in hepatocellular carcinoma (HCC). In this study, we examined the correlation between UHRF2 and acetylated H3 in HCC. Immunohistochemistry and western blot analysis demonstrated that the levels of histone H3 lysine 9 acetylation (H3K9ac) and histone H3 lysine 14 acetylation (H3K14ac) were higher in the HCC tissues and HepG2 HCC cells compared with the adjacent non-tumor tissues and L02 normal cells. The level of UHRF2 was higher in the HCC tissues compared with the adjacent non-tumor tissues, but its expression did not exhibit a significant difference between the HepG2 HCC cells and the L02 normal cells. In addition, when comparing the HCC tissues, a higher expression of UHRF2 correlated with a lower expression of H3K9ac in the HCC tissues. The overexpression of UHRF2 increased the expression of H3K9ac in L02 normal cells (P<0.01), but decreased the expression of H3K9ac in HepG2 cancer cells (P<0.05). Moreover, immunofluorescence staining and co-immunoprecipitation assay indicated that UHRF2 co-localized and interacted with H3K9ac in L02 and HepG2 cells and the plant homeodomain (PHD) finger domain was the key domain for UHRF2 directly binding to H3K9ac. Taken together, these results suggest that UHRF2 decreases the expression of H3K9ac in HepG2 HCC cells and interacts with it through the PHD domain.

Deep Feature Selection: Theory and Application to Identify Enhancers and Promoters

Sparse linear models approximate target variable(s) by a sparse linear combination of input variables. Since they are simple, fast, and able to select features, they are widely used in classification and regression. Essentially they are shallow feed-forward neural networks that have three limitations: (1) incompatibility to model nonlinearity of features, (2) inability to learn high-level features, and (3) unnatural extensions to select features in a multiclass case. Deep neural networks are models structured by multiple hidden layers with nonlinear activation functions. Compared with linear models, they have two distinctive strengths: the capability to (1) model complex systems with nonlinear structures and (2) learn high-level representation of features. Deep learning has been applied in many large and complex systems where deep models significantly outperform shallow ones. However, feature selection at the input level, which is very helpful to understand the nature of a complex system, is still not well studied. In genome research, the cis-regulatory elements in noncoding DNA sequences play a key role in the expression of genes. Since the activity of regulatory elements involves highly interactive factors, a deep tool is strongly needed to discover informative features. In order to address the above limitations of shallow and deep models for selecting features of a complex system, we propose a deep feature selection (DFS) model that (1) takes advantages of deep structures to model nonlinearity and (2) conveniently selects a subset of features right at the input level for multiclass data. Simulation experiments convince us that this model is able to correctly identify both linear and nonlinear features. We applied this model to the identification of active enhancers and promoters by integrating multiple sources of genomic information. Results show that our model outperforms elastic net in terms of size of discriminative feature subset and classification accuracy.

Transcription factor binding sites prediction based on modified nucleosomes

In computational methods, position weight matrices (PWMs) are commonly applied for transcription factor binding site (TFBS) prediction. Although these matrices are more accurate than simple consensus sequences to predict actual binding sites, they usually produce a large number of false positive (FP) predictions and so are impoverished sources of information. Several studies have employed additional sources of information such as sequence conservation or the vicinity to transcription start sites to distinguish true binding regions from random ones. Recently, the spatial distribution of modified nucleosomes has been shown to be associated with different promoter architectures. These aligned patterns can facilitate DNA accessibility for transcription factors. We hypothesize that using data from these aligned and periodic patterns can improve the performance of binding region prediction. In this study, we propose two effective features, “modified nucleosomes neighboring” and “modified nucleosomes occupancy”, to decrease FP in binding site discovery. Based on these features, we designed a logistic regression classifier which estimates the probability of a region as a TFBS. Our model learned each feature based on Sp1 binding sites on Chromosome 1 and was tested on the other chromosomes in human CD4+T cells. In this work, we investigated 21 histone modifications and found that only 8 out of 21 marks are strongly correlated with transcription factor binding regions. To prove that these features are not specific to Sp1, we combined the logistic regression classifier with the PWM, and created a new model to search TFBSs on the genome. We tested the model using transcription factors MAZ, PU.1 and ELF1 and compared the results to those using only the PWM. The results show that our model can predict Transcription factor binding regions more successfully. The relative simplicity of the model and capability of integrating other features make it a superior method for TFBS prediction.

Regulation of gene expression in the genomic context

Metazoan life is dependent on the proper temporal and spatial control of gene expression within the many cells-essentially all with the identical genome-that make up the organism. While much is understood about how individual gene regulatory elements function, many questions remain about how they interact to maintain correct regulation globally throughout the genome. In this review we summarize the basic features and functions of the crucial regulatory elements promoters, enhancers, and insulators and discuss some of the ways in which proper interactions between these elements is realized. We focus in particular on the role of core promoter sequences and propose explanations for some of the contradictory results seen in experiments aimed at understanding insulator function. We suggest that gene regulation depends on local genomic context and argue that more holistic in vivo investigations that take into account multiple local features will be necessary to understand how genome-wide gene regulation is maintained.

Characterization of the regulation and function of zinc-dependent histone deacetylases during rodent liver regeneration

The studies reported here were undertaken to define the regulation and functional importance of zinc-dependent histone deacetylase (Zn-HDAC) activity during liver regeneration using the mouse partial hepatectomy (PH) model. The results showed that hepatic HDAC activity was significantly increased in nuclear and cytoplasmic fractions following PH. Further analyses showed isoform-specific effects of PH on HDAC messenger RNA (mRNA) and protein expression, with increased expression of the class I HDACs, 1 and 8, and class II HDAC4 in regenerating liver. Hepatic expression of (class II) HDAC5 was unchanged after PH; however, HDAC5 exhibited transient nuclear accumulation in regenerating liver. These changes in hepatic HDAC expression, subcellular localization, and activity coincided with diminished histone acetylation in regenerating liver. The significance of these events was investigated by determining the effects of suberoylanilide hydroxyamic acid (SAHA, a specific inhibitor of Zn-HDAC activity) on hepatic regeneration. The results showed that SAHA treatment suppressed the effects of PH on histone deacetylation and hepatocellular bromodeoxyuridine (BrdU) incorporation. Further examination showed that SAHA blunted hepatic expression and activation of cell cycle signals downstream of induction of cyclin D1 expression in mice subjected to PH.

CONCLUSION

The data reported here demonstrate isoform-specific regulation of Zn-HDAC expression, subcellular localization, and activity in regenerating liver. These studies also indicate that HDAC activity promotes liver regeneration by regulating hepatocellular cell cycle progression at a step downstream of cyclin D1 induction.

ChIP sequencing to identify p53 targets

Mutant p53 may activate target genes through the interaction of transcription factors or through histone modifications. Chromatin immunoprecipitation (ChIP) is a method commonly used to study these types of protein interactions. In order to generate a list of target genes that may be activated through this mechanism, ChIP sequencing may be used. ChIP sequencing involves the mass parallel sequencing of ChIP DNA fragments. We describe a method by which to prepare chromatin immunoprecipitation sequencing libraries and how to analyze sequencing data. In this procedure, prepared libraries have been sent to a core facility. The results have been verified using quantitative PCR.

MicroRNA-224 is up-regulated in hepatocellular carcinoma through epigenetic mechanisms

MicroRNA-224 (miR-224) is one of the most commonly up-regulated microRNAs in hepatocellular carcinoma (HCC), which affects crucial cellular processes such as apoptosis and cell proliferation. In this study, we aim to elucidate the molecular mechanism that leads to the overexpression of miR-224 in HCC. We examined the transcript expression of miR-224 and neighboring miR-452 and genes on chromosome Xq28 in tumor and paired adjacent nontumorous tissues from 100 patients with HCC and found that miR-224 is coordinately up-regulated with its neighboring microRNA (miRNA) and genes. This coordinated up-regulation of miRNAs and genes at the Xq28 locus can be mimicked in nontransformed immortalized human liver cells by the introduction of histone deacetylase (HDAC) inhibitors, which resulted in a corresponding increase in histone H3 acetylation in this region. This miR-224-residing locus in Xq28 is reciprocally regulated by HDAC1, HDAC3, and histone acetylase protein, E1A binding protein p300 (EP300). Notably, in HCC tumors that significantly overexpress microRNA-224, EP300 is also overexpressed and displays increased binding to the Xq28 locus. In transformed HCC cells, high miR-224 expression can be attenuated through the inhibition of EP300, using either siRNA or the specific drug C646. In summary, overexpression of EP300 may account, in part, for the up-regulation of miR-224 expression in patients with HCC.

Combined chromatin and expression analysis reveals specific regulatory mechanisms within cytokine genes in the macrophage early immune response

Macrophages play a critical role in innate immunity, and the expression of early response genes orchestrate much of the initial response of the immune system. Macrophages undergo extensive transcriptional reprogramming in response to inflammatory stimuli such as Lipopolysaccharide (LPS).To identify gene transcription regulation patterns involved in early innate immune responses, we used two genome-wide approaches--gene expression profiling and chromatin immunoprecipitation-sequencing (ChIP-seq) analysis. We examined the effect of 2 hrs LPS stimulation on early gene expression and its relation to chromatin remodeling (H3 acetylation; H3Ac) and promoter binding of Sp1 and RNA polymerase II phosphorylated at serine 5 (S5P RNAPII), which is a marker for transcriptional initiation. Our results indicate novel and alternative gene regulatory mechanisms for certain proinflammatory genes. We identified two groups of up-regulated inflammatory genes with respect to chromatin modification and promoter features. One group, including highly up-regulated genes such as tumor necrosis factor (TNF), was characterized by H3Ac, high CpG content and lack of TATA boxes. The second group, containing inflammatory mediators (interleukins and CCL chemokines), was up-regulated upon LPS stimulation despite lacking H3Ac in their annotated promoters, which were low in CpG content but did contain TATA boxes. Genome-wide analysis showed that few H3Ac peaks were unique to either +/-LPS condition. However, within these, an unpacking/expansion of already existing H3Ac peaks was observed upon LPS stimulation. In contrast, a significant proportion of S5P RNAPII peaks (approx 40%) was unique to either condition. Furthermore, data indicated a large portion of previously unannotated TSSs, particularly in LPS-stimulated macrophages, where only 28% of unique S5P RNAPII peaks overlap annotated promoters. The regulation of the inflammatory response appears to occur in a very specific manner at the chromatin level for specific genes and this study highlights the level of fine-tuning that occurs in the immune response.

Systematic clustering of transcription start site landscapes

Genome-wide, high-throughput methods for transcription start site (TSS) detection have shown that most promoters have an array of neighboring TSSs where some are used more than others, forming a distribution of initiation propensities. TSS distributions (TSSDs) vary widely between promoters and earlier studies have shown that the TSSDs have biological implications in both regulation and function. However, no systematic study has been made to explore how many types of TSSDs and by extension core promoters exist and to understand which biological features distinguish them. In this study, we developed a new non-parametric dissimilarity measure and clustering approach to explore the similarities and stabilities of clusters of TSSDs. Previous studies have used arbitrary thresholds to arrive at two general classes: broad and sharp. We demonstrated that in addition to the previous broad/sharp dichotomy an additional category of promoters exists. Unlike typical TATA-driven sharp TSSDs where the TSS position can vary a few nucleotides, in this category virtually all TSSs originate from the same genomic position. These promoters lack epigenetic signatures of typical mRNA promoters and a substantial subset of them are mapping upstream of ribosomal protein pseudogenes. We present evidence that these are likely mapping errors, which have confounded earlier analyses, due to the high similarity of ribosomal gene promoters in combination with known G addition bias in the CAGE libraries. Thus, previous two-class separations of promoter based on TSS distributions are motivated, but the ultra-sharp TSS distributions will confound downstream analyses if not removed.

Tight associations between transcription promoter type and epigenetic variation in histone positioning and modification

Background

Transcription promoters are fundamental genomic cis-elements controlling gene expression. They can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. Eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. For instance, it has been recently shown that histone with H3K9 acetylation (H3K9ac) is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone H3K9 and promoter architecture.

Results

Here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. We found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. Specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. This feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. More strikingly, the overall expression levels of genes associated with broad promoters (but not peak promoters) with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones.

Conclusion

These results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture.

Transcription initiation patterns indicate divergent strategies for gene regulation at the chromatin level

The application of deep sequencing to map 5' capped transcripts has confirmed the existence of at least two distinct promoter classes in metazoans: "focused" promoters with transcription start sites (TSSs) that occur in a narrowly defined genomic span and "dispersed" promoters with TSSs that are spread over a larger window. Previous studies have explored the presence of genomic features, such as CpG islands and sequence motifs, in these promoter classes, but virtually no studies have directly investigated the relationship with chromatin features. Here, we show that promoter classes are significantly differentiated by nucleosome organization and chromatin structure. Dispersed promoters display higher associations with well-positioned nucleosomes downstream of the TSS and a more clearly defined nucleosome free region upstream, while focused promoters have a less organized nucleosome structure, yet higher presence of RNA polymerase II. These differences extend to histone variants (H2A.Z) and marks (H3K4 methylation), as well as insulator binding (such as CTCF), independent of the expression levels of affected genes. Notably, differences are conserved across mammals and flies, and they provide for a clearer separation of promoter architectures than the presence and absence of CpG islands or the occurrence of stalled RNA polymerase. Computational models support the stronger contribution of chromatin features to the definition of dispersed promoters compared to focused start sites. Our results show that promoter classes defined from 5' capped transcripts not only reflect differences in the initiation process at the core promoter but also are indicative of divergent transcriptional programs established within gene-proximal nucleosome organization.

Genetic variation in the KIAA0319 5' region as a possible contributor to dyslexia

Reading disabilities (RD) have been linked and associated with markers on chromosome 6p with results from multiple independent samples pointing to KIAA0319 as a risk gene and specifically, the 5' region of this gene. Here we focus genetic studies on a 2.3 kb region spanning the predicted promoter, the first untranslated exon, and part of the first intron, a region we identified as a region of open chromatin. Using DNA from probands with RD, we screened for genetic variants and tested select variants for association. We identified 17 DNA variants in this sample of probands, 16 of which were previously reported in public databases and one previously identified in a screen of this region. Based on the allele frequencies in the probands compared to public databases, and on possible functional consequences of the variation, we selected seven variants to test for association in a sample of families with RD, in addition to four variants which had been tested previously. We also tested two markers 5' of this region that were previously reported as associated. The strongest evidence for association was observed with alleles of the microsatellite marker located in the first untranslated exon and haplotypes of that marker. These results support previous studies indicating the 5' region of the KIAA0319 gene as the location of risk alleles contributing to RD.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5′-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP–chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.