Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein splicing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The Human Proteoform Project: Defining the human proteome

Proteins are the primary effectors of function in biology, and thus, complete knowledge of their structure and properties is fundamental to deciphering function in basic and translational research. The chemical diversity of proteins is expressed in their many proteoforms, which result from combinations of genetic polymorphisms, RNA splice variants, and posttranslational modifications. This knowledge is foundational for the biological complexes and networks that control biology yet remains largely unknown. We propose here an ambitious initiative to define the human proteome, that is, to generate a definitive reference set of the proteoforms produced from the genome. Several examples of the power and importance of proteoform-level knowledge in disease-based research are presented along with a call for improved technologies in a two-pronged strategy to the Human Proteoform Project.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Targeting the Y-box Binding Protein-1 Axis to Overcome Radiochemotherapy Resistance in Solid Tumors

Multifunctional Y-box binding protein-1 (YB-1) is highly expressed in different human solid tumors and is involved in various cellular processes. DNA damage is the major mechanism by which radiochemotherapy (RCT) induces cell death. On induction of DNA damage, a multicomponent signal transduction network, known as the DNA damage response, is activated to induce cell cycle arrest and initiate DNA repair, which protects cells against damage. YB-1 regulates nearly all cancer hallmarks described to date by participating in DNA damage response, gene transcription, mRNA splicing, translation, and tumor stemness. YB-1 lacks kinase activity, and p90 ribosomal S6 kinase and AKT are the key kinases within the RAS/mitogen-activated protein kinase and phosphoinositide 3-kinase pathways that directly activate YB-1. Thus, the molecular targeting of ribosomal S6 kinase and AKT is thought to be the most effective strategy for blocking the cellular function of YB-1 in human solid tumors. In this review, after describing the prosurvival effect of YB-1 with a focus on DNA damage repair and cancer cell stemness, clinical evidence will be provided indicating an inverse correlation between YB-1 expression and the treatment outcome of solid tumors after RCT. In the interest of being concise, YB-1 signaling cascades will be briefly discussed and the current literature on YB-1 posttranslational modifications will be summarized. Finally, the current status of targeting the YB-1 axis, especially in combination with RCT, will be highlighted.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

An integrated multi-omics approach identifies the landscape of interferon-α-mediated responses of human pancreatic beta cells

Proinflammatory cytokines are important mediators of pancreatic beta cell dysfunction and demise in the early stages of type 1 diabetes (T1D). Interferon-a (IFNa), a type I interferon member, is expressed in the islets of T1D individuals and it is expression and signaling is regulated by both genetic (T1D risk variants) and environmental factors (viral infections) associated to T1D. We presently characterized human beta cells responses to IFNa by combining ATAC-seq, RNA-seq and proteomics assays. The initial beta cell response to IFNa was characterized by major chromatin remodeling, followed by marked changes in transcriptional and translational regulation. IFNa-induced changes in alternative splicing (AS) and first exon usage increased the diversity of transcripts expressed by beta cells. This, combined with changes observed on protein modification/degradation, ER stress and MHC class I, may significantly expand the peptide repertoire presented by beta cells to the immune system. On the other hand, beta cells up-regulated checkpoint proteins, such as PDL1 and HLA-E, that may protect them against the autoimmune assault. Data mining of the present multi-omics analysis led to the identification of two compound classes that revert IFNa effects on human beta cells and may be translated to clinical trials.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1992

Evolution of EF-hand calcium-modulated proteins. II. Domains of several subfamilies have diverse evolutionary histories

In the first report in this series we described the relationships and evolution of 152 individual proteins of the EF-hand subfamilies. Here we add 66 additional proteins and define eight (CDC, TPNV, CLNB, LPS, DGK, 1F8, VIS, TCBP) new subfamilies and seven (CAL, SQUD, CDPK, EFH5, TPP, LAV, CRGP) new unique proteins, which we assume represent new subfamilies. The main focus of this study is the classification of individual EF-hand domains. Five subfamilies--calmodulin, troponin C, essential light chain, regulatory light chain, CDC31/caltractin--and three uniques--call, squidulin, and calcium-dependent protein kinase--are congruent in that all evolved from a common four-domain precursor. In contrast calpain and sarcoplasmic calcium-binding protein (SARC) each evolved from its own one-domain precursor. The remaining 19 subfamilies and uniques appear to have evolved by translocation and splicing of genes encoding the EF-hand domains that were precursors to the congruent eight and to calpain and to SARC. The rates of evolution of the EF-hand domains are slower following formation of the subfamilies and establishment of their functions. Subfamilies are not readily classified by patterns of calcium coordination, interdomain linker stability, and glycine and proline distribution. There are many homoplasies indicating that similar variants of the EF-hand evolved by independent pathways.

