Chromatin and epigenetic regulation of pre-mRNA processing

Seasonal variation in BCG-induced trained immunity

The Bacille Calmette-Guerin (BCG) vaccine is a well-established inducer of innate immune memory (also termed trained immunity), causing increased cytokine production upon heterologous secondary stimulation. Innate immune responses are known to be influenced by season, but whether seasons impact induction of trained immunity is not known. To explore the influence of season on innate immune memory induced by the BCG vaccine, we vaccinated healthy volunteers with BCG either during winter or spring. Three months later, we measured the ex vivo cytokine responses against heterologous stimuli, analyzed gene expressions and epigenetic signatures of the immune cells, and compared these with the baseline before vaccination. BCG vaccination during winter induced a stronger increase in the production of pro-inflammatory cytokines by peripheral blood mononuclear cells (PBMCs) upon stimulation with different bacterial and fungal stimuli, compared to BCG vaccination in spring. In contrast, winter BCG vaccination resulted in lower IFNγ release in PBMCs compared to spring BCG vaccination. Furthermore, NK cells of the winter-vaccinated people had a greater pro-inflammatory cytokine and IFNγ production capacity upon heterologous stimulation. BCG had only minor effects on the transcriptome of monocytes 3 months later. In contrast, we identified season-dependent epigenetic changes in monocytes and NK cells induced by vaccination, partly explaining the higher immune cell reactivity in the winter BCG vaccination group. These results suggest that BCG vaccination during winter is more prone to induce a robust trained immunity response by activating and reprogramming the immune cells, especially NK cells. (Dutch clinical trial registry no. NL58219.091.16).

N6-Methyladenosine mRNA Modification: From Modification Site Selectivity to Neurological Functions

The development of various chemical methods has enabled scientists to decipher the distribution features and biological functions of RNA modifications in the past decade. In addition to modifying noncoding RNAs such as tRNAs and rRNAs, N6-methyladenosine (m6A) has been proven to be the most abundant internal chemical modification on mRNAs in eukaryotic cells and is also the most widely studied mRNA modification to date. Extensive studies have repeatedly demonstrated the important functions of m6A in various biological conditions, ranging from embryonic organ development to adult organ function and pathogenesis. Unlike DNA methylation which is relatively stable, the reversible m6A modification on mRNA is highly dynamic and easily influenced by various internal or external factors, such as cell type, developmental stage, nutrient supply, circadian rhythm, and environmental stresses.In this Account, we review our previous findings on the site selectivity mechanisms regulating m6A formation, as well as the physiological roles of m6A modification in cerebellum development and long-term memory consolidation. In our initial efforts to profile m6A in various types of mouse and human cells, we surprisingly found that the sequence motifs surrounding m6A sites were often complementary with the seed sequences of miRNAs. By manipulating the abundance of the miRNA biogenesis enzyme Dicer or individual miRNAs or mutating miRNA sequences, we were able to reveal a new role of nucleus localized miRNAs, which is to guide the m6A methyltransferase METTL3 to bind to mRNAs and to promote m6A formation. As a result, we partially answered the question of why only a small proportion of m6A motifs within an mRNA could have m6A modification at a certain time point. We further explored the functions of m6A modification in regulating brain development and brain functions. We found that cerebellum had the most severe defects when Mettl3 was knocked out in developing mouse embryonic brain and revealed that the underlying mechanisms could be attributed to aberrant mRNA splicing and enhanced cell apoptosis under m6A deficit conditions. On the other hand, knocking out Mettl3 in postnatal hippocampus did not cause morphological defects in the mouse brain but impaired the efficacy of long-term memory consolidation. Under learning stimuli, formation of m6A modifications could be detected on transcripts encoding proteins related to dendrite growth, synapse formation, and other memory related functions. Loss of m6A modifications on these transcripts would result in translation deficiency and reduced protein production, particularly in the translation of early response genes, and therefore would compromise the efficacy of long-term memory consolidation. Interestingly, excessive training sessions or increased training intensity could overcome such m6A deficiency related memory defects, which is likely due to the longer turnover cycle and the cumulative abundance of proteins throughout the training process. In addition to revealing the roles of m6A modification in regulating long-term memory formation, our work also demonstrated an effective method for studying memory formation efficacy. As the lack of an appropriate model for studying memory formation efficacy has been a long-lasting problem in the field of neural science, our hippocampus-specific postnatal m6A knockout model could also be utilized to study other questions related to memory formation efficacy.

Associations between DNA methylation and gene regulation depend on chromatin accessibility during transgenerational plasticity

BACKGROUND

Epigenetic processes are proposed to be a mechanism regulating gene expression during phenotypic plasticity. However, environmentally induced changes in DNA methylation exhibit little-to-no association with differential gene expression in metazoans at a transcriptome-wide level. It remains unexplored whether associations between environmentally induced differential methylation and expression are contingent upon other epigenomic processes such as chromatin accessibility. We quantified methylation and gene expression in larvae of the purple sea urchin Strongylocentrotus purpuratus exposed to different ecologically relevant conditions during gametogenesis (maternal conditioning) and modeled changes in gene expression and splicing resulting from maternal conditioning as functions of differential methylation, incorporating covariates for genomic features and chromatin accessibility. We detected significant interactions between differential methylation, chromatin accessibility, and genic feature type associated with differential expression and splicing.

RESULTS

Differential gene body methylation had significantly stronger effects on expression among genes with poorly accessible transcriptional start sites while baseline transcript abundance influenced the direction of this effect. Transcriptional responses to maternal conditioning were 4-13 × more likely when accounting for interactions between methylation and chromatin accessibility, demonstrating that the relationship between differential methylation and gene regulation is partially explained by chromatin state.

CONCLUSIONS

DNA methylation likely possesses multiple associations with gene regulation during transgenerational plasticity in S. purpuratus and potentially other metazoans, but its effects are dependent on chromatin accessibility and underlying genic features.

PARP1 Regulates Circular RNA Biogenesis though Control of Transcriptional Dynamics

Circular RNAs (circRNAs) are a recently discovered class of RNAs derived from protein-coding genes that have important biological and pathological roles. They are formed through backsplicing during co-transcriptional alternative splicing; however, the unified mechanism that accounts for backsplicing decisions remains unclear. Factors that regulate the transcriptional timing and spatial organization of pre-mRNA, including RNAPII kinetics, the availability of splicing factors, and features of gene architecture, have been shown to influence backsplicing decisions. Poly (ADP-ribose) polymerase I (PARP1) regulates alternative splicing through both its presence on chromatin as well as its PARylation activity. However, no studies have investigated PARP1's possible role in regulating circRNA biogenesis. Here, we hypothesized that PARP1's role in splicing extends to circRNA biogenesis. Our results identify many unique circRNAs in PARP1 depletion and PARylation-inhibited conditions compared to the wild type. We found that while all genes producing circRNAs share gene architecture features common to circRNA host genes, genes producing circRNAs in PARP1 knockdown conditions had longer upstream introns than downstream introns, whereas flanking introns in wild type host genes were symmetrical. Interestingly, we found that the behavior of PARP1 in regulating RNAPII pausing is distinct between these two classes of host genes. We conclude that the PARP1 pausing of RNAPII works within the context of gene architecture to regulate transcriptional kinetics, and therefore circRNA biogenesis. Furthermore, this regulation of PARP1 within host genes acts to fine tune their transcriptional output with implications in gene function.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

Mechanistic aspects of reversible methylation modifications of arginine and lysine of nuclear histones and their roles in human colon cancer

