Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data

Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus

BACKGROUND

Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks.

RESULTS

In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications.

CONCLUSIONS

Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes' rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.

Epigenetic regulation of nuclear processes in fungal plant pathogens

Through the association of protein complexes to DNA, the eukaryotic nuclear genome is broadly organized into open euchromatin that is accessible for enzymes acting on DNA and condensed heterochromatin that is inaccessible. Chemical and physical alterations to chromatin may impact its organization and functionality and are therefore important regulators of nuclear processes. Studies in various fungal plant pathogens have uncovered an association between chromatin organization and expression of in planta-induced genes that are important for pathogenicity. This review discusses chromatin-based regulation mechanisms as determined in the fungal plant pathogen Verticillium dahliae and relates the importance of epigenetic transcriptional regulation and other nuclear processes more broadly in fungal plant pathogens.

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.

A novel role of Krüppel-like factor 8 as an apoptosis repressor in hepatocellular carcinoma

Background

Krüppel-like factor 8 (KLF8), a cancer-promoting factor that regulates critical gene transcription and cellular cancer-related events, has been implicated in tumor development and progression. However, the functional role of KLF8 in the pathogenesis of hepatocellular carcinoma (HCC) remains largely unknown.

Methods

The gene expression patterns and genome-wide regulatory profiles of HCC cells after KLF8 knockout were analyzed by using RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) of histone H3 lysine 27 acetylation (H3K27ac) combined with bioinformatics analysis. Transcription factor-binding motifs that recognized by KLF8 were evaluated by motif analysis. For the predicted target genes, transcriptional changes were examined by ChIP, and loss of function experiments were conducted by siRNA transfection.

Results

KLF8 functioned as a transcription repressor in HCC and mainly regulated apoptotic-related genes directly. A total of 1,816 differentially expressed genes after KLF8 knockout were identified and significantly corresponded to global changes in H3K27ac status. Furthermore, two predicted target genes, high-mobility group AT-hook 2 (HMGA2) and matrix metalloproteinase 7 (MMP7), were identified as important participants in KLF8-mediated anti-apoptotic effect in HCC. Knockout of KLF8 enhanced cell apoptosis process and caused increase in the associated H3K27ac, whereas suppression HMGA2 or MMP7 attenuated these biological effects.

Conclusions

Our work suggests a novel role and mechanism for KLF8 in the regulation of cell apoptosis in HCC and facilitates the discovery of potential therapeutic targets for HCC treatment.

AFF4 regulates osteogenic differentiation of human dental follicle cells

As a member of the AFF (AF4/FMR2) family, AFF4 is a transcription elongation factor that is a component of the super elongation complex. AFF4 serves as a scaffolding protein that connects transcription factors and promotes gene transcription through elongation and chromatin remodelling. Here, we investigated the effect of AFF4 on human dental follicle cells (DFCs) in osteogenic differentiation. In this study, we found that small interfering RNA-mediated depletion of AFF4 resulted in decreased alkaline phosphatase (ALP) activity and impaired mineralization. In addition, the expression of osteogenic-related genes (DLX5, SP7, RUNX2 and BGLAP) was significantly downregulated. In contrast, lentivirus-mediated overexpression of AFF4 significantly enhanced the osteogenic potential of human DFCs. Mechanistically, we found that both the mRNA and protein levels of ALKBH1, a critical regulator of epigenetics, changed in accordance with AFF4 expression levels. Overexpression of ALKBH1 in AFF4-depleted DFCs partially rescued the impairment of osteogenic differentiation. Our data indicated that AFF4 promoted the osteogenic differentiation of DFCs by upregulating the transcription of ALKBH1.

Distinct DNA Sequence Preference for Histone Occupancy in Primary and Transformed Cells

Genome-wide occupancy of several histone modifications in various cell types has been studied using chromatin immunoprecipitation (ChIP) sequencing. Histone occupancy depends on DNA sequence features like inter-strand symmetry of base composition and periodic occurrence of TT/AT. However, whether DNA sequence motifs act as an additional effector of histone occupancy is not known. We have analyzed the presence of DNA sequence motifs in publicly available ChIP-sequence datasets for different histone modifications. Our results show that DNA sequence motifs are associated with histone occupancy, some of which are different between primary and transformed cells. The motifs for primary and transformed cells showed different levels of GC-richness and proximity to transcription start sites (TSSs). The TSSs associated with transformed or primary cell-specific motifs showed different levels of TSS flank transcription in primary and transformed cells. Interestingly, TSSs with a motif-linked occupancy of H2AFZ, a component of positioned nucleosomes, showed a distinct pattern of RNA Polymerase II (POLR2A) occupancy and TSS flank transcription in primary and transformed cells. These results indicate that DNA sequence features dictate differential histone occupancy in primary and transformed cells, and the DNA sequence motifs affect transcription through regulation of histone occupancy.

