Nucleosomes are well positioned in exons and carry characteristic histone modifications

TET1 displays catalytic and non-catalytic functions in the adult mouse cortex

TET1/2/3 dioxygenases iteratively demethylate 5-methylcytosine, beginning with the formation of 5-hydroxymethylcytosine (5hmC). The post-mitotic brain maintains higher levels of 5hmC than most peripheral tissues, and TET1 ablation studies have underscored the critical role of TET1 in brain physiology. However, deletion of Tet1 precludes the disentangling of the catalytic and non-catalytic functions of TET1. Here, we dissect these functions of TET1 by comparing adult cortex of Tet1 wildtype (Tet1 WT), a novel Tet1 catalytically dead mutant (Tet1 HxD), and Tet1 knockout (Tet1 KO) mice. Using DNA methylation array, we uncover that Tet1 HxD and KO mutations perturb the methylation status of distinct subsets of CpG sites. Gene ontology (GO) analysis on specific differential 5hmC regions indicates that TET1's catalytic activity is linked to neuronal-specific functions. RNA-Seq further shows that Tet1 mutations predominantly impact the genes that are associated with alternative splicing. Lastly, we performed High-performance Liquid Chromatography Mass-Spectrometry lipidomics on WT and mutant cortices and uncover accumulation of lysophospholipids lysophosphatidylethanolamine and lysophosphatidylcholine in Tet1 HxD cortex. In summary, we show that Tet1 HxD does not completely phenocopy Tet1 KO, providing evidence that TET1 modulates distinct cortical functions through its catalytic and non-catalytic roles.

The Functional Relationship Between RNA Splicing and the Chromatin Landscape

Chromatin is a highly regulated and dynamic structure that has been shown to play an essential role in transcriptional and co-transcriptional regulation. In the context of RNA splicing, early evidence suggested a loose connection between the chromatin landscape and splicing. More recently, it has been shown that splicing occurs in a co-transcriptional manner, meaning that the splicing process occurs in the context of chromatin. Experimental and computational evidence have also shown that chromatin dynamics can influence the splicing process and vice versa. However, much of this evidence provides mainly correlative relationships between chromatin and splicing with just a few concrete examples providing defined molecular mechanisms by which these two processes are functionally related. Nevertheless, it is clear that chromatin and RNA splicing are tightly interconnected to one another. In this review, we highlight the current state of knowledge of the relationship between chromatin and splicing.

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

The Long Non-Coding RNA ANRIL in Cancers

ANRIL (Antisense Noncoding RNA in the INK4 Locus), a long non-coding RNA encoded in the human chromosome 9p21 region, is a critical factor for regulating gene expression by interacting with multiple proteins and miRNAs. It has been found to play important roles in various cellular processes, including cell cycle control and proliferation. Dysregulation of ANRIL has been associated with several diseases like cancers and cardiovascular diseases, for instance. Understanding the oncogenic role of ANRIL and its potential as a diagnostic and prognostic biomarker in cancer is crucial. This review provides insights into the regulatory mechanisms and oncogenic significance of the 9p21 locus and ANRIL in cancer.

Epigenetic silencing of OR and TAS2R genes expression in human orbitofrontal cortex at early stages of sporadic Alzheimer's disease

Modulation of brain olfactory (OR) and taste receptor (TASR) expression was recently reported in neurological diseases. However, there is still limited evidence of these genes' expression in the human brain and the transcriptional regulation mechanisms involved remain elusive. We explored the possible expression and regulation of selected OR and TASR in the human orbitofrontal cortex (OFC) of sporadic Alzheimer's disease (AD) and non-demented control specimens using quantitative real-time RT-PCR and ELISA. Global H3K9me3 amounts were measured on OFC total histone extracts, and H3K9me3 binding at each chemoreceptor locus was examined through native chromatin immunoprecipitation. To investigate the potential interactome of the repressive histone mark H3K9me3 in OFC specimens, native nuclear complex co-immunoprecipitation (Co-IP) was combined with reverse phase-liquid chromatography coupled to mass spectrometry analysis. Interaction between H3K9me3 and MeCP2 was validated by reciprocal Co-IP, and global MeCP2 levels were quantitated. We found that OR and TAS2R genes are expressed and markedly downregulated in OFC at early stages of sporadic AD, preceding the progressive reduction in their protein levels and the appearance of AD-associated neuropathology. The expression pattern did not follow disease progression suggesting transcriptional regulation through epigenetic mechanisms. We discovered an increase of OFC global H3K9me3 levels and a substantial enrichment of this repressive signature at ORs and TAS2Rs proximal promoter at early stages of AD, ultimately lost at advanced stages. We revealed the interaction between H3K9me3 and MeCP2 at early stages and found that MeCP2 protein is increased in sporadic AD. Findings suggest MeCP2 might be implicated in OR and TAS2R transcriptional regulation through interaction with H3K9me3, and as an early event, it may uncover a novel etiopathogenetic mechanism of sporadic AD.

The TUDOR domain of SMN is an H3K79me1 histone mark reader