Developmental proceedings and maintenance of cellular homeostasis are regulated by the precise orchestration of a series of epigenetic events that eventually control gene expression. DNA methylation and post-translational modifications (PTMs) of histones are well-characterized epigenetic events responsible for fine-tuning gene expression. PTMs of histones bear molecular logic of gene expression at chromosomal territory and have become a fascinating field of epigenetics. Nowadays, reversible methylation on histone arginine and lysine is gaining increasing attention as a significant PTM related to reorganizing local nucleosomal structure, chromatin dynamics, and transcriptional regulation. It is now well-accepted and reported that histone marks play crucial roles in colon cancer initiation and progression by encouraging abnormal epigenomic reprogramming. It is becoming increasingly clear that multiple PTM marks at the N-terminal tails of the core histones cross-talk with one another to intricately regulate DNA-templated biological processes such as replication, transcription, recombination, and damage repair in several malignancies, including colon cancer. These functional cross-talks provide an additional layer of message, which spatiotemporally fine-tunes the overall gene expression regulation. Nowadays, it is evident that several PTMs instigate colon cancer development. How colon cancer-specific PTM patterns or codes are generated and how they affect downstream molecular events are uncovered to some extent. Future studies would address more about epigenetic communication, and the relationship between histone modification marks to define cellular functions in depth. This chapter will comprehensively highlight the importance of histone arginine and lysine-based methylation modifications and their functional cross-talk with other histone marks from the perspective of colon cancer development.

Association of DNA methylation with steroidogenic enzymes in Cushing's adenoma

DNA methylation and demethylation regulate the transcription of genes. DNA methylation-associated gene expression of adrenal steroidogenic enzymes may regulate cortisol production in cortisol-producing adenoma (CPA). We aimed to determine the DNA methylation levels of all genes encoding steroidogenic enzymes involved in CPA. Additionally, the aims were to clarify the DNA methylation-associated gene expression and evaluate the difference of CPA genotype from others using DNA methylation data. Twenty-five adrenal CPA and six nonfunctioning adrenocortical adenoma (NFA) samples were analyzed. RNA sequencing and DNA methylation array were performed. The methylation levels at 118 methylation sites of the genes were investigated, and their methylation and mRNA levels were subsequently integrated. Among all the steroidogenic enzyme genes studied, CYP17A1 gene was mainly found to be hypomethylated in CPA compared to that in NFA, and the Benjamini-Hochberg procedure demonstrated that methylation levels at two sites in the CYP17A1 gene body were statistically significant. PRKACA mutant CPAs predominantly exhibited hypomethylation of CYP17A1 gene compared with the GNAS mutant CPAs. Inverse associations between CYP17A1 methylation in three regions of the gene body and its mRNA levels were observed in the NFAs and CPAs. In applying clustering analysis using CYP17A1 methylation and mRNA levels, CPAs with PRKACA mutation were differentiated from NFAs and CPAs with a GNAS mutation. We demonstrated that CPAs exhibited hypomethylation of the CYP17A1 gene body in CPA, especially in the PRKACA mutant CPAs. Methylation of CYP17A1 gene may influence its transcription levels.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.

Genetic variant rs10251977 (G>A) in EGFR-AS1 modulates the expression of EGFR isoforms A and D

Tyrosine kinase inhibitor is an effective chemo-therapeutic drug against tumors with deregulated EGFR pathway. Recently, a genetic variant rs10251977 (G>A) in exon 20 of EGFR reported to act as a prognostic marker for HNSCC. Genotyping of this polymorphism in oral cancer patients showed a similar frequency in cases and controls. EGFR-AS1 expressed significantly high level in tumors and EGFR-A isoform expression showed significant positive correlation (r = 0.6464, p < 0.0001) with reference to EGFR-AS1 expression levels, consistent with larger TCGA HNSCC tumor dataset. Our bioinformatic analysis showed enrichment of alternative splicing marks H3K36me3 and presence of intronic polyA sites spanning around exon 15a and 15b of EGFR facilitates skipping of exon 15b, thereby promoting the splicing of EGFR-A isoform. In addition, high level expression of PTBP1 and its binding site in EGFR and EGFR-AS1 enhances the expression of EGFR-A isoform (r = 0.7404, p < 0.0001) suggesting that EGFR-AS1 expression modulates the EGFR-A and D isoforms through alternative splicing. In addition, this polymorphism creates a binding site for miR-891b in EGFR-AS1 and may negatively regulate the EGFR-A. Collectively, our results suggested the presence of genetic variant in EGFR-AS1 modulates the expression of EGFR-D and A isoforms.

Methylation multiplicity and its clinical values in cancer

Methylation at DNA, RNA and protein levels plays critical roles in many cellular processes and is associated with diverse differentiation events, physiological activities and human diseases. To aid in the diagnostic and therapeutic design for cancer treatment utilising methylation, this review provides a boutique yet comprehensive overview on methylation at different levels including the mechanisms, cross-talking and clinical implications with a particular focus on cancers. We conclude that DNA methylation is the sole type of methylation that has been largely translated into clinics and used for, mostly, early diagnosis. Translating the onco-therapeutic and prognostic values of RNA and protein methylations into clinical use deserves intensive efforts. Simultaneous examination of methylations at multiple levels or together with other forms of molecular markers represents an interesting research direction with profound clinical translational potential.

Crosstalk Between mRNA 3'-End Processing and Epigenetics

The majority of eukaryotic genes produce multiple mRNA isoforms by using alternative poly(A) sites in a process called alternative polyadenylation (APA). APA is a dynamic process that is highly regulated in development and in response to extrinsic or intrinsic stimuli. Mis-regulation of APA has been linked to a wide variety of diseases, including cancer, neurological and immunological disorders. Since the first example of APA was described 40 years ago, the regulatory mechanisms of APA have been actively investigated. Conventionally, research in this area has focused primarily on the roles of regulatory cis-elements and trans-acting RNA-binding proteins. Recent studies, however, have revealed important functions for epigenetic mechanisms, including DNA and histone modifications and higher-order chromatin structures, in APA regulation. Here we will discuss these recent findings and their implications for our understanding of the crosstalk between epigenetics and mRNA 3'-end processing.

Statistical inference of the rate of RNA polymerase II elongation by total RNA sequencing

Motivation

Sequencing total RNA without poly-A selection enables us to obtain a transcriptomic profile of nascent RNAs undergoing transcription with co-transcriptional splicing. In general, the RNA-seq reads exhibit a sawtooth pattern in a gene, which is characterized by a monotonically decreasing gradient across introns in the 5’–3’ direction, and by substantially higher levels of RNA-seq reads present in exonic regions. Such patterns result from the process of underlying transcription elongation by RNA polymerase II, which traverses the DNA strand in a 5’–3’ direction as it performs a complex series of mRNA synthesis and processing. Therefore, data of sequenced total RNAs could be utilized to infer the rate of transcription elongation by solving the inverse problem.

Results

Though solving the inverse problem in total RNA-seq has the great potential, statistical methods have not yet been fully developed. We demonstrate what extent the newly developed method can be useful. The objective is to reconstruct the spatial distribution of transcription elongation rates in a gene from a given noisy, sawtooth-like profile. It is necessary to recover the signal source of the elongation rates separately from several types of nuisance factors, such as unobserved modes of co-transcriptionally occurring mRNA splicing, which exert significant influences on the sawtooth shape. The present method was tested using published total RNA-seq data derived from mouse embryonic stem cells. We investigated the spatial characteristics of the estimated elongation rates, focusing especially on the relation to promoter-proximal pausing of RNA polymerase II, nucleosome occupancy and histone modification patterns.

