Translation of Circular RNAs: Functions of Translated Products and Related Bioinformatics Approaches.

Over the past two decades, studies have discovered a special form of alternative splicing (AS) that produces a circular form of RNA. This stands in contrast to normal AS, which produces a linear form of RNA. Although these circRNAs have garnered considerable attention in the scientific community for their biogenesis and functions, the focus of these studies has been on the regulatory role of circRNAs with the assumption that circRNAs are non-coding. As non-coding RNAs, they may regulate mRNA transcription, tumor initiation, and translation by sponging miRNAs and RNA-binding proteins (RBPs). In addition to these regulatory roles of circRNAs, however, recent studies have provided strong evidence for their translation. The translation of circRNAs is expected to have an important role in promoting cancer cell growth and activating molecular pathways related to cancer development. In some cases, the translation of circRNAs is shown to be efficiently driven by an internal ribosome entry site (IRES). The development of a computational tool for identifying and characterizing the translation of circRNAs using high-throughput sequencing and IRES increases identifiable proteins translated from circRNAs. In turn, it has a substantial impact on helping researchers understand the functional role of proteins derived from circRNAs. New web resources for aggregating, cataloging, and visualizing translational information of circRNAs derived from previous studies have been developed. In this paper, general concepts of circRNA, circRNA biogenesis, translation of circRNA, and existing circRNA tools and databases are summarized to provide new insight into circRNA studies.

An extracellular vesicular mutant KRAS-associated protein complex promotes lung inflammation and tumor growth.

Extracellular vesicles (EVs) contain more than 100 proteins. Whether there are EVs proteins that act as an 'organiser' of protein networks to generate a new or different biological effect from that identified in EV-producing cells has never been demonstrated. Here, as a proof-of-concept, we demonstrate that EV-G12D-mutant KRAS serves as a leader that forms a protein complex and promotes lung inflammation and tumour growth via the Fn1/IL-17A/FGF21 axis. Mechanistically, in contrast to cytosol derived G12D-mutant KRAS complex from EVs-producing cells, EV-G12D-mutant KRAS interacts with a group of extracellular vesicular factors via fibronectin-1 (Fn1), which drives the activation of the IL-17A/FGF21 inflammation pathway in EV recipient cells. We show that: (i), depletion of EV-Fn1 leads to a reduction of a number of inflammatory cytokines including IL-17A; (ii) induction of IL-17A promotes lung inflammation, which in turn leads to IL-17A mediated induction of FGF21 in the lung; and (iii) EV-G12D-mutant KRAS complex mediated lung inflammation is abrogated in IL-17 receptor KO mice. These findings establish a new concept in EV function with potential implications for novel therapeutic interventions in EV-mediated disease processes.

Gut bacterial isoamylamine promotes age-related cognitive dysfunction by promoting microglial cell death.

The intestinal microbiome releases a plethora of small molecules. Here, we show that the Ruminococcaceae metabolite isoamylamine (IAA) is enriched in aged mice and elderly people, whereas Ruminococcaceae phages, belonging to the Myoviridae family, are reduced. Young mice orally administered IAA show cognitive decline, whereas Myoviridae phage administration reduces IAA levels. Mechanistically, IAA promotes apoptosis of microglial cells by recruiting the transcriptional regulator p53 to the S100A8 promoter region. Specifically, IAA recognizes and binds the S100A8 promoter region to facilitate the unwinding of its self-complementary hairpin structure, thereby subsequently enabling p53 to access the S100A8 promoter and enhance S100A8 expression. Thus, our findings provide evidence that small molecules released from the gut microbiome can directly bind genomic DNA and act as transcriptional coregulators by recruiting transcription factors. These findings further unveil a molecular mechanism that connects gut metabolism to gene expression in the brain with implications for disease development.

Ginger nanoparticles mediated induction of Foxa2 prevents high-fat diet-induced insulin resistance.

