Quantitative prediction of variant effects on alternative splicing in MAPT using endogenous pre-messenger RNA structure probing
- University of North Carolina at Chapel Hill, United States
- Clemson University, United States
- Baylor College of Medicine, United States
- University of Rochester, United States
Abstract
Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor messenger RNA (mRNA) structure and downstream function is particularly challenging. Here we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on precursor mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the Bact complex with precursor mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a b-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and six newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.
Data availability
Sequencing data have been deposited in SRA under BioProject ID PRJNA762079 and PRJNA812003.DMS Reactivities are available as SNRNASMs at https://bit.ly/2WaDw6FAll data generated or analyzed during this study are included in the manuscript and supporting files; Source Data files have been provided for Figures 1,2,4,5 and 6.Modeling and feature generation code is uploaded at https://git.io/JuSW8
Article and author information
Author details
Funding
National Institutes of Health (R01 HL111527)
- Alain Laederach
National Institutes of Health (R35 GM 140844)
- Alain Laederach
National Institutes of Health (R01 GM076485)
- David H Mathews
National Institutes of Health (R35 GM122532)
- Kevin Weeks
Cancer Prevention and Research Institute of Texas (CPRIT Scholar)
- Anthony M Mustoe
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Kumar et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,906
- views
-
- 438
- downloads
-
- 6
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Quantitative prediction of variant effects on alternative splicing in MAPT using endogenous pre-messenger RNA structure probing
eLife 11:e73888.
https://doi.org/10.7554/eLife.73888
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B-cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution, and dynamics. We present HILARy (high-precision inference of lineages in antibody repertoires), an efficient, fast, and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and dN/dS ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
Antimicrobial resistance is responsible for an alarming number of deaths, estimated at 5 million per year. To combat priority pathogens, like Helicobacter pylori, the development of novel therapies is of utmost importance. Understanding the molecular alterations induced by medications is critical for the design of multi-targeting treatments capable of eradicating the infection and mitigating its pathogenicity. However, the application of bulk omics approaches for unraveling drug molecular mechanisms of action is limited by their inability to discriminate between target-specific modifications and off-target effects. This study introduces a multi-omics method to overcome the existing limitation. For the first time, the Proteome Integral Solubility Alteration (PISA) assay is utilized in bacteria in the PISA-Express format to link proteome solubility with different and potentially immediate responses to drug treatment, enabling us the resolution to understand target-specific modifications and off-target effects. This study introduces a comprehensive method for understanding drug mechanisms and optimizing the development of multi-targeting antimicrobial therapies.