H3K4 demethylase KDM5B regulates global dynamics of transcription elongation and alternative splicing in embryonic stem cells

Epigenetic regulation of chromatin plays a critical role in controlling embryonic stem (ES) cell self-renewal and pluripotency. However, the roles of histone demethylases and activating histone modifications such as trimethylated histone 3 lysine 4 (H3K4me3) in transcriptional events such as RNA polymerase II (RNAPII) elongation and alternative splicing are largely unknown. In this study, we show that KDM5B, which demethylates H3K4me3, plays an integral role in regulating RNAPII occupancy, transcriptional initiation and elongation, and alternative splicing events in ES cells. Depletion of KDM5B leads to altered RNAPII promoter occupancy, and decreased RNAPII initiation and elongation rates at active genes and at genes marked with broad H3K4me3 domains. Moreover, our results demonstrate that spreading of H3K4me3 from promoter to gene body regions, which is mediated by depletion of KDM5B, modulates RNAPII elongation rates and RNA splicing in ES cells. We further show that KDM5B is enriched nearby alternatively spliced exons, and depletion of KDM5B leads to altered levels of H3K4 methylation in alternatively spliced exon regions, which is accompanied by differential expression of these alternatively splice exons. Altogether, our data indicate an epigenetic role for KDM5B in regulating RNAPII elongation and alternative splicing, which may support the diverse mRNA repertoire in ES cells.

Decrypting ENCODEd epigenetic marks of human tRN-A-RS genes in normal, stem and cancer cell lines

Screening large-scale ENCODE data of 625 cytoplasmic transfer RNA (tRNAs) and 37 aminoacyl tRNA synthetase (AARSs) human genes, we deconstruct the array of relations between 10 histone marks affecting 15 chromatin states; their tissue specificity and variations and interchange amongst normal, cancerous and stem cells. The histone marks of RNA Pol II transcribed AARS genes share, but also contrast with that on RNA Pol III transcribed tRNA genes. tRNAs with identical/similar sequences may be in significantly varying states even within the same cell line; the chromatin scaffold, where the tRNA gene resides, is the key determinant. Hepatocellular carcinoma cell line has dominant H3K27me3, and singular clustering of other marks. Leukaemic cell line has hyperactive genes. The quiescence of the stem cells is encoded in the markers. Leaving aside the important exceptions in stem cells and elsewhere, tRNAs with cove scores above 50 have active markers and precise sets of transcription factors, and are usually well conserved compared to the low-scoring ones. Pseudo tRNAs are in heterochromatin/repressed state with anomalous exceptions in cancer cells. We motivate that Epigenetic-Phishing hacks the translation apparatus through the chromatin states governed by the histone marks of tRNA and AARS genes, and speculate on their therapeutic implications in cancer and on stem cells.

Baseline Chromatin Modification Levels May Predict Interindividual Variability in Ozone-Induced Gene Expression

Traditional toxicological paradigms have relied on factors such as age, genotype, and disease status to explain variability in responsiveness to toxicant exposure; however, these are neither sufficient to faithfully identify differentially responsive individuals nor are they modifiable factors that can be leveraged to mitigate the exposure effects. Unlike these factors, the epigenome is dynamic and shaped by an individual's environment. We sought to determine whether baseline levels of specific chromatin modifications correlated with the interindividual variability in their ozone (O3)-mediated induction in an air-liquid interface model using primary human bronchial epithelial cells from a panel of 11 donors. We characterized the relationship between the baseline abundance of 6 epigenetic markers with established roles as key regulators of gene expression-histone H3 lysine 4 trimethylation (H3K4me3), H3K27 acetylation (H3K27ac), pan-acetyl H4 (H4ac), histone H3K27 di/trimethylation (H3K27me2/3), unmodified H3, and 5-hydroxymethylcytosine (5-hmC)-and the variability in the O3-induced expression of IL-8, IL-6, COX2, and HMOX1. Baseline levels of H3K4me3, H3K27me2/3, and 5-hmC, but not H3K27ac, H4ac, and total H3, correlated with the interindividual variability in O3-mediated induction of HMOX1 and COX2. In contrast, none of the chromatin modifications that we examined correlated with the induction of IL-8 and IL-6. From these findings, we propose an "epigenetic seed and soil" model in which chromatin modification states between individuals differ in the relative abundance of specific modifications (the "soil") that govern how receptive the gene is to toxicant-mediated cellular signals (the "seed") and thus regulate the magnitude of exposure-related gene induction.

