Structure-guided isoform identification for the human transcriptome
- Johns Hopkins University, United States
- Seoul National University, Republic of Korea
Abstract
Recently-developed methods to predict three-dimensional protein structure with high accuracy have opened new avenues for genome and proteome research. We explore a new hypothesis in genome annotation, namely whether computationally predicted structures can help to identify which of multiple possible gene isoforms represents a functional protein product. Guided by protein structure predictions, we evaluated over 230,000 isoforms of human protein-coding genes assembled from over 10,000 RNA sequencing experiments across many human tissues. From this set of assembled transcripts, we identified hundreds of isoforms with more confidently predicted structure and potentially superior function in comparison to canonical isoforms in the latest human gene database. We illustrate our new method with examples where structure provides a guide to function in combination with expression and evolutionary evidence. Additionally, we provide the complete set of structures as a resource to better understand the function of human genes and their isoforms. These results demonstrate the promise of protein structure prediction as a genome annotation tool, allowing us to refine even the most highly-curated catalog of human proteins. More generally we demonstrate a practical, structure-guided approach that can be used to enhance the annotation of any genome.
Data availability
Gene identifiers for all predicted protein isoforms as well as pLDDT scores and evolutionary conservation data from mouse can be found in table S1. Predicted scores and GTEx expression data for all isoforms overlapping a MANE locus can be found in table S2. Data for the 401 alternate isoforms with evidence of relatively superior structure, and possibly superior function, can be found in table S3. Additionally, all data can be downloaded from the project website, isoform.io.
Article and author information
Author details
Funding
National Institutes of Health (R01-HG006677)
- Steven L Salzberg
National Institutes of Health (R35-GM130151)
- Steven L Salzberg
National Research Foundation of Korea (2019R1-A6A1-A10073437)
- Martin Steinegger
National Research Foundation of Korea (2020M3-A9G7-103933)
- Martin Steinegger
National Research Foundation of Korea (2021-R1C1-C102065)
- Martin Steinegger
National Research Foundation of Korea (2021-M3A9-I4021220)
- Martin Steinegger
Seoul National University (Creative-Pioneering Researchers Program)
- Martin Steinegger
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Sommer 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
-
- 5,008
- views
-
- 527
- downloads
-
- 23
- 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)
Structure-guided isoform identification for the human transcriptome
eLife 11:e82556.
https://doi.org/10.7554/eLife.82556
Further reading
-
- Computational and Systems Biology
- Genetics and Genomics
Untranslated regions (UTRs) contain crucial regulatory elements for RNA stability, translation and localization, so their integrity is indispensable for gene expression. Approximately 3.7% of genetic variants associated with diseases occur in UTRs, yet a comprehensive understanding of UTR variant functions remains limited due to inefficient experimental and computational assessment methods. To systematically evaluate the effects of UTR variants on RNA stability, we established a massively parallel reporter assay on 6555 UTR variants reported in human disease databases. We examined the RNA degradation patterns mediated by the UTR library in two cell lines, and then applied LASSO regression to model the influential regulators of RNA stability. We found that UA dinucleotides and UA-rich motifs are the most prominent destabilizing element. Gain of UA dinucleotide outlined mutant UTRs with reduced stability. Studies on endogenous transcripts indicate that high UA-dinucleotide ratios in UTRs promote RNA degradation. Conversely, elevated GC content and protein binding on UA dinucleotides protect high-UA RNA from degradation. Further analysis reveals polarized roles of UA-dinucleotide-binding proteins in RNA protection and degradation. Furthermore, the UA-dinucleotide ratio of both UTRs is a common characteristic of genes in innate immune response pathways, implying a coordinated stability regulation through UTRs at the transcriptomic level. We also demonstrate that stability-altering UTRs are associated with changes in biobank-based health indices, underscoring the importance of precise UTR regulation for wellness. Our study highlights the importance of RNA stability regulation through UTR primary sequences, paving the way for further exploration of their implications in gene networks and precision medicine.
-
- Genetics and Genomics
Flavin-containing monooxygenases (FMOs) are a conserved family of xenobiotic enzymes upregulated in multiple longevity interventions, including nematode and mouse models. Previous work supports that C. elegans fmo-2 promotes longevity, stress resistance, and healthspan by rewiring endogenous metabolism. However, there are five C. elegans FMOs and five mammalian FMOs, and it is not known whether promoting longevity and health benefits is a conserved role of this gene family. Here, we report that expression of C. elegans fmo-4 promotes lifespan extension and paraquat stress resistance downstream of both dietary restriction and inhibition of mTOR. We find that overexpression of fmo-4 in just the hypodermis is sufficient for these benefits, and that this expression significantly modifies the transcriptome. By analyzing changes in gene expression, we find that genes related to calcium signaling are significantly altered downstream of fmo-4 expression. Highlighting the importance of calcium homeostasis in this pathway, fmo-4 overexpressing animals are sensitive to thapsigargin, an ER stressor that inhibits calcium flux from the cytosol to the ER lumen. This calcium/fmo-4 interaction is solidified by data showing that modulating intracellular calcium with either small molecules or genetics can change expression of fmo-4 and/or interact with fmo-4 to affect lifespan and stress resistance. Further analysis supports a pathway where fmo-4 modulates calcium homeostasis downstream of activating transcription factor-6 (atf-6), whose knockdown induces and requires fmo-4 expression. Together, our data identify fmo-4 as a longevity-promoting gene whose actions interact with known longevity pathways and calcium homeostasis.
