Abstract: Although the function of DNA methylation in gene promoter regions is well established in transcriptional repression, the function of the evolutionarily conserved widespread distribution of DNA methylation in gene body regions remains incompletely understood. Here, we show that DNA methylation is enriched in included alternatively spliced exons (ASEs), and that inhibition of DNA methylation results in aberrant splicing of ASEs. The methyl-CpG-binding protein MeCP2 is enriched in included ASEs, particularly those that are also highly methylated, and inhibition of DNA methylation disrupts specific targeting of MeCP2 to exons. Interestingly, ablation of MeCP2 results in increased histone acetylation and aberrant ASE-skipping events. We further show that inhibition of histone deacetylase (HDAC) activity leads to exon skipping that shows a highly significant degree of overlap with that caused by MeCP2 knockdown. Together, our data indicate that intragenic DNA methylation operates in exon definition to modulate alternative RNA splicing and can enhance exon recognition via recruitment of the multifunctional protein MeCP2, which thereby maintains local histone hypoacetylation through the subsequent recruitment of HDACs.

Abstract: Huntington disease (HD) is a devastating, late-onset, inherited neurodegenerative disorder that manifests with personality changes, movement disorders, and cognitive decline. It is caused by a CAG repeat expansion in exon 1 of the HTT gene that translates to a polyglutamine tract in the huntingtin protein (HTT). The formation of HTT fragments has been implicated as an essential step in the molecular pathogenesis of HD and several proteases that cleave HTT have been identified. However, the importance of smaller N-terminal fragments has been highlighted by their presence in HD postmortem brains and by the fact that nuclear inclusions are only detected by antibodies to the N terminus of HTT. Despite an intense research effort, the precise length of these fragments and the mechanism by which they are generated remains unknown. Here we show that CAG repeat length–dependent aberrant splicing of exon 1 HTT results in a short polyadenylated mRNA that is translated into an exon 1 HTT protein. Given that mutant exon 1 HTT proteins have consistently been shown to be highly pathogenic in HD mouse models, the aberrant splicing of HTT mRNA provides a mechanistic basis for the molecular pathogenesis of HD. RNA-targeted therapeutic strategies designed to lower the levels of HTT are under development. Many of these approaches would not prevent the production of exon 1 HTT and should be reviewed in light of our findings.

Abstract: Transactivating response region DNA binding protein (TDP-43) is the major protein component of ubiquitinated inclusions found in amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) with ubiquitinated inclusions. Two ALS-causing mutants (TDP-43(Q331K) and TDP-43(M337V)), but not wild-type human TDP-43, are shown here to provoke age-dependent, mutant-dependent, progressive motor axon degeneration and motor neuron death when expressed in mice at levels and in a cell type-selective pattern similar to endogenous TDP-43. Mutant TDP-43-dependent degeneration of lower motor neurons occurs without: (i) loss of TDP-43 from the corresponding nuclei, (ii) accumulation of TDP-43 aggregates, and (iii) accumulation of insoluble TDP-43. Computational analysis using splicing-sensitive microarrays demonstrates alterations of endogenous TDP-43-dependent alternative splicing events conferred by both human wild-type and mutant TDP-43(Q331K), but with high levels of mutant TDP-43 preferentially enhancing exon exclusion of some target pre-mRNAs affecting genes involved in neurological transmission and function. Comparison with splicing alterations following TDP-43 depletion demonstrates that TDP-43(Q331K) enhances normal TDP-43 splicing function for some RNA targets but loss-of-function for others. Thus, adult-onset motor neuron disease does not require aggregation or loss of nuclear TDP-43, with ALS-linked mutants producing loss and gain of splicing function of selected RNA targets at an early disease stage.

Abstract: Alternative cleavage and polyadenylation (APA) generates diverse mRNA isoforms. We developed 3' region extraction and deep sequencing (3'READS) to address mispriming issues that commonly plague poly(A) site (pA) identification, and we used the method to comprehensively map pAs in the mouse genome. Thorough annotation of gene 3' ends revealed over 5,000 previously overlooked pAs (∼8% of total) flanked by A-rich sequences, underscoring the necessity of using an accurate tool for pA mapping. About 79% of mRNA genes and 66% of long noncoding RNA genes undergo APA, but these two gene types have distinct usage patterns for pAs in introns and upstream exons. Quantitative analysis of APA isoforms by 3'READS indicated that promoter-distal pAs, regardless of intron or exon locations, become more abundant during embryonic development and cell differentiation and that upregulated isoforms have stronger pAs, suggesting global modulation of the 3' end-processing activity in development and differentiation.

