Histone modifications are associated with transcript isoform diversity in normal and cancer cells

Age-dependent regulation of ELP1 exon 20 splicing in Familial Dysautonomia by RNA Polymerase II kinetics and chromatin structure

Familial Dysautonomia (FD) is a rare disease caused by ELP1 exon 20 skipping. Here we clarify the role of RNA Polymerase II (RNAPII) and chromatin on this splicing event. A slow RNAPII mutant and chromatin-modifying chemicals that reduce the rate of RNAPII elongation induce exon skipping whereas chemicals that create a more relaxed chromatin exon inclusion. In the brain of a mouse transgenic for the human FD-ELP1 we observed on this gene an age-dependent decrease in the RNAPII density profile that was most pronounced on the alternative exon, a robust increase in the repressive marks H3K27me3 and H3K9me3 and a decrease of H3K27Ac, together with a progressive reduction in ELP1 exon 20 inclusion level. In HEK 293T cells, selective drug-induced demethylation of H3K27 increased RNAPII elongation on ELP1 and SMN2, promoted the inclusion of the corresponding alternative exons, and, by RNA-sequencing analysis, induced changes in several alternative splicing events. These data suggest a co-transcriptional model of splicing regulation in which age-dependent changes in H3K27me3/Ac modify the rate of RNAPII elongation and affect processing of ELP1 alternative exon 20.

Differential exon usage of developmental genes is associated with deregulated epigenetic marks

Alternative exon usage is known to affect a large portion of genes in mammalian genomes. Importantly, different splice isoforms sometimes possess distinctly different protein functions. Here, we analyzed data from the Human Epigenome Atlas for 11 different human adult tissues and for 8 cultured cells that mimic early developmental stages. We found a significant enrichment of cases where differential usage of exons in various developmental stages of human cells and tissues is associated with differential epigenetic modifications in the flanking regions of individual exons. Many of the genes that were differentially regulated at the exon level and showed deregulated histone marks at the respective exon flanks are functionally associated with development and metabolism.

Histone Marks-Dependent Effect on Alternative Splicing: New Perspectives for Targeted Splicing Modulation in Cancer?

Alternative splicing (AS) is a tightly regulated mechanism that generates the complex human proteome from a small number of genes. Cis-regulatory RNA motifs in exons and introns control AS, recruiting positive and negative trans-acting splicing regulators. At a higher level, chromatin affects splicing events. Growing evidence indicates that the popular histone code hypothesis can be extended to RNA-level processes, such as AS. In addition to nucleosome positioning, which can generate transcriptional barriers to shape the final splicing outcome, histone post-translational modifications can contribute to the detailed regulation of single exon inclusion/exclusion. A histone-based system can identify alternatively spliced chromatin stretches, affecting RNAPII elongation locally or recruiting splicing components via adaptor complexes. In tumor cells, several mechanisms trigger misregulated AS events and produce cancer-associated transcripts. On a genome-wide level, aberrant AS can be the consequence of dysfunctional epigenetic splicing code, including altered enrichment in histone post-translational modifications. This review describes the main findings related to the effect of histone modifications and variants on splicing outcome and how a dysfunctional epigenetic splicing code triggers aberrant AS in cancer. In addition, it highlights recent advances in programmable DNA-targeting technologies and their possible application for AS targeted epigenetic modulation.

Splicing-associated chromatin signatures: a combinatorial and position-dependent role for histone marks in splicing definition

Alternative splicing relies on the combinatorial recruitment of splicing regulators to specific RNA binding sites. Chromatin has been shown to impact this recruitment. However, a limited number of histone marks have been studied at a global level. In this work, a machine learning approach, applied to extensive epigenomics datasets in human H1 embryonic stem cells and IMR90 foetal fibroblasts, has identified eleven chromatin modifications that differentially mark alternatively spliced exons depending on the level of exon inclusion. These marks act in a combinatorial and position-dependent way, creating characteristic splicing-associated chromatin signatures (SACS). In support of a functional role for SACS in coordinating splicing regulation, changes in the alternative splicing of SACS-marked exons between ten different cell lines correlate with changes in SACS enrichment levels and recruitment of the splicing regulators predicted by RNA motif search analysis. We propose the dynamic nature of chromatin modifications as a mechanism to rapidly fine-tune alternative splicing when necessary.

Chromatin is known to regulate splicing by modulating recruitment of splicing factors. Using machine learning approaches, the authors have underlined a chromatin code for alternative splicing regulation that is conserved amongst cell lines.