Spinal muscular atrophy is the leading genetic cause of infant mortality and results from depleted levels of functional survival of motor neuron (SMN) protein by either deletion or mutation of the SMN1 gene. SMN is characterized by a central TUDOR domain, which mediates the association of SMN with arginine methylated (Rme) partners, such as coilin, fibrillarin, and RNA pol II (RNA polymerase II). Herein, we biochemically demonstrate that SMN also associates with histone H3 monomethylated on lysine 79 (H3K79me1), defining SMN as not only the first protein known to associate with the H3K79me1 histone modification but also the first histone mark reader to recognize both methylated arginine and lysine residues. Mutational analyzes provide evidence that SMNTUDOR associates with H3 via an aromatic cage. Importantly, most SMNTUDOR mutants found in spinal muscular atrophy patients fail to associate with H3K79me1.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Regulation of pre-mRNA splicing: roles in physiology and disease, and therapeutic prospects

The removal of introns from mRNA precursors and its regulation by alternative splicing are key for eukaryotic gene expression and cellular function, as evidenced by the numerous pathologies induced or modified by splicing alterations. Major recent advances have been made in understanding the structures and functions of the splicing machinery, in the description and classification of physiological and pathological isoforms and in the development of the first therapies for genetic diseases based on modulation of splicing. Here, we review this progress and discuss important remaining challenges, including predicting splice sites from genomic sequences, understanding the variety of molecular mechanisms and logic of splicing regulation, and harnessing this knowledge for probing gene function and disease aetiology and for the design of novel therapeutic approaches.

The sound of silence: mechanisms and implications of HUSH complex function

The vertebrate genome is under constant threat of invasion by genetic parasites. Whether the host can immediately recognize and respond to invading elements has been unclear. The discovery of the human silencing hub (HUSH) complex, and the finding that it provides immediate protection from genome invasion by silencing products of reverse transcription, have important implications for mammalian genome evolution. In this review, we summarize recent insights into HUSH function and describe how cellular introns provide a novel means of self-nonself discrimination, allowing HUSH to recognize and transcriptionally repress a broad range of intronless genetic elements. We discuss how HUSH contributes to genome evolution, and highlight studies reporting the critical role of HUSH in development and implicating HUSH in the control of immune signaling and cancer progression.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

SET domain containing 2 (SETD2) influences metabolism and alternative splicing during myogenesis

Epigenetic regulatory mechanisms are increasingly recognized as crucial determinants of cellular specification and differentiation. During muscle cell differentiation (myogenesis), extensive remodelling of histone acetylation and methylation occurs. Several of these histone modifications aid in the expression of muscle-specific genes and the silencing of genes that block lineage commitment. Therefore, the identification of new epigenetic regulatory mechanisms is of high interest. Still, the functional relevance of numerous histone modifications during myogenesis remain completely uncertain. In this study, we focus on the function of H3K36me3 and its epigenetic writer, SET domain containing 2 (SETD2), in the context of muscle cell differentiation. We first observed that SETD2 expression increases during myogenesis. Targeted depletion of SETD2 in undifferentiated (myoblasts) and differentiated (myotubes) muscle cells reduced H3K36me3 levels and induced profound changes in gene expression and slight alterations in alternative splicing, as determined by deep RNA-sequencing analysis. Enzymes that function in metabolic pathways were upregulated in response to SETD2 depletion. Furthermore, we demonstrated that upregulation of several glycolytic enzymes was associated with an increase in intracellular pyruvate levels in SETD2-depleted cells, indicating a novel role for SETD2 in metabolic programming during myogenesis. Together, our results provide new insight into the signalling pathways controlled by chromatin-modifying enzymes and their associated histone modifications during muscle cell differentiation.

PARP1′s Involvement in RNA Polymerase II Elongation: Pausing and Releasing Regulation through the Integrator and Super Elongation Complex

RNA polymerase elongation along the gene body is tightly regulated to ensure proper transcription and alternative splicing events. Understanding the mechanism and factors critical in regulating the rate of RNA polymerase II elongation and processivity is clearly important. Recently we showed that PARP1, a well-known DNA repair protein, when bound to chromatin, regulates RNA polymerase II elongation. However, the mechanism by which it does so is not known. In the current study, we aimed to tease out how PARP1 regulates RNAPII elongation. We show, both in vivo and in vitro, that PARP1 binds directly to the Integrator subunit 3 (IntS3), a member of the elongation Integrator complex. The association between the two proteins is mediated via the C-terminal domain of PARP1 to the C-terminal domain of IntS3. Interestingly, the occupancy of IntS3 along two PARP1 target genes mimicked that of PARP1, suggesting a role in its recruitment/assembly of elongation factors. Indeed, the knockdown of PARP1 resulted in differential chromatin association and gene occupancy of IntS3 and other key elongation factors. Most of these PARP1-mediated effects were due to the physical presence of PARP1 rather than its PARylation activity. These studies argue that PARP1 controls the progressive RNAPII elongation complexes. In summary, we present a platform to begin to decipher PARP1′s role in recruiting/scaffolding elongation factors along the gene body regions during RNA polymerase II elongation and gene regulation.

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.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