Tissue-Specific Expression of Estrogen Receptor 1 Is Regulated by DNA Methylation in a T-DMR

The mechanism controlling tissue-specific expression of estrogen receptor 1 (ESR1) is unclear. In other genes, DNA methylation of a region called the tissue-dependent and differentially methylated region (T-DMR) has been associated with tissue-specific gene expression. This study investigated whether human ESR1 has a T-DMR and whether DNA methylation of the T-DMR regulates its expression. ESR1 expression was tissue-specific, being high in the endometrium and mammary gland and low/nil in the placenta and skin. Therefore, DNA methylation profiles of the promoter of ESR1 were analyzed in these tissues and in breast cancer tissues. In all of the normal tissues, the proximal promoter regions were unmethylated. On the other hand, the distal regions (T-DMR) were unmethylated in the endometrium and mammary gland, but were moderately methylated and hypermethylated in the placenta and skin, respectively. T-DMR-methylated reporter assay was performed to examine whether DNA methylation at the T-DMR suppresses ESR1 transcription. T-DMR, but not the promoter region, had transcriptional activities and DNA methylation of the T-DMR suppressed ESR1 transcription. Early growth response protein 1 was shown to be a possible transcription factor to bind the T-DMR and up-regulate ESR1 expression. ESR1 has several upstream exons, and each upstream exon, Exon-A/Exon-B/Exon-C, had its own T-DMR. In some breast cancer cases and breast cancer cell lines, ESR1 expression was not regulated by DNA methylation at T-DMR as it is in normal tissues. In conclusion, ESR1 has a T-DMR. DNA methylation status at the T-DMR is involved in tissue-specific ESR1 expression in normal tissues but not always in breast cancer.

Tet1 and Tet2 Protect DNA Methylation Canyons against Hypermethylation

DNA methylation is a dynamic epigenetic modification with an important role in cell fate specification and reprogramming. The Ten eleven translocation (Tet) family of enzymes converts 5-methylcytosine to 5-hydroxymethylcytosine, which promotes passive DNA demethylation and functions as an intermediate in an active DNA demethylation process. Tet1/Tet2 double-knockout mice are characterized by developmental defects and epigenetic instability, suggesting a requirement for Tet-mediated DNA demethylation for the proper regulation of gene expression during differentiation. Here, we used whole-genome bisulfite and transcriptome sequencing to characterize the underlying mechanisms. Our results uncover the hypermethylation of DNA methylation canyons as the genomic key feature of Tet1/Tet2 double-knockout mouse embryonic fibroblasts. Canyon hypermethylation coincided with disturbed regulation of associated genes, suggesting a mechanistic explanation for the observed Tet-dependent differentiation defects. Based on these results, we propose an important regulatory role of Tet-dependent DNA demethylation for the maintenance of DNA methylation canyons, which prevents invasive DNA methylation and allows functional regulation of canyon-associated genes.

Epigenetic modifications are associated with inter-species gene expression variation in primates

BACKGROUND

Changes in gene regulation have long been thought to play an important role in evolution and speciation, especially in primates. Over the past decade, comparative genomic studies have revealed extensive inter-species differences in gene expression levels, yet we know much less about the extent to which regulatory mechanisms differ between species.

RESULTS

To begin addressing this gap, we perform a comparative epigenetic study in primate lymphoblastoid cell lines, to query the contribution of RNA polymerase II and four histone modifications, H3K4me1, H3K4me3, H3K27ac, and H3K27me3, to inter-species variation in gene expression levels. We find that inter-species differences in mark enrichment near transcription start sites are significantly more often associated with inter-species differences in the corresponding gene expression level than expected by chance alone. Interestingly, we also find that first-order interactions among the five marks, as well as chromatin states, do not markedly contribute to the degree of association between the marks and inter-species variation in gene expression levels, suggesting that the marginal effects of the five marks dominate this contribution.

