Splice-Junction-Based Mapping of Alternative Isoforms in the Human Proteome

A Ratiometric Catalog of Protein Isoform Shifts in the Cardiac Fetal Gene Program

The fetal genetic program orchestrates cardiac development and the re-expression of fetal genes is thought to underlie cardiac disease and adaptation. Here, a proteomics ratio test using mass spectrometry is applied to find protein isoforms with statistically significant usage differences in the fetal vs. postnatal mouse heart. Changes in isoform usage ratios are pervasive at the protein level, with 104 significant events observed, including 88 paralog-derived isoform switching events and 16 splicing-derived isoform switching events between fetal and postnatal hearts. The ratiometric proteomic comparisons rediscovered hallmark fetal gene signatures including a postnatal switch from fetal β (MYH7) toward ɑ (MYH6) myosin heavy chains and from slow skeletal muscle (TNNI1) toward cardiac (TNNI3) troponin I. Altered usages in metabolic proteins are prominent, including a platelet to muscle phosphofructokinase (PFKP - PFKM), enolase 1 to 3 (ENO1 - ENO3), and alternative splicing of pyruvate kinase M2 toward M1 (PKM2 - PKM1) isoforms in glycolysis. The data also revealed a parallel change in mitochondrial proteins in cardiac development, suggesting the shift toward aerobic respiration involves also a remodeling of the mitochondrial protein isoform proportion. Finally, a number of glycolytic protein isoforms revert toward their fetal forms in adult hearts under pathological cardiac hypertrophy, suggesting their functional roles in adaptive or maladaptive response, but this reversal is partial. In summary, this work presents a catalog of ratiometric protein markers of the fetal genetic program of the mouse heart, including previously unreported splice isoform markers.

Proteomics applications in next generation induced pluripotent stem cell models

INTRODUCTION

Induced pluripotent stem (iPS) cell technology has transformed biomedical research. New opportunities now exist to create new organoids, microtissues, and body-on-a-chip systems for basic biology investigations and clinical translations.

AREAS COVERED

We discuss the utility of proteomics for attaining an unbiased view into protein expression changes during iPS cell differentiation, cell maturation, and tissue generation. The ability to discover cell-type specific protein markers during the differentiation and maturation of iPS-derived cells has led to new strategies to improve cell production yield and fidelity. In parallel, proteomic characterization of iPS-derived organoids is helping to realize the goal of bridging in vitro and in vivo systems.

EXPERT OPINIONS

We discuss some current challenges of proteomics in iPS cell research and future directions, including the integration of proteomic and transcriptomic data for systems-level analysis.

rMATS-turbo: an efficient and flexible computational tool for alternative splicing analysis of large-scale RNA-seq data

Pre-mRNA alternative splicing is a prevalent mechanism for diversifying eukaryotic transcriptomes and proteomes. Regulated alternative splicing plays a role in many biological processes, and dysregulated alternative splicing is a feature of many human diseases. Short-read RNA sequencing (RNA-seq) is now the standard approach for transcriptome-wide analysis of alternative splicing. Since 2011, our laboratory has developed and maintained Replicate Multivariate Analysis of Transcript Splicing (rMATS), a computational tool for discovering and quantifying alternative splicing events from RNA-seq data. Here we provide a protocol for the contemporary version of rMATS, rMATS-turbo, a fast and scalable re-implementation that maintains the statistical framework and user interface of the original rMATS software, while incorporating a revamped computational workflow with a substantial improvement in speed and data storage efficiency. The rMATS-turbo software scales up to massive RNA-seq datasets with tens of thousands of samples. To illustrate the utility of rMATS-turbo, we describe two representative application scenarios. First, we describe a broadly applicable two-group comparison to identify differential alternative splicing events between two sample groups, including both annotated and novel alternative splicing events. Second, we describe a quantitative analysis of alternative splicing in a large-scale RNA-seq dataset (~1,000 samples), including the discovery of alternative splicing events associated with distinct cell states. We detail the workflow and features of rMATS-turbo that enable efficient parallel processing and analysis of large-scale RNA-seq datasets on a compute cluster. We anticipate that this protocol will help the broad user base of rMATS-turbo make the best use of this software for studying alternative splicing in diverse biological systems.

IS-PRM-based peptide targeting informed by long-read sequencing for alternative proteome detection

Alternative splicing is a major contributor of transcriptomic complexity, but the extent to which transcript isoforms are translated into stable, functional protein isoforms is unclear. Furthermore, detection of relatively scarce isoform-specific peptides is challenging, with many protein isoforms remaining uncharted due to technical limitations. Recently, a family of advanced targeted MS strategies, termed internal standard parallel reaction monitoring (IS-PRM), have demonstrated multiplexed, sensitive detection of pre-defined peptides of interest. Such approaches have not yet been used to confirm existence of novel peptides. Here, we present a targeted proteogenomic approach that leverages sample-matched long-read RNA sequencing (LR RNAseq) data to predict potential protein isoforms with prior transcript evidence. Predicted tryptic isoform-specific peptides, which are specific to individual gene product isoforms, serve as "triggers" and "targets" in the IS-PRM method, Tomahto. Using the model human stem cell line WTC11, LR RNAseq data were generated and used to inform the generation of synthetic standards for 192 isoform-specific peptides (114 isoforms from 55 genes). These synthetic "trigger" peptides were labeled with super heavy tandem mass tags (TMT) and spiked into TMT-labeled WTC11 tryptic digest, predicted to contain corresponding endogenous "target" peptides. Compared to DDA mode, Tomahto increased detectability of isoforms by 3.6-fold, resulting in the identification of five previously unannotated isoforms. Our method detected protein isoform expression for 43 out of 55 genes corresponding to 54 resolved isoforms. This LR RNA seq-informed Tomahto targeted approach, called LRP-IS-PRM, is a new modality for generating protein-level evidence of alternative isoforms - a critical first step in designing functional studies and eventually clinical assays.