Modeling circRNA expression pattern with integrated sequence and epigenetic features demonstrates the potential involvement of H3K79me2 in circRNA expression

MOTIVATION

CircRNAs are an abundant class of non-coding RNAs with widespread, cell-/tissue-specific patterns. Previous work suggested that epigenetic features might be related to circRNA expression. However, the contribution of epigenetic changes to circRNA expression has not been investigated systematically. Here, we built a machine learning framework named CIRCScan, to predict circRNA expression in various cell lines based on the sequence and epigenetic features.

RESULTS

The predicted accuracy of the expression status models was high with area under the curve of receiver operating characteristic (ROC) values of 0.89-0.92 and the false-positive rates of 0.17-0.25. Predicted expressed circRNAs were further validated by RNA-seq data. The performance of expression-level prediction models was also good with normalized root-mean-square errors of 0.28-0.30 and Pearson's correlation coefficient r over 0.4 in all cell lines, along with Spearman's correlation coefficient ρ of 0.33-0.46. Noteworthy, H3K79me2 was highly ranked in modeling both circRNA expression status and levels across different cells. Further analysis in additional nine cell lines demonstrated a significant enrichment of H3K79me2 in circRNA flanking intron regions, supporting the potential involvement of H3K79me2 in circRNA expression regulation.

AVAILABILITY AND IMPLEMENTATION

The CIRCScan assembler is freely available online for academic use at https://github.com/johnlcd/CIRCScan.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Specific histone modifications associate with alternative exon selection during mammalian development

Alternative splicing (AS) is frequent during early mouse embryonic development. Specific histone post-translational modifications (hPTMs) have been shown to regulate exon splicing by either directly recruiting splice machinery or indirectly modulating transcriptional elongation. In this study, we hypothesized that hPTMs regulate expression of alternatively spliced genes for specific processes during differentiation. To address this notion, we applied an innovative machine learning approach to relate global hPTM enrichment to AS regulation during mammalian tissue development. We found that specific hPTMs, H3K36me3 and H3K4me1, play a role in skipped exon selection among all the tissues and developmental time points examined. In addition, we used iterative random forest model and found that interactions of multiple hPTMs most strongly predicted splicing when they included H3K36me3 and H3K4me1. Collectively, our data demonstrated a link between hPTMs and alternative splicing which will drive further experimental studies on the functional relevance of these modifications to alternative splicing.

Alternative splicing in human cancer cells is modulated by the amiloride derivative 3,5-diamino-6-chloro-N-(N-(2,6-dichlorobenzoyl)carbamimidoyl)pyrazine-2-carboxide

Alternative splicing (AS) is a process that enables the generation of multiple protein isoforms with different biological properties from a single mRNA. Cancer cells often use the maneuverability conferred by AS to produce proteins that contribute to growth and survival. In our previous studies, we identified that amiloride modulates AS in cancer cells. However, the effective concentration of amiloride required to modulate AS is too high for use in cancer treatment. In this study, we used computational algorithms to screen potential amiloride derivatives for their ability to regulate AS in cancer cells. We found that 3,5-diamino-6-chloro-N-(N-(2,6-dichlorobenzoyl)carbamimidoyl)pyrazine-2-carboxamide (BS008) can regulate AS of apoptotic gene transcripts, including HIPK3, SMAC, and BCL-X, at a lower concentration than amiloride. This splicing regulation involved various splicing factors, and it was accompanied by a change in the phosphorylation state of serine/arginine-rich proteins (SR proteins). RNA sequencing was performed to reveal that AS of many other apoptotic gene transcripts, such as AATF, ATM, AIFM1, NFKB1, and API5, was also modulated by BS008. In vivo experiments further indicated that treatment of tumor-bearing mice with BS008 resulted in a marked decrease in tumor size. BS008 also had inhibitory effects in vitro, either alone or in a synergistic combination with the cytotoxic chemotherapeutic agents sorafenib and nilotinib. BS008 enabled sorafenib dose reduction without compromising antitumor activity. These findings suggest that BS008 may possess therapeutic potential for cancer treatment.

Bayesian nonparametric discovery of isoforms and individual specific quantification