CONCLUSIONS

Our observations suggest that epigenetic modifications are substantially associated with changes in gene expression levels among primates and may represent important molecular mechanisms in primate evolution.

Nonparametric Tests for Differential Histone Enrichment with ChIP-Seq Data

Chromatin immunoprecipitation sequencing (ChIP-seq) is a powerful method for analyzing protein interactions with DNA. It can be applied to identify the binding sites of transcription factors (TFs) and genomic landscape of histone modification marks (HMs). Previous research has largely focused on developing peak-calling procedures to detect the binding sites for TFs. However, these procedures may fail when applied to ChIP-seq data of HMs, which have diffuse signals and multiple local peaks. In addition, it is important to identify genes with differential histone enrichment regions between two experimental conditions, such as different cellular states or different time points. Parametric methods based on Poisson/negative binomial distribution have been proposed to address this differential enrichment problem and most of these methods require biological replications. However, many ChIP-seq data usually have a few or even no replicates. We propose a nonparametric method to identify the genes with differential histone enrichment regions even without replicates. Our method is based on nonparametric hypothesis testing and kernel smoothing in order to capture the spatial differences in histone-enriched profiles. We demonstrate the method using ChIP-seq data on a comparative epigenomic profiling of adipogenesis of murine adipose stromal cells and the Encyclopedia of DNA Elements (ENCODE) ChIP-seq data. Our method identifies many genes with differential H3K27ac histone enrichment profiles at gene promoter regions between proliferating preadipocytes and mature adipocytes in murine 3T3-L1 cells. The test statistics also correlate with the gene expression changes well and are predictive to gene expression changes, indicating that the identified differentially enriched regions are indeed biologically meaningful.

Promoter or enhancer, what's the difference? Deconstruction of established distinctions and presentation of a unifying model

Gene transcription is strictly controlled by the interplay of regulatory events at gene promoters and gene-distal regulatory elements called enhancers. Despite extensive studies of enhancers, we still have a very limited understanding of their mechanisms of action and their restricted spatio-temporal activities. A better understanding would ultimately lead to fundamental insights into the control of gene transcription and the action of regulatory genetic variants involved in disease. Here, I review and discuss pros and cons of state-of-the-art genomics methods to localize and infer the activity of enhancers. Among the different approaches, profiling of enhancer RNAs yields the highest specificity and may be superior in detecting in vivo activity. I discuss their apparent similarities to promoters, which challenge the established view of enhancers and promoters as distinct entities, and present a unifying model of regulatory elements in transcriptional regulation, in which activity, transcriptional output and regulatory function is context specific.

Distinct promoters, subjected to epigenetic regulation, drive the expression of two clusterin mRNAs in prostate cancer cells

The human clusterin (CLU) gene codes for several mRNAs characterized by different sequences at their 5' end. We investigated the expression of two CLU mRNAs, called CLU 1 and CLU 2, in immortalized (PNT1a) and tumorigenic (PC3 and DU145) prostate epithelial cells, as well as in normal fetal fibroblasts (WI38) following the administration of the epigenetic drugs 5-aza-2'-deoxycytidine (AZDC) and trichostatin A (TSA) given either as single or combined treatment (AZDC-TSA). Our experimental evidences show that: a) CLU 1 is the most abundant transcript variant. b) CLU 2 is expressed at a low level in normal fibroblasts and virtually absent in prostate cancer cells. c) CLU 1, and to a greater extent CLU 2 expression, increased by AZDC-TSA treatment in prostate cancer cells. d) Both CLU 1 and CLU 2 encode for secreted CLU. e) P2, a novel promoter that overlaps the CLU 2 Transcription Start Site (TSS), drives CLU 2 expression. f) A CpG island, methylated in prostate cancer cells and not in normal fibroblasts, is responsible for long-term heritable regulation of CLU 1 expression. g) ChIP assay of histone tail modifications at CLU promoters (P1 and P2) shows that treatment of prostate cancer cells with AZDC-TSA causes enrichment of Histone3(Lys9)acetylated (H3K9ac) and reduction of Histone3(Lys27)trimethylated (H3K27me3), inducing active transcription of both CLU variants. In conclusion, we show for the first time that the expression of CLU 2 mRNA is driven by a novel promoter, P2, whose activity responds to epigenetic drugs treatment through changes in histone modifications.