Sarcomeric tropomyosin expression during human iPSC differentiation into cardiomyocytes

Tropomyosin (TPM) is an essential sarcomeric component, stabilizing the thin filament and facilitating actin's interaction with myosin. In mammals, including humans, there are four TPM genes (TPM1, TPM2, TPM3, and TPM4) each of which generates a multitude of TPM isoforms via alternative splicing and using different promoters. In this study, we have examined the expression of transcripts as well as proteins of various sarcomeric TPM isoforms during human inducible pluripotent stem cell differentiation into cardiomyocytes. During the differentiation time course, we harvested cells on Days 0, 5, 10, 15, and 20 to analyze for various sarcomeric TPM transcripts by qRT-PCR and for sarcomeric TPM proteins using two-dimensional Western blot with sarcomeric TPM-specific CH1 monoclonal antibody followed by mass spectra analyses. Our results show increasing levels of total TPM transcripts and proteins during the period of differentiation, but varying levels of specific TPM isoforms during the same period. By Day 20, the rank order of TPM transcripts was TPM1α > TPM1κ > TPM2α > TPM1μ > TPM3α > TPM4α. TPM1α was the dominant protein produced with some TPM2 and much less TPM1κ and μ. Interestingly, small amounts of two lower molecular weight TPM3 isoforms were detected on Day 15. To the best of our knowledge this is the first demonstration of TPM1μ non-muscle isoform protein expression before and during cardiac differentiation.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched with intrinsically disordered regions. Moreover, over two-thirds of such regions are predicted to function in protein binding and RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

Integrative Proteogenomics for Differential Expression and Splicing Variation in a DM1 Mouse Model

Dysregulated mRNA splicing is involved in the pathogenesis of many diseases including cancer, neurodegenerative diseases, and muscular dystrophies such as myotonic dystrophy type 1 (DM1). Comprehensive assessment of dysregulated splicing on the transcriptome and proteome level has been methodologically challenging, and thus investigations have often been targeting only few genes. Here, we performed a large-scale coordinated transcriptomic and proteomic analysis to characterize a DM1 mouse model (HSALR) in comparison to wild type. Our integrative proteogenomics approach comprised gene- and splicing-level assessments for mRNAs and proteins. It recapitulated many known instances of aberrant mRNA splicing in DM1 and identified new ones. It enabled the design and targeting of splicing-specific peptides and confirmed the translation of known instances of aberrantly spliced disease-related genes (e.g., Atp2a1, Bin1, Ryr1), complemented by novel findings (Flnc and Ywhae). Comparative analysis of large-scale mRNA and protein expression data showed quantitative agreement of differentially expressed genes and splicing patterns between disease and wild type. We hence propose this work as a suitable blueprint for a robust and scalable integrative proteogenomic strategy geared toward advancing our understanding of splicing-based disorders. With such a strategy, splicing-based biomarker candidates emerge as an attractive and accessible option, as they can be efficiently asserted on the mRNA and protein level in coordinated fashion.

Expression and splicing mediate distinct biological signals

BACKGROUND

Through alternative splicing, most human genes produce multiple isoforms in a cell-, tissue-, and disease-specific manner. Numerous studies show that alternative splicing is essential for development, diseases, and their treatments. Despite these important examples, the extent and biological relevance of splicing are currently unknown.

RESULTS

To solve this problem, we developed pairedGSEA and used it to profile transcriptional changes in 100 representative RNA-seq datasets. Our systematic analysis demonstrates that changes in splicing, on average, contribute to 48.1% of the biological signal in expression analyses. Gene-set enrichment analysis furthermore indicates that expression and splicing both convey shared and distinct biological signals.

CONCLUSIONS

These findings establish alternative splicing as a major regulator of the human condition and suggest that most contemporary RNA-seq studies likely miss out on critical biological insights. We anticipate our results will contribute to the transition from a gene-centric to an isoform-centric research paradigm.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched in intrinsically disordered regions, and over two-thirds of such regions are predicted to function in protein binding and/or RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

SEPepQuant enhances the detection of possible isoform regulations in shotgun proteomics

Shotgun proteomics is essential for protein identification and quantification in biomedical research, but protein isoform characterization is challenging due to the extensive number of peptides shared across proteins, hindering our understanding of protein isoform regulation and their roles in normal and disease biology. We systematically assess the challenge and opportunities of shotgun proteomics-based protein isoform characterization using in silico and experimental data, and then present SEPepQuant, a graph theory-based approach to maximize isoform characterization. Using published data from one induced pluripotent stem cell study and two human hepatocellular carcinoma studies, we demonstrate the ability of SEPepQuant in addressing the key limitations of existing methods, providing more comprehensive isoform-level characterization, identifying hundreds of isoform-level regulation events, and facilitating streamlined cross-study comparisons. Our analysis provides solid evidence to support a widespread role of protein isoform regulation in normal and disease processes, and SEPepQuant has broad applications to biological and translational research.

Nanopore sequencing unveils the complexity of the cold-activated murine brown adipose tissue transcriptome

Alternative transcription increases transcriptome complexity by expression of multiple transcripts per gene. Annotation and quantification of transcripts using short-read sequencing is non-trivial. Long-read sequencing aims at overcoming these problems by sequencing full-length transcripts. Activation of brown adipose tissue (BAT) thermogenesis involves major transcriptomic remodeling and positively affects metabolism via increased energy expenditure. We benchmark Oxford Nanopore Technology (ONT) long-read sequencing protocols to Illumina short-read sequencing assessing alignment characteristics, gene and transcript detection and quantification, differential gene and transcript expression, transcriptome reannotation, and differential transcript usage (DTU). We find ONT sequencing is superior to Illumina for transcriptome reassembly, reducing the risk of false-positive events by unambiguously mapping reads to transcripts. We identified novel isoforms of genes undergoing DTU in cold-activated BAT including Cars2, Adtrp, Acsl5, Scp2, Aldoa, and Pde4d, validated by real-time PCR. The reannotated murine BAT transcriptome established here provides a framework for future investigations into the regulation of BAT.