Availability and implementation

A C implementation of PolSter and sample data are available at https://github.com/yoshida-lab/PolSter.

Supplementary information

Supplementary data are available at Bioinformatics online.

DAPK3 participates in the mRNA processing of immediate early genes in chronic lymphocytic leukaemia

Cross‐linking of the B‐cell receptor (BCR) induces transcriptional activation of immediate early genes (IEGs) including EGR1 and DUSP2 in chronic lymphocytic leukaemia (CLL). Here, we have shown that this transcriptional activation correlated with histone H3 threonine 6 and 11 phosphorylation. Both transcription and histone post‐translational modifications are repressed by ibrutinib, a small molecule inhibitor used in CLL treatment. Moreover, we have identified the death‐associated protein kinase 3 (DAPK3), as the kinase mediating these histone phosphorylation marks in response to activation of the BCR signalling pathway with this kinase being recruited to RNA polymerase II in an anti‐IgM‐dependent manner. DAPK inhibition mimics ibrutinib‐induced repression of both IEG mRNA and histone H3 phosphorylation and has anti‐proliferative effect comparable to ibrutinib in CLL in vitro. DAPK inhibitor does not repress transcription itself but impacts on mRNA processing and has a broader anti‐tumour effect than ibrutinib, by repressing both anti‐IgM‐ and CD40L‐dependent activation.

EWS-FLI1 modulated alternative splicing of ARID1A reveals novel oncogenic function through the BAF complex

Connections between epigenetic reprogramming and transcription or splicing create novel mechanistic networks that can be targeted with tailored therapies. Multiple subunits of the chromatin remodeling BAF complex, including ARID1A, play a role in oncogenesis, either as tumor suppressors or oncogenes. Recent work demonstrated that EWS-FLI1, the oncogenic driver of Ewing sarcoma (ES), plays a role in chromatin regulation through interactions with the BAF complex. However, the specific BAF subunits that interact with EWS-FLI1 and the precise role of the BAF complex in ES oncogenesis remain unknown. In addition to regulating transcription, EWS-FLI1 also alters the splicing of many mRNA isoforms, but the role of splicing modulation in ES oncogenesis is not well understood. We have identified a direct connection between the EWS-FLI1 protein and ARID1A isoform protein variant ARID1A-L. We demonstrate here that ARID1A-L is critical for ES maintenance and supports oncogenic transformation. We further report a novel feed-forward cycle in which EWS-FLI1 leads to preferential splicing of ARID1A-L, promoting ES growth, and ARID1A-L reciprocally promotes EWS-FLI1 protein stability. Dissecting this interaction may lead to improved cancer-specific drug targeting.

Coordinated Control of mRNA and rRNA Processing Controls Embryonic Stem Cell Pluripotency and Differentiation

Stem cell-specific transcriptional networks are well known to control pluripotency, but constitutive cellular processes such as mRNA splicing and protein synthesis can add complex layers of regulation with poorly understood effects on cell-fate decisions. Here, we show that the RNA binding protein HTATSF1 controls embryonic stem cell differentiation by regulating multiple aspects of RNA processing during ribosome biogenesis. HTATSF1, in a complex with splicing factor SF3B1, controls intron removal from ribosomal protein transcripts and regulates ribosomal RNA transcription and processing, thereby controlling 60S ribosomal abundance and protein synthesis. HTATSF1-dependent protein synthesis is essential for naive pre-implantation epiblast to transition into post-implantation epiblast, a stage with transiently low protein synthesis, and further differentiation toward neuroectoderm. Together, these results identify coordinated regulation of ribosomal RNA and protein synthesis by HTATSF1 and show that this essential mechanism controls protein synthesis during early mammalian embryogenesis.

Increased intron retention is a post-transcriptional signature associated with progressive aging and Alzheimer's disease

Intron retention (IR) by alternative splicing is a conserved regulatory mechanism that can affect gene expression and protein function during adult development and age-onset diseases. However, it remains unclear whether IR undergoes spatial or temporal changes during different stages of aging or neurodegeneration like Alzheimer's disease (AD). By profiling the transcriptome of Drosophila head cells at different ages, we observed a significant increase in IR events for many genes during aging. Differential IR affects distinct biological functions at different ages and occurs at several AD-associated genes in older adults. The increased nucleosome occupancy at the differentially retained introns in young animals suggests that it may regulate the level of IR during aging. Notably, an increase in the number of IR events was also observed in healthy older mouse and human brain tissues, as well as in the cerebellum and frontal cortex from independent AD cohorts. Genes with differential IR shared many common features, including shorter intron length, no perturbation in their mRNA level, and enrichment for biological functions that are associated with mRNA processing and proteostasis. The differentially retained introns identified in AD frontal cortex have higher GC content, with many of their mRNA transcripts showing an altered level of protein expression compared to control samples. Taken together, our results suggest that an increased IR is an conserved signature that is associated with aging. By affecting pathways involved in mRNA and protein homeostasis, changes of IR pattern during aging may regulate the transition from healthy to pathological state in late-onset sporadic AD.

The potential role of a retrotransposed gene and a long noncoding RNA in regulating an X-linked chromatin gene (KDM5C): Novel epigenetic mechanism in autism