Non-NASA Center↗
DOE OSTI2025

Spatial top-down proteomics for the functional characterization of human kidney

Background: The Human Proteome Project has credibly detected nearly 93% of the roughly 20,000 proteins which are predicted by the human genome. However, the proteome is enigmatic, where alterations in amino acid sequences from polymorphisms and alternative splicing, errors in translation, and post-translational modifications result in a proteome depth estimated at several million unique proteoforms. Recently mass spectrometry has been demonstrated in several landmark efforts mapping the human proteoform landscape in bulk analyses. Herein, we developed an integrated workflow for characterizing proteoforms from human tissue in a spatially resolved manner by coupling laser capture microdissection, nanoliter-scale sample preparation, and mass spectrometry imaging. Results: Using healthy human kidney sections as the case study, we focused our analyses on the major functional tissue units including glomeruli, tubules, and medullary rays. After laser capture microdissection, these isolated functional tissue units were processed with microPOTS (microdroplet processing in one-pot for trace samples) for sensitive top-down proteomics measurement. This provided a quantitative database of 616 proteoforms that was further leveraged as a library for mass spectrometry imaging with near-cellular spatial resolution over the entire section. Notably, several mitochondrial proteoforms were found to be differentially abundant between glomeruli and convoluted tubules, and further spatial contextualization was provided by mass spectrometry imaging confirming unique differences identified by microPOTS, and further expanding the field-of-view for unique distributions such as enhanced abundance of a truncated form (1-74) of ubiquitin within cortical regions. Conclusions: We developed an integrated workflow to directly identify proteoforms and reveal their spatial distributions. Where of the 20 differentially abundant proteoforms identified as discriminate between tubules and glomeruli by microPOTS, the vast majority of tubular proteoforms were of mitochondrial origin (8 of 10) where discriminate proteoforms in glomeruli were primarily hemoglobin subunits (9 of 10). These trends were also identified within ion images demonstrating spatially resolved characterization of proteoforms that has the potential to reshape discovery-based proteomics because the proteoforms are the ultimate effector of cellular functions. Applications of this technology have the potential to unravel etiology and pathophysiology of disease states, informing on biologically active proteoforms, which remodel the proteomic landscape in chronic and acute disorders.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

mTOR regulates aerobic glycolysis through NEAT1 and nuclear paraspeckle-mediated mechanism in hepatocellular carcinoma