Modulation of RNA Splicing by Oligonucleotides: Mechanisms of Action and Therapeutic Implications

Dysregulation of RNA splicing causes many diseases and disorders. Several therapeutic approaches have been developed to correct aberrant alternative splicing events for the treatment of cancers and hereditary diseases, including gene therapy and redirecting splicing, using small molecules or splice switching oligonucleotides (SSO). Significant advances in the chemistry and pharmacology of nucleic acid have led to the development of clinically approved SSO drugs for the treatment of spinal muscular dystrophy and Duchenne muscular dystrophy (DMD). In this review, we discuss the mechanisms of SSO action with emphasis on "less common" approaches to modulate alternative splicing, including bipartite and bifunctional SSO, oligonucleotide decoys for splice factors and SSO-mediated mRNA degradation via AS-NMD and NGD pathways. We briefly discuss the current progress and future perspectives of SSO therapy for rare and ultrarare diseases.

Antibiofilm Activity of Allicin and Quercetin in Treating Biofilm-Associated Orthopaedics Infection

Biofilms formed by bacteria are the group of sessile microbial cells that remain encompassed by self-secreted polymeric substances and have resulted in great health-care concern. The extracellular polymeric substances (EPS) prevent the penetration of antibiotics and other drugs, thereby resulting in the development of multi-drug resistance or antibiotic resistance. The biofilm-associated prosthetics being places at the joins of bone injury are the common sites for the development of biofilm-associated infection. This often spreads and results in the development of orthopaedic infections. Most of the infections are associated with musculoskeletal system and originate from non-living surfaces. The biofilm prevents the penetration of drugs, thereby resulting in the development of antibiotic resistance or multi-drug resistance. The minimum inhibitory concentration (MIC) for allicin and quercetin was found to be 80 µg/mL for quercetin and 100 µg/mL for amoxicillin against the sessile communities of Pseudomonas aeruginosa associated with the orthopaedic infection. The role of quercetin and allicin in reduction of protein, carbohydrate and eDNA content of the exopolysaccharides (EPS) was tested. The anti-quorum sensing activity of quercetin and allicin was confirmed both by biochemical and by photomicrographic studies. The antibiofilm and antimicrobial activities of quercetin and allicin were determined both by in vitro and in silico studies on P. aeruginosa bacterial strain from biofilm-associated orthopaedic infection.

Mapping cis-regulatory elements in the midgestation mouse placenta

The placenta is a temporary organ that provides the developing fetus with nutrients, oxygen, and protection in utero. Defects in its development, which may be caused by misregulated gene expression, can lead to devastating outcomes for the mother and fetus. In mouse, placental defects during midgestation commonly lead to embryonic lethality. However, the regulatory mechanisms controlling expression of genes during this period have not been thoroughly investigated. Therefore, we generated and analyzed ChIP-seq data for multiple histone modifications known to mark cis-regulatory regions. We annotated active and poised promoters and enhancers, as well as regions generally associated with repressed gene expression. We found that poised promoters were associated with neuronal development genes, while active promoters were largely associated with housekeeping genes. Active and poised enhancers were associated with placental development genes, though only active enhancers were associated with genes that have placenta-specific expression. Motif analysis within active enhancers identified a large network of transcription factors, including those that have not been previously studied in the placenta and are candidates for future studies. The data generated and genomic regions annotated provide researchers with a foundation for future studies, aimed at understanding how specific genes in the midgestation mouse placenta are regulated.

High frequency of intron retention and clustered H3K4me3-marked nucleosomes in short first introns of human long non-coding RNAs

BACKGROUND

It is established that protein-coding exons are preferentially localized in nucleosomes. To examine whether the same is true for non-coding exons, we analysed nucleosome occupancy in and adjacent to internal exons in genes encoding long non-coding RNAs (lncRNAs) in human CD4+ T cells and K562 cells.

RESULTS

We confirmed that internal exons in lncRNAs are preferentially associated with nucleosomes, but also observed an elevated signal from H3K4me3-marked nucleosomes in the sequences upstream of these exons. Examination of 200 genomic lncRNA loci chosen at random across all chromosomes showed that high-density regions of H3K4me3-marked nucleosomes, which we term 'slabs', are associated with genomic regions exhibiting intron retention. These retained introns occur in over 50% of lncRNAs examined and are mostly first introns with an average length of just 354 bp, compared to the average length of all human introns of 6355 and 7987 bp in mRNAs and lncRNAs, respectively. Removal of short introns from the dataset abrogated the high upstream H3K4me3 signal, confirming that the association of slabs and short lncRNA introns with intron retention holds genome-wide. The high upstream H3K4me3 signal is also associated with alternatively spliced exons, known to be prominent in lncRNAs. This phenomenon was not observed with mRNAs.

CONCLUSIONS

There is widespread intron retention and clustered H3K4me3-marked nucleosomes in short first introns of human long non-coding RNAs, which raises intriguing questions about the relationship of IR to lncRNA function and chromatin organization.

NucHMM: a method for quantitative modeling of nucleosome organization identifying functional nucleosome states distinctly associated with splicing potentiality