Abstract: In contrast to other primary epidermal growth factor receptor (EGFR) mutations in lung adenocarcinomas, insertions in exon 20 of EGFR have been generally associated with resistance to EGFR-tyrosine kinase inhibitors. Their molecular spectrum, clinicopathologic characteristics, and prevalence are not well established. Tumors harboring EGFR exon 20 insertions were identified through an algorithmic screen of 1,500 lung adenocarcinomas. Cases were first tested for common mutations in EGFR (exons 19 and 21) and KRAS (exon 2) and, if negative, further analyzed for EGFR exon 20 insertions. All samples underwent extended genotyping for other driver mutations in EGFR, KRAS, BRAF, ERBB2/HER2, NRAS, PIK3CA, MEK1, and AKT by mass spectrometry; a subset was evaluated for ALK rearrangements. We identified 33 EGFR exon 20 insertion cases [2.2%, 95% confidence interval (CI), 1.6-3.1], all mutually exclusive with mutations in the other genes tested (except PIK3CA). They were more common among never-smokers (P < 0.0001). There was no association with age, sex, race, or stage. Morphologically, tumors were similar to those with common EGFR mutations but with frequent solid histology. Insertions were highly variable in position and size, ranging from 3 to 12 bp, resulting in 13 different insertions, which, by molecular modeling, are predicted to have potentially different effects on erlotinib binding. EGFR exon 20 insertion testing identifies a distinct subset of lung adenocarcinomas, accounting for at least 9% of all EGFR-mutated cases, representing the third most common type of EGFR mutation after exon 19 deletions and L858R. Insertions are structurally heterogeneous with potential implications for response to EGFR inhibitors.

Abstract: Alternative splicing (AS) enables programmed diversity of gene expression across tissues and development. We show here that binding in distal intronic regions (>500 nucleotides (nt) from any exon) by Rbfox splicing factors important in development is extensive and is an active mode of splicing regulation. Similarly to exon-proximal sites, distal sites contain evolutionarily conserved GCATG sequences and are associated with AS activation and repression upon modulation of Rbfox abundance in human and mouse experimental systems. As a proof of principle, we validated the activity of two specific Rbfox enhancers in KIF21A and ENAH distal introns and showed that a conserved long-range RNA-RNA base-pairing interaction (an RNA bridge) is necessary for Rbfox-mediated exon inclusion in the ENAH gene. Thus we demonstrate a previously unknown RNA-mediated mechanism for AS control by distally bound RNA-binding proteins.