Hepatocellular Carcinoma (HCC) is a major form of liver cancer and a leading cause of cancer-related death worldwide. New insights into HCC pathobiology and mechanism of drug actions are urgently needed to improve patient outcomes. HCC undergoes metabolic reprogramming of glucose metabolism from respiration to aerobic glycolysis, a phenomenon known as the ‘Warburg Effect’ that supports rapid cancer cell growth, survival, and invasion. mTOR is known to promote Warburg Effect, but the underlying mechanism(s) remains poorly defined. The aim of this study is to understand the mechanism(s) and significance of mTOR regulation of aerobic glycolysis in HCC. We profiled mTORC1-dependent long non-coding RNAs (lncRNAs) by RNA-seq of HCC cells treated with rapamycin. Chromatin immunoprecipitation (ChIP) and luciferase reporter assays were used to explore the transcriptional regulation of NEAT1 by mTORC1. [U- 13 C]-glucose labeling and metabolomic analysis, extracellular acidification Rate (ECAR) by Seahorse XF Analyzer, and glucose uptake assay were used to investigate the role of mTOR-NEAT1-NONO signaling in the regulation of aerobic glycolysis. RNA immunoprecipitation (RIP) and NONO-binding motif scanning were performed to identify the regulatory mechanism of pre-mRNA splicing by mTOR-NEAT1. Myristoylated AKT1 (mAKT1)/NRASV 12 -driven HCC model developed by hydrodynamic transfection (HDT) was employed to explore the significance of mTOR-NEAT1 signaling in HCC tumorigenesis and mTOR-targeted therapy. mTOR regulates lncRNA transcriptome in HCC and that NEAT1 is a major mTOR transcriptional target. Interestingly, although both NEAT1_1 and NEAT1_2 are down-regulated in HCC, only NEAT1_2 is significantly correlated with poor overall survival of HCC patients. NEAT1_2 is the organizer of nuclear paraspeckles that sequester the RNA-binding proteins NONO and SFPQ. We show that upon oncogenic activation, mTORC1 suppresses NEAT1_2 expression and paraspeckle biogenesis, liberating NONO/SFPQ, which in turn, binds to U5 within the spliceosome, stimulating mRNA splicing and expression of key glycolytic enzymes. This series of actions lead to enhanced glucose transport, aerobic glycolytic flux, lactate production, and HCC growth both in vitro and in vivo. Furthermore, the paraspeckle-mediated mechanism is important for the anticancer action of US FDA-approved drugs rapamycin/temsirolimus. These findings reveal a molecular mechanism by which mTOR promotes the ‘Warburg Effect’, which is important for the metabolism and development of HCC, and anticancer response of mTOR-targeted therapy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Structural complementarity facilitates E7820-mediated degradation of RBM39 by DCAF15

The investigational drugs E7820, indisulam and tasisulam (aryl-sulfonamides) promote the degradation of the splicing factor RBM39 in a proteasome-dependent mechanism. While the activity critically depends on the cullin RING ligase substrate receptor DCAF15, the molecular details remain elusive. Here we present the cryo-EM structure of the DDB1–DCAF15–DDA1 core ligase complex bound to RBM39 and E7820 at a resolution of 4.4 Å, together with crystal structures of engineered subcomplexes. We show that DCAF15 adopts a new fold stabilized by DDA1, and that extensive protein–protein contacts between the ligase and substrate mitigate low affinity interactions between aryl-sulfonamides and DCAF15. Here, our data demonstrate how aryl-sulfonamides neo-functionalize a shallow, non-conserved pocket on DCAF15 to selectively bind and degrade RBM39 and the closely related splicing factor RBM23 without the requirement for a high-affinity ligand, which has broad implications for the de novo discovery of molecular glue degraders.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

Left-handed Z-DNA: structure and function

Z-DNA is a high energy conformer of B-DNA that forms in vivo during transcription as a result of torsional strain generated by a moving polymerase. An understanding of the biological role of Z-DNA has advanced with the discovery that the RNA editing enzyme double-stranded RNA adenosine deaminase type I (ADAR1) has motifs specific for the Z-DNA conformation. Editing by ADAR1 requires a double-stranded RNA substrate. In the cases known, the substrate is formed by folding an intron back onto the exon that is targeted for modification. The use of introns to direct processing of exons requires that editing occurs before splicing. Recognition of Z-DNA by ADAR1 may allow editing of nascent transcripts to be initiated immediately after transcription, ensuring that editing and splicing are performed in the correct sequence. Structural characterization of the Z-DNA binding domain indicates that it belongs to the winged helix-turn-helix class of proteins and is similar to the globular domain of histone-H5.

Review↗
DOE OSTI2020

Hemiptera phylogenomic resources: Tree‐based orthology prediction and conserved exon identification

Abstract High‐throughput sequencing of transcriptomes and targeted genomic regions are advancing our knowledge of The Tree of Life. Building phylogenies with regions of the genome requires 1‐to‐1 orthologue resources of genes and noncoding loci. One organismal group that has received little attention in this area is the Hemiptera, the fifth largest insect order represented by ~103,590 named species. Here, we present a set of 3,872 Hemiptera 1‐to‐1 orthogroups based on tree‐based orthology inference of eight Hemiptera species with publicly available genome sequences. We also estimate a set of 406 orthologous exons with similar mRNA splice sites that can be used for Sanger sequencing and develop enrichment probes for targeted genome sequencing for phylogenomic inference. We show this novel set of orthologues is informative at the protein, coding sequence and exon molecular levels and provides robust branch support in both gene tree–species tree methods and concatenated sequence phylogenies. In addition, we demonstrate the utility of these loci to resolve relationships in whiteflies, Bemisia tabaci , a large species complex with few phylogenomic resources. Last, we compare our Hemiptera phylogeny with previously published phylogenies and other orthologue databases, while providing suggestions on further improvement to this phylogenomic resource.