A growing body of evidence supports the potential role of the circadian system and chromatin remodeling genes in autism. Considering the heterogeneity and gender discrepancy in autism, and the complex nature of the epigenetic landscape, identification of biologically relevant epigenetic factors requires reducing heterogeneity using proper subtyping. For this study, we used X chromosome inactivation (XCI) status in females with autism as an epigenetic marker for subtyping and examined the expression level of members of KDM5, a chromatin remodeling gene family. KDM5 are histone demethylases involved in the circadian molecular machinery. We used human blood samples to characterize alternatively spliced KDM5 isoforms and noticed that KDM5C undergoes a complex splicing process. We also identified a KDM5C isoform (KDM5C-3'UTR-lncRNA) containing a novel 3'UTR originated from a retrotransposed gene (retro-SUV39H2) of an autosomal methyltransferase (SUV39H2). This 3'UTR shows 84% sequence homology with long ncRNAs (lncRNAs) and is located 32 kb downstream of KDM5C. The KDM5C-3'UTR-lncRNA isoform was differentially expressed in autistic females with XCI skewness compared with controls. KDM5C plays a crucial role in balancing histone H3K4 methylation states. The identified retro-SUV39H2 originated lncRNA also shows H3K4 marks. By assessing the expression level of alternatively spliced Kdm5 isoforms at different circadian time-points, we showed that some isoforms follow a circadian oscillation pattern in wild type mouse brain.This study provides the first evidence and a suggestive model for the potential role of retrotransposed elements in autism through linking methylases and demethylases, two functionally complementary components of chromatin remodeling, which may collectively contribute to disease etiology through lncRNAs. Autism Res 2019, 12: 1007-1021. © 2019 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: Genes do not function in isolated conditions and their proper expression level also depends on a mechanism called gene regulation. An example of gene regulation is when changes outside DNA sequences influence the function of autism susceptibility genes. Alternative splicing is one type of gene regulation, which produces several versions of a gene (called variants) that may slightly differ from each other and be expressed at different levels in response to environmental changes. The circadian clock is an essential timing mechanism that enables organisms to maintain internal processes in sync with the dynamic environment brought about by the day-night cycle. The goal of this study was to assess if a subset of females with autism with certain genetic marker had a unique pattern of alternative splicing of three circadian genes. We identified a novel variant that is differentially expressed in this subset. Our study provides a novel subject stratification strategy, and a suggestive model of how biologically relevant components of a gene regulatory process may be linked and, possibly, collectively contribute to the etiology of autism.

Potential epigenomic co-management in rare diseases and epigenetic therapy

The purpose of this review is to highlight the impact of the alternative splicing process on human disease. Epigenetic regulation determines not only what parts of the genome are expressed but also how they are spliced. The recent progress in the field of epigenetics has important implications for the study of rare diseases. The role of epigenetics in rare diseases is a key issue in molecular physiology and medicine because not only rare diseases can benefit from epigenetic research, but can also provide useful principles for other common and complex disorders such as cancer, cardiovascular, type 2 diabetes, obesity, and neurological diseases. Predominantly, epigenetic modifications include DNA methylation, histone modification, and RNA-associated silencing. These modifications in the genome regulate numerous cellular activities. Disruption of epigenetic regulation process can contribute to the etiology of numerous diseases during both prenatal and postnatal life. Here, I discuss current knowledge about this matter including some current epigenetic therapies and future directions in the field by emphasizing on the RNA-based therapy via antisense oligonucleotides to correct splicing defects.

Splicing and neurodegeneration: Insights and mechanisms

Splicing is the global cellular process whereby intervening sequences (introns) in precursor messenger RNA (pre-mRNA) are removed and expressed regions (exons) are ligated together, resulting in a mature mRNA transcript that is exported and translated in the cytoplasm. The tightly regulated splicing cycle is also flexible allowing for the inclusion or exclusion of some sequences depending on the specific cellular context. Alternative splicing allows for the generation of many transcripts from a single gene, thereby expanding the proteome. Although all cells require the function of the spliceosome, neurons are highly sensitive to splicing perturbations with numerous neurological diseases linked to splicing defects. The sensitivity of neurons to splicing alterations is largely due to the complex neuronal cell types and functions in the nervous system that require specific splice isoforms to maintain cellular homeostasis. In the past several years, the relationship between RNA splicing and the nervous system has been the source of significant investigation. Here, we review the current knowledge on RNA splicing in neurobiology and discuss its potential role and impact in neurodegenerative diseases. We will examine the impact of alternative splicing and the role of splicing regulatory proteins on neurodegeneration, highlighting novel animal models including mouse and zebrafish. We will also examine emerging technologies and therapeutic interventions that aim to "drug" the spliceosome. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Development.

The chromatin remodeler ZmCHB101 impacts alternative splicing contexts in response to osmotic stress

Maize SWI3-type chromatin remodeler impacts alternative splicing contexts in response to osmotic stress by altering nucleosome density and affecting transcriptional elongation rate. Alternative splicing (AS) is commonly found in higher eukaryotes and is an important posttranscriptional regulatory mechanism to generate transcript diversity. AS has been widely accepted as playing essential roles in different biological processes including growth, development, signal transduction and responses to biotic and abiotic stresses in plants. However, whether and how chromatin remodeling complex functions in AS in plant under osmotic stress remains unknown. Here, we show that a maize SWI3D protein, ZmCHB101, impacts AS contexts in response to osmotic stress. Genome-wide analysis of mRNA contexts in response to osmotic stress using ZmCHB101-RNAi lines reveals that ZmCHB101 impacts alternative splicing contexts of a subset of osmotic stress-responsive genes. Intriguingly, ZmCHB101-mediated regulation of gene expression and AS is largely uncoupled, pointing to diverse molecular functions of ZmCHB101 in transcriptional and posttranscriptional regulation. We further found ZmCHB101 impacts the alternative splicing contexts by influencing alteration of chromatin and histone modification status as well as transcriptional elongation rates mediated by RNA polymerase II. Taken together, our findings suggest a novel insight of how plant chromatin remodeling complex impacts AS under osmotic stress .

StereoGene: rapid estimation of genome-wide correlation of continuous or interval feature data

Motivation

Genomics features with similar genome-wide distributions are generally hypothesized to be functionally related, for example, colocalization of histones and transcription start sites indicate chromatin regulation of transcription factor activity. Therefore, statistical algorithms to perform spatial, genome-wide correlation among genomic features are required.

Results

Here, we propose a method, StereoGene, that rapidly estimates genome-wide correlation among pairs of genomic features. These features may represent high-throughput data mapped to reference genome or sets of genomic annotations in that reference genome. StereoGene enables correlation of continuous data directly, avoiding the data binarization and subsequent data loss. Correlations are computed among neighboring genomic positions using kernel correlation. Representing the correlation as a function of the genome position, StereoGene outputs the local correlation track as part of the analysis. StereoGene also accounts for confounders such as input DNA by partial correlation. We apply our method to numerous comparisons of ChIP-Seq datasets from the Human Epigenome Atlas and FANTOM CAGE to demonstrate its wide applicability. We observe the changes in the correlation between epigenomic features across developmental trajectories of several tissue types consistent with known biology and find a novel spatial correlation of CAGE clusters with donor splice sites and with poly(A) sites. These analyses provide examples for the broad applicability of StereoGene for regulatory genomics.

Availability and implementation

The StereoGene C ++ source code, program documentation, Galaxy integration scripts and examples are available from the project homepage http://stereogene.bioinf.fbb.msu.ru/.

Contact

favorov@sensi.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

A reassessment of DNA-immunoprecipitation-based genomic profiling

DNA immunoprecipitation followed by sequencing (DIP-seq) is a common enrichment method for profiling DNA modifications in mammalian genomes. However, the results of independent DIP-seq studies often show considerable variation between profiles of the same genome and between profiles obtained by alternative methods. Here we show that these differences are primarily due to the intrinsic affinity of IgG for short unmodified DNA repeats. This pervasive experimental error accounts for 50-99% of regions identified as 'enriched' for DNA modifications in DIP-seq data. Correction of this error profoundly altered DNA-modification profiles for numerous cell types, including mouse embryonic stem cells, and subsequently revealed novel associations among DNA modifications, chromatin modifications and biological processes. We conclude that both matched input and IgG controls are essential in order for the results of DIP-based assays to be interpreted correctly, and that complementary, non-antibody-based techniques should be used to validate DIP-based findings to avoid further misinterpretation of genome-wide profiling data.

Epigenetics of Circadian Rhythms in Imprinted Neurodevelopmental Disorders

DNA sequence information alone cannot account for the immense variability between chromosomal alleles within diverse cell types in the brain, whether these differences are observed across time, cell type, or parental origin. The complex control and maintenance of gene expression and modulation are regulated by a multitude of molecular and cellular mechanisms that layer on top of the genetic code. The integration of genetic and environmental signals required for regulating brain development and function is achieved in part by a dynamic epigenetic landscape that includes DNA methylation, histone modifications, and noncoding RNAs. These epigenetic mechanisms establish and maintain core biological processes, including genomic imprinting and entrainment of circadian rhythms. This chapter will focus on how the epigenetic layers of DNA methylation and long, noncoding RNAs interact with circadian rhythms at specific imprinted chromosomal loci associated with the human neurodevelopmental disorders Prader-Willi, Angelman, Kagami-Ogata, and Temple syndromes.

Sleep Deprivation and the Epigenome

Sleep deprivation disrupts the lives of millions of people every day and has a profound impact on the molecular biology of the brain. These effects begin as changes within a neuron, at the DNA and RNA level, and result in alterations in neuronal plasticity and dysregulation of many cognitive functions including learning and memory. The epigenome plays a critical role in regulating gene expression in the context of memory storage. In this review article, we begin by describing the effects of epigenetic alterations on the regulation of gene expression, focusing on the most common epigenetic mechanisms: (i) DNA methylation; (ii) histone modifications; and (iii) non-coding RNAs. We then discuss evidence suggesting that sleep loss impacts the epigenome and that these epigenetic alterations might mediate the changes in cognition seen following disruption of sleep. The link between sleep and the epigenome is only beginning to be elucidated, but clear evidence exists that epigenetic alterations occur following sleep deprivation. In the future, these changes to the epigenome could be utilized as biomarkers of sleep loss or as therapeutic targets for sleep-related disorders.

DNA methylation in epigenetic inheritance of metabolic diseases through the male germ line

The global rise in metabolic diseases can be attributed to a complex interplay between biology, behavior and environmental factors. This article reviews the current literature concerning DNA methylation-based epigenetic inheritance (intergenerational and transgenerational) of metabolic diseases through the male germ line. Included are a presentation of the basic principles for DNA methylation in developmental programming, and a description of windows of susceptibility for the inheritance of environmentally induced aberrations in DNA methylation and their associated metabolic disease phenotypes. To this end, escapees, genomic regions with the intrinsic potential to transmit acquired paternal epigenetic information across generations by escaping the extensive programmed DNA demethylation that occurs during gametogenesis and in the zygote, are described. The ongoing descriptive and functional examinations of DNA methylation in the relevant biological samples, in conjugation with analyses of non-coding RNA and histone modifications, hold promise for improved delineation of the effect size and mechanistic background for epigenetic inheritance of metabolic diseases.

The Expanding Landscape of Alternative Splicing Variation in Human Populations

Alternative splicing is a tightly regulated biological process by which the number of gene products for any given gene can be greatly expanded. Genomic variants in splicing regulatory sequences can disrupt splicing and cause disease. Recent developments in sequencing technologies and computational biology have allowed researchers to investigate alternative splicing at an unprecedented scale and resolution. Population-scale transcriptome studies have revealed many naturally occurring genetic variants that modulate alternative splicing and consequently influence phenotypic variability and disease susceptibility in human populations. Innovations in experimental and computational tools such as massively parallel reporter assays and deep learning have enabled the rapid screening of genomic variants for their causal impacts on splicing. In this review, we describe technological advances that have greatly increased the speed and scale at which discoveries are made about the genetic variation of alternative splicing. We summarize major findings from population transcriptomic studies of alternative splicing and discuss the implications of these findings for human genetics and medicine.

Microdroplet PCR for Highly Multiplexed Targeted Bisulfite Sequencing

Many methods exist for examining CpG DNA methylation. However, many of these are qualitative, laborious to apply to a large number of genes simultaneously, or are not easy to target to specific regions of interest. Microdroplet PCR-based bisulfite sequencing allows for quantitative single base resolution analysis of investigator selected regions of interest. Following bisulfite conversion of genomic DNA, targeted microdroplet PCR is conducted with custom primer libraries. Samples are then fragmented, concatenated, and sequenced by high-throughput sequencing. The most recent technology allows for this method to be conducted with as little as 250 ng of bisulfite-converted DNA. The primary advantage of this method is the ability to hand-select the targeted regions covered by up to 10,000 amplicons of 500-600 bp. Moreover, the nature of microdroplet PCR virtually eliminates PCR bias and allows for the amplification of all targets simultaneously in a single tube.

Dependence-induced increase of alcohol self-administration and compulsive drinking mediated by the histone methyltransferase PRDM2

Epigenetic processes have been implicated in the pathophysiology of alcohol dependence, but the specific molecular mechanisms mediating dependence-induced neuroadaptations remain largely unknown. Here, we found that a history of alcohol dependence persistently decreased the expression of Prdm2, a histone methyltransferase that monomethylates histone 3 at the lysine 9 residue (H3K9me1), in the rat dorsomedial prefrontal cortex (dmPFC). Downregulation of Prdm2 was associated with decreased H3K9me1, supporting that changes in Prdm2 mRNA levels affected its activity. Chromatin immunoprecipitation followed by massively parallel DNA sequencing showed that genes involved in synaptic communication are epigenetically regulated by H3K9me1 in dependent rats. In non-dependent rats, viral-vector-mediated knockdown of Prdm2 in the dmPFC resulted in expression changes similar to those observed following a history of alcohol dependence. Prdm2 knockdown resulted in increased alcohol self-administration, increased aversion-resistant alcohol intake and enhanced stress-induced relapse to alcohol seeking, a phenocopy of postdependent rats. Collectively, these results identify a novel epigenetic mechanism that contributes to the development of alcohol-seeking behavior following a history of dependence.

The histone variant H2A.Z promotes efficient cotranscriptional splicing in S. cerevisiae

In eukaryotes, a dynamic ribonucleic protein machine known as the spliceosome catalyzes the removal of introns from premessenger RNA (pre-mRNA). Recent studies show the processes of RNA synthesis and RNA processing to be spatio-temporally coordinated, indicating that RNA splicing takes place in the context of chromatin. H2A.Z is a highly conserved histone variant of the canonical histone H2A. In Saccharomyces cerevisiae, H2A.Z is deposited into chromatin by the SWR-C complex, is found near the 5' ends of protein-coding genes, and has been implicated in transcription regulation. Here we show that splicing of intron-containing genes in cells lacking H2A.Z is impaired, particularly under suboptimal splicing conditions. Cells lacking H2A.Z are especially dependent on a functional U2 snRNP (small nuclear RNA [snRNA] plus associated proteins), as H2A.Z shows extensive genetic interactions with U2 snRNP-associated proteins, and RNA sequencing (RNA-seq) reveals that introns with nonconsensus branch points are particularly sensitive to H2A.Z loss. Consistently, H2A.Z promotes efficient spliceosomal rearrangements involving the U2 snRNP, as H2A.Z loss results in persistent U2 snRNP association and decreased recruitment of downstream snRNPs to nascent RNA. H2A.Z impairs transcription elongation, suggesting that spliceosome rearrangements are tied to H2A.Z's role in elongation. Depletion of disassembly factor Prp43 suppresses H2A.Z-mediated splice defects, indicating that, in the absence of H2A.Z, stalled spliceosomes are disassembled, and unspliced RNAs are released. Together, these data demonstrate that H2A.Z is required for efficient pre-mRNA splicing and indicate a role for H2A.Z in coordinating the kinetics of transcription elongation and splicing.

Histone H3 lysine 36 methylation affects temperature-induced alternative splicing and flowering in plants

Background

Global warming severely affects flowering time and reproductive success of plants. Alternative splicing of pre-messenger RNA (mRNA) is an important mechanism underlying ambient temperature-controlled responses in plants, yet its regulation is poorly understood. An increase in temperature promotes changes in plant morphology as well as the transition from the vegetative to the reproductive phase in Arabidopsis thaliana via changes in splicing of key regulatory genes. Here we investigate whether a particular histone modification affects ambient temperature-induced alternative splicing and flowering time.

Results

We use a genome-wide approach and perform RNA-sequencing (RNA-seq) analyses and histone H3 lysine 36 tri-methylation (H3K36me3) chromatin immunoprecipitation sequencing (ChIP-seq) in plants exposed to different ambient temperatures. Analysis and comparison of these datasets reveal that temperature-induced differentially spliced genes are enriched in H3K36me3. Moreover, we find that reduction of H3K36me3 deposition causes alteration in temperature-induced alternative splicing. We also show that plants with mutations in H3K36me3 writers, eraser, or readers have altered high ambient temperature-induced flowering.

Conclusions

Our results show a key role for the histone mark H3K36me3 in splicing regulation and plant plasticity to fluctuating ambient temperature. Our findings open new perspectives for the breeding of crops that can better cope with environmental changes due to climate change.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1235-x) contains supplementary material, which is available to authorized users.