Rationale: The obesity epidemic has expanded globally, due in large part to the increased consumption of high-fat diets (HFD), and has increased the risk of major chronic diseases, including type 2 diabetes. Diet manipulation is the foundation of prevention and treatment of obesity and diabetes. The molecular mechanisms that mediate the diet-based prevention of insulin resistance, however, remain to be identified. Here, we report that treatment with orally administered ginger-derived nanoparticles (GDNP) prevents insulin resistance by restoring homeostasis in gut epithelial Foxa2 mediated signaling in mice fed a high-fat diet (HFD). Methods: Ginger-derived nanoparticles (GDNP) were added into drinking water to treat high-fat diet fed mice for at least one year or throughout their life span. A micro array profile of intestinal, liver and fat tissue of GDNP treated mice was used to analyze their gene expression profile. Genes associated with metabolism or insulin signaling were further quantified using the real time polymerase chain reaction (RT-PCR). Surface plasmon resonance (SPR) was used for determining the interaction between Foxa2 protein and phosphatic acid lipid nanoparticles. Results: HFD-feeding inhibited the expression of Foxa2; the GDNPs increased the expression of Foxa2 and protected Foxa2 against Akt-1 mediated phosphorylation and subsequent inactivation of Foxa2. Increasing expression of Foxa2 leads to altering the composition of intestinal epithelial cell (IEC) exosomes of mice fed a HFD and prevents IEC exosome mediated insulin resistance. Collectively, oral administration of GDNP prevents insulin resistance in HFD mice. Interestingly, oral administration of GDNP also extended the life span of the mice and inhibited skin inflammation. Conclusion: Our findings showed that GDNP treatment can prevent HFD-induced obesity and insulin resistance via protecting the Foxa2 from Akt-1 mediated phosphorylation. GDNP treatment provides an alternative approach based on diet manipulation for the development of therapeutic interventions for obesity.

Plant-derived exosomal microRNAs inhibit lung inflammation induced by exosomes SARS-CoV-2 Nsp12.