We develop a novel computational method, NucHMM, to identify functional nucleosome states associated with cell type-specific combinatorial histone marks and nucleosome organization features such as phasing, spacing and positioning. We test it on publicly available MNase-seq and ChIP-seq data in MCF7, H1, and IMR90 cells and identify 11 distinct functional nucleosome states. We demonstrate these nucleosome states are distinctly associated with the splicing potentiality of skipping exons. This advances our understanding of the chromatin function at the nucleosome level and offers insights into the interplay between nucleosome organization and splicing processes.

CTCF as a regulator of alternative splicing: new tricks for an old player

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the 'master weaver of the genome'. Given its multiple connections, further characterization of CTCF's emerging role in splicing regulation might extend its functional repertoire towards a 'conductor of the splicing orchestra'.

RNA polymerase II speed: a key player in controlling and adapting transcriptome composition

RNA polymerase II (RNA Pol II) speed or elongation rate, i.e., the number of nucleotides synthesized per unit of time, is a major determinant of transcriptome composition. It controls co-transcriptional processes such as splicing, polyadenylation, and transcription termination, thus regulating the production of alternative splice variants, circular RNAs, alternatively polyadenylated transcripts, or read-through transcripts. RNA Pol II speed itself is regulated in response to intra- and extra-cellular stimuli and can in turn affect the transcriptome composition in response to these stimuli. Evidence points to a potentially important role of transcriptome composition modification through RNA Pol II speed regulation for adaptation of cells to a changing environment, thus pointing to a function of RNA Pol II speed regulation in cellular physiology. Analyzing RNA Pol II speed dynamics may therefore be central to fully understand the regulation of physiological processes, such as the development of multicellular organisms. Recent findings also raise the possibility that RNA Pol II speed deregulation can be detrimental and participate in disease progression. Here, we review initial and current approaches to measure RNA Pol II speed, as well as providing an overview of the factors controlling speed and the co-transcriptional processes which are affected. Finally, we discuss the role of RNA Pol II speed regulation in cell physiology.

Tissue-specific cell-free DNA degradation quantifies circulating tumor DNA burden

Profiling of circulating tumor DNA (ctDNA) may offer a non-invasive approach to monitor disease progression. Here, we develop a quantitative method, exploiting local tissue-specific cell-free DNA (cfDNA) degradation patterns, that accurately estimates ctDNA burden independent of genomic aberrations. Nucleosome-dependent cfDNA degradation at promoters and first exon-intron junctions is strongly associated with differential transcriptional activity in tumors and blood. A quantitative model, based on just 6 regulatory regions, could accurately predict ctDNA levels in colorectal cancer patients. Strikingly, a model restricted to blood-specific regulatory regions could predict ctDNA levels across both colorectal and breast cancer patients. Using compact targeted sequencing (<25 kb) of predictive regions, we demonstrate how the approach could enable quantitative low-cost tracking of ctDNA dynamics and disease progression.

Circulating tumour DNA (ctDNA) represents a non-invasive option to monitor cancer progression. Here, the authors perform deep sequencing of plasma cell-free DNA, and find that nucleosome-dependent cfDNA degradation at 6 specific regulatory regions is predictive of ctDNA burden.

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.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

BRG1 knockdown inhibits proliferation through multiple cellular pathways in prostate cancer

Background

BRG1 (encoded by SMARCA4) is a catalytic component of the SWI/SNF chromatin remodelling complex, with key roles in modulating DNA accessibility. Dysregulation of BRG1 is observed, but functionally uncharacterised, in a wide range of malignancies. We have probed the functions of BRG1 on a background of prostate cancer to investigate how BRG1 controls gene expression programmes and cancer cell behaviour.

Results

Our investigation of SMARCA4 revealed that BRG1 is over-expressed in the majority of the 486 tumours from The Cancer Genome Atlas prostate cohort, as well as in a complementary panel of 21 prostate cell lines. Next, we utilised a temporal model of BRG1 depletion to investigate the molecular effects on global transcription programmes. Depleting BRG1 had no impact on alternative splicing and conferred only modest effect on global expression. However, of the transcriptional changes that occurred, most manifested as down-regulated expression. Deeper examination found the common thread linking down-regulated genes was involvement in proliferation, including several known to increase prostate cancer proliferation (KLK2, PCAT1 and VAV3). Interestingly, the promoters of genes driving proliferation were bound by BRG1 as well as the transcription factors, AR and FOXA1. We also noted that BRG1 depletion repressed genes involved in cell cycle progression and DNA replication, but intriguingly, these pathways operated independently of AR and FOXA1. In agreement with transcriptional changes, depleting BRG1 conferred G1 arrest.

Conclusions

Our data have revealed that BRG1 promotes cell cycle progression and DNA replication, consistent with the increased cell proliferation associated with oncogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01023-7.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

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.

Splicing to Keep Cycling: The Importance of Pre-mRNA Splicing during the Cell Cycle