Differences in expression of genes in the MyD88 and TRIF signalling pathways and methylation of TLR4 and TRIF in Tibetan chickens and DaHeng S03 chickens infected with Salmonella enterica serovar enteritidis

Salmonella enterica serovar (S. enteritidis) is a pathogenic bacterium that can cause symptoms of food poisoning, leading to death of poultry, resulting in serious economic losses. The MyD88 and TRIF signalling pathways play important roles in activating innate and adaptive immunity in chickens infected with S. enteritidis. The objective of the present study was to characterize in vivo mRNA expressions, protein levels and methylation levels of genes in the above two pathways in both Tibetan chickens and DaHeng S03 chickens infected with S. enteritidis. MyD88-dependent and TRIF-dependent signalling pathway were activated by infection, and the MyD88 signalling pathway induced cytokines LITAF and IL-8 played important roles in fighting against the S. enteritidis infection in vivo. The TLR4 methylation might alter expression of genes involved in the MyD88 signalling pathway, and thus different breeds of chickens might show differences in susceptibility to the S. enteritidis. The increased expression of INF β was activated by S. enteritidis, but its expressions were different in levels of mRNA and protein in Tibetan chickens and DaHeng chickens, suggesting its functions on the resistance to S. enteritidis infection in chickens. This study contributes to the understanding of two pathways activated in response to S. enteritidis infection, and gives indications on the mechanisms underlying resistance of Tibetan chickens and DaHeng chickens to S. enteritidis.