Owen, Christopher L.↗
NASA NTRS1993

An intron within the 16S ribosomal RNA gene of the archaeon Pyrobaculum aerophilum

The 16S rRNA genes of Pyrobaculum aerophilum and Pyrobaculum islandicum were amplified by the polymerase chain reaction, and the resulting products were sequenced directly. The two organisms are closely related by this measure (over 98% similar). However, they differ in that the (lone) 16S rRNA gene of Pyrobaculum aerophilum contains a 713-bp intron not seen in the corresponding gene of Pyrobaculum islandicum. To our knowledge, this is the only intron so far reported in the small subunit rRNA gene of a prokaryote. Upon excision the intron is circularized. A secondary structure model of the intron-containing rRNA suggests a splicing mechanism of the same type as that invoked for the tRNA introns of the Archaea and Eucarya and 23S rRNAs of the Archaea. The intron contains an open reading frame whose protein translation shows no certain homology with any known protein sequence.

NASA Discipline Exobiology↗
DOE OSTI2022

SARS-CoV-2 Diverges from Other Betacoronaviruses in Only Partially Activating the IRE1α/XBP1 Endoplasmic Reticulum Stress Pathway in Human Lung-Derived Cells

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed over 6 million individuals worldwide and continues to spread in countries where vaccines are not yet widely available or its citizens are hesitant to become vaccinated. Therefore, it is critical to unravel the molecular mechanisms that allow SARS-CoV-2 and other coronaviruses to infect and overtake the host machinery of human cells. Coronavirus replication triggers endoplasmic reticulum (ER) stress and activation of the unfolded protein response (UPR), a key host cell pathway widely believed to be essential for viral replication. We examined the master UPR sensor IRE1α kinase/RNase and its downstream transcription factor effector XBP1s, which is processed through an IRE1α-mediated mRNA splicing event, in human lung-derived cells infected with betacoronaviruses. We found that human respiratory coronavirus OC43 (HCoV-OC43), Middle East respiratory syndrome coronavirus (MERS-CoV), and murine coronavirus (MHV) all induce ER stress and strongly trigger the kinase and RNase activities of IRE1α as well as XBP1 splicing. In contrast, SARS-CoV-2 only partially activates IRE1α through autophosphorylation, but its RNase activity fails to splice XBP1. Moreover, while IRE1α was dispensable for replication in human cells for all coronaviruses tested, it was required for maximal expression of genes associated with several key cellular functions, including the interferon signaling pathway, during SARS-CoV-2 infection. Our data suggest that SARS-CoV-2 actively inhibits the RNase of autophosphorylated IRE1α, perhaps as a strategy to eliminate detection by the host immune system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Top-down proteomics

Proteoforms arising from posttranslational modifications, genetic polymorphisms, and RNA splice variants, play a pivotal role as the key drivers in biology. Thus, a comprehensive understanding of proteoforms is essential for unraveling the intricacies of biological systems and bridging the gap between genotype and phenotype. By analyzing whole proteins without digestion, top-down proteomics (TDP) provides a holistic view of the proteome and presents a next-generation approach for deciphering protein function, uncovering disease mechanisms, and advancing precision medicine. This Primer embarks on a journey into the world of TDP by encapsulating its historical context, underlying principles, recent advances, and an outlook on the future of TDP. The experimental section navigates instrumentation, sample preparation, intact protein separation, tandem mass spectrometry techniques, and data collection. Results decipher raw data, visualize intact protein spectra, unravel data analysis, and explain proteoform identification, characterization, and quantitation, as well as statistical analysis. Various applications of TDP spanning the human proteoform project, biomedical, biopharmaceutical, and clinical applications are described. These are complemented by discussions on measurement reproducibility, limitations, and a forward-looking perspective outlining uncharted waters where the field can advance, and potential exciting future applications of TDP.

