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.

Regulation of alternative splicing through coupling with transcription and chromatin structure.

Alternative precursor messenger RNA (pre-mRNA) splicing plays a pivotal role in the flow of genetic information from DNA to proteins by expanding the coding capacity of genomes. Regulation of alternative splicing is as important as regulation of transcription to determine cell- and tissue-specific features, normal cell functioning, and responses of eukaryotic cells to external cues. Its importance is confirmed by the evolutionary conservation and diversification of alternative splicing and the fact that its deregulation causes hereditary disease and cancer. This review discusses the multiple layers of cotranscriptional regulation of alternative splicing in which chromatin structure, DNA methylation, histone marks, and nucleosome positioning play a fundamental role in providing a dynamic scaffold for interactions between the splicing and transcription machineries. We focus on evidence for how the kinetics of RNA polymerase II (RNAPII) elongation and the recruitment of splicing factors and adaptor proteins to chromatin components act in coordination to regulate alternative splicing.

Gene architecture directs splicing outcome in separate nuclear spatial regions.

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

Bi-allelic variants in neuronal cell adhesion molecule cause a neurodevelopmental disorder characterized by developmental delay, hypotonia, neuropathy/spasticity.

Cell adhesion molecules are membrane-bound proteins predominantly expressed in the central nervous system along principal axonal pathways with key roles in nervous system development, neural cell differentiation and migration, axonal growth and guidance, myelination, and synapse formation. Here, we describe ten affected individuals with bi-allelic variants in the neuronal cell adhesion molecule NRCAM that lead to a neurodevelopmental syndrome of varying severity; the individuals are from eight families. This syndrome is characterized by developmental delay/intellectual disability, hypotonia, peripheral neuropathy, and/or spasticity. Computational analyses of NRCAM variants, many of which cluster in the third fibronectin type III (Fn-III) domain, strongly suggest a deleterious effect on NRCAM structure and function, including possible disruption of its interactions with other proteins. These findings are corroborated by previous in vitro studies of murine Nrcam-deficient cells, revealing abnormal neurite outgrowth, synaptogenesis, and formation of nodes of Ranvier on myelinated axons. Our studies on zebrafish nrcamaΔ mutants lacking the third Fn-III domain revealed that mutant larvae displayed significantly altered swimming behavior compared to wild-type larvae (p < 0.03). Moreover, nrcamaΔ mutants displayed a trend toward increased amounts of α-tubulin fibers in the dorsal telencephalon, demonstrating an alteration in white matter tracts and projections. Taken together, our study provides evidence that NRCAM disruption causes a variable form of a neurodevelopmental disorder and broadens the knowledge on the growing role of the cell adhesion molecule family in the nervous system.

Alternative splicing regulates biogenesis of miRNAs located across exon-intron junctions.

The initial step in microRNA (miRNA) biogenesis requires processing of the precursor miRNA (pre-miRNA) from a longer primary transcript. Many pre-miRNAs originate from introns, and both a mature miRNA and a spliced RNA can be generated from the same transcription unit. We have identified a mechanism in which RNA splicing negatively regulates the processing of pre-miRNAs that overlap exon-intron junctions. Computational analysis identified dozens of such pre-miRNAs, and experimental validation demonstrated competitive interaction between the Microprocessor complex and the splicing machinery. Tissue-specific alternative splicing regulates maturation of one such miRNA, miR-412, resulting in effects on its targets that code a protein network involved in neuronal cell death processes. This mode of regulation specifically controls maturation of splice-site-overlapping pre-miRNAs but not pre-miRNAs located completely within introns or exons of the same transcript. Our data present a biological role of alternative splicing in regulation of miRNA biogenesis.

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.

The importance of DNA methylation of exons on alternative splicing.