Modulation the alternative splicing of GLA (IVS4+919G>A) in Fabry disease

While a base substitution in intron 4 of GLA (IVS4+919G>A) that causes aberrant alternative splicing resulting in Fabry disease has been reported, its molecular mechanism remains unclear. Here we reported that upon IVS4+919G>A transversion, H3K36me3 was enriched across the alternatively spliced region. PSIP1, an adapter of H3K36me3, together with Hsp70 and NONO were recruited and formed a complex with SF2/ASF and SRp20, which further promoted GLA splicing. Amiloride, a splicing regulator in cancer cells, could reverse aberrant histone modification patterns and disrupt the association of splicing complex with GLA. It could also reverse aberrant GLA splicing in a PP1-dependant manner. Our findings revealed the alternative splicing mechanism of GLA (IVS4+919G>A), and a potential treatment for this specific genetic type of Fabry disease by amiloride in the future.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

Intragenic CpG islands play important roles in bivalent chromatin assembly of developmental genes

CpG, 5'-C-phosphate-G-3', islands (CGIs) have long been known for their association with enhancers, silencers, and promoters, and for their epigenetic signatures. They are maintained in embryonic stem cells (ESCs) in a poised but inactive state via the formation of bivalent chromatin containing both active and repressive marks. CGIs also occur within coding sequences, where their functional role has remained obscure. Intragenic CGIs (iCGIs) are largely absent from housekeeping genes, but they are found in all genes associated with organ development and cell lineage control. In this paper, we investigated the epigenetic status of iCGIs and found that they too reside in bivalent chromatin in ESCs. Cell type-specific DNA methylation of iCGIs in differentiated cells was linked to the loss of both the H3K4me3 and H3K27me3 marks, and disruption of physical interaction with promoter regions, resulting in transcriptional activation of key regulators of differentiation such as PAXs, HOXs, and WNTs. The differential epigenetic modification of iCGIs appears to be mediated by cell type-specific transcription factors distinct from those bound by promoter, and these transcription factors may be involved in the hypermethylation of iCGIs upon cell differentiation. iCGIs thus play a key role in the cell type-specific regulation of transcription.

Methylation of glucocorticoid receptor gene promoter modulates morphine dependence and accompanied hypothalamus-pituitary-adrenal axis dysfunction

Previous studies demonstrated that dysfunction of the hypothalamus-pituitary-adrenal (HPA) axis played an important role in morphine dependence. Nonetheless, the molecular mechanism underlying morphine-induced HPA axis dysfunction and morphine dependence remains unclear. In the current study, 5'-aza-2'-deoxycytidine (5-aza), an inhibitor of DNA methyltransferases (DNMTs), was used to examine the effects of glucocorticoid receptor (GR) promoter 1₇ methylation on chronic morphine-induced HPA axis dysfunction and behavioral changes in rats and the underlying mechanism. Our results showed that chronic but not acute morphine downregulated the expression of nuclear GR protein and GR exon 1₇ variant mRNA, and upregulated the methylation of GR 1₇ exon promoter in the hippocampus of rats. Meanwhile, 5-aza per se had no effect on observed molecular and behavior change. In contrast, pretreatment of 5-aza into rat hippocampus reversed chronic morphine-induced hypermethylation of GR 1₇ promoter and decrease in GR expression. Moreover, pretreatment of 5-aza attenuated chronic morphine-enhanced HPA axis reactivity and the naloxone-precipitated somatic signs in morphine-dependent rats. Our results suggest that chronic morphine induced hypermethylation of GR 1₇ promoter, which then downregulated the expression of hippocampal GR, and was thus involved in chronic morphine-induced dysfunction of the HPA axis and the modulation of morphine dependence. Moreover, chronic morphine-induced hypermethylation of GR 1₇ promoter may be at least partially due to the increase in hippocampal DNMT 1 expression and its binding at GR 1₇ promoter in the rat hippocampus. © 2016 Wiley Periodicals, Inc.

Genomic aberrations frequently alter chromatin regulatory genes in chordoma

Chordoma is a rare primary bone neoplasm that is resistant to standard chemotherapies. Despite aggressive surgical management, local recurrence and metastasis is not uncommon. To identify the specific genetic aberrations that play key roles in chordoma pathogenesis, we utilized a genome-wide high-resolution SNP-array and next generation sequencing (NGS)-based molecular profiling platform to study 24 patient samples with typical histopathologic features of chordoma. Matching normal tissues were available for 16 samples. SNP-array analysis revealed nonrandom copy number losses across the genome, frequently involving 3, 9p, 1p, 14, 10, and 13. In contrast, copy number gain is uncommon in chordomas. Two minimum deleted regions were observed on 3p within a ∼8 Mb segment at 3p21.1-p21.31, which overlaps SETD2, BAP1 and PBRM1. The minimum deleted region on 9p was mapped to CDKN2A locus at 9p21.3, and homozygous deletion of CDKN2A was detected in 5/22 chordomas (∼23%). NGS-based molecular profiling demonstrated an extremely low level of mutation rate in chordomas, with an average of 0.5 mutations per sample for the 16 cases with matched normal. When the mutated genes were grouped based on molecular functions, many of the mutation events (∼40%) were found in chromatin regulatory genes. The combined copy number and mutation profiling revealed that SETD2 is the single gene affected most frequently in chordomas, either by deletion or by mutations. Our study demonstrated that chordoma belongs to the C-class (copy number changes) tumors whose oncogenic signature is non-random multiple copy number losses across the genome and genomic aberrations frequently alter chromatin regulatory genes. © 2016 Wiley Periodicals, Inc.

TINTIN, at the interface of chromatin, transcription elongation, and mRNA processing