Multi-omic stratification of the missense variant cysteinome

Cancer genomes are rife with genetic variants; one key outcome of this variation is gain-ofcysteine, which is the most frequently acquired amino acid due to missense variants in COSMIC. Acquired cysteines are both driver mutations and sites targeted by precision therapies. However, despite their ubiquity, nearly all acquired cysteines remain uncharacterized. Here, we pair cysteine chemoproteomics-a technique that enables proteome-wide pinpointing of functional, redox sensitive, and potentially druggable residues-with genomics to reveal the hidden landscape of cysteine acquisition. For both cancer and healthy genomes, we find that cysteine acquisition is a ubiquitous consequence of genetic variation that is further elevated in the context of decreased DNA repair. Our chemoproteogenomics platform integrates chemoproteomic, whole exome, and RNA-seq data, with a customized 2-stage false discovery rate (FDR) error controlled proteomic search, further enhanced with a user-friendly FragPipe interface. Integration of CADD predictions of deleteriousness revealed marked enrichment for likely damaging variants that result in acquisition of cysteine. By deploying chemoproteogenomics across eleven cell lines, we identify 116 gain-of-cysteines, of which 10 were liganded by electrophilic druglike molecules. Reference cysteines proximal to missense variants were also found to be pervasive, 791 in total, supporting heretofore untapped opportunities for proteoform-specific chemical probe development campaigns. As chemoproteogenomics is further distinguished by sample-matched combinatorial variant databases and compatible with redox proteomics and small molecule screening, we expect widespread utility in guiding proteoform-specific biology and therapeutic discovery.

Defining the Sarcomeric Proteoform Landscape in Ischemic Cardiomyopathy by Top-Down Proteomics

Ischemic cardiomyopathy (ICM) is a prominent form of heart failure, but the molecular mechanisms underlying ICM remain relatively understudied due to marked phenotypic heterogeneity. Alterations in post-translational modifications (PTMs) and isoform switches in sarcomeric proteins play important roles in cardiac pathophysiology. Thus, it is essential to define sarcomeric proteoform landscape to better understand ICM. Herein, we have implemented a top-down liquid chromatography (LC)-mass spectrometry (MS)-based proteomics method for the identification and quantification of sarcomeric proteoforms in the myocardia of donors without heart diseases (n = 16) compared to end-stage ICM patients (n = 16). Importantly, quantification of post-translational modifications (PTMs) and expression reveal significant changes in various sarcomeric proteins extracted from ICM tissues. Changes include altered phosphorylation and expression of cardiac troponin I (cTnI) and enigma homologue 2 (ENH2) as well as an increase in muscle LIM protein (MLP) and calsarcin-1 (Cal-1) phosphorylation in ICM hearts. Our results imply that the contractile apparatus of the sarcomere is severely dysregulated during ICM. Thus, this is the first study to uncover significant molecular changes to multiple sarcomeric proteins in the LV myocardia of the end-stage ICM patients using liquid chromatography-mass spectrometry (LC-MS)-based top-down proteomics. Raw data are available via the PRIDE repository with identifier PXD038066.

Mass spectrometry-based targeted proteomics for analysis of protein mutations

Cancers are caused by accumulated DNA mutations. This recognition of the central role of mutations in cancer and recent advances in next-generation sequencing, has initiated the massive screening of clinical samples and the identification of 1000s of cancer-associated gene mutations. However, proteomic analysis of the expressed mutation products lags far behind genomic (transcriptomic) analysis. With comprehensive global proteomics analysis, only a small percentage of single nucleotide variants detected by DNA and RNA sequencing have been observed as single amino acid variants due to current technical limitations. Proteomic analysis of mutations is important with the potential to advance cancer biomarker development and the discovery of new therapeutic targets for more effective disease treatment. Targeted proteomics using selected reaction monitoring (also known as multiple reaction monitoring) and parallel reaction monitoring, has emerged as a powerful tool with significant advantages over global proteomics for analysis of protein mutations in terms of detection sensitivity, quantitation accuracy and overall practicality (e.g., reliable identification and the scale of quantification). Herein we review recent advances in the targeted proteomics technology for enhancing detection sensitivity and multiplexing capability and highlight its broad biomedical applications for analysis of protein mutations in human bodily fluids, tissues, and cell lines. Furthermore, we review recent applications of top-down proteomics for analysis of protein mutations. Unlike the commonly used bottom-up proteomics which requires digestion of proteins into peptides, top-down proteomics directly analyzes intact proteins for more precise characterization of mutation isoforms. Finally, general perspectives on the potential of achieving both high sensitivity and high sample throughput for large-scale targeted detection and quantification of important protein mutations are discussed.

Transcriptomic and Proteomic Analyses of Celery Cytoplasmic Male Sterile Line and Its Maintainer Line

Male sterility is a common phenomenon in the plant kingdom and based on the organelles harboring the male-sterility genes, it can be classified into the genic male sterility (GMS) and the cytoplasmic male sterility (CMS). In every generation, CMS can generate 100% male-sterile population, which is very important for the breeders to take advantage of the heterosis and for the seed producers to guarantee the seed purity. Celery is a cross-pollinated plant with the compound umbel type of inflorescence which carries hundreds of small flowers. These characteristics make CMS the only option to produce the commercial hybrid celery seeds. In this study, transcriptomic and proteomic analyses were performed to identify genes and proteins that are associated with celery CMS. A total of 1255 differentially expressed genes (DEGs) and 89 differentially expressed proteins (DEPs) were identified between the CMS and its maintainer line, then 25 genes were found to differentially expressed at both the transcript and protein levels. Ten DEGs involved in the fleece layer and outer pollen wall development were identified by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, most of which were down-regulated in the sterile line W99A. These DEGs and DEPs were mainly enriched in the pathways of "phenylpropanoid/sporopollenin synthesis/metabolism", "energy metabolism", "redox enzyme activity" and "redox processes". Results obtained in this study laid a foundation for the future investigation of mechanisms of pollen development as well as the reasons for the CMS in celery.