Control of VEGF-A transcriptional programs by pausing and genomic compartmentalization

Vascular endothelial growth factor A (VEGF-A) is a master regulator of angiogenesis, vascular development and function. In this study we investigated the transcriptional regulation of VEGF-A-responsive genes in primary human aortic endothelial cells (HAECs) and human umbilical vein endothelial cells (HUVECs) using genome-wide global run-on sequencing (GRO-Seq). We demonstrate that half of VEGF-A-regulated gene promoters are characterized by a transcriptionally competent paused RNA polymerase II (Pol II). We show that transition into productive elongation is a major mechanism of gene activation of virtually all VEGF-regulated genes, whereas only ∼40% of the genes are induced at the level of initiation. In addition, we report a comprehensive chromatin interaction map generated in HUVECs using tethered conformation capture (TCC) and characterize chromatin interactions in relation to transcriptional activity. We demonstrate that sites of active transcription are more likely to engage in chromatin looping and cell type-specific transcriptional activity reflects the boundaries of chromatin interactions. Furthermore, we identify large chromatin compartments with a tendency to be coordinately transcribed upon VEGF-A stimulation. We provide evidence that these compartments are enriched for clusters of regulatory regions such as super-enhancers and for disease-associated single nucleotide polymorphisms (SNPs). Collectively, these findings provide new insights into mechanisms behind VEGF-A-regulated transcriptional programs in endothelial cells.

A wavelet approach to detect enriched regions and explore epigenomic landscapes

Epigenetic landscapes represent how cells regulate gene activity. To understand their effect on gene regulation, it is important to detect their occupancy in the genome. Unlike transcription factors whose binding regions are limited to narrow regions, histone modification marks are enriched over broader areas. The stochastic characteristics unique to each mark make it hard to detect their enrichment. Classically, a predefined window has been used to detect their enrichment. However, these approaches heavily rely on the predetermined parameters. Also, the window-based approaches cannot handle the enrichment of multiple marks. We propose a novel algorithm, called SeqW, to detect enrichment of multiple histone modification marks. SeqW applies a zooming approach to detect a broadly enriched domain. The zooming approach helps domain detection by increasing signal-to-noise ratio. The borders of the domains are detected by studying the characteristics of signals in the wavelet domain. We show that SeqW outperformed previous predictors in detecting broad peaks. Also, we applied SeqW in studying spatial combinations of histone modification patterns.

Genome-wide analysis of histone modifications in human endometrial stromal cells

Dramatic changes of gene expressions occur in human endometrial stromal cells (ESCs) during decidualization. The changes in gene expression are associated with changes of chromatin structure, which are regulated by histone modifications. Here we investigated genome-wide changes in histone modifications associated with decidualization in human ESCs using chromatin immunoprecipitation combined with next-generation sequencing. ESCs were incubated with estradiol and medroxyprogesterone acetate for 14 days to induce decidualization. The chromatin immunoprecipitation-sequence data showed that induction of decidualization increased H3K27ac and H3K4me3 signals in many genomic regions but decreased in only a few regions. Most of the H3K27ac-increased regions (80%) and half of the H3K4me3-increased regions were located in the distal promoter regions (more than 3 kb upstream or downstream of the transcription start site). RNA sequence showed that induction of decidualization up-regulated 881 genes, 223 of which had H3K27ac- or H3K4me3-increased regions in the proximal and distal promoter regions. Induction of decidualization increased the mRNA levels of these genes more than it increased the mRNA levels of genes without H3K27ac- or H3K4me3-increased regions. Pathway analysis revealed that up-regulated genes with the H3K27ac- or H3K4me3-increased regions were associated with the insulin signaling, which may be involved in glucose uptake that is necessary for ESCs to undergo decidualization. These results show that histone modification statuses on a genome-wide basis change in human ESCs during decidualization. The main changes of histone modifications are increases of H3K27ac and H3K4me3 in both the proximal and distal promoter regions, which are involved in the up-regulation of gene expression that occurs during decidualization.

CBX8, a novel DNA repair protein, promotes tumorigenesis in human esophageal carcinoma