Pre-mRNA splicing is a fundamental process in mammalian gene expression, and alternative splicing plays an extensive role in generating protein diversity. Because the majority of genes undergo pre-mRNA splicing, most cellular processes depend on proper spliceosome function. We focus on the cell cycle and describe its dependence on pre-mRNA splicing and accurate alternative splicing. We outline the key cell-cycle factors and their known alternative splicing isoforms. We discuss different levels of pre-mRNA splicing regulation such as post-translational modifications and changes in the expression of splicing factors. We describe the effect of chromatin dynamics on pre-mRNA splicing during the cell cycle. In addition, we focus on spliceosome component SF3B1, which is mutated in many types of cancer, and describe the link between SF3B1 and its inhibitors and the cell cycle.

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.

Transcription and splicing: A two-way street

RNA synthesis by RNA polymerase II and RNA processing are closely coupled during the transcription cycle of protein-coding genes. This coupling affords opportunities for quality control and regulation of gene expression and the effects can go in both directions. For example, polymerase speed can affect splice site selection and splicing can increase transcription and affect the chromatin landscape. Here we review the many ways that transcription and splicing influence one another, including how splicing "talks back" to transcription. We will also place the connections between transcription and splicing in the context of other RNA processing events that define the exons that will make up the final mRNA. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

H3K36 Methylation and the Chromodomain Protein Eaf3 Are Required for Proper Cotranscriptional Spliceosome Assembly

In the eukaryotic cell, spliceosomes assemble onto pre-mRNA cotranscriptionally. Spliceosome assembly takes place in the context of the chromatin environment, suggesting that the state of the chromatin may affect splicing. The molecular details and mechanisms through which chromatin affects splicing, however, are still unclear. Here, we show a role for the histone methyltransferase Set2 and its histone modification, H3K36 methylation, in pre-mRNA splicing through high-throughput sequencing. Moreover, the effect of H3K36 methylation on pre-mRNA splicing is mediated through the chromodomain protein Eaf3. We find that Eaf3 is recruited to intron-containing genes and that Eaf3 interacts with the splicing factor Prp45. Eaf3 acts with Prp45 and Prp19 after formation of the precatalytic B complex around the time of splicing activation, thus revealing the step in splicing that is regulated by H3K36 methylation. These studies support a model whereby H3K36 facilitates recruitment of an "adapter protein" to support efficient, constitutive splicing.

Epigenome-based splicing prediction using a recurrent neural network

Alternative RNA splicing provides an important means to expand metazoan transcriptome diversity. Contrary to what was accepted previously, splicing is now thought to predominantly take place during transcription. Motivated by emerging data showing the physical proximity of the spliceosome to Pol II, we surveyed the effect of epigenetic context on co-transcriptional splicing. In particular, we observed that splicing factors were not necessarily enriched at exon junctions and that most epigenetic signatures had a distinctly asymmetric profile around known splice sites. Given this, we tried to build an interpretable model that mimics the physical layout of splicing regulation where the chromatin context progressively changes as the Pol II moves along the guide DNA. We used a recurrent-neural-network architecture to predict the inclusion of a spliced exon based on adjacent epigenetic signals, and we showed that distinct spatio-temporal features of these signals were key determinants of model outcome, in addition to the actual nucleotide sequence of the guide DNA strand. After the model had been trained and tested (with >80% precision-recall curve metric), we explored the derived weights of the latent factors, finding they highlight the importance of the asymmetric time-direction of chromatin context during transcription.

The DNA methylation landscape in cancer

As one of the most abundant and well-studied epigenetic modifications, DNA methylation plays an essential role in normal development and cellular biology. Global alterations to the DNA methylation landscape contribute to alterations in the transcriptome and deregulation of cellular pathways. Indeed, improved methods to study DNA methylation patterning and dynamics at base pair resolution and across individual DNA molecules on a genome-wide scale has highlighted the scope of change to the DNA methylation landscape in disease states, particularly during tumorigenesis. More recently has been the development of DNA hydroxymethylation profiling techniques, which allows differentiation between 5mC and 5hmC profiles and provides further insights into DNA methylation dynamics and remodeling in tumorigenesis. In this review, we describe the distribution of DNA methylation and DNA hydroxymethylation in different genomic contexts, first in normal cells, and how this is altered in cancer. Finally, we discuss DNA methylation profiling technologies and the most recent advances in single-cell methods, bisulfite-free approaches and ultra-long read sequencing techniques.

Native elongation transcript sequencing reveals temperature dependent dynamics of nascent RNAPII transcription in Arabidopsis

Temperature profoundly affects the kinetics of biochemical reactions, yet how large molecular complexes such as the transcription machinery accommodate changing temperatures to maintain cellular function is poorly understood. Here, we developed plant native elongating transcripts sequencing (plaNET-seq) to profile genome-wide nascent RNA polymerase II (RNAPII) transcription during the cold-response of Arabidopsis thaliana with single-nucleotide resolution. Combined with temporal resolution, these data revealed transient genome-wide reprogramming of nascent RNAPII transcription during cold, including characteristics of RNAPII elongation and thousands of non-coding transcripts connected to gene expression. Our results suggest a role for promoter-proximal RNAPII stalling in predisposing genes for transcriptional activation during plant-environment interactions. At gene 3'-ends, cold initially facilitated transcriptional termination by limiting the distance of read-through transcription. Within gene bodies, cold reduced the kinetics of co-transcriptional splicing leading to increased intragenic stalling. Our data resolved multiple distinct mechanisms by which temperature transiently altered the dynamics of nascent RNAPII transcription and associated RNA processing, illustrating potential biotechnological solutions and future focus areas to promote food security in the context of a changing climate.