Alternative splicing (AS) contributes to proteome diversity. As splicing occurs cotranscriptionally, epigenetic determinants such as DNA methylation likely play a part in regulation of AS. Previously, we have shown that DNA methylation marks exons and that a loss of DNA methylation alters splicing patterns in a genome-wide manner. To investigate the influence of DNA methylation on splicing of individual genes, we developed a method to manipulate DNA methylation in vivo in a site-specific manner using the deactivated endonuclease Cas9 fused to enzymes that methylate or demethylate DNA. We used this system to directly change the DNA methylation pattern of selected exons and introns. We demonstrated that changes in the methylation pattern of alternatively spliced exons, but not constitutively spliced exons or introns, altered inclusion levels. This is the first direct demonstration that DNA methylation of exon-encoding regions is directly involved in regulation of AS.

How Are Short Exons Flanked by Long Introns Defined and Committed to Splicing?

The splice sites (SSs) delimiting an intron are brought together in the earliest step of spliceosome assembly yet it remains obscure how SS pairing occurs, especially when introns are thousands of nucleotides long. Splicing occurs in vivo in mammals within minutes regardless of intron length, implying that SS pairing can instantly follow transcription. Also, factors required for SS pairing, such as the U1 small nuclear ribonucleoprotein (snRNP) and U2AF65, associate with RNA polymerase II (RNAPII), while nucleosomes preferentially bind exonic sequences and associate with U2 snRNP. Based on recent publications, we assume that the 5' SS-bound U1 snRNP can remain tethered to RNAPII until complete synthesis of the downstream intron and exon. An additional U1 snRNP then binds the downstream 5' SS, whereas the RNAPII-associated U2AF65 binds the upstream 3' SS to facilitate SS pairing along with exon definition. Next, the nucleosome-associated U2 snRNP binds the branch site to advance splicing complex assembly. This may explain how RNAPII and chromatin are involved in spliceosome assembly and how introns lengthened during evolution with a relatively minimal compromise in splicing.

A network-based analysis of colon cancer splicing changes reveals a tumorigenesis-favoring regulatory pathway emanating from ELK1.

Splicing aberrations are prominent drivers of cancer, yet the regulatory pathways controlling them are mostly unknown. Here we develop a method that integrates physical interaction, gene expression, and alternative splicing data to construct the largest map of transcriptomic and proteomic interactions leading to cancerous splicing aberrations defined to date, and identify driver pathways therein. We apply our method to colon adenocarcinoma and non-small-cell lung carcinoma. By focusing on colon cancer, we reveal a novel tumor-favoring regulatory pathway involving the induction of the transcription factor MYC by the transcription factor ELK1, as well as the subsequent induction of the alternative splicing factor PTBP1 by both. We show that PTBP1 promotes specific RAC1,NUMB, and PKM splicing isoforms that are major triggers of colon tumorigenesis. By testing the pathway's activity in patient tumor samples, we find ELK1,MYC, and PTBP1 to be overexpressed in conjunction with oncogenic KRAS mutations, and show that these mutations increase ELK1 levels via the RAS-MAPK pathway. We thus illuminate, for the first time, a full regulatory pathway connecting prevalent cancerous mutations to functional tumor-inducing splicing aberrations. Our results demonstrate our method is applicable to different cancers to reveal regulatory pathways promoting splicing aberrations.

Phosphatidylserine Ameliorates Neurodegenerative Symptoms and Enhances Axonal Transport in a Mouse Model of Familial Dysautonomia.

Familial Dysautonomia (FD) is a neurodegenerative disease in which aberrant tissue-specific splicing of IKBKAP exon 20 leads to reduction of IKAP protein levels in neuronal tissues. Here we generated a conditional knockout (CKO) mouse in which exon 20 of IKBKAP is deleted in the nervous system. The CKO FD mice exhibit developmental delays, sensory abnormalities, and less organized dorsal root ganglia (DRGs) with attenuated axons compared to wild-type mice. Furthermore, the CKO FD DRGs show elevated HDAC6 levels, reduced acetylated α-tubulin, unstable microtubules, and impairment of axonal retrograde transport of nerve growth factor (NGF). These abnormalities in DRG properties underlie neuronal degeneration and FD symptoms. Phosphatidylserine treatment decreased HDAC6 levels and thus increased acetylation of α-tubulin. Further PS treatment resulted in recovery of axonal outgrowth and enhanced retrograde axonal transport by decreasing histone deacetylase 6 (HDAC6) levels and thus increasing acetylation of α-tubulin levels. Thus, we have identified the molecular pathway that leads to neurodegeneration in FD and have demonstrated that phosphatidylserine treatment has the potential to slow progression of neurodegeneration.