Identification of Alternative Splicing in Proteomes of Human Melanoma Cell Lines without RNA Sequencing Data

Alternative splicing is one of the main regulation pathways in living cells beyond simple changes in the level of protein expression. Most of the approaches proposed in proteomics for the identification of specific splicing isoforms require a preliminary deep transcriptomic analysis of the sample under study, which is not always available, especially in the case of the re-analysis of previously acquired data. Herein, we developed new algorithms for the identification and validation of protein splice isoforms in proteomic data in the absence of RNA sequencing of the samples under study. The bioinformatic approaches were tested on the results of proteome analysis of human melanoma cell lines, obtained earlier by high-resolution liquid chromatography and mass spectrometry (LC-MS). A search for alternative splicing events for each of the cell lines studied was performed against the database generated from all known transcripts (RefSeq) and the one composed of peptide sequences, which included all biologically possible combinations of exons. The identifications were filtered using the prediction of both retention times and relative intensities of fragment ions in the corresponding mass spectra. The fragmentation mass spectra corresponding to the discovered alternative splicing events were additionally examined for artifacts. Selected splicing events were further validated at the mRNA level by quantitative PCR.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

Multi-omics approach to identifying isoform variants as therapeutic targets in cancer patients

Cancer-specific alternatively spliced events (ASE) play a role in cancer pathogenesis and can be targeted by immunotherapy, oligonucleotide therapy, and small molecule inhibition. However, identifying actionable ASE targets remains challenging due to the uncertainty of its protein product, structure impact, and proteoform (protein isoform) function. Here we argue that an integrated multi-omics profiling strategy can overcome these challenges, allowing us to mine this untapped source of targets for therapeutic development. In this review, we will provide an overview of current multi-omics strategies in characterizing ASEs by utilizing the transcriptome, proteome, and state-of-art algorithms for protein structure prediction. We will discuss limitations and knowledge gaps associated with each technology and informatics analytics. Finally, we will discuss future directions that will enable the full integration of multi-omics data for ASE target discovery.

Proteogenomics reveals sex-biased aging genes and coordinated splicing in cardiac aging

The risks of heart diseases are significantly modulated by age and sex, but how these factors influence baseline cardiac gene expression remains incompletely understood. Here, we used RNA sequencing and mass spectrometry to compare gene expression in female and male young adult (4 mo) and early aging (20 mo) mouse hearts, identifying thousands of age- and sex-dependent gene expression signatures. Sexually dimorphic cardiac genes are broadly distributed, functioning in mitochondrial metabolism, translation, and other processes. In parallel, we found over 800 genes with differential aging response between male and female, including genes in cAMP and PKA signaling. Analysis of the sex-adjusted aging cardiac transcriptome revealed a widespread remodeling of exon usage patterns that is largely independent from differential gene expression, concomitant with upstream changes in RNA-binding protein and splice factor transcripts. To evaluate the impact of the splicing events on cardiac proteoform composition, we applied an RNA-guided proteomics computational pipeline to analyze the mass spectrometry data and detected hundreds of putative splice variant proteins that have the potential to rewire the cardiac proteome. Taken together, the results here suggest that cardiac aging is associated with 1) widespread sex-biased aging genes and 2) a rewiring of RNA splicing programs, including sex- and age-dependent changes in exon usages and splice patterns that have the potential to influence cardiac protein structure and function. These changes contribute to the emerging evidence for considerable sexual dimorphism in the cardiac aging process that should be considered in the search for disease mechanisms.NEW & NOTEWORTHY Han et al. used proteogenomics to compare male and female mouse hearts at 4 and 20 mo. Sex-biased cardiac genes function in mitochondrial metabolism, translation, autophagy, and other processes. Hundreds of cardiac genes show sex-by-age interactions, that is, sex-biased aging genes. Cardiac aging is accompanied with a remodeling of exon usage in functionally coordinated genes, concomitant with differential expression of RNA-binding proteins and splice factors. These features represent an underinvestigated aspect of cardiac aging that may be relevant to the search for disease mechanisms.

Harnessing gene fusion-derived neoantigens for 'cold' breast and prostate tumor immunotherapy

Breast and prostate cancers are generally considered immunologically 'cold' tumors due to multiple mechanisms rendering them unresponsive to immune checkpoint blockade therapies. With little success in garnering positive outcomes in modern immunotherapeutic clinical trials, it is prudent to re-examine the role of immunogenic neoantigens in these cold tumors. Gene fusions are driver mutations in hormone-driven cancers that can result in alternative mutation-specific neoantigens to promote immunotherapy sensitivity. This review focuses on 1) gene fusion formation mechanisms in neoantigen generation; 2) gene fusion neoantigens in cancer immunotherapeutic strategies and associated clinical trials; and 3) challenges and opportunities in computational and liquid biopsy technologies. This review is anticipated to initiate further research into gene fusion neoantigens of cold tumors for further experimental validation.

Proteotranscriptomics - A facilitator in omics research

Applications in omics research, such as comparative transcriptomics and proteomics, require the knowledge of the species-specific gene sequence and benefit from a comprehensive high-quality annotation of the coding genes to achieve high coverage. While protein-coding genes can in simple cases be detected by scanning the genome for open reading frames, in more complex genomes exonic sequences are separated by introns. Despite advances in sequencing technologies that allow for ever-growing numbers of genomes, the quality of many of the provided genome assemblies do not reach reference quality. These non-contiguous assemblies with gaps and the necessity to predict splice sites limit accurate gene annotation from solely genomic data. In contrast, the transcriptome only contains transcribed gene regions, is devoid of introns and thus provides the optimal basis for the identification of open reading frames. The additional integration of proteomics data to validate predicted protein-coding genes further enriches for accurate gene models. This review outlines the principles of the proteotranscriptomics approach, discusses common challenges and suggests methods for improvement.