Epigenetic biomarkers and preterm birth

Preterm birth (PTB) is a major public health challenge, and novel, sensitive approaches to predict PTB are still evolving. Epigenomic markers are being explored as biomarkers of PTB because of their molecular stability compared to gene expression. This approach is also relatively new compared to gene-based diagnostics, which relies on mutations or single nucleotide polymorphisms. The fundamental principle of epigenome diagnostics is that epigenetic reprogramming in the target tissue (e.g. placental tissue) might be captured by more accessible surrogate tissue (e.g. blood) using biochemical epigenome assays on circulating DNA that incorporate methylation, histone modifications, nucleosome positioning, and/or chromatin accessibility. Epigenomic-based biomarkers may hold great potential for early identification of the majority of PTBs that are not associated with genetic variants or mutations. In this review, we discuss recent advances made in the development of epigenome assays focusing on its potential exploration for association and prediction of PTB. We also summarize population-level cohort studies conducted in the USA and globally that provide opportunities for genetic and epigenetic marker development for PTB. In addition, we summarize publicly available epigenome resources and published PTB studies. We particularly focus on ongoing genome-wide DNA methylation and epigenome-wide association studies. Finally, we review the limitations of current research, the importance of establishing a comprehensive biobank, and possible directions for future studies in identifying effective epigenome biomarkers to enhance health outcomes for pregnant women at risk of PTB and their infants.

Nascent RNA analyses: tracking transcription and its regulation

The programmes that direct an organism's development and maintenance are encoded in its genome. Decoding of this information begins with regulated transcription of genomic DNA into RNA. Although transcription and its control can be tracked indirectly by measuring stable RNAs, it is only by directly measuring nascent RNAs that the immediate regulatory changes in response to developmental, environmental, disease and metabolic signals are revealed. Multiple complementary methods have been developed to quantitatively track nascent transcription genome-wide at nucleotide resolution, all of which have contributed novel insights into the mechanisms of gene regulation and transcription-coupled RNA processing. Here we critically evaluate the array of strategies used for investigating nascent transcription and discuss the recent conceptual advances they have provided.

Histone H1.5 binds over splice sites in chromatin and regulates alternative splicing

Chromatin organization and epigenetic markers influence splicing, though the magnitudes of these effects and the mechanisms are largely unknown. Here, we demonstrate that linker histone H1.5 influences mRNA splicing. We observed that linker histone H1.5 binds DNA over splice sites of short exons in human lung fibroblasts (IMR90 cells). We found that association of H1.5 with these splice sites correlated with the level of inclusion of alternatively spliced exons. Exons marked by H1.5 had more RNA polymerase II (RNAP II) stalling near the 3' splice site than did exons not associated with H1.5. In cells depleted of H1.5, we showed that the inclusion of five exons evaluated decreased and that RNAP II levels over these exons were also reduced. Our findings indicate that H1.5 is involved in regulation of splice site selection and alternative splicing, a function not previously demonstrated for linker histones.

LMX1B mRNA expression and its gene body CpG methylation are valuable prognostic biomarkers for laryngeal squamous cell carcinoma