Phosphatidylserine enhances IKBKAP transcription by activating the MAPK/ERK signaling pathway.

Familial dysautonomia (FD) is a genetic disorder manifested due to abnormal development and progressive degeneration of the sensory and autonomic nervous system. FD is caused by a point mutation in the IKBKAP gene encoding the IKAP protein, resulting in decreased protein levels. A promising potential treatment for FD is phosphatidylserine (PS); however, the manner by which PS elevates IKAP levels has yet to be identified. Analysis of ChIP-seq results of the IKBKAP promoter region revealed binding of the transcription factors CREB and ELK1, which are regulated by the mitogen-activated protein kinase (MAPK)/extracellular-regulated kinase (ERK) signaling pathway. We show that PS treatment enhanced ERK phosphorylation in cells derived from FD patients. ERK activation resulted in elevated IKBKAP transcription and IKAP protein levels, whereas pretreatment with the MAPK inhibitor U0126 blocked elevation of the IKAP protein level. Overexpression of either ELK1 or CREB activated the IKBKAP promoter, whereas downregulation of these transcription factors resulted in a decrease of the IKAP protein. Additionally, we show that PS improves cell migration, known to be enhanced by MAPK/ERK activation and abrogated in FD cells. In conclusion, our results demonstrate that PS activates the MAPK/ERK signaling pathway, resulting in activation of transcription factors that bind the promoter region of IKBKAP and thus enhancing its transcription. Therefore, compounds that activate the MAPK/ERK signaling pathway could constitute potential treatments for FD.

The alternative role of DNA methylation in splicing regulation.

Although DNA methylation was originally thought to only affect transcription, emerging evidence shows that it also regulates alternative splicing. Exons, and especially splice sites, have higher levels of DNA methylation than flanking introns, and the splicing of about 22% of alternative exons is regulated by DNA methylation. Two different mechanisms convey DNA methylation information into the regulation of alternative splicing. The first involves modulation of the elongation rate of RNA polymerase II (Pol II) by CCCTC-binding factor (CTCF) and methyl-CpG binding protein 2 (MeCP2); the second involves the formation of a protein bridge by heterochromatin protein 1 (HP1) that recruits splicing factors onto transcribed alternative exons. These two mechanisms, however, regulate only a fraction of such events, implying that more underlying mechanisms remain to be found.

Cotranscriptional histone H2B monoubiquitylation is tightly coupled with RNA polymerase II elongation rate.

Various histone modifications decorate nucleosomes within transcribed genes. Among these, monoubiquitylation of histone H2B (H2Bub1) and methylation of histone H3 on lysines 36 (H3K36me2/3) and 79 (H3K79me2/3) correlate positively with gene expression. By measuring the progression of the transcriptional machinery along genes within live cells, we now report that H2B monoubiquitylation occurs cotranscriptionally and accurately reflects the advance of RNA polymerase II (Pol II). In contrast, H3K36me3 and H3K79me2 are less dynamic and represent Pol II movement less faithfully. High-resolution ChIP-seq reveals that H2Bub1 levels are selectively reduced at exons and decrease in an exon-dependent stepwise manner toward the 3' end of genes. Exonic depletion of H2Bub1 in gene bodies is highly correlated with Pol II pausing at exons, suggesting elongation rate changes associated with intron-exon structure. In support of this notion, H2Bub1 levels were found to be significantly correlated with transcription elongation rates measured in various cell lines. Overall, our data shed light on the organization of H2Bub1 within transcribed genes and single out H2Bub1 as a reliable marker for ongoing transcription elongation.

DNA-methylation effect on cotranscriptional splicing is dependent on GC architecture of the exon-intron structure.