DNA damage response and repair are carried out by certain proteins following damage by environmental clastogens, such as ionizing radiation and reactive oxygen species. It has been reported that many carcinomas that are characterized by resistance to chemotherapy and poor outcomes show dysfunction of these proteins. Chromobox homologue 8 (CBX8), a member of the polycomb group of proteins, has been identified as a factor that protects tumor cells from the detrimental effects of ionizing radiation (IR) or hydrogen peroxide (H2O2). In this study, we found that CBX8 was up-regulated in esophageal carcinoma tissues compared with adjacent non-cancerous tissues (P<0.01) and correlated with TNM stage in esophageal squamous cell carcinoma patients. Depletion of CBX8 decreased cell proliferation both in vitro and in vivo and increased the phosphorylation levels of p21, Wee1, and CHK1, which result in cyclin-dependent kinase inhibition and cell-cycle delay. CBX8 depletion also led to accumulation of spontaneous DNA damage and raised the sensitivity of tumor cells to IR or H2O2. We also found that the total level of CBX8 in the cells was increased after treating tumor cells with clastogens. In addition, our data showed that decreased CBX8 expression was accompanied by the reduction of EZH2 and EED, which have been reported to participate in DNA damage repair. Collectively, CBX8 might emerge as an oncogene for promoting the proliferation of tumor cells and raising the resistance of neoplasms to chemotherapy.

Combinatorial epigenetic patterns as quantitative predictors of chromatin biology

Background

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) is the most widely used method for characterizing the epigenetic states of chromatin on a genomic scale. With the recent availability of large genome-wide data sets, often comprising several epigenetic marks, novel approaches are required to explore functionally relevant interactions between histone modifications. Computational discovery of "chromatin states" defined by such combinatorial interactions enabled descriptive annotations of genomes, but more quantitative approaches are needed to progress towards predictive models.

Results

We propose non-negative matrix factorization (NMF) as a new unsupervised method to discover combinatorial patterns of epigenetic marks that frequently co-occur in subsets of genomic regions. We show that this small set of combinatorial "codes" can be effectively displayed and interpreted. NMF codes enable dimensionality reduction and have desirable statistical properties for regression and classification tasks. We demonstrate the utility of codes in the quantitative prediction of Pol2-binding and the discrimination between Pol2-bound promoters and enhancers. Finally, we show that specific codes can be linked to molecular pathways and targets of pluripotency genes during differentiation.

Conclusions

We have introduced and evaluated a new computational approach to represent combinatorial patterns of epigenetic marks as quantitative variables suitable for predictive modeling and supervised machine learning. To foster widespread adoption of this method we make it available as an open-source software-package – epicode at https://github.com/mcieslik-mctp/epicode.

Epigenetic regulation of gene responsiveness in Arabidopsis

The regulation of chromatin structure is inevitable for proper transcriptional response in eukaryotes. Recent reports in Arabidopsis have suggested that gene responsiveness is modulated by particular chromatin status. One such feature is H2A.Z, a histone variant conserved among eukaryotes. In Arabidopsis, H2A.Z is enriched within gene bodies of transcriptionally variable genes, which is in contrast to genic DNA methylation found within constitutive genes. In the absence of H2A.Z, the genes normally harboring H2A.Z within gene bodies are transcriptionally misregulated, while DNA methylation is unaffected. Therefore, H2A.Z may promote variability of gene expression without affecting genic DNA methylation. Another epigenetic information that could be important for gene responsiveness is trimethylation of histone H3 lysine 4 (H3K4me3). The level of H3K4me3 increases when stress responsive genes are transcriptionally activated, and it decreases after recovery from the stress. Even after the recovery, however, H3K4me3 is kept at some atypical levels, suggesting possible role of H3K4me3 for a stress memory. In this review, we summarize and discuss the growing evidences connecting chromatin features and gene responsiveness.

Importance of C/EBPβ binding and histone acetylation status in the promoter regions for induction of IGFBP-1, PRL, and Mn-SOD by cAMP in human endometrial stromal cells