Most human protein-coding genes can be transcribed into multiple distinct mRNA isoforms. These alternative splicing patterns encourage molecular diversity, and dysregulation of isoform expression plays an important role in disease etiology. However, isoforms are difficult to characterize from short-read RNA-seq data because they share identical subsequences and occur in different frequencies across tissues and samples. Here, we develop biisq, a Bayesian nonparametric model for isoform discovery and individual specific quantification from short-read RNA-seq data. biisq does not require isoform reference sequences but instead estimates an isoform catalog shared across samples. We use stochastic variational inference for efficient posterior estimates and demonstrate superior precision and recall for simulations compared to state-of-the-art isoform reconstruction methods. biisq shows the most gains for low abundance isoforms, with 36% more isoforms correctly inferred at low coverage versus a multi-sample method and 170% more versus single-sample methods. We estimate isoforms in the GEUVADIS RNA-seq data and validate inferred isoforms by associating genetic variants with isoform ratios.

Alternative splicing leads to transcript isoform diversity. Here, Aguiar et al. develop biisq, a Bayesian nonparametric approach to discover and quantify isoforms from RNA-seq data.

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.

Aryl hydrocarbon receptor nuclear translocator (ARNT) isoforms control lymphoid cancer cell proliferation through differentially regulating tumor suppressor p53 activity

The aryl hydrocarbon receptor nuclear translocator (ARNT) is involved in xenobiotic and hypoxic responses, and we previously showed that ARNT also regulates nuclear factor-κB (NF-κB) signaling by altering the DNA binding activity of the RelB subunit. However, our initial study of ARNT-mediated RelB modulation was based on simultaneous suppression of the two ARNT isoforms, isoform 1 and 3, and precluded the examination of their individual functions. We find here that while normal lymphocytes harbor equal levels of isoform 1 and 3, lymphoid malignancies exhibit a shift to higher levels of ARNT isoform 1. These elevated levels of ARNT isoform 1 are critical to the proliferation of these cancerous cells, as suppression of isoform 1 in a human multiple myeloma (MM) cell line, and an anaplastic large cell lymphoma (ALCL) cell line, triggered S-phase cell cycle arrest, spontaneous apoptosis, and sensitized cells to doxorubicin treatment. Furthermore, co-suppression of RelB or p53 with ARNT isoform 1 prevented cell cycle arrest and blocked doxorubicin induced apoptosis. Together our findings reveal that certain blood cancers rely on ARNT isoform 1 to potentiate proliferation by antagonizing RelB and p53-dependent cell cycle arrest and apoptosis. Significantly, our results identify ARNT isoform 1 as a potential target for anticancer therapies.

Reduced H3K27me3 Expression Is Common in Nodular Melanomas of Childhood Associated With Congenital Melanocytic Nevi But Not in Proliferative Nodules

The formation of a nodule within a congenital melanocytic nevus (CMN) raises concerns about possible melanoma. Most new nodular growths that develop during childhood, however, are benign proliferative nodules (PN); melanoma is very rare. The distinction of melanoma from PN can at times be difficult clinically and histopathologically, requiring ancillary molecular tests for diagnosis. Although the application of molecular methods has revealed new insights into the mutational and genomic landscape of childhood melanomas, little is known about epigenetic events that may drive the growth of a melanoma or PN in a CMN. In this study we compared the expression of H3K27me3, a key regulator in chromatin remodelling-controlled transcription, in PNs and pediatric nodular melanomas arising within medium-sized to large CMN by immunohistochemistry. Significant loss of H3K27me3 expression was seen in 4 of 5 melanomas, but not in any of the 20 PNs. This observation suggests that epigenetic events likely play a role in the pathogenesis of melanoma developing in the dermis or subcutis of CMN. Furthermore, assessing for H3K27me3 expression by immunohistochemistry may be diagnostically useful for problematic cases.

Histone posttranslational modifications predict specific alternative exon subtypes in mammalian brain

A compelling body of literature, based on next generation chromatin immunoprecipitation and RNA sequencing of reward brain regions indicates that the regulation of the epigenetic landscape likely underlies chronic drug abuse and addiction. It is now critical to develop highly innovative computational strategies to reveal the relevant regulatory transcriptional mechanisms that may underlie neuropsychiatric disease. We have analyzed chromatin regulation of alternative splicing, which is implicated in cocaine exposure in mice. Recent literature has described chromatin-regulated alternative splicing, suggesting a novel function for drug-induced neuroepigenetic remodeling. However, the extent of the genome-wide association between particular histone modifications and alternative splicing remains unexplored. To address this, we have developed novel computational approaches to model the association between alternative splicing and histone posttranslational modifications in the nucleus accumbens (NAc), a brain reward region. Using classical statistical methods and machine learning to combine ChIP-Seq and RNA-Seq data, we found that specific histone modifications are strongly associated with various aspects of differential splicing. H3K36me3 and H3K4me1 have the strongest association with splicing indicating they play a significant role in alternative splicing in brain reward tissue.