This study aimed to explore the prognostic value of LMX1B mRNA expression and the methylation of its CpG sites in patients with laryngeal squamous cell carcinoma (LSCC). An in-silicon analysis was performed using data from the cancer genome atlas (TCGA)-Head and Neck Squamous Carcinoma (HNSC). After screening, 112 LSCC and 10 adjacent normal tissues were identified as eligible samples for analysis. Results showed that LMX1B expression was significantly upregulated in the cancer tissues (p < 0.01) and was an independent prognostic indicator in terms of OS (HR: 1.233, 95%CI: 1.082-1.405, p = 0.002) and RFS (HR: 1.200, 95%CI: 1.002-1.438, p = 0.048). By examining the methylation profile of 55 CpG sites in LMX1B locus, we found that the promoter methylation status was irrelevant to LMX1B expression. In comparison, LMX1B expression was generally positively correlated with gene body methylation. Among the gene body CpG sites, cg13600622 methylation showed a better predictive value than LMX1B expression in terms of OS (HR: 12.363, 95%CI: 1.076-142.033, p = 0.043), while cg14204784 methylation was a better marker of shorter RFS (HR: 12.363, 95%CI: 1.076-142.033, p = 0.043). Among the known downstream genes of LMX1B, only NR4A2 expression showed a moderately negative correlation (Pearson's r = -0.54) with it in LSCC tissues. However, this correlation was inconsistent with previous publications those reported a positive correlation between them. Based on these findings, we infer that upregulated LMX1B mRNA expression had an independent prognostic value in LSCC patients. Increased gene body methylation might be an important mechanism of its upregulation. Among the gene body CpG sites, cg13600622 and cg14204784 methylation level might be better prognostic markers than LMX1B mRNA expression in terms of OS and RFS respectively.

Regulation of alternative splicing by p300-mediated acetylation of splicing factors

Splicing of precursor mRNA (pre-mRNA) is an important regulatory step in gene expression. Recent evidence points to a regulatory role of chromatin-related proteins in alternative splicing regulation. Using an unbiased approach, we have identified the acetyltransferase p300 as a key chromatin-related regulator of alternative splicing. p300 promotes genome-wide exon inclusion in both a transcription-dependent and -independent manner. Using CD44 as a paradigm, we found that p300 regulates alternative splicing by modulating the binding of splicing factors to pre-mRNA. Using a tethering strategy, we found that binding of p300 to the CD44 promoter region promotes CD44v exon inclusion independently of RNAPII transcriptional elongation rate. Promoter-bound p300 regulates alternative splicing by acetylating splicing factors, leading to exclusion of hnRNP M from CD44 pre-mRNA and activation of Sam68. p300-mediated CD44 alternative splicing reduces cell motility and promotes epithelial features. Our findings reveal a chromatin-related mechanism of alternative splicing regulation and demonstrate its impact on cellular function.

Entanglement of Genetics and Epigenetics in Parkinson's Disease

Parkinson disease (PD) is a common neurodegenerative disorder that progresses with age, with an increasing number of symptoms. Some of the efforts to understand PD progression have been focusing on the regulation of epigenetic mechanisms, that generally include small molecular modifications to the DNA and histones that are essential for regulating gene activity. Here, we have pointed out difficulties to untangle genetic and epigenetic mechanisms, and reviewed several studies that have aimed for untangling. Some of those have enabled more solid claims on independent roles for epigenetic mechanisms. Hereby, evidence that specific DNA hydroxymethylation, global hyperacetylation, and histone deacetylase (HDAC) dependent regulation of SNCA, one of the hallmark genes involved in PD, have become more prominent from the current perspective, than mechanisms that directly involve DNA methylation. In the absence of current epigenetic clinical targets to counteract PD progression, we also hypothesize how several mechanisms may affect local and global epigenetics in PD neurons, including inflammation, oxidative stress, autophagy and DNA repair mechanisms which may lead to future therapeutic targets.

Insights into the biogenesis and potential functions of exonic circular RNA

Circular RNAs (circRNAs) exhibit unique properties due to their covalently closed nature. Models of circRNAs synthesis and function are emerging but much remains undefined about this surprisingly prevalent class of RNA. Here, we identified exonic circRNAs from human and mouse RNA-sequencing datasets, documenting multiple new examples. Addressing function, we found that many circRNAs co-sediment with ribosomes, indicative of their translation potential. By contrast, circRNAs with potential to act as microRNA sponges were scarce, with some support for a collective sponge function by groups of circRNAs. Addressing circRNA biogenesis, we delineated several features commonly associated with circRNA occurrence. CircRNA-producing genes tend to be longer and to contain more exons than average. Back-splice acceptor exons are strongly enriched at ordinal position 2 within genes, and circRNAs typically have a short exon span with two exons being the most prevalent. The flanking introns either side of circRNA loci are exceptionally long. Of note also, single-exon circRNAs derive from unusually long exons while multi-exon circRNAs are mostly generated from exons of regular length. These findings independently validate and extend similar observations made in a number of prior studies. Furthermore, we analysed high-resolution RNA polymerase II occupancy data from two separate human cell lines to reveal distinctive transcription dynamics at circRNA-producing genes. Specifically, RNA polymerase II traverses the introns of these genes at above average speed concomitant with an accentuated slow-down at exons. Collectively, these features indicate how a perturbed balance between transcription and linear splicing creates important preconditions for circRNA production. We speculate that these preconditions need to be in place so that looping interactions between flanking introns can promote back-splicing to raise circRNA production to appreciable levels.

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 .

Structure and function of an ectopic Polycomb chromatin domain

Polycomb group proteins (PcGs) drive target gene repression and form large chromatin domains. In Drosophila, DNA elements known as Polycomb group response elements (PREs) recruit PcGs to the DNA. We have shown that, within the invected-engrailed (inv-en) Polycomb domain, strong, constitutive PREs are dispensable for Polycomb domain structure and function. We suggest that the endogenous chromosomal location imparts stability to this Polycomb domain. To test this possibility, a 79-kb en transgene was inserted into other chromosomal locations. This transgene is functional and forms a Polycomb domain. The spreading of the H3K27me3 repressive mark, characteristic of PcG domains, varies depending on the chromatin context of the transgene. Unlike at the endogenous locus, deletion of the strong, constitutive PREs from the transgene leads to both loss- and gain-of function phenotypes, demonstrating the important role of these regulatory elements. Our data show that chromatin context plays an important role in Polycomb domain structure and function.

Human exonization through differential nucleosome occupancy

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, whereas the roles of nucleosome binding in shaping changes through evolution remain to be addressed. Here we performed a comparative study to clarify the roles of nucleosome occupancy in exon origination. By profiling a high-resolution, cross-species mononucleosome landscape for mammalian tissues, we found nucleosome occupancy profiles are conserved across tissues and species. Further, through a phylogenetic approach, we found that the feature of differential nucleosome occupancy appears prior to the origination of new exons and, presumably, facilitates the origin of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content, which suggests the function of nucleosome binding in exonization.

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, but the roles of nucleosome occupancy in shaping changes through evolution remain to be addressed. Here we present high-resolution nucleosome occupancy profiles for multiple tissues derived from human, macaque, tree shrew, mouse, and pig. Genome-wide comparison reveals conserved nucleosome occupancy profiles across both different species and tissue types. Notably, we found significantly higher levels of nucleosome occupancy in exons than in introns, a pattern correlated with the different exon–intron GC content. We then determined whether this biased occupancy may play roles in the origination of new exons through evolution, rather than being a downstream effect of exonization, through a comparative approach to sequentially trace the order of the exonization and biased nucleosome binding. By identifying recently evolved exons in human but not in macaque using matched RNA sequencing, we found that higher exonic nucleosome occupancy also existed in macaque regions orthologous to these exons. Presumably, such biased nucleosome occupancy facilitates the origination of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content. These data thus support a model that sites bound by nucleosomes are more likely to evolve into exons, which we term the “nucleosome-first” model.

Nucleosome Positioning of Intronless Genes in the Human Genome

Nucleosomes, the basic units of chromatin, are involved in transcription regulation and DNA replication. Intronless genes, which constitute 3 percent of the human genome, differ from intron-containing genes in evolution and function. Our analysis reveals that nucleosome positioning shows a distinct pattern in intronless and intron-containing genes. The nucleosome occupancy upstream of transcription start sites of intronless genes is lower than that of intron-containing genes. In contrast, high occupancy and well positioned nucleosomes are observed along the gene body of intronless genes, which is perfectly consistent with the barrier nucleosome model. Intronless genes have a significantly lower expression level than intron-containing genes and most of them are not expressed in CD4+ T cell lines and GM12878 cell lines, which results from their tissue specificity. However, the highly expressed genes are at the same expression level between the two types of genes. The highly expressed intronless genes require a higher density of RNA Pol II in an elongating state to compensate for the lack of introns. Additionally, 5' and 3' nucleosome depleted regions of highly expressed intronless genes are deeper than those of highly expressed intron-containing genes.

Integrative analysis reveals functional and regulatory roles of H3K79me2 in mediating alternative splicing

Background

Accumulating evidence suggests alternative splicing (AS) is a co-transcriptional splicing process not only controlled by RNA-binding splicing factors, but also mediated by epigenetic regulators, such as chromatin structure, nucleosome density, and histone modification. Aberrant AS plays an important role in regulating various diseases, including cancers.

Methods

In this study, we integrated AS events derived from RNA-seq with H3K79me2 ChIP-seq data across 34 different normal and cancer cell types and found the higher enrichment of H3K79me2 in two AS types, skipping exon (SE) and alternative 3′ splice site (A3SS).

Results

Interestingly, by applying self-organizing map (SOM) clustering, we unveiled two clusters mainly comprised of blood cancer cell types with a strong correlation between H3K79me2 and SE. Remarkably, the expression of transcripts associated with SE was not significantly different from that of those not associated with SE, indicating the involvement of H3K79me2 in splicing has little impact on full mRNA transcription. We further showed that the deletion of DOT1L1, the sole H3K79 methyltransferase, impeded leukemia cell proliferation as well as switched exon skipping to the inclusion isoform in two MLL-rearranged acute myeloid leukemia cell lines. Our data demonstrate H3K79me2 was involved in mediating SE processing, which might in turn influence transformation and disease progression in leukemias.

Conclusions

Collectively, our work for the first time reveals that H3K79me2 plays functional and regulatory roles through a co-transcriptional splicing mechanism.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0538-1) contains supplementary material, which is available to authorized users.

Histone and RNA-binding protein interaction creates crosstalk network for regulation of alternative splicing

Alternative splicing is an essential process in eukaryotes, as it increases the complexity of gene expression by generating multiple proteins from a single pre-mRNA. However, information on the regulatory mechanisms for alternative splicing is lacking, because splicing occurs over a short period via the transient interactions of proteins within functional complexes of the spliceosome. Here, we investigated in detail the molecular mechanisms connecting alternative splicing with epigenetic mechanisms. We identified interactions between histone proteins and splicing factors such as Rbfox2, Rbfox3, and splicing factor proline and glutamine rich protein (SFPQ) by in vivo crosslinking and immunoprecipitation. Furthermore, we confirmed that splicing factors were bound to specific modified residues of histone proteins. Additionally, changes in histone methylation due to histone methyltransferase inhibitor treatment notably affected alternative splicing in selected genes. Therefore, we suggested that there may be crosstalk mechanisms connecting histone modifications and RNA-binding proteins that increase the local concentration of RNA-binding proteins in alternative exon loci of nucleosomes by binding specific modified histone proteins, leading to alternative splicing. This crosstalk mechanism may play a major role in epigenetic processes such as histone modification and the regulation of alternative splicing.