Roberts, David S.↗
DOE OSTI2025

Sorted-cell proteomics reveals an AT1-associated epithelial cornification phenotype and suggests endothelial redox imbalance in human bronchopulmonary dysplasia

Bronchopulmonary dysplasia (BPD) is a neonatal lung disease characterized by inflammation and scarring leading to long-term tissue damage. Previous whole tissue proteomics identified BPD-specific proteome changes and cell type shifts. Little is known about the proteome-level changes within specific cell populations in disease. Here, we sorted epithelial (EPI) and endothelial (ENDO) cell populations based on their differential surface markers from normal and BPD human lungs. Using a low-input compatible sample preparation method (MicroPOT), proteins were extracted and digested into peptides and subjected to liquid chromatography-tandem mass spectrometry (LC-MS/MS) proteome analysis. Of the 4,970 proteins detected, 293 were modulated in abundance or detection in the EPI population and 422 were modulated in ENDO cells. Modulation of proteins associated with actin-cytoskeletal function, such as SCEL, LMO7, and TBA1B was observed in the BPD EPIs. Using confocal imaging and analysis, we validated the presence of aberrant multilayer-like structures comprising SCEL and LMO7, known to be associated with epidermal cornification, in the human BPD lung. This is the first report of the accumulation of cornification-associated proteins in BPD. Their localization in the alveolar parenchyma, primarily associated with alveolar type 1 (AT1) cells, suggests a role in the BPD postinjury response. In the ENDOs, redox balance and mitochondrial function pathways were modulated. Alternative mRNA splicing and cell proliferative functions were elevated in both populations, suggesting potential dysregulation of cell progenitor fate. This study characterized the proteome of epithelial and endothelial cells from the BPD lung for the first time, identifying population-specific changes in BPD pathogenesis.

BPD↗
DOE OSTI2020

Global analysis of SBP gene family in Brachypodium distachyon reveals its association with spike development

SQUAMOSA-promoter binding like proteins (SBPs/SPLs) are plant specific transcription factors targeted by miR156 and involved in various biological pathways, playing multi-faceted developmental roles. This gene family is not well characterized in Brachypodium. We identified a total of 18 SBP genes in B. distachyon genome. Phylogenetic analysis revealed that SBP gene family in Brachypodium expanded through large scale duplication. A total of 10 BdSBP genes were identified as targets of miR156. Transcript cleavage analysis of selected BdSBPs by miR156 confirmed their antagonistic connection. Alternative splicing was observed playing an important role in BdSBPs and miR156 interaction. Characterization of T-DNA Bdsbp9 mutant showed reduced plant growth and spike length, reflecting its involvement in the spike development. Expression of a majority of BdSBPs elevated during spikelet initiation. Specifically, BdSBP1 and BdSBP3 differentially expressed in response to vernalization. Differential transcript abundance of BdSBP1, BdSBP3, BdSBP8, BdSBP9, BdSBP14, BdSBP18 and BdSBP23 genes was observed during the spike development under high temperature. Co-expression network, protein–protein interaction and biological pathway analysis indicate that BdSBP genes mainly regulate transcription, hormone, RNA and transport pathways. Our work reveals the multi-layered control of SBP genes and demonstrates their association with spike development and temperature sensitivity in Brachypodium.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

A complex ligase ribozyme evolved in vitro from a group I ribozyme domain

Like most proteins, complex RNA molecules often are modular objects made up of distinct structural and functional domains. The component domains of a protein can associate in alternative combinations to form molecules with different functions. These observations raise the possibility that complex RNAs also can be assembled from preexisting structural and functional domains. To test this hypothesis, an in vitro evolution procedure was used to isolate a previously undescribed class of complex ligase ribozymes, starting from a pool of 10(16) different RNA molecules that contained a constant region derived from a large structural domain that occurs within self-splicing group I ribozymes. Attached to this constant region were three hypervariable regions, totaling 85 nucleotides, that gave rise to the catalytic motif within the evolved catalysts. The ligase ribozymes catalyze formation of a 3',5'-phosphodiester linkage between adjacent template-bound oligonucleotides, one bearing a 3' hydroxyl and the other a 5' triphosphate. Ligation occurs in the context of a Watson-Crick duplex, with a catalytic rate of 0.26 min(-1) under optimal conditions. The constant region is essential for catalytic activity and appears to retain the tertiary structure of the group I ribozyme. This work demonstrates that complex RNA molecules, like their protein counterparts, can share common structural domains while exhibiting distinct catalytic functions.