DNA methylation is known to regulate transcription and was recently found to be involved in exon recognition via cotranscriptional splicing. We recently observed that exon-intron architectures can be grouped into two classes: one with higher GC content in exons compared to the flanking introns, and the other with similar GC content in exons and introns. The first group has higher nucleosome occupancy on exons than introns, whereas the second group exhibits weak nucleosome marking of exons, suggesting another type of epigenetic marker distinguishes exons from introns when GC content is similar. We find different and specific patterns of DNA methylation in each of the GC architectures; yet in both groups, DNA methylation clearly marks the exons. Exons of the leveled GC architecture exhibit a significantly stronger DNA methylation signal in relation to their flanking introns compared to exons of the differential GC architecture. This is accentuated by a reduction of the DNA methylation level in the intronic sequences in proximity to the splice sites and shows that different epigenetic modifications mark the location of exons already at the DNA level. Also, lower levels of methylated CpGs on alternative exons can successfully distinguish alternative exons from constitutive ones. Three positions at the splice sites show high CpG abundance and accompany elevated nucleosome occupancy in a leveled GC architecture. Overall, these results suggest that DNA methylation affects exon recognition and is influenced by the GC architecture of the exon and flanking introns.

Alternative splicing and evolution: diversification, exon definition and function.

Over the past decade, it has been shown that alternative splicing (AS) is a major mechanism for the enhancement of transcriptome and proteome diversity, particularly in mammals. Splicing can be found in species from bacteria to humans, but its prevalence and characteristics vary considerably. Evolutionary studies are helping to address questions that are fundamental to understanding this important process: how and when did AS evolve? Which AS events are functional? What are the evolutionary forces that shaped, and continue to shape, AS? And what determines whether an exon is spliced in a constitutive or alternative manner? In this Review, we summarize the current knowledge of AS and evolution and provide insights into some of these unresolved questions.

Phosphatidylserine increases IKBKAP levels in a humanized knock-in IKBKAP mouse model.

Familial dysautonomia (FD) is a severe neurodegenerative genetic disorder restricted to the Ashkenazi Jewish population. The most common mutation in FD patients is a T-to-C transition at position 6 of intron 20 of the IKBKAP gene. This mutation causes aberrant skipping of exon 20 in a tissue-specific manner, leading to reduction of the IκB kinase complex-associated protein (IKAP) protein in the nervous system. We established a homozygous humanized mouse strain carrying human exon 20 and its two flanking introns; the 3' intron has the transition observed in the IKBKAP gene of FD patients. Although our FD humanized mouse does not display FD symptoms, the unique, tissue-specific splicing pattern of the IKBKAP in these mice allowed us to evaluate the effect of therapies on gene expression and exon 20 splicing. The FD mice were supplemented with phosphatidylserine (PS), a safe food supplement that increases mRNA and protein levels of IKBKAP in cell lines generated from FD patients. Here we demonstrated that PS treatment increases IKBAKP mRNA and IKAP protein levels in various tissues of FD mice without affecting exon 20 inclusion levels. We also observed that genes associated with transcription regulation and developmental processes were up-regulated in the cerebrum of PS-treated mice. Thus, PS holds promise for the treatment of FD.

Changes in exon-intron structure during vertebrate evolution affect the splicing pattern of exons.

Exon-intron architecture is one of the major features directing the splicing machinery to the short exons that are located within long flanking introns. However, the evolutionary dynamics of exon-intron architecture and its impact on splicing is largely unknown. Using a comparative genomic approach, we analyzed 17 vertebrate genomes and reconstructed the ancestral motifs of both 3' and 5' splice sites, as also the ancestral length of exons and introns. Our analyses suggest that vertebrate introns increased in length from the shortest ancestral introns to the longest primate introns. An evolutionary analysis of splice sites revealed that weak splice sites act as a restrictive force keeping introns short. In contrast, strong splice sites allow recognition of exons flanked by long introns. Reconstruction of the ancestral state suggests these phenomena were not prevalent in the vertebrate ancestor, but appeared during vertebrate evolution. By calculating evolutionary rate shifts in exons, we identified cis-acting regulatory sequences that became fixed during the transition from early vertebrates to mammals. Experimental validations performed on a selection of these hexamers confirmed their regulatory function. We additionally revealed many features of exons that can discriminate alternative from constitutive exons. These features were integrated into a machine-learning approach to predict whether an exon is alternative. Our algorithm obtains very high predictive power (AUC of 0.91), and using these predictions we have identified and successfully validated novel alternatively spliced exons. Overall, we provide novel insights regarding the evolutionary constraints acting upon exons and their recognition by the splicing machinery.