Abstract: Alternative splicing is widely recognized for its roles in regulating genes and creating gene diversity. However, despite many efforts, the repertoire of gene splicing variation is still incompletely characterized, even in humans. Here we describe a new computational system, ASprofile, and its application to RNA-seq data from Illumina’s Human Body Map project (>2.5 billion reads). Using the system, we identified putative alternative splicing events in 16 different human tissues, which provide a dynamic picture of splicing variation across the tissues. We detected 26,989 potential exon skipping events representing differences in splicing patterns among the tissues. A large proportion of the events (>60%) were novel, involving new exons (~3000), new introns (~16000), or both. When tracing these events across the sixteen tissues, only a small number (4-7%) appeared to be differentially expressed (‘switched’) between two tissues, while 30-45% showed little variation, and the remaining 50-65% were not present in one or both tissues compared. Novel exon skipping events appeared to be slightly less variable than known events, but were more tissue-specific. Our study represents the first effort to build a comprehensive catalog of alternative splicing in normal human tissues from RNA-seq data, while providing insights into the role of alternative splicing in shaping tissue transcriptome differences. The catalog of events and the ASprofile software are freely available from the Zenodo repository ( http://zenodo.org/record/7068 ; doi: 10.5281/zenodo.7068 ) and from our web site http://ccb.jhu.edu/software/ASprofile .

Abstract: TDP-43 encodes an alternative-splicing regulator with tandem RNA-recognition motifs (RRMs) The protein regulates cystic fibrosis transmembrane regulator (CFTR) exon 9 splicing through binding to long UG-rich RNA sequences and is found in cytoplasmic inclusions of several neurodegenerative diseases We solved the solution structure of the TDP-43 RRMs in complex with UG-rich RNA Ten nucleotides are bound by both RRMs, and six are recognized sequence specifically Among these, a central G interacts with both RRMs and stabilizes a new tandem RRM arrangement Mutations that eliminate recognition of this key nucleotide or crucial inter-RRM interactions disrupt RNA binding and TDP-43-dependent splicing regulation In contrast, point mutations that affect base-specific recognition in either RRM have weaker effects Our findings reveal not only how TDP-43 recognizes UG repeats but also how RNA binding-dependent inter-RRM interactions are crucial for TDP-43 function

Abstract: One of the fundamental issues in RNA splicing research is represented by understanding how the spliceosome can successfully define exons and introns in a huge variety of pre-mRNA molecules with nucleotide-precision. Since its first description, researchers in this field have identified and characterized many fundamental elements and players capable of affecting the splicing process, both in a negative and positive manner. Indeed, it can be argued that today we know a great deal about the forces that make an exon, an exon and an intron, an intron. As will be discussed in this review, these decisions are a result of a complex combinatorial control resulting from many different factors/influences. Most importantly, these influences act across several levels of complexity starting from the relatively simple interaction between two consensus 5' and 3' splice sites to much more complex factors: such as the interplay between silencer or enhancer sequences, transcriptional processivity, genomic milieu, nucleosome positioning, and histone modifications at the chromatin level. Depending on local contexts, all these factors will act either antagonistically or synergistically to decide the exon/intron fate of any given RNA sequence. At present, however, what we still lack is a precise understanding of how all these processes add up to help the spliceosome reach a decision. Therefore, it is expected that future challenges in splicing research will be the careful characterization of all these influences to improve our ability to predict splicing choices in different organisms or in specific contexts.

Abstract: This study identifies MBNL proteins as negative regulators of alternative splicing events that are differentially regulated between ES cells and other cell types; several lines of evidence show that these proteins repress an ES cell alternative splicing program and the reprogramming of somatic cells to induced pluripotent stem cells. Ben Blencowe and colleagues identify the muscleblind-like RNA binding proteins MBNL1 and MBNL2 as negative regulators of alternative splicing events that are differentially regulated between embryonic stem cells and other cell types. Several lines of evidence show that they are involved in the regulation of embryonic-stem-cell-like alternative splicing patterns. The authors also identify a regulatory role during the reprogramming of fibroblasts to induced pluripotent stem (iPS) cells. Previous investigations of the core gene regulatory circuitry that controls the pluripotency of embryonic stem (ES) cells have largely focused on the roles of transcription, chromatin and non-coding RNA regulators1,2,3. Alternative splicing represents a widely acting mode of gene regulation4,5,6,7,8, yet its role in regulating ES-cell pluripotency and differentiation is poorly understood. Here we identify the muscleblind-like RNA binding proteins, MBNL1 and MBNL2, as conserved and direct negative regulators of a large program of cassette exon alternative splicing events that are differentially regulated between ES cells and other cell types. Knockdown of MBNL proteins in differentiated cells causes switching to an ES-cell-like alternative splicing pattern for approximately half of these events, whereas overexpression of MBNL proteins in ES cells promotes differentiated-cell-like alternative splicing patterns. Among the MBNL-regulated events is an ES-cell-specific alternative splicing switch in the forkhead family transcription factor FOXP1 that controls pluripotency9. Consistent with a central and negative regulatory role for MBNL proteins in pluripotency, their knockdown significantly enhances the expression of key pluripotency genes and the formation of induced pluripotent stem cells during somatic cell reprogramming.

Abstract: CLP1 was the first mammalian RNA kinase to be identified. However, determining its in vivo function has been elusive. Here we generated kinase-dead Clp1 (Clp1(K/K)) mice that show a progressive loss of spinal motor neurons associated with axonal degeneration in the peripheral nerves and denervation of neuromuscular junctions, resulting in impaired motor function, muscle weakness, paralysis and fatal respiratory failure. Transgenic rescue experiments show that CLP1 functions in motor neurons. Mechanistically, loss of CLP1 activity results in accumulation of a novel set of small RNA fragments, derived from aberrant processing of tyrosine pre-transfer RNA. These tRNA fragments sensitize cells to oxidative-stress-induced p53 (also known as TRP53) activation and p53-dependent cell death. Genetic inactivation of p53 rescues Clp1(K/K) mice from the motor neuron loss, muscle denervation and respiratory failure. Our experiments uncover a mechanistic link between tRNA processing, formation of a new RNA species and progressive loss of lower motor neurons regulated by p53.

Abstract: Pollen grains of Arabidopsis (Arabidopsis thaliana) contain two haploid sperm cells enclosed in a haploid vegetative cell. Upon germination, the vegetative cell extrudes a pollen tube that carries the sperm to an ovule for fertilization. Knowing the identity, relative abundance, and splicing patterns of pollen transcripts will improve our understanding of pollen and allow investigation of tissue-specific splicing in plants. Most Arabidopsis pollen transcriptome studies have used the ATH1 microarray, which does not assay splice variants and lacks specific probe sets for many genes. To investigate the pollen transcriptome, we performed high-throughput sequencing (RNA-Seq) of Arabidopsis pollen and seedlings for comparison. Gene expression was more diverse in seedling, and genes involved in cell wall biogenesis were highly expressed in pollen. RNA-Seq detected at least 4,172 protein-coding genes expressed in pollen, including 289 assayed only by nonspecific probe sets. Additional exons and previously unannotated 5′ and 3′ untranslated regions for pollen-expressed genes were revealed. We detected regions in the genome not previously annotated as expressed; 14 were tested and 12 were confirmed by polymerase chain reaction. Gapped read alignments revealed 1,908 high-confidence new splicing events supported by 10 or more spliced read alignments. Alternative splicing patterns in pollen and seedling were highly correlated. For most alternatively spliced genes, the ratio of variants in pollen and seedling was similar, except for some encoding proteins involved in RNA splicing. This study highlights the robustness of splicing patterns in plants and the importance of ongoing annotation and visualization of RNA-Seq data using interactive tools such as Integrated Genome Browser.

Abstract: Alternative splicing is regulated by splicing factors that modulate splice site selection. In some cases, however, splicing factors show antagonistic activities by either activating or repressing splicing. Here, we show that these opposing outcomes are based on their binding location relative to regulated 59 splice sites. SR proteins enhance splicing only when they are recruited to the exon. However, they interfere with splicing by simply relocating them to the opposite intronic side of the splice site. hnRNP splicing factors display analogous opposing activities, but in a reversed position dependence. Activation by SR or hnRNP proteins increases splice site recognition at the earliest steps of exon definition, whereas splicing repression promotes the assembly of nonproductive complexes that arrest spliceosome assembly prior to splice site pairing. Thus, SR and hnRNP splicing factors exploit similar mechanisms to positively or negatively influence splice site selection.

Abstract: 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.

Abstract: The incidence of lung diseases and cancer caused by cigarette smoke is increasing. The molecular mechanisms of gene regulation induced by cigarette smoke that ultimately lead to cancer remain unclear. This report describes a novel long noncoding RNA (lncRNA) that is induced by cigarette smoke extract (CSE) both in vitro and in vivo and is elevated in numerous lung cancer cell lines. We have termed this lncRNA the smoke and cancer–associated lncRNA–1 (SCAL1). This lncRNA is located in chromosome 5, and initial sequencing analysis reveals a transcript with four exons and three introns. The expression of SCAL1 is regulated transcriptionally by nuclear factor erythroid 2–related factor (NRF2), as determined by the small, interfering RNA (siRNA) knockdown of NRF2 and kelch-like ECH-associated protein 1 (KEAP1). A nuclear factor erythroid-derived 2 (NF-E2) motif was identified in the promoter region that shows binding to NRF2 after its activation. Functionally, the siRNA knockdown of SCAL1 in human bronchial epithelial cells shows a significant potentiation of cytotoxicity induced by CSE in vitro. Altogether, these results identify a novel and intriguing new noncoding RNA that may act downstream of NRF2 to regulate gene expression and mediate oxidative stress protection in airway epithelial cells.

Abstract: While splicing differences between tissues, sexes and species are well documented, little is known about the extent and the nature of splicing changes that take place during human or mammalian development and aging. Here, using high-throughput transcriptome sequencing, we have characterized splicing changes that take place during whole human lifespan in two brain regions: prefrontal cortex and cerebellum. Identified changes were confirmed using independent human and rhesus macaque RNA-seq data sets, exon arrays and PCR, and were detected at the protein level using mass spectrometry. Splicing changes across lifespan were abundant in both of the brain regions studied, affecting more than a third of the genes expressed in the human brain. Approximately 15% of these changes differed between the two brain regions. Across lifespan, splicing changes followed discrete patterns that could be linked to neural functions, and associated with the expression profiles of the corresponding splicing factors. More than 60% of all splicing changes represented a single splicing pattern reflecting preferential inclusion of gene segments potentially targeting transcripts for nonsense-mediated decay in infants and elderly.

Abstract: Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

Abstract: The gene encoding a DNA/RNA binding protein FUS/TLS is frequently mutated in amyotrophic lateral sclerosis (ALS). Mutations commonly affect its carboxy-terminal nuclear localization signal, resulting in varying deficiencies of FUS nuclear localization and abnormal cytoplasmic accumulation. Increasing evidence suggests deficiencies in FUS nuclear function may contribute to neuron degeneration. Here we report a novel FUS autoregulatory mechanism and its deficiency in ALS-associated mutants. Using FUS CLIP-seq, we identified significant FUS binding to a highly conserved region of exon 7 and the flanking introns of its own pre-mRNAs. We demonstrated that FUS is a repressor of exon 7 splicing and that the exon 7-skipped splice variant is subject to nonsense-mediated decay (NMD). Overexpression of FUS led to the repression of exon 7 splicing and a reduction of endogenous FUS protein. Conversely, the repression of exon 7 was reduced by knockdown of FUS protein, and moreover, it was rescued by expression of EGFP-FUS. This dynamic regulation of alternative splicing describes a novel mechanism of FUS autoregulation. Given that ALS-associated FUS mutants are deficient in nuclear localization, we examined whether cells expressing these mutants would be deficient in repressing exon 7 splicing. We showed that FUS harbouring R521G, R522G or ΔExon15 mutation (minor, moderate or severe cytoplasmic localization, respectively) directly correlated with respectively increasing deficiencies in both exon 7 repression and autoregulation of its own protein levels. These data suggest that compromised FUS autoregulation can directly exacerbate the pathogenic accumulation of cytoplasmic FUS protein in ALS. We showed that exon 7 skipping can be induced by antisense oligonucleotides targeting its flanking splice sites, indicating the potential to alleviate abnormal cytoplasmic FUS accumulation in ALS. Taken together, FUS autoregulation by alternative splicing provides insight into a molecular mechanism by which FUS-regulated pre-mRNA processing can impact a significant number of targets important to neurodegeneration.

Abstract: Temperature is one of key environmental parameters that affect the whole life of fishes and an increasing number of studies have been directed towards understanding the mechanisms of cold acclimation in fish. However, the adaptation of larvae to cold stress and the cold-specific transcriptional alterations in fish larvae remain largely unknown. In this study, we characterized the development of cold-tolerance in zebrafish larvae and investigated the transcriptional profiles under cold stress using RNA-seq. Pre-exposure of 96 hpf zebrafish larvae to cold stress (16°C) for 24 h significantly increased their survival rates under severe cold stress (12°C). RNA-seq generated 272 million raw reads from six sequencing libraries and about 92% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 1,431 up- and 399 down-regulated genes. Gene ontology enrichment analysis of cold-induced genes revealed that RNA splicing, ribosome biogenesis and protein catabolic process were the most highly overrepresented biological processes. Spliceosome, proteasome, eukaryotic ribosome biogenesis and RNA transport were the most highly enriched pathways for genes up-regulated by cold stress. Moreover, alternative splicing of 197 genes and promoter switching of 64 genes were found to be regulated by cold stress. A shorter isoform of stk16 that lacks 67 amino acids at the N-terminus was specifically generated by skipping the second exon in cold-treated larvae. Alternative promoter usage was detected for per3 gene under cold stress, which leading to a highly up-regulated transcript encoding a truncated protein lacking the C-terminal domains. These findings indicate that zebrafish larvae possess the ability to build cold-tolerance under mild low temperature and transcriptional and post-transcriptional regulations are extensively involved in this acclimation process.

Abstract: SHORT syndrome is a rare, multisystem disease characterized by short stature, anterior-chamber eye anomalies, characteristic facial features, lipodystrophy, hernias, hyperextensibility, and delayed dentition. As part of the FORGE (Finding of Rare Disease Genes) Canada Consortium, we studied individuals with clinical features of SHORT syndrome to identify the genetic etiology of this rare disease. Whole-exome sequencing in a family trio of an affected child and unaffected parents identified a de novo frameshift insertion, c.1906_1907insC (p.Asn636Thrfs∗18), in exon 14 of PIK3R1. Heterozygous mutations in exon 14 of PIK3R1 were subsequently identified by Sanger sequencing in three additional affected individuals and two affected family members. One of these mutations, c.1945C>T (p.Arg649Trp), was confirmed to be a de novo mutation in one affected individual and was also identified and shown to segregate with the phenotype in an unrelated family. The other mutation, a de novo truncating mutation (c.1971T>G [p.Tyr657∗]), was identified in another affected individual. PIK3R1 is involved in the phosphatidylinositol 3 kinase (PI3K) signaling cascade and, as such, plays an important role in cell growth, proliferation, and survival. Functional studies on lymphoblastoid cells with the PIK3R1 c.1906_1907insC mutation showed decreased phosphorylation of the downstream S6 target of the PI3K-AKT-mTOR pathway. Our findings show that PIK3R1 mutations are the major cause of SHORT syndrome and suggest that the molecular mechanism of disease might involve downregulation of the PI3K-AKT-mTOR pathway.

Abstract: The RNA transcriptome varies in response to cellular differentiation as well as environmental factors, and can be characterized by the diversity and abundance of transcript isoforms. Differential transcription analysis, the detection of differences between the transcriptomes of different cells, may improve understanding of cell differentiation and development and enable the identification of biomarkers that classify disease types. The availability of high-throughput short-read RNA sequencing technologies provides in-depth sampling of the transcriptome, making it possible to accurately detect the differences between transcriptomes. In this article, we present a new method for the detection and visualization of differential transcription. Our approach does not depend on transcript or gene annotations. It also circumvents the need for full transcript inference and quantification, which is a challenging problem because of short read lengths, as well as various sampling biases. Instead, our method takes a divide-and-conquer approach to localize the difference between transcriptomes in the form of alternative splicing modules (ASMs), where transcript isoforms diverge. Our approach starts with the identification of ASMs from the splice graph, constructed directly from the exons and introns predicted from RNA-seq read alignments. The abundance of alternative splicing isoforms residing in each ASM is estimated for each sample and is compared across sample groups. A non-parametric statistical test is applied to each ASM to detect significant differential transcription with a controlled false discovery rate. The sensitivity and specificity of the method have been assessed using simulated data sets and compared with other state-of-the-art approaches. Experimental validation using qRT-PCR confirmed a selected set of genes that are differentially expressed in a lung differentiation study and a breast cancer data set, demonstrating the utility of the approach applied on experimental biological data sets. The software of DiffSplice is available at http://www.netlab.uky.edu/p/bioinfo/DiffSplice.

Abstract: Most human pre-mRNA transcripts are alternatively spliced, but the significance and fine-tuning of alternative splicing in different biological processes is only starting to be understood. SRSF3 (SRp20) is a member of a highly conserved family of splicing factors that have critical roles in key biological processes, including tumor progression. Here, we show that SRSF3 regulates cellular senescence, a p53-mediated process to suppress tumorigenesis, through TP53 alternative splicing. Downregulation of SRSF3 was observed in normal human fibroblasts undergoing replicative senescence, and was associated with the upregulation of p53β, an alternatively spliced isoform of p53 that promotes p53-mediated senescence. Knockdown of SRSF3 by short interfering RNA (siRNA) in early-passage fibroblasts induced senescence, which was associated with elevated expression of p53β at mRNA and protein levels. Knockdown of p53 partially rescued SRSF3-knockdown-induced senescence, suggesting that SRSF3 acts on p53-mediated cellular senescence. RNA pulldown assays demonstrated that SRSF3 binds to an alternatively spliced exon uniquely included in p53β mRNA through the consensus SRSF3-binding sequences. RNA crosslinking and immunoprecipitation assays (CLIP) also showed that SRSF3 in vivo binds to endogenous p53 pre-mRNA at the region containing the p53β-unique exon. Splicing assays using a transfected TP53 minigene in combination with siRNA knockdown of SRSF3 showed that SRSF3 functions to inhibit the inclusion of the p53β-unique exon in splicing of p53 pre-mRNA. These data suggest that downregulation of SRSF3 represents an endogenous mechanism for cellular senescence that directly regulates the TP53 alternative splicing to generate p53β. This study uncovers the role for general splicing machinery in tumorigenesis, and suggests that SRSF3 is a direct regulator of p53.

Abstract: Lewy body diseases include Parkinson disease and dementia with Lewy bodies and are characterized by the widespread distribution of Lewy bodies in virtually every brain area. The main component of Lewy bodies is alpha-synuclein (AS). Accumulating evidence suggests that AS oligomerization and aggregation are strongly associated with the pathogenesis of Lewy body diseases. AS is a small soluble protein with aggregation-prone properties under certain conditions. These properties are enhanced by posttranslational modifications such as phosphorylation, ubiquitination, nitration, and truncation. Accordingly, Lewy bodies contain abundant phosphorylated, nitrated, and monoubiquitinated AS. However, alternative splicing of the AS gene is also known to modify AS aggregation propensities. Splicing gives rise to four related forms of the protein, the main transcript and those that lack exon 4, exon 6, or both. Since AS structure and properties have been extensively studied, it is possible to predict the consequences of the splicing out of the two aforesaid exons. The present review discusses the latest insights on the mechanisms of AS posttranslational modifications and intends to depict their role in the pathogenesis of Lewy body diseases. The implications of deregulated alternative splicing are examined as well, and a hypothesis for the development of the pure form of dementia with Lewy bodies is proposed.

Abstract: DICER1 is an endoribonuclease central to the generation of microRNAs (miRNAs) and short interfering RNAs (siRNAs). Germline mutations in DICER1 have been associated with a pleiotropic tumour predisposition syndrome and Wilms tumour (WT) is a rare manifestation of this syndrome. Three WTs, each in a child with a deleterious germline DICER1 mutation, were screened for somatic DICER1 mutations and were found to bear specific mutations in either the RNase IIIa ( n= 1) or the RNase IIIb domain ( n= 2). In the two latter cases, we demonstrate that the germline and somatic DICER1 mutations were in trans, suggesting that the two-hit hypothesis of tumour formation applies for these examples of WT. Among 191 apparently sporadic WTs, we identified five different missense or deletion somatic DICER1 mutations (2.6%) in four individual WTs; one tumour had two very likely deleterious somatic mutations in trans in the RNase IIIb domain (c.5438A>G and c.5452G>A). In vitro studies of two somatic single-base substitutions (c.5429A>G and c.5438A>G) demonstrated exon 25 skipping from the transcript, a phenomenon not previously reported in DICER1. Further we show that DICER1 transcripts lacking exon 25 can be translated in vitro. This study has demonstrated that a subset of WTs exhibits two ‘hits’ in DICER1, suggesting that these mutations could be key events in the pathogenesis of these tumours. Copyright  2013 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Abstract: Endoplasmic reticulum aminopeptidase 1 (ERAP1) trims peptides for MHC class I presentation, influencing the degree and specificity of CD8+ T cell responses. Single-nucleotide polymorphisms within the exons encoding ERAP1 are associated with autoimmune diseases and cervical carcinoma, but it is not known whether they act independently or as disease-associated haplotypes. We sequenced ERAP1 from 20 individuals and show that single-nucleotide polymorphisms occur as distinct haplotypes in the human population and that these haplotypes encode functionally distinct ERAP1 alleles. Using a wide range of substrates, we are able to demonstrate that for any given substrate distinct ERAP1 alleles can be “normal,” “hypofunctional,” or “hyperfunctional” and that each allele has a trend bias toward one of these three activities. Thus, the repertoire of peptides presented at the cell surface for recognition by CTL is likely to depend on the precise combination of both MHC class I and ERAP1 alleles expressed within an individual, and has important implications for predisposition to disease.