Connecting the dots: chromatin and alternative splicing in EMT

Nature has devised sophisticated cellular machinery to process mRNA transcripts produced by RNA Polymerase II, removing intronic regions and connecting exons together, to produce mature RNAs. This process, known as splicing, is very closely linked to transcription. Alternative splicing, or the ability to produce different combinations of exons that are spliced together from the same genomic template, is a fundamental means of regulating protein complexity. Similar to transcription, both constitutive and alternative splicing can be regulated by chromatin and its associated factors in response to various signal transduction pathways activated by external stimuli. This regulation can vary between different cell types, and interference with these pathways can lead to changes in splicing, often resulting in aberrant cellular states and disease. The epithelial to mesenchymal transition (EMT), which leads to cancer metastasis, is influenced by alternative splicing events of chromatin remodelers and epigenetic factors such as DNA methylation and non-coding RNAs. In this review, we will discuss the role of epigenetic factors including chromatin, chromatin remodelers, DNA methyltransferases, and microRNAs in the context of alternative splicing, and discuss their potential involvement in alternative splicing during the EMT process.

Global variability in gene expression and alternative splicing is modulated by mitochondrial content

Noise in gene expression is a main determinant of phenotypic variability. Increasing experimental evidence suggests that genome-wide cellular constraints largely contribute to the heterogeneity observed in gene products. It is still unclear, however, which global factors affect gene expression noise and to what extent. Since eukaryotic gene expression is an energy demanding process, differences in the energy budget of each cell could determine gene expression differences. Here, we quantify the contribution of mitochondrial variability (a natural source of ATP variation) to global variability in gene expression. We find that changes in mitochondrial content can account for ∼50% of the variability observed in protein levels. This is the combined result of the effect of mitochondria dosage on transcription and translation apparatus content and activities. Moreover, we find that mitochondrial levels have a large impact on alternative splicing, thus modulating both the abundance and type of mRNAs. A simple mathematical model in which mitochondrial content simultaneously affects transcription rate and splicing site choice can explain the alternative splicing data. The results of this study show that mitochondrial content (and/or probably function) influences mRNA abundance, translation, and alternative splicing, which ultimately affects cellular phenotype.

Effects of deranged metabolism on epigenetic changes in cancer

The concept of epigenetics is now providing the mechanisms by which cells transfer their new environmental-change-induced phenotypes to their daughter cells. However, how extracellular or cytoplasmic environmental cues are connected to the nuclear epigenome remains incompletely understood. Recently emerging evidence suggests that epigenetic changes are correlated with metabolic changes via chromatin remodeling. As many human complex diseases including cancer harbor both epigenetic changes and metabolic dysregulation, understanding the molecular processes linking them has huge implications for disease pathogenesis and therapeutic intervention. In this review, the impacts of metabolic changes on cancer epigenetics are discussed, along with the current knowledge on cancer metabolism and epigenetics.

Histone modifications involved in cassette exon inclusions: a quantitative and interpretable analysis

Background

Chromatin structure and epigenetic modifications have been shown to involve in the co-transcriptional splicing of RNA precursors. In particular, some studies have suggested that some types of histone modifications (HMs) may participate in the alternative splicing and function as exon marks. However, most existing studies pay attention to the qualitative relationship between epigenetic modifications and exon inclusion. The quantitative analysis that reveals to what extent each type of epigenetic modification is responsible for exon inclusion is very helpful for us to understand the splicing process.

Results

In this paper, we focus on the quantitative analysis of HMs’ influence on the inclusion of cassette exons (CEs) into mature RNAs. With the high-throughput ChIP-seq and RNA-seq data obtained from ENCODE website, we modeled the association of HMs with CE inclusions by logistic regression whose coefficients are meaningful and interpretable for us to reveal the effect of each type of HM. Three type of HMs, H3K36me3, H3K9me3 and H4K20me1, were found to play major role in CE inclusions. HMs’ effect on CE inclusions is conservative across cell types, and does not depend on the expression levels of the genes hosting CEs. HMs located in the flanking regions of CEs were also taken into account in our analysis, and HMs within bounded flanking regions were shown to affect moderately CE inclusions. Moreover, we also found that HMs on CEs whose length is approximately close to nucleosomal-DNA length affect greatly on CE inclusion.

Conclusions

We suggested that a few types of HMs correlate closely to alternative splicing and perhaps function jointly with splicing machinery to regulate the inclusion level of exons. Our findings are helpful to understand HMs’ effect on exon definition, as well as the mechanism of co-transcriptional splicing.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1148) contains supplementary material, which is available to authorized users.