Recent work including high-resolution genome-wide analysis uncovered a new trimeric complex involved in transcription elongation, both as an integral part of the NuA4 histone acetyltransferase and as an independent functional entity. The complex is conserved in eukaryotes and is named TINTIN, for Trimer Independent of NuA4 for transcription Interactions with Nucleosomes. This point of view covers the current knowledge regarding TINTIN's function in modulating chromatin structure and influencing transcription elongation in eukaryotes. It also points to several physical and functional links to co-transcriptional processes, including interactions with the mRNA splicing machinery and the nuclear exosome.

RNA-seq analysis of impact of PNN on gene expression and alternative splicing in corneal epithelial cells

PURPOSE

The specialized corneal epithelium requires differentiated properties, specific for its role at the anterior surface of the eye. Thus, tight maintenance of the differentiated qualities of the corneal epithelial is essential. Pinin (PNN) is an exon junction component (EJC) that has dramatic implications for corneal epithelial cell differentiation and may act as a stabilizer of the corneal epithelial cell phenotype. Our studies revealed that PNN is involved in transcriptional repression complexes and spliceosomal complexes, placing PNN at the fulcrum between chromatin and mRNA splicing. Transcriptome analysis of PNN-knockdown cells revealed clear and reproducible alterations in transcript profiles and splicing patterns of a subset of genes that would significantly impact the epithelial cell phenotype. We further investigated PNN's role in the regulation of gene expression and alternative splicing (AS) in a corneal epithelial context.

METHODS

Human corneal epithelial (HCET) cells that carry the doxycycline-inducible PNN-knockdown shRNA vector were used to perform RNA-seq to determine differential gene expression and differential AS events.

RESULTS

Multiple genes and AS events were identified as differentially expressed between PNN-knockdown and control cells. Genes upregulated by PNN knockdown included a large proportion of genes that are associated with enhanced cell migration and ECM remodeling processes, such as MMPs, ADAMs, HAS2, LAMA3, CXCRs, and UNC5C. Genes downregulated in response to PNN depletion included IGFBP5, FGD3, FGFR2, PAX6, RARG, and SOX10. AS events in PNN-knockdown cells compared to control cells were also more likely to be detected, and upregulated. In particular, 60% of exon-skipping events, detected in only one condition, were detected in PNN-knockdown cells and of the shared exon-skipping events, 92% of those differentially expressed were more frequent in the PNN knockdown.

CONCLUSIONS

These data suggest that lowering of PNN levels in epithelial cells results in dramatic transformation in the number and composition of splicing variants and that PNN plays a crucial role in the selection of which RNA isoforms differentiating cells produce. Many of the genes affected by PNN knockdown are known to affect the epithelial phenotype. This window into the complexity of RNA splicing in the corneal epithelium implies that PNN exerts broad influence over the regulation and maintenance of the epithelial cell phenotype.

Global analysis of physical and functional RNA targets of hnRNP L reveals distinct sequence and epigenetic features of repressed and enhanced exons

HnRNP L is a ubiquitous splicing-regulatory protein that is critical for the development and function of mammalian T cells. Previous work has identified a few targets of hnRNP L-dependent alternative splicing in T cells and has described transcriptome-wide association of hnRNP L with RNA. However, a comprehensive analysis of the impact of hnRNP L on mRNA expression remains lacking. Here we use next-generation sequencing to identify transcriptome changes upon depletion of hnRNP L in a model T-cell line. We demonstrate that hnRNP L primarily regulates cassette-type alternative splicing, with minimal impact of hnRNP L depletion on transcript abundance, intron retention, or other modes of alternative splicing. Strikingly, we find that binding of hnRNP L within or flanking an exon largely correlates with exon repression by hnRNP L. In contrast, exons that are enhanced by hnRNP L generally lack proximal hnRNP L binding. Notably, these hnRNP L-enhanced exons share sequence and context features that correlate with poor nucleosome positioning, suggesting that hnRNP may enhance inclusion of a subset of exons via a cotranscriptional or epigenetic mechanism. Our data demonstrate that hnRNP L controls inclusion of a broad spectrum of alternative cassette exons in T cells and suggest both direct RNA regulation as well as indirect mechanisms sensitive to the epigenetic landscape.

The nucleosome landscape of Plasmodium falciparum reveals chromatin architecture and dynamics of regulatory sequences

In eukaryotes, the chromatin architecture has a pivotal role in regulating all DNA-associated processes and it is central to the control of gene expression. For Plasmodium falciparum, a causative agent of human malaria, the nucleosome positioning profile of regulatory regions deserves particular attention because of their extreme AT-content. With the aid of a highly controlled MNase-seq procedure we reveal how positioning of nucleosomes provides a structural and regulatory framework to the transcriptional unit by demarcating landmark sites (transcription/translation start and end sites). In addition, our analysis provides strong indications for the function of positioned nucleosomes in splice site recognition. Transcription start sites (TSSs) are bordered by a small nucleosome-depleted region, but lack the stereotypic downstream nucleosome arrays, highlighting a key difference in chromatin organization compared to model organisms. Furthermore, we observe transcription-coupled eviction of nucleosomes on strong TSSs during intraerythrocytic development and demonstrate that nucleosome positioning and dynamics can be predictive for the functionality of regulatory DNA elements. Collectively, the strong nucleosome positioning over splice sites and surrounding putative transcription factor binding sites highlights the regulatory capacity of the nucleosome landscape in this deadly human pathogen.

Histone demethylase JARID1C inactivation triggers genomic instability in sporadic renal cancer

Mutations in genes encoding chromatin-remodeling proteins are often identified in a variety of cancers. For example, the histone demethylase JARID1C is frequently inactivated in patients with clear cell renal cell carcinoma (ccRCC); however, it is largely unknown how JARID1C dysfunction promotes cancer. Here, we determined that JARID1C binds broadly to chromatin domains characterized by the trimethylation of lysine 9 (H3K9me3), which is a histone mark enriched in heterochromatin. Moreover, we found that JARID1C localizes on heterochromatin, is required for heterochromatin replication, and forms a complex with established players of heterochromatin assembly, including SUV39H1 and HP1α, as well as with proteins not previously associated with heterochromatin assembly, such as the cullin 4 (CUL4) complex adaptor protein DDB1. Transcription on heterochromatin is tightly suppressed to safeguard the genome, and in ccRCC cells, JARID1C inactivation led to the unrestrained expression of heterochromatic noncoding RNAs (ncRNAs) that in turn triggered genomic instability. Moreover, ccRCC patients harboring JARID1C mutations exhibited aberrant ncRNA expression and increased genomic rearrangements compared with ccRCC patients with tumors endowed with other genetic lesions. Together, these data suggest that inactivation of JARID1C in renal cancer leads to heterochromatin disruption, genomic rearrangement, and aggressive ccRCCs. Moreover, our results shed light on a mechanism that underlies genomic instability in sporadic cancers.

Convergence of Acquired Mutations and Alternative Splicing of CD19 Enables Resistance to CART-19 Immunotherapy

The CD19 antigen, expressed on most B-cell acute lymphoblastic leukemias (B-ALL), can be targeted with chimeric antigen receptor-armed T cells (CART-19), but relapses with epitope loss occur in 10% to 20% of pediatric responders. We detected hemizygous deletions spanning the CD19 locus and de novo frameshift and missense mutations in exon 2 of CD19 in some relapse samples. However, we also discovered alternatively spliced CD19 mRNA species, including one lacking exon 2. Pull-down/siRNA experiments identified SRSF3 as a splicing factor involved in exon 2 retention, and its levels were lower in relapsed B-ALL. Using genome editing, we demonstrated that exon 2 skipping bypasses exon 2 mutations in B-ALL cells and allows expression of the N-terminally truncated CD19 variant, which fails to trigger killing by CART-19 but partly rescues defects associated with CD19 loss. Thus, this mechanism of resistance is based on a combination of deleterious mutations and ensuing selection for alternatively spliced RNA isoforms.

