ABSTRACT
The cochlea's ability to discriminate sound frequencies is facilitated by a special topography along its longitudinal axis known as tonotopy. Auditory hair cells located at the base of the cochlea respond to high-frequency sounds, whereas hair cells at the apex respond to lower frequencies. Gradual changes in morphological and physiological features along the length of the cochlea determine each region's frequency selectivity, but it remains unclear how tonotopy is established during cochlear development. Recently, sonic hedgehog (SHH) was proposed to initiate the establishment of tonotopy by conferring regional identity to the primordial cochlea. Here, using mouse genetics, we provide in vivo evidence that regional identity in the embryonic cochlea acts as a framework upon which tonotopy-specific properties essential for frequency selectivity in the mature cochlea develop. We found that follistatin (FST) is required for the maintenance of apical cochlear identity, but dispensable for its initial induction. In a fate-mapping analysis, we found that FST promotes expansion of apical cochlear cells, contributing to the formation of the apical cochlear domain. SHH, in contrast, is required both for the induction and maintenance of apical identity. In the absence of FST or SHH, mice produce a short cochlea lacking its apical domain. This results in the loss of apex-specific anatomical and molecular properties and low-frequency-specific hearing loss.
Subject(s)
Follistatin , Hedgehog Proteins , Animals , Mice , Follistatin/genetics , Hedgehog Proteins/genetics , Hedgehog Proteins/metabolism , Cochlea/physiology , Hearing/physiology , Mammals/metabolismABSTRACT
Gut microbiota function has numerous effects on humans and the diet humans consume has emerged as a pivotal determinant of gut microbiota function. Here, a new concept that gut microbiota can be trained by diet-derived exosome-like nanoparticles (ELNs) to release healthy outer membrane vesicles (OMVs) is introduced. Specifically, OMVs released from garlic ELN (GaELNs) trained human gut Akkermansia muciniphila (A. muciniphila) can reverse high-fat diet-induced type 2 diabetes (T2DM) in mice. Oral administration of OMVs released from GaELNs trained A. muciniphila can traffick to the brain where they are taken up by microglial cells, resulting in inhibition of high-fat diet-induced brain inflammation. GaELNs treatment increases the levels of OMV Amuc-1100, P9, and phosphatidylcholines. Increasing the levels of Amuc-1100 and P9 leads to increasing the GLP-1 plasma level. Increasing the levels of phosphatidylcholines is required for inhibition of cGas and STING-mediated inflammation and GLP-1R crosstalk with the insulin pathway that leads to increasing expression of Insulin Receptor Substrate (IRS1 and IRS2) on OMV targeted cells. These findings reveal a molecular mechanism whereby OMVs from plant nanoparticle-trained gut bacteria regulate genes expressed in the brain, and have implications for the treatment of brain dysfunction caused by a metabolic syndrome.
Subject(s)
Brain-Gut Axis , Diabetes Mellitus, Type 2 , Exosomes , Garlic , Gastrointestinal Microbiome , Nanoparticles , Diabetes Mellitus, Type 2/metabolism , Garlic/chemistry , Animals , Nanoparticles/chemistry , Exosomes/metabolism , Mice , Akkermansia , Humans , Male , Diet, High-Fat , Mice, Inbred C57BL , Brain/metabolism , Brain/pathologyABSTRACT
Bark protects the tree against environmental insults. Here, we analyzed whether this defensive strategy could be utilized to broadly enhance protection against colitis. As a proof of concept, we show that exosome-like nanoparticles (MBELNs) derived from edible mulberry bark confer protection against colitis in a mouse model by promoting heat shock protein family A (Hsp70) member 8 (HSPA8)-mediated activation of the AhR signaling pathway. Activation of this pathway in intestinal epithelial cells leads to the induction of COP9 Constitutive Photomorphogenic Homolog Subunit 8 (COPS8). Utilizing a gut epithelium-specific knockout of COPS8, we demonstrate that COPS8 acts downstream of the AhR pathway and is required for the protective effect of MBELNs by inducing an array of anti-microbial peptides. Our results indicate that MBELNs represent an undescribed mode of inter-kingdom communication in the mammalian intestine through an AhR-COPS8-mediated anti-inflammatory pathway. These data suggest that inflammatory pathways in a microbiota-enriched intestinal environment are regulated by COPS8 and that edible plant-derived ELNs may hold the potential as new agents for the prevention and treatment of gut-related inflammatory disease.
Subject(s)
Colitis , Exosomes , Morus , Nanoparticles , Animals , Colitis/chemically induced , Colitis/metabolism , Colitis/prevention & control , Disease Models, Animal , Exosomes/metabolism , Mice , Mice, Inbred C57BL , Plant Bark/metabolismABSTRACT
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.
Subject(s)
COVID-19/metabolism , Exosomes/metabolism , MicroRNAs/metabolism , Plants/metabolism , Pneumonia/metabolism , A549 Cells , Animals , Cell Line , Cell Line, Tumor , Chlorocebus aethiops , Cytokines/metabolism , Epithelial Cells/metabolism , Humans , Interleukin-6/metabolism , Macrophages, Alveolar/metabolism , Male , Mice , Mice, Inbred C57BL , NF-kappa B/metabolism , SARS-CoV-2/pathogenicity , Tumor Necrosis Factor-alpha/metabolism , U937 Cells , Vero CellsABSTRACT
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.
Subject(s)
Alternative Splicing , Carrier Proteins/metabolism , Chromatin/metabolism , Histones/metabolism , RNA Precursors/genetics , RNA, Messenger/genetics , Base Sequence , Carrier Proteins/genetics , Cell Cycle Proteins , Cell Line, Tumor , Chromatin/genetics , Co-Repressor Proteins , DNA Methylation/genetics , DNA-Binding Proteins , HeLa Cells , Histone-Lysine N-Methyltransferase/genetics , Histone-Lysine N-Methyltransferase/metabolism , Histones/genetics , Humans , Introns/genetics , Lysine/genetics , Lysine/metabolism , Peptide Elongation Factors/antagonists & inhibitors , Peptide Elongation Factors/genetics , Peptide Elongation Factors/metabolism , Protein Binding , Protein Structure, Tertiary , RNA Interference , RNA Processing, Post-Transcriptional/genetics , RNA, Small Interfering , Ribonucleoprotein, U5 Small Nuclear/antagonists & inhibitors , Ribonucleoprotein, U5 Small Nuclear/genetics , Ribonucleoprotein, U5 Small Nuclear/metabolism , Sequence Analysis, RNA , Spliceosomes/geneticsABSTRACT
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.
Subject(s)
Computational Biology , Deep Learning , RNA, Circular , RNA, Long Noncoding , Computational Biology/methods , RNA, Circular/classification , RNA, Circular/genetics , RNA, Long Noncoding/genetics , Reproducibility of ResultsABSTRACT
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.
Subject(s)
Gene Expression Profiling/methods , High-Throughput Nucleotide Sequencing/methods , RNA, Circular , Sequence Analysis, RNA/methods , Transcriptome/genetics , Animals , Computational Biology , Databases, Genetic , Humans , RNA, Circular/genetics , RNA, Circular/metabolism , Rats , SoftwareABSTRACT
Superantigens (SAgs) produced by Staphylococcus aureus at high concentrations induce proliferation of T cells bearing specific TCR Vß sequences and massive cytokinemia that cause toxic shock syndrome. However, the biological relevance of SAgs produced at very low concentrations during asymptomatic colonization or chronic infections is not understood. In this study, we demonstrate that suboptimal stimulation of human PBMCs with a low concentration (1 ng/ml) of staphylococcal enterotoxin C1, at which half-maximal T cell proliferation was observed, induced CD8+CD25+ T cells expressing markers related to regulatory T cells (Tregs), such as IFN-γ, IL-10, TGF-ß, FOXP3, CD28, CTLA4, TNFR2, CD45RO, and HLA-DR. Importantly, these CD8+CD25+ T cells suppressed responder cell proliferation mediated in contact-dependent and soluble factor-dependent manners, involving galectin-1 and granzymes, respectively. In contrast, optimal stimulation of human PBMCs with a high concentration (1 µg/ml) of staphylococcal enterotoxin C1, at which maximal T cell proliferation was observed, also induced similar expression of markers related to Tregs, including FOXP3 in CD8+CD25+ cells, but these T cells were not functionally immunosuppressive. We further demonstrated that SAg-induced TCR Vß-restricted and MHC class II-restricted expansion of immunosuppressive CD8+CD25+ T cells is independent of CD4+ T cells. Our results suggest that the concentration of SAg strongly affects the functional characteristics of activated T cells, and low concentrations of SAg produced during asymptomatic colonization or chronic S. aureus infection induce immunosuppressive CD8+ Tregs, potentially promoting colonization, propagation, and invasion of S. aureus in the host.
Subject(s)
CD8-Positive T-Lymphocytes/immunology , Enterotoxins/immunology , Immunomodulation , Staphylococcus aureus/immunology , T-Lymphocyte Subsets/immunology , T-Lymphocytes, Regulatory/immunology , Adolescent , Adult , Biomarkers , CD8-Positive T-Lymphocytes/metabolism , Cytokines/genetics , Cytokines/metabolism , Forkhead Transcription Factors/metabolism , Genes, T-Cell Receptor beta/genetics , Humans , Immunization , Immunophenotyping , Interleukin-2 Receptor alpha Subunit/metabolism , Phenotype , T-Lymphocyte Subsets/metabolism , T-Lymphocytes, Regulatory/metabolism , Young AdultABSTRACT
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.
Subject(s)
Heterogeneous-Nuclear Ribonucleoprotein Group C/metabolism , Osteocytes/metabolism , RNA Splicing , Transcription, Genetic , Vitamin D/metabolism , Alternative Splicing , Cell Line, Tumor , Cluster Analysis , Gene Expression Profiling , Gene Expression Regulation , Gene Knockdown Techniques , Heterogeneous-Nuclear Ribonucleoprotein Group C/genetics , High-Throughput Nucleotide Sequencing , Humans , Osteocytes/drug effects , Promoter Regions, Genetic , Protein Binding , RNA Interference , RNA Precursors , Vitamin D3 24-Hydroxylase/genetics , Vitamin D3 24-Hydroxylase/metabolismABSTRACT
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.
Subject(s)
Heterogeneous-Nuclear Ribonucleoproteins/genetics , Huntington Disease/genetics , Motor Cortex/metabolism , Polypyrimidine Tract-Binding Protein/genetics , Transcriptome/genetics , Adult , Aged , Alternative Splicing/genetics , Animals , Autopsy , Female , Gene Expression Regulation , Genetic Predisposition to Disease , Heterogeneous-Nuclear Ribonucleoproteins/biosynthesis , High-Throughput Nucleotide Sequencing , Humans , Huntingtin Protein/genetics , Huntington Disease/physiopathology , Male , Mice , Middle Aged , Motor Cortex/pathology , Polypyrimidine Tract-Binding Protein/biosynthesisABSTRACT
RNA-binding proteins (RBPs) play a critical role in the regulation of alternative splicing (AS), a prevalent mechanism for generating transcriptomic and proteomic diversity in eukaryotic cells. Studies have shown that AS can be regulated by RBPs in a binding-site-position dependent manner. Depending on where RBPs bind, splicing of an alternative exon can be enhanced or suppressed. Therefore, spatial analyses of RBP motifs and binding sites around alternative exons will help elucidate splicing regulation by RBPs. The development of high-throughput sequencing technologies has allowed transcriptome-wide analyses of AS and RBP-RNA interactions. Given a set of differentially regulated alternative exons obtained from RNA sequencing (RNA-seq) experiments, the rMAPS web server (http://rmaps.cecsresearch.org) performs motif analyses of RBPs in the vicinity of alternatively spliced exons and creates RNA maps that depict the spatial patterns of RBP motifs. Similarly, rMAPS can also perform spatial analyses of RBP-RNA binding sites identified by cross-linking immunoprecipitation sequencing (CLIP-seq) experiments. We anticipate rMAPS will be a useful tool for elucidating RBP regulation of alternative exon splicing using high-throughput sequencing data.
Subject(s)
Exons , RNA-Binding Proteins/chemistry , RNA/chemistry , Transcriptome , User-Computer Interface , Alternative Splicing , Binding Sites , Computer Graphics , Eukaryotic Cells/metabolism , High-Throughput Nucleotide Sequencing , Internet , Introns , Protein Binding , RNA/genetics , RNA/metabolism , RNA-Binding Proteins/genetics , RNA-Binding Proteins/metabolismABSTRACT
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.
Subject(s)
Alternative Splicing , Heterogeneous-Nuclear Ribonucleoproteins/metabolism , Nuclear Proteins/metabolism , Pyrimidines/metabolism , RNA-Binding Proteins/metabolism , Ribonucleoprotein, U2 Small Nuclear/metabolism , Ribonucleoproteins/metabolism , Transcriptome , Adaptor Proteins, Signal Transducing/genetics , Adaptor Proteins, Signal Transducing/metabolism , Base Sequence , Binding Sites , Cyclin-Dependent Kinase 2/genetics , Cyclin-Dependent Kinase 2/metabolism , Cytosine/metabolism , DNA-Binding Proteins , Exons , Gene Expression Regulation , Heterogeneous-Nuclear Ribonucleoproteins/genetics , Humans , Introns , K562 Cells , Molecular Sequence Data , Nuclear Proteins/genetics , Polymers/metabolism , Protein Binding , RNA-Binding Proteins/genetics , Ribonucleoprotein, U2 Small Nuclear/genetics , Ribonucleoproteins/genetics , Sequence Analysis, RNA , Splicing Factor U2AFABSTRACT
RNA-seq has become a popular technology for studying genetic variation of pre-mRNA alternative splicing. Commonly used RNA-seq aligners rely on the consensus splice site dinucleotide motifs to map reads across splice junctions. Consequently, genomic variants that create novel splice site dinucleotides may produce splice junction RNA-seq reads that cannot be mapped to the reference genome. We developed and evaluated an approach to identify 'hidden' splicing variations in personal transcriptomes, by mapping personal RNA-seq data to personal genomes. Computational analysis and experimental validation indicate that this approach identifies personal specific splice junctions at a low false positive rate. Applying this approach to an RNA-seq data set of 75 individuals, we identified 506 personal specific splice junctions, among which 437 were novel splice junctions not documented in current human transcript annotations. 94 splice junctions had splice site SNPs associated with GWAS signals of human traits and diseases. These involve genes whose splicing variations have been implicated in diseases (such as OAS1), as well as novel associations between alternative splicing and diseases (such as ICA1). Collectively, our work demonstrates that the personal genome approach to RNA-seq read alignment enables the discovery of a large but previously unknown catalog of splicing variations in human populations.
Subject(s)
Alternative Splicing , Gene Expression Profiling/methods , Genome, Human , Polymorphism, Single Nucleotide , RNA Splice Sites , Sequence Analysis, RNA/methods , Disease/genetics , Genome-Wide Association Study , Humans , TranscriptomeABSTRACT
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.
Subject(s)
Alternative Splicing , Computer Simulation , High-Throughput Nucleotide Sequencing/methods , RNA/genetics , Cell Line, Tumor , Humans , Male , Prostatic Neoplasms/genetics , Prostatic Neoplasms/pathologyABSTRACT
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.
Subject(s)
RNA-Binding Proteins/metabolism , Ureter/cytology , Ureter/metabolism , Animals , Epithelial Cells/cytology , Epithelial Cells/metabolism , Female , Kidney/cytology , Kidney/embryology , Male , Mice , Mice, Knockout , Nephrons/cytology , Nephrons/metabolism , RNA Splicing/genetics , RNA Splicing/physiology , RNA-Binding Proteins/genetics , Real-Time Polymerase Chain Reaction , Reverse Transcriptase Polymerase Chain Reaction , Signal Transduction/genetics , Signal Transduction/physiologyABSTRACT
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.
Subject(s)
Alternative Splicing , High-Throughput Nucleotide Sequencing , Sequence Analysis, RNA , Bayes Theorem , Brain/metabolism , Cell Line, Tumor , Humans , Multivariate Analysis , RNA-Binding Proteins/metabolism , Reverse Transcriptase Polymerase Chain ReactionABSTRACT
Background: Multiple Sclerosis (MS) is an autoimmune neurodegenerative disease affecting approximately 3 million people globally. Despite rigorous research on MS, aspects of its development and progression remain unclear. Understanding molecular mechanisms underlying MS is crucial to providing insights into disease pathways, identifying potential biomarkers for early diagnosis, and revealing novel therapeutic targets for improved patient outcomes. Methods: We utilized publicly available RNA-seq data (GSE138614) from post-mortem white matter tissues of five donors without any neurological disorder and ten MS patient donors. This data was interrogated for differential gene expression, alternative splicing and single nucleotide variants as well as for functional enrichments in the resulting datasets. Results: A comparison of non-MS white matter (WM) to MS samples yielded differentially expressed genes involved in adaptive immune response, cell communication, and developmental processes. Genes with expression changes positively correlated with tissue inflammation were enriched in the immune system and receptor interaction pathways. Negatively correlated genes were enriched in neurogenesis, nervous system development, and metabolic pathways. Alternatively spliced transcripts between WM and MS lesions included genes that play roles in neurogenesis, myelination, and oligodendrocyte differentiation, such as brain enriched myelin associated protein (BCAS1), discs large MAGUK scaffold protein 1 (DLG1), KH domain containing RNA binding (QKI), and myelin basic protein (MBP). Our approach to comparing normal appearing WM (NAWM) and active lesion (AL) from one donor and NAWM and chronic active (CA) tissues from two donors, showed that different IgH and IgK gene subfamilies were differentially expressed. We also identified pathways involved in white matter injury repair and remyelination in these tissues. Differentially spliced genes between these lesions were involved in axon and dendrite structure stability. We also identified exon skipping events and spontaneous single nucleotide polymorphisms in membrane associated ring-CH-type finger 1 (MARCHF1), UDP glycosyltransferase 8 (UGT8), and other genes important in autoimmunity and neurodegeneration. Conclusion: Overall, we identified unique genes, pathways, and novel splicing events affecting disease progression that can be further investigated as potential novel drug targets for MS treatment.
ABSTRACT
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.
ABSTRACT
Essential amino acids (EAA) and microRNAs (miRs) control biological activity of a cell. Whether EAA regulates the activity of miR has never been demonstrated. Here, as proof-of-concept, a tryptophan (Trp, an EAA) complex containing Argonaute 2 (Ago2) and miRs including miR-193a (Trp/Ago2/miR-193a) is identified. Trp binds miR-193a-3p and interacts with Ago2. Trp/Ago2/miR-193a increases miR-193a-3p activity via enhancing Argonaute 2 (Ago2) RNase activity. Other miRs including miR-103 and miR-107 in the Trp complex enhance miR-193a activity by targeting the same genes. Mechanistically, the Trp/Ago2/miR-193a complex interacts with Trp-binding pockets of the PIWI domain of Ago2 to enhance Ago2 mediated miR activity. This newly formed Ago2/Trp/miR-193a-3p complex is more efficient than miR-193a-3p alone in inhibiting the expression of targeted genes and inhibiting colon cancer liver metastasis. The findings show that Trp regulates miR activity through communication with the RNA-induced silencing complexes (RISC), which provides the basis for tryptophan based miR therapy.
Subject(s)
Argonaute Proteins , Colonic Neoplasms , Liver Neoplasms , MicroRNAs , RNA-Induced Silencing Complex , Tryptophan , Tryptophan/metabolism , Liver Neoplasms/genetics , Liver Neoplasms/metabolism , Liver Neoplasms/secondary , Humans , Argonaute Proteins/metabolism , Argonaute Proteins/genetics , RNA-Induced Silencing Complex/metabolism , RNA-Induced Silencing Complex/genetics , Colonic Neoplasms/genetics , Colonic Neoplasms/pathology , Colonic Neoplasms/metabolism , MicroRNAs/genetics , MicroRNAs/metabolism , Mice , Animals , Cell Line, Tumor , Disease Models, AnimalABSTRACT
MicroRNAs (miRNAs) are small non-coding RNAs that are known for their role in the post-transcriptional regulation of target genes. Typically, their functions are predicted by first identifying their target genes and then finding biological processes enriched in these targets. Current tools for miRNA functional analysis use only genes with physical binding sites as their targets and exclude other genes that are indirectly targeted transcriptionally through transcription factors. Here, we introduce a method to predict gene ontology (GO) annotations indirectly targeted by microRNAs. The proposed method resulted in better performance in predicting known miRNA-GO term associations compared to the canonical approach. To facilitate miRNA GO enrichment analysis, we developed an R Shiny application, miRinGO, that is freely available online at GitHub.