Lung inflammation is a hallmark of coronavirus disease 2019 (COVID-19). In this study, we show that mice develop inflamed lung tissue after being administered exosomes released from the lung epithelial cells exposed to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Nsp12 and Nsp13 (exosomesNsp12Nsp13). Mechanistically, we show that exosomesNsp12Nsp13 are taken up by lung macrophages, leading to activation of nuclear factor κB (NF-κB) and the subsequent induction of an array of inflammatory cytokines. Induction of tumor necrosis factor (TNF)-α, interleukin (IL)-6, and IL-1ß from exosomesNsp12Nsp13-activated lung macrophages contributes to inducing apoptosis in lung epithelial cells. Induction of exosomesNsp12Nsp13-mediated lung inflammation was abolished with ginger exosome-like nanoparticle (GELN) microRNA (miRNA aly-miR396a-5p. The role of GELNs in inhibition of the SARS-CoV-2-induced cytopathic effect (CPE) was further demonstrated via GELN aly-miR396a-5p- and rlcv-miR-rL1-28-3p-mediated inhibition of expression of Nsp12 and spike genes, respectively. Taken together, our results reveal exosomesNsp12Nsp13 as potentially important contributors to the development of lung inflammation, and GELNs are a potential therapeutic agent to treat COVID-19.

High-fat diet-induced upregulation of exosomal phosphatidylcholine contributes to insulin resistance.

High-fat diet (HFD) decreases insulin sensitivity. How high-fat diet causes insulin resistance is largely unknown. Here, we show that lean mice become insulin resistant after being administered exosomes isolated from the feces of obese mice fed a HFD or from patients with type II diabetes. HFD altered the lipid composition of exosomes from predominantly phosphatidylethanolamine (PE) in exosomes from lean animals (L-Exo) to phosphatidylcholine (PC) in exosomes from obese animals (H-Exo). Mechanistically, we show that intestinal H-Exo is taken up by macrophages and hepatocytes, leading to inhibition of the insulin signaling pathway. Moreover, exosome-derived PC binds to and activates AhR, leading to inhibition of the expression of genes essential for activation of the insulin signaling pathway, including IRS-2, and its downstream genes PI3K and Akt. Together, our results reveal HFD-induced exosomes as potential contributors to the development of insulin resistance. Intestinal exosomes thus have potential as broad therapeutic targets.

seekCRIT: Detecting and characterizing differentially expressed circular RNAs using high-throughput sequencing data.

Over the past two decades, researchers have discovered a special form of alternative splicing that produces a circular form of RNA. Although these circular RNAs (circRNAs) have garnered considerable attention in the scientific community for their biogenesis and functions, the focus of current studies has been on the tissue-specific circRNAs that exist only in one tissue but not in other tissues or on the disease-specific circRNAs that exist in certain disease conditions, such as cancer, but not under normal conditions. This approach was conducted in the relative absence of methods that analyze a group of common circRNAs that exist in both conditions, but are more abundant in one condition relative to another (differentially expressed). Studies of differentially expressed circRNAs (DECs) between two conditions would serve as a significant first step in filling this void. Here, we introduce a novel computational tool, seekCRIT (seek for differentially expressed CircRNAs In Transcriptome), that identifies the DECs between two conditions from high-throughput sequencing data. Using rat retina RNA-seq data from ischemic and normal conditions, we show that over 74% of identifiable circRNAs are expressed in both conditions and over 40 circRNAs are differentially expressed between two conditions. We also obtain a high qPCR validation rate of 90% for DECs with a FDR of < 5%. Our results demonstrate that seekCRIT is a novel and efficient approach to detect DECs using rRNA depleted RNA-seq data. seekCRIT is freely downloadable at https://github.com/UofLBioinformatics/seekCRIT. The source code is licensed under the MIT License. seekCRIT is developed and tested on Linux CentOS-7.

circDeep: deep learning approach for circular RNA classification from other long non-coding RNA.

MOTIVATION: Over the past two decades, a circular form of RNA (circular RNA), produced through alternative splicing, has become the focus of scientific studies due to its major role as a microRNA (miRNA) activity modulator and its association with various diseases including cancer. Therefore, the detection of circular RNAs is vital to understanding their biogenesis and purpose. Prediction of circular RNA can be achieved in three steps: distinguishing non-coding RNAs from protein coding gene transcripts, separating short and long non-coding RNAs and predicting circular RNAs from other long non-coding RNAs (lncRNAs). However, the available tools are less than 80 percent accurate for distinguishing circular RNAs from other lncRNAs due to difficulty of classification. Therefore, the availability of a more accurate and fast machine learning method for the identification of circular RNAs, which considers the specific features of circular RNA, is essential to the development of systematic annotation. RESULTS: Here we present an End-to-End deep learning framework, circDeep, to classify circular RNA from other lncRNA. circDeep fuses an RCM descriptor, ACNN-BLSTM sequence descriptor and a conservation descriptor into high level abstraction descriptors, where the shared representations across different modalities are integrated. The experiments show that circDeep is not only faster than existing tools but also performs at an unprecedented level of accuracy by achieving a 12 percent increase in accuracy over the other tools. AVAILABILITY AND IMPLEMENTATION: https://github.com/UofLBioinformatics/circDeep. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

AF1q inhibited T cell attachment to breast cancer cell by attenuating Intracellular Adhesion Molecule-1 expression.

AIM: To investigate whether AF1q, overexpressed in metastatic cells compared with the primary tumor cells, plays a pivotal role in breast cancer metastasis. METHODS: To investigate whether AF1q has a responsibility in the acquisition of a metastatic phenotype, we performed RNA-sequencing (RNA-Seq) to identify the gene signature and applied the Metacore direct interactions network building algorithm with the top 40 amplicons of RNA-Seq. RESULTS: Most genes were directly linked with intercellular adhesion molecule-1 (ICAM-1). Likewise, we identified that ICAM-1 expression is attenuated in metastatic cells compared to primary tumor cells. Moreover, overexpression of AF1q attenuated ICAM-1 expression, whereas suppression of AF1q elicited the opposite effect. AF1q had an effect on ICAM-1 promoter region and regulated its transcription. Decreased ICAM-1 expression affected the attachment of T cells to a breast cancer cell monolayer. We confirmed the finding by performing the analysis on Burkitt's lymphoma. CONCLUSION: Attenuation of ICAM-1 by AF1q on tumor cells disadvantages host anti-tumor defenses through the trafficking of lymphocytes, which affects tumor progression and metastasis.

Concerted effects of heterogeneous nuclear ribonucleoprotein C1/C2 to control vitamin D-directed gene transcription and RNA splicing in human bone cells.

Traditionally recognized as an RNA splicing regulator, heterogeneous nuclear ribonucleoprotein C1/C2 (hnRNPC1/C2) can also bind to double-stranded DNA and function in trans as a vitamin D response element (VDRE)-binding protein. As such, hnRNPC1/C2 may couple transcription induced by the active form of vitamin D, 1,25-dihydroxyvitamin D (1,25(OH)2D) with subsequent RNA splicing. In MG63 osteoblastic cells, increased expression of the 1,25(OH)2D target gene CYP24A1 involved immunoprecipitation of hnRNPC1/C2 with CYP24A1 chromatin and RNA. Knockdown of hnRNPC1/C2 suppressed expression of CYP24A1, but also increased expression of an exon 10-skipped CYP24A1 splice variant; in a minigene model the latter was attenuated by a functional VDRE in the CYP24A1 promoter. In genome-wide analyses, knockdown of hnRNPC1/C2 resulted in 3500 differentially expressed genes and 2232 differentially spliced genes, with significant commonality between groups. 1,25(OH)2D induced 324 differentially expressed genes, with 187 also observed following hnRNPC1/C2 knockdown, and a further 168 unique to hnRNPC1/C2 knockdown. However, 1,25(OH)2D induced only 10 differentially spliced genes, with no overlap with differentially expressed genes. These data indicate that hnRNPC1/C2 binds to both DNA and RNA and influences both gene expression and RNA splicing, but these actions do not appear to be linked through 1,25(OH)2D-mediated induction of transcription.

Ablation of the epithelial-specific splicing factor Esrp1 results in ureteric branching defects and reduced nephron number.

BACKGROUND: Abnormalities in ureteric bud (UB) branching morphogenesis lead to congenital anomalies of the kidney and reduced nephron numbers associated with chronic kidney disease (CKD) and hypertension. Previous studies showed that the epithelial fibroblast growth factor receptor 2 (Fgfr2) IIIb splice variant supports ureteric morphogenesis in response to ligands from the metanephric mesenchyme during renal organogenesis. The epithelial-specific splicing regulator Esrp1 is required for expression of Fgfr2-IIIb and other epithelial-specific splice variants. Our objective was to determine whether Esrp1 is required for normal kidney development. RESULTS: Ablation of Esrp1 in mice, alone or together with its paralog Esrp2, was associated with reduced kidney size and increased incidence of renal aplasia. Three-dimensional imaging showed that embryonic Esrp1 knockout (KO) kidneys had fewer ureteric tips and reduced nephron numbers. Analysis of alternative splicing in Esrp-null ureteric epithelial cells by RNA-Seq confirmed a splicing switch in Fgfr2 as well as numerous other transcripts. CONCLUSIONS: Our findings reveal that Esrp1-regulated splicing in ureteric epithelial cells plays an important role in renal development. Defects in Esrp1 KO kidneys likely reflect reduced and/or absent ureteric branching, leading to decreased nephron induction secondary to incorrect Fgfr2 splicing and other splicing alterations. Developmental Dynamics 245:991-1000, 2016. © 2016 The Authors. Developmental Dynamics published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.

Transcriptome sequencing reveals aberrant alternative splicing in Huntington's disease.

Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by a CAG expansion in the gene-encoding Huntingtin (HTT). Transcriptome dysregulation is a major feature of HD pathogenesis, as revealed by a large body of work on gene expression profiling of tissues from human HD patients and mouse models. These studies were primarily focused on transcriptional changes affecting steady-state overall gene expression levels using microarray based approaches. A major missing component, however, has been the study of transcriptome changes at the post-transcriptional level, such as alternative splicing. Alternative splicing is a critical mechanism for expanding regulatory and functional diversity from a limited number of genes, and is particularly complex in the mammalian brain. Here we carried out a deep RNA-seq analysis of the BA4 (Brodmann area 4) motor cortex from seven human HD brains and seven controls to systematically discover aberrant alternative splicing events and characterize potential associated splicing factors in HD. We identified 593 differential alternative splicing events between HD and control brains. Using two expanded panels with a total of 108 BA4 tissues from patients and controls, we identified four splicing factors exhibiting significantly altered expression levels in HD patient brains. Moreover, follow-up molecular analyses of one splicing factor PTBP1 revealed its impact on disease-associated splicing patterns in HD. Collectively, our data provide genomic evidence for widespread splicing dysregulation in HD brains, and suggest the role of aberrant alternative splicing in the pathogenesis of HD.

Determination of a Comprehensive Alternative Splicing Regulatory Network and Combinatorial Regulation by Key Factors during the Epithelial-to-Mesenchymal Transition.

The epithelial-to-mesenchymal transition (EMT) is an essential biological process during embryonic development that is also implicated in cancer metastasis. While the transcriptional regulation of EMT has been well studied, the role of alternative splicing (AS) regulation in EMT remains relatively uncharacterized. We previously showed that the epithelial cell-type-specific proteins epithelial splicing regulatory proteins 1 (ESRP1) and ESRP2 are important for the regulation of many AS events that are altered during EMT. However, the contributions of the ESRPs and other splicing regulators to the AS regulatory network in EMT require further investigation. Here, we used a robust in vitro EMT model to comprehensively characterize splicing switches during EMT in a temporal manner. These investigations revealed that the ESRPs are the major regulators of some but not all AS events during EMT. We determined that the splicing factor RBM47 is downregulated during EMT and also regulates numerous transcripts that switch splicing during EMT. We also determined that Quaking (QKI) broadly promotes mesenchymal splicing patterns. Our study highlights the broad role of posttranscriptional regulation during the EMT and the important role of combinatorial regulation by different splicing factors to fine tune gene expression programs during these physiological and developmental transitions.

αCP binding to a cytosine-rich subset of polypyrimidine tracts drives a novel pathway of cassette exon splicing in the mammalian transcriptome.

Alternative splicing (AS) is a robust generator of mammalian transcriptome complexity. Splice site specification is controlled by interactions of cis-acting determinants on a transcript with specific RNA binding proteins. These interactions are frequently localized to the intronic U-rich polypyrimidine tracts (PPT) located 5' to the majority of splice acceptor junctions. αCPs (also referred to as polyC-binding proteins (PCBPs) and hnRNPEs) comprise a subset of KH-domain proteins with high affinity and specificity for C-rich polypyrimidine motifs. Here, we demonstrate that αCPs promote the splicing of a defined subset of cassette exons via binding to a C-rich subset of polypyrimidine tracts located 5' to the αCP-enhanced exonic segments. This enhancement of splice acceptor activity is linked to interactions of αCPs with the U2 snRNP complex and may be mediated by cooperative interactions with the canonical polypyrimidine tract binding protein, U2AF65. Analysis of αCP-targeted exons predicts a substantial impact on fundamental cell functions. These findings lead us to conclude that the αCPs play a direct and global role in modulating the splicing activity and inclusion of an array of cassette exons, thus driving a novel pathway of splice site regulation within the mammalian transcriptome.

The contribution of Alu exons to the human proteome.

BACKGROUND: Alu elements are major contributors to lineage-specific new exons in primate and human genomes. Recent studies indicate that some Alu exons have high transcript inclusion levels or tissue-specific splicing profiles, and may play important regulatory roles in modulating mRNA degradation or translational efficiency. However, the contribution of Alu exons to the human proteome remains unclear and controversial. The prevailing view is that exons derived from young repetitive elements, such as Alu elements, are restricted to regulatory functions and have not had adequate evolutionary time to be incorporated into stable, functional proteins. RESULTS: We adopt a proteotranscriptomics approach to systematically assess the contribution of Alu exons to the human proteome. Using RNA sequencing, ribosome profiling, and proteomics data from human tissues and cell lines, we provide evidence for the translational activities of Alu exons and the presence of Alu exon derived peptides in human proteins. These Alu exon peptides represent species-specific protein differences between primates and other mammals, and in certain instances between humans and closely related primates. In the case of the RNA editing enzyme ADARB1, which contains an Alu exon peptide in its catalytic domain, RNA sequencing analyses of A-to-I editing demonstrate that both the Alu exon skipping and inclusion isoforms encode active enzymes. The Alu exon derived peptide may fine tune the overall editing activity and, in limited cases, the site selectivity of ADARB1 protein products. CONCLUSIONS: Our data indicate that Alu elements have contributed to the acquisition of novel protein sequences during primate and human evolution.

The splicing regulators Esrp1 and Esrp2 direct an epithelial splicing program essential for mammalian development.

Tissue- and cell-type-specific regulators of alternative splicing (AS) are essential components of posttranscriptional gene regulation, necessary for normal cellular function, patterning, and development. Mice with ablation of Epithelial splicing regulatory protein (Esrp1) develop cleft lip and palate. Loss of both Esrp1 and its paralog Esrp2 results in widespread developmental defects with broad implications to human disease. Deletion of the Esrps in the epidermis revealed their requirement for establishing a proper skin barrier, a primary function of epithelial cells comprising the epidermis. We profiled the global Esrp-mediated splicing regulatory program in epidermis, which revealed large-scale programs of epithelial cell-type-specific splicing required for epithelial cell functions. These mice represent a valuable model for evaluating the essential role for AS in development and function of epithelial cells, which play essential roles in tissue homeostasis in numerous organs, and provide a genetic tool to evaluate important functional properties of epithelial-specific splice variants in vivo.

Transcriptome-wide landscape of pre-mRNA alternative splicing associated with metastatic colonization.

UNLABELLED: Metastatic colonization is an ominous feature of cancer progression. Recent studies have established the importance of pre-mRNA alternative splicing (AS) in cancer biology. However, little is known about the transcriptome-wide landscape of AS associated with metastatic colonization. Both in vitro and in vivo models of metastatic colonization were utilized to study AS regulation associated with cancer metastasis. Transcriptome profiling of prostate cancer cells and derivatives crossing in vitro or in vivo barriers of metastasis revealed splicing factors with significant gene expression changes associated with metastatic colonization. These include splicing factors known to be differentially regulated in epithelial-mesenchymal transition (ESRP1, ESRP2, and RBFOX2), a cellular process critical for cancer metastasis, as well as novel findings (NOVA1 and MBNL3). Finally, RNA-seq indicated a large network of AS events regulated by multiple splicing factors with altered gene expression or protein activity. These AS events are enriched for pathways important for cell motility and signaling, and affect key regulators of the invasive phenotype such as CD44 and GRHL1. IMPLICATIONS: Transcriptome-wide remodeling of AS is an integral regulatory process underlying metastatic colonization, and AS events affect the metastatic behavior of cancer cells.

rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data.

Ultra-deep RNA sequencing (RNA-Seq) has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We previously developed multivariate analysis of transcript splicing (MATS), a statistical method for detecting differential alternative splicing between two RNA-Seq samples. Here we describe a new statistical model and computer program, replicate MATS (rMATS), designed for detection of differential alternative splicing from replicate RNA-Seq data. rMATS uses a hierarchical model to simultaneously account for sampling uncertainty in individual replicates and variability among replicates. In addition to the analysis of unpaired replicates, rMATS also includes a model specifically designed for paired replicates between sample groups. The hypothesis-testing framework of rMATS is flexible and can assess the statistical significance over any user-defined magnitude of splicing change. The performance of rMATS is evaluated by the analysis of simulated and real RNA-Seq data. rMATS outperformed two existing methods for replicate RNA-Seq data in all simulation settings, and RT-PCR yielded a high validation rate (94%) in an RNA-Seq dataset of prostate cancer cell lines. Our data also provide guiding principles for designing RNA-Seq studies of alternative splicing. We demonstrate that it is essential to incorporate biological replicates in the study design. Of note, pooling RNAs or merging RNA-Seq data from multiple replicates is not an effective approach to account for variability, and the result is particularly sensitive to outliers. The rMATS source code is freely available at rnaseq-mats.sourceforge.net/. As the popularity of RNA-Seq continues to grow, we expect rMATS will be useful for studies of alternative splicing in diverse RNA-Seq projects.

BS69/ZMYND11 reads and connects histone H3.3 lysine 36 trimethylation-decorated chromatin to regulated pre-mRNA processing.

BS69 (also called ZMYND11) contains tandemly arranged PHD, BROMO, and PWWP domains, which are chromatin recognition modalities. Here, we show that BS69 selectively recognizes histone variant H3.3 lysine 36 trimethylation (H3.3K36me3) via its chromatin-binding domains. We further identify BS69 association with RNA splicing regulators, including the U5 snRNP components of the spliceosome, such as EFTUD2. Remarkably, RNA sequencing shows that BS69 mainly regulates intron retention (IR), which is the least understood RNA alternative splicing event in mammalian cells. Biochemical and genetic experiments demonstrate that BS69 promotes IR by antagonizing EFTUD2 through physical interactions. We further show that regulation of IR by BS69 also depends on its binding to H3K36me3-decorated chromatin. Taken together, our study identifies an H3.3K36me3-specific reader and a regulator of IR and reveals that BS69 connects histone H3.3K36me3 to regulated RNA splicing, providing significant, important insights into chromatin regulation of pre-mRNA processing.

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data.

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT-PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT-PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.