SIGNIFICANCE

CART-19 yield 70% response rates in patients with B-ALL, but also produce escape variants. We discovered that the underlying mechanism is the selection for preexisting alternatively spliced CD19 isoforms with the compromised CART-19 epitope. This mechanism suggests a possibility of targeting alternative CD19 ectodomains, which could improve survival of patients with B-cell neoplasms.

Transcriptome profile analysis reflects rat liver and kidney damage following chronic ultra-low dose Roundup exposure

BACKGROUND

Glyphosate-based herbicides (GBH) are the major pesticides used worldwide. Converging evidence suggests that GBH, such as Roundup, pose a particular health risk to liver and kidneys although low environmentally relevant doses have not been examined. To address this issue, a 2-year study in rats administering 0.1 ppb Roundup (50 ng/L glyphosate equivalent) via drinking water (giving a daily intake of 4 ng/kg bw/day of glyphosate) was conducted. A marked increased incidence of anatomorphological and blood/urine biochemical changes was indicative of liver and kidney structure and functional pathology. In order to confirm these findings we have conducted a transcriptome microarray analysis of the liver and kidneys from these same animals.

RESULTS

The expression of 4224 and 4447 transcript clusters (a group of probes corresponding to a known or putative gene) were found to be altered respectively in liver and kidney (p < 0.01, q < 0.08). Changes in gene expression varied from -3.5 to 3.7 fold in liver and from -4.3 to 5.3 in kidneys. Among the 1319 transcript clusters whose expression was altered in both tissues, ontological enrichment in 3 functional categories among 868 genes were found. First, genes involved in mRNA splicing and small nucleolar RNA were mostly upregulated, suggesting disruption of normal spliceosome activity. Electron microscopic analysis of hepatocytes confirmed nucleolar structural disruption. Second, genes controlling chromatin structure (especially histone-lysine N-methyltransferases) were mostly upregulated. Third, genes related to respiratory chain complex I and the tricarboxylic acid cycle were mostly downregulated. Pathway analysis suggests a modulation of the mTOR and phosphatidylinositol signalling pathways. Gene disturbances associated with the chronic administration of ultra-low dose Roundup reflect a liver and kidney lipotoxic condition and increased cellular growth that may be linked with regeneration in response to toxic effects causing damage to tissues. Observed alterations in gene expression were consistent with fibrosis, necrosis, phospholipidosis, mitochondrial membrane dysfunction and ischemia, which correlate with and thus confirm observations of pathology made at an anatomical, histological and biochemical level.

CONCLUSION

Our results suggest that chronic exposure to a GBH in an established laboratory animal toxicity model system at an ultra-low, environmental dose can result in liver and kidney damage with potential significant health implications for animal and human populations.

Regulation of RNA metabolism is important for in vitro dedifferentiation of plant cells

The characteristic high regenerative ability of plants has been exploited to develop in vitro plant regeneration techniques, which are usually initiated by an in vitro dedifferentiation step induced by artificial phytohormone treatment. Recent advances in plant molecular biological and genetic technologies have revealed the importance of the regulation of RNA metabolism, including the control of rRNA biosynthesis, pre-mRNA splicing, and miRNA-based RNA decay, in successful in vitro dedifferentiation. This review provides a brief overview of current knowledge of the roles of RNA metabolism in the dedifferentiation of plant cells in vitro. In addition, the possibility that plant-specific aspects of RNA metabolism regulation are linked closely to their high regenerative ability is discussed.

Acetylation of intragenic histones on HPV16 correlates with enhanced HPV16 gene expression

We report that many histone modifications are unevenly distributed over the HPV16 genome in cervical cancer cells as well as in HPV16-immortalized keratinocytes. For example, H3K36me3 and H3K9Ac that are common in highly expressed cellular genes and over exons, were more common in the early than in the late region of the HPV16 genome. In contrast, H3K9me3, H4K20me3, H2BK5me1 and H4K16Ac were more frequent in the HPV16 late region. Furthermore, a region encompassing the HPV16 early polyadenylation signal pAE displayed high levels of histone H3 acetylation. Histone deacetylase (HDAC) inhibitors caused a 2- to 8-fold induction of HPV16 early and late mRNAs in cervical cancer cells and in immortalized keratinocytes, while at the same time increasing the levels of acetylated histones in the cells and on the HPV16 genome specifically. We concluded that increased histone acetylation on the HPV16 genome correlates with increased HPV16 gene expression.

Co-transcriptional mRNP formation is coordinated within a molecular mRNP packaging station in S. cerevisiae

In eukaryotes, the messenger RNA (mRNA), the blueprint of a protein‐coding gene, is processed and packaged into a messenger ribonucleoprotein particle (mRNP) by mRNA‐binding proteins in the nucleus. The steps of mRNP formation – transcription, processing, packaging, and the orchestrated release of the export‐competent mRNP from the site of transcription for nuclear mRNA export – are tightly coupled to ensure a highly efficient and regulated process. The importance of highly accurate nuclear mRNP formation is illustrated by the fact that mutations in components of this pathway lead to cellular inviability or to severe diseases in metazoans. We hypothesize that efficient mRNP formation is realized by a molecular mRNP packaging station, which is built by several recruitment platforms and coordinates the individual steps of mRNP formation.

DNMT3A in haematological malignancies

DNA methylation patterns are disrupted in various malignancies, suggesting a role in the development of cancer, but genetic aberrations directly linking the DNA methylation machinery to malignancies were rarely observed, so this association remained largely correlative. Recently, however, mutations in the gene encoding DNA methyltransferase 3A (DNMT3A) were reported in patients with acute myeloid leukaemia (AML), and subsequently in patients with various other haematological malignancies, pointing to DNMT3A as a critically important new tumour suppressor. Here, we review the clinical findings related to DNMT3A, tie these data to insights from basic science studies conducted over the past 20 years and present a roadmap for future research that should advance the agenda for new therapeutic strategies.

Unexpected selection to retain high GC content and splicing enhancers within exons of multiexonic lncRNA loci

If sequencing was possible only for genomes, and not for RNAs or proteins, then functional protein-coding exons would be recognizable by their unusual patterns of nucleotide composition, specifically a high GC content across the body of exons, and an unusual nucleotide content near their edges. RNAs and proteins can, of course, be sequenced but the extent of functionality of intergenic long noncoding RNAs (lncRNAs) remains under question owing to their low nucleotide conservation. Inspired by the nucleotide composition patterns of protein-coding exons, we sought evidence for functionality across lncRNA loci from diverse species. We found that such patterns across multiexonic lncRNA loci mirror those of proteincoding genes, although to a lesser degree: Specifically, compared with introns, lncRNA exons are GC rich. Additionally we report evidence for the action of purifying selection to preserve exonic splicing enhancers within human multiexonic lncRNAs and nucleotide composition in fruit fly lncRNAs. Our findings provide evidence for selection for more efficient rates of transcription and splicing within lncRNA loci. Despite only a minor proportion of their RNA bases being constrained, multiexonic intergenic lncRNAs appear to require accurate splicing of their exons to transact their function.