Dynamic changes of gene expressions occur in human endometrial stromal cells (ESCs) during decidualization. CCAAT/enhancer-binding proteinβ (C/EBPβ) regulates the expression of a number of decidualization-related genes. In addition to transcription factors, it is important to know the role of epigenetic mechanisms, such as histone modifications in the regulation of decidualization-related genes. This study investigated the molecular and epigenetic mechanisms by which cAMP up-regulates the expression of IGF-binding protein-1 (IGFBP-1), prolactin (PRL), and manganese superoxide dismutase (Mn-SOD) in ESC. ESCs isolated from proliferative phase endometrium were incubated with cAMP to induce decidualization. IGFBP-1, PRL, and Mn-SOD mRNA expressions were determined by real-time RT-PCR. The C/EBPβ binding and histone modification status (acetylation of histone-H3 lysine-27 [H3K27ac]) in the promoter were examined by chromatin immunoprecipitation assay. Knockdowns of C/EBPβ were performed using the small interfering RNA method. cAMP induced mRNA expressions of IGFBP-1 and PRL accompanied by the increases in both C/EBPβ binding activities and H3K27ac levels in the promoters. The stimulatory effects of cAMP on mRNA levels and H3K27ac levels were completely abolished by C/EBPβ knockdown. cAMP increased Mn-SOD mRNA levels and C/EBPβ binding activities in the enhancer region. C/EBPβ knockdown inhibited Mn-SOD mRNA levels. The H3K27ac levels in the enhancer were high before cAMP stimulus but were not further increased by cAMP and were not inhibited by C/EBPβ knockdown. These results show that C/EBPβ regulates the expression of IGFBP-1 and PRL by altering the histone acetylation status of their promoters but differently regulates Mn-SOD gene expression in human ESC during decidualization.

MYC-induced epigenetic activation of GATA4 in lung adenocarcinoma

Human lung cancer is a disease with high incidence and accounts for most cancer-related deaths in both men and women. Metastasis is a common event in non-small cell lung carcinoma (NSCLC), diminishing the survival chance of the patients with this type of tumor. It has been shown that MYC is involved in the development of metastasis from NSCLC, but the mechanism underlying this switch remained to be identified. Here, we focus on GATA4 as a MYC target in the development of metastasis with origin in lung adenocarcinoma, the most common type of NSCLC. Epigenetic alterations at the GATA4 promoter level were observed after MYC expression in lung adenocarcinoma in vivo and in vitro. Such alterations include site-specific demethylation that accompanies the displacement of the MYC-associated zinc finger protein (MAZ) from the GATA4 promoter, which leads to GATA4 expression. Histone modification analysis of the GATA4 promoter revealed a switch from repressive histone marks to active histone marks after MYC binding, which corresponds to active GATA4 expression. Our results thus identify a novel epigenetic mechanism by which MYC activates GATA4 leading to metastasis in lung adenocarcinoma, suggesting novel potential targets for the development of antimetastatic therapy.

DNA hypermethylation of alternatively spliced and repeat sequences in humans

DNA methylation is presently accepted as a tentative regulatory parameter in splicing. Recently, we reported significant methylation differences among various exonic splicing-enhancing elements and alternative splicing events, based on CpG methylation data from the Human Epigenome Project for chromosomes 6, 20 and 22. Presently, using a different computational approach and the same database, we report: (a) significant increase of hypermethylation in intronic and exonic sequences close to acceptor sites, relative to overall introns and exons, respectively (1,973 CpGs examined); (b) frequent CpGs, mostly hypomethylated, in donors and infrequent CpGs mostly hypermethylated, in acceptors; and (c) hypermethylation in cassette exons which are occasionally spliced and have weaker average splicing potential, relative to constitutive exons (p < 0.0001). CpGs are hypomethylated in non-coding exons (only 16 % hypermethylation). Single-exon genes, similarly to first exons, frequently contain hypomethylated CpGs, while in internal and last exons CpGs are more frequently hypermethylated. Methylation is also more frequent in strange introns and splice sites processed by the minor spliceosome, e.g., ATAC, (p < 0.0001 in all cases), but not in sites of incomplete processing, e.g., retained introns or bleeding exons, (p = 0.706 and p = 0.313, respectively). Most Alus, which are known to contribute to transcript presentation, are heavily methylated, in contrast with other Alus, e.g., AluJo and mammalian interspersed repetitive elements which have been previously associated with alternative expression. These results elucidate the role of intragenic methylation in association with alternative splicing and facilitate the evaluation of genomic variations/polymorphisms and the development of tools for the prediction of alternative splicing events.

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.