NASA Discipline Exobiology↗
DOE OSTI2022

An auto-inhibited state of protein kinase G and implications for selective activation

Cyclic GMP-dependent protein kinases (PKGs) are key mediators of the nitric oxide/cyclic guanosine monophosphate (cGMP) signaling pathway that regulates biological functions as diverse as smooth muscle contraction, cardiac function, and axon guidance. Understanding how cGMP differentially triggers mammalian PKG isoforms could lead to new therapeutics that inhibit or activate PKGs, complementing drugs that target nitric oxide synthases and cyclic nucleotide phosphodiesterases in this signaling axis. Alternate splicing of PRKG1 transcripts confers distinct leucine zippers, linkers, and auto-inhibitory (AI) pseudo-substrate sequences to PKG Iα and Iβ that result in isoform-specific activation properties, but the mechanism of enzyme auto-inhibition and its alleviation by cGMP is not well understood. Here, we present a crystal structure of PKG Iβ in which the AI sequence and the cyclic nucleotide-binding (CNB) domains are bound to the catalytic domain, providing a snapshot of the auto-inhibited state. Specific contacts between the PKG Iβ AI sequence and the enzyme active site help explain isoform-specific activation constants and the effects of phosphorylation in the linker. We also present a crystal structure of a PKG I CNB domain with an activating mutation linked to Thoracic Aortic Aneurysms and Dissections. Similarity of this structure to wildtype cGMP-bound domains and differences with the auto-inhibited enzyme provide a mechanistic basis for constitutive activation. We show that PKG Iβ auto-inhibition is mediated by contacts within each monomer of the native full-length dimeric protein, and using the available structural and biochemical data we develop a model for the regulation and cooperative activation of PKGs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Transcriptomic and proteomic signatures of stemness and differentiation in the colon crypt

The intestinal epithelium is replaced weekly by non-quiescent stem cells with kinetics that rely on a rapid loss of stemness and choice for secretory or absorptive lineage differentiation. To determine how the cellular transcriptome and proteome changes during these transitions, we developed a new cell sorting method to purify stem cells, secretory and absorptive progenitor cells, and mature, differentiated cells. Transcriptome analyses revealed that as stem cells transition to the progenitor stage, alternative mRNA splicing and polyadenylation dominate changes in the transcriptome. In contrast, as progenitors differentiate into mature cell types, alterations in gene expression and mRNA levels drive the changes. RNA processing targets mRNAs encoding regulators of cell cycle, RNA regulators, cell adhesion, SUMOylation, and Wnt and Notch signaling. Additionally, global proteome profiling of sorted cell populations detected >2,800 proteins and revealed RNA:protein patterns of abundance and correlation. Paired together, these data highlight new potentials for autocrine and feedback regulation and provide new insights into cell state transitions in the crypt.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

An expectation–maximization framework for comprehensive prediction of isoform-specific functions

Advances in RNA sequencing technologies have achieved an unprecedented accuracy in the quantification of mRNA isoforms, but our knowledge of isoform-specific functions has lagged behind. There is a need to understand the functional consequences of differential splicing, which could be supported by the generation of accurate and comprehensive isoform-specific gene ontology annotations. We present isoform interpretation, a method that uses expectation–maximization to infer isoform-specific functions based on the relationship between sequence and functional isoform similarity. We predicted isoform-specific functional annotations for 85 617 isoforms of 17 900 protein-coding human genes spanning a range of 17 430 distinct gene ontology terms. Comparison with a gold-standard corpus of manually annotated human isoform functions showed that isoform interpretation significantly outperforms state-of-the-art competing methods. We provide experimental evidence that functionally related isoforms predicted by isoform interpretation show a higher degree of domain sharing and expression correlation than functionally related genes. We also show that isoform sequence similarity correlates better with inferred isoform function than with gene-level function.

59 BASIC BIOLOGICAL SCIENCES↗