Chromatin density and splicing destiny: on the cross-talk between chromatin structure and splicing.

How are short exonic sequences recognized within the vast intronic oceans in which they reside? Despite decades of research, this remains one of the most fundamental, yet enigmatic, questions in the field of pre-mRNA splicing research. For many years, studies aiming to shed light on this process were focused at the RNA level, characterizing the manner by which splicing factors and auxiliary proteins interact with splicing signals, thereby enabling, facilitating and regulating splicing. However, we increasingly understand that splicing is not an isolated process; rather it occurs co-transcriptionally and is presumably also regulated by transcription-related processes. In fact, studies by our group and others over the past year suggest that DNA structure in terms of nucleosome positioning and specific histone modifications, which have a well established role in transcription, may also have a role in splicing. In this review we discuss evidence for the coupling between transcription and splicing, focusing on recent findings suggesting a link between chromatin structure and splicing, and highlighting challenges this emerging field is facing.

A network-based method for associating genes with autism spectrum disorder.

Autism spectrum disorder (ASD) is a highly heritable complex disease that affects 1% of the population, yet its underlying molecular mechanisms are largely unknown. Here we study the problem of predicting causal genes for ASD by combining genome-scale data with a network propagation approach. We construct a predictor that integrates multiple omic data sets that assess genomic, transcriptomic, proteomic, and phosphoproteomic associations with ASD. In cross validation our predictor yields mean area under the ROC curve of 0.87 and area under the precision-recall curve of 0.89. We further show that it outperforms previous gene-level predictors of autism association. Finally, we show that we can use the model to predict genes associated with Schizophrenia which is known to share genetic components with ASD.

IKAP/Elp1 involvement in cytoskeleton regulation and implication for familial dysautonomia.

Deficiency in the IKAP/Elp1 protein leads to the recessive sensory autosomal congenital neuropathy which is called familial dysautonomia (FD). This protein was originally identified as a role player in transcriptional elongation being a subunit of the RNAPII transcriptional Elongator multi-protein complex. Subsequently, IKAP/Elp1 was shown to play various functions in the cytoplasm. Here, we describe experiments performed with IKAP/Elp1 downregulated cell lines and FD-derived cells and tissues. Immunostaining of the cytoskeleton component α-tubulin in IKAP/Elp1 downregulated cells revealed disorganization of the microtubules (MTs) that was reflected in aberrant cell shape and process formation. In contrast to a recent report on the decrease in α-tubulin acetylation in IKAP/Elp1 downregulated cells, we were unable to observe any effect of IKAP/Elp1 deficiency on α-tubulin acetylation in the FD cerebrum and in a variety of IKAP/Elp1 downregulated cell lines. To explore possible candidates involved in the observed aberrations in MTs, we focused on superior cervical ganglion-10 protein (SCG10), also called STMN2, which is known to be an MT destabilizing protein. We have found that SCG10 is upregulated in the IKAP/Elp1-deficient FD cerebrum, FD fibroblasts and in IKAP/Elp1 downregulated neuroblastoma cell line. To better understand the effect of IKAP/Elp1 deficiency on SCG10 expression, we investigated the possible involvement of RE-1-silencing transcription factor (REST), a known repressor of the SCG10 gene. Indeed, REST was downregulated in the IKAP/Elp1-deficient FD cerebrum and IKAP/Elp1 downregulated neuroblastoma cell line. These results could shed light on a possible link between IKAP/Elp1 deficiency and cytoskeleton destabilization.