Transcriptome innovations in primates revealed by single-molecule long-read sequencing

Transcriptomic diversity greatly contributes to the fundamentals of disease, lineage-specific biology, and environmental adaptation. However, much of the actual isoform repertoire contributing to shaping primate evolution remains unknown. Here, we combined deep long- and short-read sequencing complemented with mass spectrometry proteomics in a panel of lymphoblastoid cell lines (LCLs) from human, three other great apes, and rhesus macaque, producing the largest full-length isoform catalog in primates to date. Around half of the captured isoforms are not annotated in their reference genomes, significantly expanding the gene models in primates. Furthermore, our comparative analyses unveil hundreds of transcriptomic innovations and isoform usage changes related to immune function and immunological disorders. The confluence of these evolutionary innovations with signals of positive selection and their limited impact in the proteome points to changes in alternative splicing in genes involved in immune response as an important target of recent regulatory divergence in primates.

Identification of Protein Isoforms Using Reference Databases Built from Long and Short Read RNA-Sequencing

Alternative splicing can lead to distinct protein isoforms. These can have different functions in specific cells and tissues or in different developmental stages. In this study, we explored whether transcripts assembled from long read, nanopore-based, direct RNA-sequencing (RNA-seq) could improve the identification of protein isoforms in human K562 cells. By comparing with Illumina-based short read RNA-seq, we showed that a large proportion of Ensembl transcripts (5949/14,326) and genes expressing alternatively spliced transcripts (486/2981) identified with long direct reads were missed by short paired-end reads. By co-analyzing proteomic and transcriptomic data, we also showed that some peptides (826/35,976), proteins (262/3215), and protein isoforms arising from distinct transcript variants (574/1212) identified with isoform-specific peptides via custom long-read-based databases were missed in Illumina-derived databases. Finally, we generated unequivocal peptide evidence for a set of protein isoforms and showed that long read, direct RNA-seq allows the discovery of novel protein isoforms not already in reference databases or custom databases built from short read RNA-seq data. Our analysis highlights the benefits of long read RNA-seq data in the generation of reference databases to increase tandem mass spectrometry (MS/MS) identification of protein isoforms.

Abnormal global alternative RNA splicing in COVID-19 patients

Viral infections can alter host transcriptomes by manipulating host splicing machinery. Despite intensive transcriptomic studies on SARS-CoV-2, a systematic analysis of alternative splicing (AS) in severe COVID-19 patients remains largely elusive. Here we integrated proteomic and transcriptomic sequencing data to study AS changes in COVID-19 patients. We discovered that RNA splicing is among the major down-regulated proteomic signatures in COVID-19 patients. The transcriptome analysis showed that SARS-CoV-2 infection induces widespread dysregulation of transcript usage and expression, affecting blood coagulation, neutrophil activation, and cytokine production. Notably, CD74 and LRRFIP1 had increased skipping of an exon in COVID-19 patients that disrupts a functional domain, which correlated with reduced antiviral immunity. Furthermore, the dysregulation of transcripts was strongly correlated with clinical severity of COVID-19, and splice-variants may contribute to unexpected therapeutic activity. In summary, our data highlight that a better understanding of the AS landscape may aid in COVID-19 diagnosis and therapy.

Despite intensive studies on the transcriptional signatures of COVID-19 patients, how SARS-CoV-2 affects AS landscape and the contribution of AS to the pathogenesis of COVID-19 remain largely elusive. By profiling the lung transcriptome and lung proteome of nine patients who died of COVID-19 during the first wave of the pandemic in Wuhan, China, we obtained molecular insights into the AS of cellular transcripts upon SARS-CoV-2 infection. Interestingly, SARS-CoV-2 proteins directly engage host spliceosome to dysregulate essential steps of mature mRNA production and result in widespread dysregulation of cellular function. Taken together, our findings shed light on COVID-19 molecular mechanism and offer potential therapeutic targets for severe COVID-19 disease.

Enhanced protein isoform characterization through long-read proteogenomics

BACKGROUND

The detection of physiologically relevant protein isoforms encoded by the human genome is critical to biomedicine. Mass spectrometry (MS)-based proteomics is the preeminent method for protein detection, but isoform-resolved proteomic analysis relies on accurate reference databases that match the sample; neither a subset nor a superset database is ideal. Long-read RNA sequencing (e.g., PacBio or Oxford Nanopore) provides full-length transcripts which can be used to predict full-length protein isoforms.

RESULTS

We describe here a long-read proteogenomics approach for integrating sample-matched long-read RNA-seq and MS-based proteomics data to enhance isoform characterization. We introduce a classification scheme for protein isoforms, discover novel protein isoforms, and present the first protein inference algorithm for the direct incorporation of long-read transcriptome data to enable detection of protein isoforms previously intractable to MS-based detection. We have released an open-source Nextflow pipeline that integrates long-read sequencing in a proteomic workflow for isoform-resolved analysis.

CONCLUSIONS

Our work suggests that the incorporation of long-read sequencing and proteomic data can facilitate improved characterization of human protein isoform diversity. Our first-generation pipeline provides a strong foundation for future development of long-read proteogenomics and its adoption for both basic and translational research.

Transcription factor protein interactomes reveal genetic determinants in heart disease

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

Splice variants denote differences between a cancer stem cell side population of EWSR1‑ERG‑based Ewing sarcoma cells, its main population and EWSR1‑FLI‑based cells

Ewing sarcoma is a challenging cancer entity, which, besides the characteristic presence of a fusion gene, is driven by multiple alternative splicing events. So far, splice variants in Ewing sarcoma cells were mainly analyzed for EWSR1‑FLI1. The present study provided a comprehensive alternative splicing study on CADO‑ES1, an Ewing model cell line for an EWSR1‑ERG fusion gene. Based on a well‑-characterized RNA‑sequencing dataset with extensive control mechanisms across all levels of analysis, the differential spliced genes in Ewing cancer stem cells were ATP13A3 and EPB41, while the main population was defined by ACADVL, NOP58 and TSPAN3. All alternatively spliced genes were further characterized by their Gene Ontology (GO) terms and by their membership in known protein complexes. These results confirm and extend previous studies towards a systematic whole‑transcriptome analysis. A highlight is the striking segregation of GO terms associated with five basic splice events. This mechanistic insight, together with a coherent integration of all observations with prior knowledge, indicates that EWSR1‑ERG is truly a close twin to EWSR1‑FLI1, but still exhibits certain individuality. Thus, the present study provided a measure of variability in Ewing sarcoma, whose understanding is essential both for clinical procedures and basic mechanistic insight.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high‐throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function.

This article is categorized under:

RNA Processing > Splicing Regulation/Alternative Splicing

Translation > Regulation

RNA Evolution and Genomics > Computational Analyses of RNA

Characterization of protein isoform diversity in human umbilical vein endothelial cells via long-read proteogenomics

Endothelial cells (ECs) comprise the lumenal lining of all blood vessels and are critical for the functioning of the cardiovascular system. Their phenotypes can be modulated by alternative splicing of RNA to produce distinct protein isoforms. To characterize the RNA and protein isoform landscape within ECs, we applied a long read proteogenomics approach to analyse human umbilical vein endothelial cells (HUVECs). Transcripts delineated from PacBio sequencing serve as the basis for a sample-specific protein database used for downstream mass-spectrometry (MS) analysis to infer protein isoform expression. We detected 53,863 transcript isoforms from 10,426 genes, with 22,195 of those transcripts being novel. Furthermore, the predominant isoform in HUVECs does not correspond with the accepted "reference isoform" 25% of the time, with vascular pathway-related genes among this group. We found 2,597 protein isoforms supported through unique peptides, with an additional 2,280 isoforms nominated upon incorporation of long-read transcript evidence. We characterized a novel alternative acceptor for endothelial-related gene CDH5, suggesting potential changes in its associated signalling pathways. Finally, we identified novel protein isoforms arising from a diversity of RNA splicing mechanisms supported by uniquely mapped novel peptides. Our results represent a high-resolution atlas of known and novel isoforms of potential relevance to endothelial phenotypes and function.[Figure: see text].

Transcriptome features of striated muscle aging and predictability of protein level changes

We performed total RNA sequencing and multi-omics analysis comparing skeletal muscle and cardiac muscle in young adult (4 months) vs. early aging (20 months) mice to examine the molecular mechanisms of striated muscle aging. We observed that aging cardiac and skeletal muscles both invoke transcriptomic changes in innate immune system and mitochondria pathways but diverge in extracellular matrix processes. On an individual gene level, we identified 611 age-associated signatures in skeletal and cardiac muscles, including a number of myokine and cardiokine encoding genes. Because RNA and protein levels correlate only partially, we reason that differentially expressed transcripts that accurately reflect their protein counterparts will be more valuable proxies for proteomic changes and by extension physiological states. We applied a computational data analysis workflow to estimate which transcriptomic changes are more likely relevant to protein-level regulation using large proteogenomics data sets. We estimate about 48% of the aging-associated transcripts predict protein levels well (r ≥ 0.5). In parallel, a comparison of the identified aging-regulated genes with public human transcriptomics data showed that only 35-45% of the identified genes show an age-dependent expression in corresponding human tissues. Thus, integrating both RNA-protein correlation and human conservation across data sources, we nominate 134 prioritized aging striated muscle signatures that are predicted to correlate strongly with protein levels and that show age-dependent expression in humans. The results here reveal new details into how aging reshapes gene expression in striated muscles at the transcript and protein levels.

Bottom-up proteomic analysis of human adult cardiac tissue and isolated cardiomyocytes

The heart is composed of multiple cell types, each with a specific function. Cell-type-specific approaches are necessary for defining the intricate molecular mechanisms underlying cardiac development, homeostasis, and pathology. While single-cell RNA-seq studies are beginning to define the chamber-specific cellular composition of the heart, our views of the proteome are more limited because most proteomics studies have utilized homogenized human cardiac tissue. To promote future cell-type specific analyses of the human heart, we describe the first method for cardiomyocyte isolation from cryopreserved human cardiac tissue followed by flow cytometry for purity assessment. We also describe a facile method for preparing isolated cardiomyocytes and whole cardiac tissue homogenate for bottom-up proteomic analyses. Prior experience in dissociating cardiac tissue or proteomics is not required to execute these methods. We compare different sample preparation workflows and analysis methods to demonstrate how these can impact the depth of proteome coverage achieved. We expect this how-to guide will serve as a starting point for investigators interested in general and cell-type-specific views of the cardiac proteome.

A novel XBP1 variant is highly enriched in cancer tissues and is specifically required for cancer cell survival

XBP1 is a basic leucine zipper (bZIP) transcription factor and a key mediator of the endoplasmic reticulum (ER) stress-activated unfolded protein response (UPR). XBP1-mediated transcription facilitates cell adaptation to ER stress and also promotes tumor progression, while suppressing anti-tumor immunity. Here we report a novel XBP1 variant, namely XBP1 variant 1 (XBP1v1, Xv1 for short), that is specifically required for survival of cancer cells. Xv1 contains a cryptic first exon that is conserved only in humans and great apes. Comparing to XBP1, Xv1 encodes a protein with a different N-terminal sequence containing 25 amino acids. Analysis of RNAseq database reveals that Xv1 is broadly expressed across cancer types but almost none in normal tissues. Elevated Xv1 expression is associated with poor survival of patients with several types of cancer. Knockdown of Xv1 induces death of multiple cancer cell lines but has little effect on non-cancerous cells in vitro. Moreover, knockdown of Xv1 also inhibits growth of a xenograft breast tumor in mice. Together, our results indicate that Xv1 is essential for survival of cancer cells.

Importance of evaluating protein glycosylation in pluripotent stem cell-derived cardiomyocytes for research and clinical applications

Proper protein glycosylation is critical to normal cardiomyocyte physiology. Aberrant glycosylation can alter protein localization, structure, drug interactions, and cellular function. The in vitro differentiation of human pluripotent stem cells into cardiomyocytes (hPSC-CM) has become increasingly important to the study of protein function and to the fields of cardiac disease modeling, drug testing, drug discovery, and regenerative medicine. Here, we offer our perspective on the importance of protein glycosylation in hPSC-CM. Protein glycosylation is dynamic in hPSC-CM, but the timing and extent of glycosylation are still poorly defined. We provide new data highlighting how observed changes in hPSC-CM glycosylation may be caused by underlying differences in the protein or transcript abundance of enzymes involved in building and trimming the glycan structures or glycoprotein gene products. We also provide evidence that alternative splicing results in altered sites of glycosylation within the protein sequence. Our findings suggest the need to precisely define protein glycosylation events that may have a critical impact on the function and maturation state of hPSC-CM. Finally, we provide an overview of analytical strategies available for studying protein glycosylation and identify opportunities for the development of new bioinformatic approaches to integrate diverse protein glycosylation data types. We predict that these tools will promote the accurate assessment of protein glycosylation in future studies of hPSC-CM that will ultimately be of significant experimental and clinical benefit.

Computation-assisted targeted proteomics of alternative splicing protein isoforms in the human heart

Alternative splicing is prevalent in the heart and implicated in many cardiovascular diseases, but not every alternative transcript is translated and detecting non-canonical isoforms at the protein level remains challenging. Here we show the use of a computation-assisted targeted proteomics workflow to detect protein alternative isoforms in the human heart. We build on a recent strategy to integrate deep RNA-seq and large-scale mass spectrometry data to identify candidate translated isoform peptides. A machine learning approach is then applied to predict their fragmentation patterns and design protein isoform-specific parallel reaction monitoring detection (PRM) assays. As proof-of-principle, we built PRM assays for 29 non-canonical isoform peptides and detected 22 peptides in a human heart lysate. The predictions-aided PRM assays closely mirrored synthetic peptide standards for non-canonical sequences. This approach may be useful for validating non-canonical protein identification and discovering functionally relevant isoforms in the heart.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Peptide Correlation Analysis (PeCorA) Reveals Differential Proteoform Regulation

Shotgun proteomics techniques infer the presence and quantity of proteins using peptide proxies produced by cleavage of the proteome with a protease. Most protein quantitation strategies assume that multiple peptides derived from a protein will behave quantitatively similar across treatment groups, but this assumption may be false due to (1) heterogeneous proteoforms and (2) technical artifacts. Here we describe a strategy called peptide correlation analysis (PeCorA) that detects quantitative disagreements between peptides mapped to the same protein. PeCorA fits linear models to assess whether a peptide's change across treatment groups differs from all other peptides assigned to the same protein. PeCorA revealed that ∼15% of proteins in a mouse microglia stress data set contain at least one discordant peptide. Inspection of the discordant peptides shows the utility of PeCorA for the direct and indirect detection of regulated post-translational modifications (PTMs) and also for the discovery of poorly quantified peptides. The exclusion of poorly quantified peptides before protein quantity summarization decreased false-positives in a benchmark data set. Finally, PeCorA suggests that the inactive isoform of prothrombin, a coagulation cascade protease, is more abundant in plasma from COVID-19 patients relative to non-COVID-19 controls. PeCorA is freely available as an R package that works with arbitrary tables of quantified peptides.

Deep Learning in Proteomics

Proteomics, the study of all the proteins in biological systems, is becoming a data-rich science. Protein sequences and structures are comprehensively catalogued in online databases. With recent advancements in tandem mass spectrometry (MS) technology, protein expression and post-translational modifications (PTMs) can be studied in a variety of biological systems at the global scale. Sophisticated computational algorithms are needed to translate the vast amount of data into novel biological insights. Deep learning automatically extracts data representations at high levels of abstraction from data, and it thrives in data-rich scientific research domains. Here, a comprehensive overview of deep learning applications in proteomics, including retention time prediction, MS/MS spectrum prediction, de novo peptide sequencing, PTM prediction, major histocompatibility complex-peptide binding prediction, and protein structure prediction, is provided. Limitations and the future directions of deep learning in proteomics are also discussed. This review will provide readers an overview of deep learning and how it can be used to analyze proteomics data.

Proteomic signatures of acute oxidative stress response to paraquat in the mouse heart

The heart is sensitive to oxidative damage but a global view on how the cardiac proteome responds to oxidative stressors has yet to fully emerge. Using quantitative tandem mass spectrometry, we assessed the effects of acute exposure of the oxidative stress inducer paraquat on protein expression in mouse hearts. We observed widespread protein expression changes in the paraquat-exposed heart especially in organelle-containing subcellular fractions. During cardiac response to acute oxidative stress, proteome changes are consistent with a rapid reduction of mitochondrial metabolism, coupled with activation of multiple antioxidant proteins, reduction of protein synthesis and remediation of proteostasis. In addition to differential expression, we saw evidence of spatial reorganizations of the cardiac proteome including the translocation of hexokinase 2 to more soluble fractions. Treatment with the antioxidants Tempol and MitoTEMPO reversed many proteomic signatures of paraquat but this reversal was incomplete. We also identified a number of proteins with unknown function in the heart to be triggered by paraquat, suggesting they may have functions in oxidative stress response. Surprisingly, protein expression changes in the heart correlate poorly with those in the lung, consistent with differential sensitivity or stress response in these two organs. The results and data set here could provide insights into oxidative stress responses in the heart and avail the search for new therapeutic targets.

Determining Alternative Protein Isoform Expression Using RNA Sequencing and Mass Spectrometry

Alternative splicing greatly expands the coding capacity of the human genome, but how many alternative transcripts are translated as proteins or carry functional importance remains unknown and awaits experimental verification. Here, we describe a protocol that combines transcriptomics (RNA-seq) and proteomics (mass spectrometry [MS]) analyses to identify alternative isoforms in proteomes. This workflow is applicable to custom-generated RNA-seq and MS data from matching samples, as well as the reanalysis of existing transcriptomics and proteomics datasets in public repositories. For complete details on the use and execution of this protocol, please refer to Lau et al. (2019).

A high-stringency blueprint of the human proteome

The Human Proteome Organization (HUPO) launched the Human Proteome Project (HPP) in 2010, creating an international framework for global collaboration, data sharing, quality assurance and enhancing accurate annotation of the genome-encoded proteome. During the subsequent decade, the HPP established collaborations, developed guidelines and metrics, and undertook reanalysis of previously deposited community data, continuously increasing the coverage of the human proteome. On the occasion of the HPP's tenth anniversary, we here report a 90.4% complete high-stringency human proteome blueprint. This knowledge is essential for discerning molecular processes in health and disease, as we demonstrate by highlighting potential roles the human proteome plays in our understanding, diagnosis and treatment of cancers, cardiovascular and infectious diseases.

An analysis of tissue-specific alternative splicing at the protein level

The role of alternative splicing is one of the great unanswered questions in cellular biology. There is strong evidence for alternative splicing at the transcript level, and transcriptomics experiments show that many splice events are tissue specific. It has been suggested that alternative splicing evolved in order to remodel tissue-specific protein-protein networks. Here we investigated the evidence for tissue-specific splicing among splice isoforms detected in a large-scale proteomics analysis. Although the data supporting alternative splicing is limited at the protein level, clear patterns emerged among the small numbers of alternative splice events that we could detect in the proteomics data. More than a third of these splice events were tissue-specific and most were ancient: over 95% of splice events that were tissue-specific in both proteomics and RNAseq analyses evolved prior to the ancestors of lobe-finned fish, at least 400 million years ago. By way of contrast, three in four alternative exons in the human gene set arose in the primate lineage, so our results cannot be extrapolated to the whole genome. Tissue-specific alternative protein forms in the proteomics analysis were particularly abundant in nervous and muscle tissues and their genes had roles related to the cytoskeleton and either the structure of muscle fibres or cell-cell connections. Our results suggest that this conserved tissue-specific alternative splicing may have played a role in the development of the vertebrate brain and heart.

We manually curated a set of 255 splice events detected in a large-scale tissue-based proteomics experiment and found that more than a third had evidence of significant tissue-specific differences. Events that were significantly tissue-specific at the protein level were highly conserved; almost 75% evolved over 400 million years ago. The tissues in which we found most evidence for tissue-specific splicing were nervous tissues and cardiac tissues. Genes with tissue-specific events in these two tissues had functions related to important cellular structures in brain and heart tissues. These splice events may have been essential for the development of vertebrate heart and muscle. However, our data set may not be representative of alternative exons as a whole. We found that most tissue specific splicing was strongly conserved, but just 5% of annotated alternative exons in the human gene set are ancient. More than three quarters of alternative exons are primate-derived. Although the analysis does not provide a definitive answer to the question of the functional role of alternative splicing, our results do indicate that alternative splice variants may have played a significant part in the evolution of brain and heart tissues in vertebrates.

Distinct hypertrophic cardiomyopathy genotypes result in convergent sarcomeric proteoform profiles revealed by top-down proteomics

Hypertrophic cardiomyopathy (HCM) is the most common heritable heart disease. Although the genetic cause of HCM has been linked to mutations in genes encoding sarcomeric proteins, the ability to predict clinical outcomes based on specific mutations in HCM patients is limited. Moreover, how mutations in different sarcomeric proteins can result in highly similar clinical phenotypes remains unknown. Posttranslational modifications (PTMs) and alternative splicing regulate the function of sarcomeric proteins; hence, it is critical to study HCM at the level of proteoforms to gain insights into the mechanisms underlying HCM. Herein, we employed high-resolution mass spectrometry-based top-down proteomics to comprehensively characterize sarcomeric proteoforms in septal myectomy tissues from HCM patients exhibiting severe outflow track obstruction (n = 16) compared to nonfailing donor hearts (n = 16). We observed a complex landscape of sarcomeric proteoforms arising from combinatorial PTMs, alternative splicing, and genetic variation in HCM. A coordinated decrease of phosphorylation in important myofilament and Z-disk proteins with a linear correlation suggests PTM cross-talk in the sarcomere and dysregulation of protein kinase A pathways in HCM. Strikingly, we discovered that the sarcomeric proteoform alterations in the myocardium of HCM patients undergoing septal myectomy were remarkably consistent, regardless of the underlying HCM-causing mutations. This study suggests that the manifestation of severe HCM coalesces at the proteoform level despite distinct genotype, which underscores the importance of molecular characterization of HCM phenotype and presents an opportunity to identify broad-spectrum treatments to mitigate the most severe manifestations of this genetically heterogenous disease.

Identification of modified peptides using localization-aware open search

Identification of post-translationally or chemically modified peptides in mass spectrometry-based proteomics experiments is a crucial yet challenging task. We have recently introduced a fragment ion indexing method and the MSFragger search engine to empower an open search strategy for comprehensive analysis of modified peptides. However, this strategy does not consider fragment ions shifted by unknown modifications, preventing modification localization and limiting the sensitivity of the search. Here we present a localization-aware open search method, in which both modification-containing (shifted) and regular fragment ions are indexed and used in scoring. We also implement a fast mass calibration and optimization method, allowing optimization of the mass tolerances and other key search parameters. We demonstrate that MSFragger with mass calibration and localization-aware open search identifies modified peptides with significantly higher sensitivity and accuracy. Comparing MSFragger to other modification-focused tools (pFind3, MetaMorpheus, and TagGraph) shows that MSFragger remains an excellent option for fast, comprehensive, and sensitive searches for modified peptides in shotgun proteomics data.