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.

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.

Simultaneous proteome localization and turnover analysis reveals spatiotemporal features of protein homeostasis disruptions.

The spatial and temporal distributions of proteins are critical to protein function, but cannot be directly assessed by measuring protein bundance. Here we describe a mass spectrometry-based proteomics strategy, Simultaneous Proteome Localization and Turnover (SPLAT), to measure concurrently protein turnover rates and subcellular localization in the same experiment. Applying the method, we find that unfolded protein response (UPR) has different effects on protein turnover dependent on their subcellular location in human AC16 cells, with proteome-wide slowdown but acceleration among stress response proteins in the ER and Golgi. In parallel, UPR triggers broad differential localization of proteins including RNA-binding proteins and amino acid transporters. Moreover, we observe newly synthesized proteins including EGFR that show a differential localization under stress than the existing protein pools, reminiscent of protein trafficking disruptions. We next applied SPLAT to an induced pluripotent stem cell derived cardiomyocyte (iPSC-CM) model of cancer drug cardiotoxicity upon treatment with the proteasome inhibitor carfilzomib. Paradoxically, carfilzomib has little effect on global average protein half-life, but may instead selectively disrupt sarcomere protein homeostasis. This study provides a view into the interactions of protein spatial and temporal dynamics and demonstrates a method to examine protein homeostasis regulations in stress and drug response.

Simultaneous proteome localization and turnover analysis reveals spatiotemporal features of protein homeostasis disruptions.

The functions of proteins depend on their spatial and temporal distributions, which are not directly measured by static protein abundance. Under endoplasmic reticulum (ER) stress, the unfolded protein response (UPR) pathway remediates proteostasis in part by altering the turnover kinetics and spatial distribution of proteins. A global view of these spatiotemporal changes has yet to emerge and it is unknown how they affect different cellular compartments and pathways. Here we describe a mass spectrometry-based proteomics strategy and data analysis pipeline, termed Simultaneous Proteome Localization and Turnover (SPLAT), to measure concurrently the changes in protein turnover and subcellular distribution in the same experiment. Investigating two common UPR models of thapsigargin and tunicamycin challenge in human AC16 cells, we find that the changes in protein turnover kinetics during UPR varies across subcellular localizations, with overall slowdown but an acceleration in endoplasmic reticulum and Golgi proteins involved in stress response. In parallel, the spatial proteomics component of the experiment revealed an externalization of amino acid transporters and ion channels under UPR, as well as the migration of RNA-binding proteins toward an endosome co-sedimenting compartment. The SPLAT experimental design classifies heavy and light SILAC labeled proteins separately, allowing the observation of differential localization of new and old protein pools and capturing a partition of newly synthesized EGFR and ITGAV to the ER under stress that suggests protein trafficking disruptions. Finally, application of SPLAT toward human induced pluripotent stem cell derived cardiomyocytes (iPSC-CM) exposed to the cancer drug carfilzomib, identified a selective disruption of proteostasis in sarcomeric proteins as a potential mechanism of carfilzomib-mediated cardiotoxicity. Taken together, this study provides a global view into the spatiotemporal dynamics of human cardiac cells and demonstrates a method for inferring the coordinations between spatial and temporal proteome regulations in stress and drug response.

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.

Senescence in human AC16 cardiac cells is associated with thymidine kinase induction and histone loss.

AC16 cells are a transformed human cardiac cell line commonly used to study cardiomyocyte biology. We show that reduced proliferation and senescence markers can be robustly induced in AC16 cells cultured in low serum condition and treated with (i) low-dose doxorubicin, (ii) UV 254 nm, or (iii) H 2 O 2 exposure for up to 48 hours. Increased p21 (CDKN1A) and H2A.X variant histone (H2AX) levels serve as reliable molecular markers upon all three treatment conditions, but the up-regulation of another common senescence marker, p16 (CDKN2A) was not observed. A proteomics screen further shows that the loss of histones and the increased expression of thymidine kinases (TK1) are prominent features of AC16 cells under doxorubicin induced senescence.

Co-culture with normal and senescent AC16 cardiomyocytes modulates fibrotic response in primary human ventricular fibroblasts.

Transwell co-culture with human AC16 cardiomyocyte-like cells modifies the response of primary human ventricular fibroblasts to TGF-ß stimulation. Fibrotic response markers including collagen I (COL1A1) and ɑ-smooth muscle actin (ACTA2) are amplified in the presence of AC16 cells, whereas others including periostin (POSTN) and fibronectin (FN1) are suppressed. Similar modulation is observed when the ventricular fibroblasts are co-cultured with AC16 cells under baseline and induced senescence conditions. Given that the response to TGF-ß stimulation is commonly measured to study fibrotic signaling and drug treatments in vitro, the results here suggest that the effect of cellular crosstalk should be more broadly considered.

Protein prediction models support widespread post-transcriptional regulation of protein abundance by interacting partners.

Protein and mRNA levels correlate only moderately. The availability of proteogenomics data sets with protein and transcript measurements from matching samples is providing new opportunities to assess the degree to which protein levels in a system can be predicted from mRNA information. Here we examined the contributions of input features in protein abundance prediction models. Using large proteogenomics data from 8 cancer types within the Clinical Proteomic Tumor Analysis Consortium (CPTAC) data set, we trained models to predict the abundance of over 13,000 proteins using matching transcriptome data from up to 958 tumor or normal adjacent tissue samples each, and compared predictive performances across algorithms, data set sizes, and input features. Over one-third of proteins (4,648) showed relatively poor predictability (elastic net r ≤ 0.3) from their cognate transcripts. Moreover, we found widespread occurrences where the abundance of a protein is considerably less well explained by its own cognate transcript level than that of one or more trans locus transcripts. The incorporation of additional trans-locus transcript abundance data as input features increasingly improved the ability to predict sample protein abundance. Transcripts that contribute to non-cognate protein abundance primarily involve those encoding known or predicted interaction partners of the protein of interest, including not only large multi-protein complexes as previously shown, but also small stable complexes in the proteome with only one or few stable interacting partners. Network analysis further shows a complex proteome-wide interdependency of protein abundance on the transcript levels of multiple interacting partners. The predictive model analysis here therefore supports that protein-protein interaction including in small protein complexes exert post-transcriptional influence on proteome compositions more broadly than previously recognized. Moreover, the results suggest mRNA and protein co-expression analysis may have utility for finding gene interactions and predicting expression changes in biological systems.

Sex-specific responses to slow progressive pressure overload in a large animal model of HFpEF.

Approximately 50% of all heart failure (HF) diagnoses can be classified as HF with preserved ejection fraction (HFpEF). HFpEF is more prevalent in females compared with males, but the underlying mechanisms are unknown. We previously showed that pressure overload (PO) in male felines induces a cardiopulmonary phenotype with essential features of human HFpEF. The goal of this study was to determine if slow progressive PO induces distinct cardiopulmonary phenotypes in females and males in the absence of other pathological stressors. Female and male felines underwent aortic constriction (banding) or sham surgery after baseline echocardiography, pulmonary function testing, and blood sampling. These assessments were repeated at 2 and 4 mo postsurgery to document the effects of slow progressive pressure overload. At 4 mo, invasive hemodynamic studies were also performed. Left ventricle (LV) tissue was collected for histology, myofibril mechanics, extracellular matrix (ECM) mass spectrometry, and single-nucleus RNA sequencing (snRNAseq). The induced pressure overload (PO) was not different between sexes. PO also induced comparable changes in LV wall thickness and myocyte cross-sectional area in both sexes. Both sexes had preserved ejection fraction, but males had a slightly more robust phenotype in hemodynamic and pulmonary parameters. There was no difference in LV fibrosis and ECM composition between banded male and female animals. LV snRNAseq revealed changes in gene programs of individual cell types unique to males and females after PO. Based on these results, both sexes develop cardiopulmonary dysfunction but the phenotype is somewhat less advanced in females.NEW & NOTEWORTHY We performed a comprehensive assessment to evaluate the effects of slow progressive pressure overload on cardiopulmonary function in a large animal model of heart failure with preserved ejection fraction (HFpEF) in males and females. Functional and structural assessments were performed at the organ, tissue, cellular, protein, and transcriptional levels. This is the first study to compare snRNAseq and ECM mass spectrometry of HFpEF myocardium from males and females. The results broaden our understanding of the pathophysiological response of both sexes to pressure overload. Both sexes developed a robust cardiopulmonary phenotype, but the phenotype was equal or a bit less robust in females.

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.

Shotgun Proteomics Sample Processing Automated by an Open-Source Lab Robot.

Mass spectrometry-based shotgun proteomics experiments require multiple sample preparation steps, including enzymatic protein digestion and clean-up, which can take up significant person-hours of bench labor and present a source of batch-to-batch variability. Lab automation with pipetting robots can reduce manual work, maximize throughput, and increase research reproducibility. Still, the steep starting prices of standard automation stations make them unaffordable for many academic laboratories. This article describes a proteomics sample preparation workflow using an affordable, open-source automation system (The Opentrons OT-2), including instructions for setting up semi-automated protein reduction, alkylation, digestion, and clean-up steps; as well as accompanying open-source Python scripts to program the OT-2 system through its application programming interface.

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.

HDAC Inhibition Reverses Preexisting Diastolic Dysfunction and Blocks Covert Extracellular Matrix Remodeling.

BACKGROUND: Diastolic dysfunction (DD) is associated with the development of heart failure and contributes to the pathogenesis of other cardiac maladies, including atrial fibrillation. Inhibition of histone deacetylases (HDACs) has been shown to prevent DD by enhancing myofibril relaxation. We addressed the therapeutic potential of HDAC inhibition in a model of established DD with preserved ejection fraction. METHODS: Four weeks after uninephrectomy and implantation with deoxycorticosterone acetate pellets, when DD was clearly evident, 1 cohort of mice was administered the clinical-stage HDAC inhibitor ITF2357/Givinostat. Echocardiography, blood pressure measurements, and end point invasive hemodynamic analyses were performed. Myofibril mechanics and intact cardiomyocyte relaxation were assessed ex vivo. Cardiac fibrosis was evaluated by picrosirius red staining and second harmonic generation microscopy of left ventricle (LV) sections, RNA sequencing of LV mRNA, mass spectrometry-based evaluation of decellularized LV biopsies, and atomic force microscopy determination of LV stiffness. Mechanistic studies were performed with primary rat and human cardiac fibroblasts. RESULTS: HDAC inhibition normalized DD without lowering blood pressure in this model of systemic hypertension. In contrast to previous models, myofibril relaxation was unimpaired in uninephrectomy/deoxycorticosterone acetate mice. Furthermore, cardiac fibrosis was not evident in any mouse cohort on the basis of picrosirius red staining or second harmonic generation microscopy. However, mass spectrometry revealed induction in the expression of >100 extracellular matrix proteins in LVs of uninephrectomy/deoxycorticosterone acetate mice, which correlated with profound tissue stiffening based on atomic force microscopy. ITF2357/Givinostat treatment blocked extracellular matrix expansion and LV stiffening. The HDAC inhibitor was subsequently shown to suppress cardiac fibroblast activation, at least in part, by blunting recruitment of the profibrotic chromatin reader protein BRD4 (bromodomain-containing protein 4) to key gene regulatory elements. CONCLUSIONS: These findings demonstrate the potential of HDAC inhibition as a therapeutic intervention to reverse existing DD and establish blockade of extracellular matrix remodeling as a second mechanism by which HDAC inhibitors improve ventricular filling. Our data reveal the existence of pathophysiologically relevant covert or hidden cardiac fibrosis that is below the limit of detection of histochemical stains such as picrosirius red, highlighting the need to evaluate fibrosis of the heart using diverse methodologies.

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.

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.

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.

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

The protein-level translational status and function of many alternative splicing events remain poorly understood. We use an RNA sequencing (RNA-seq)-guided proteomics method to identify protein alternative splicing isoforms in the human proteome by constructing tissue-specific protein databases that prioritize transcript splice junction pairs with high translational potential. Using the custom databases to reanalyze ∼80 million mass spectra in public proteomics datasets, we identify more than 1,500 noncanonical protein isoforms across 12 human tissues, including ∼400 sequences undocumented on TrEMBL and RefSeq databases. We apply the method to original quantitative mass spectrometry experiments and observe widespread isoform regulation during human induced pluripotent stem cell cardiomyocyte differentiation. On a proteome scale, alternative isoform regions overlap frequently with disordered sequences and post-translational modification sites, suggesting that alternative splicing may regulate protein function through modulating intrinsically disordered regions. The described approach may help elucidate functional consequences of alternative splicing and expand the scope of proteomics investigations in various systems.

Mining the Proteome Associated with Rheumatic and Autoimmune Diseases.

A steady increase in the incidence of osteoarthritis and other rheumatic diseases has been observed in recent decades, including autoimmune conditions such as rheumatoid arthritis, spondyloarthropathies, systemic lupus erythematosus, systemic sclerosis, and Sjögren's syndrome. Rheumatic and autoimmune diseases (RADs) are characterized by the inflammation of joints, muscles, or other connective tissues. In addition to often experiencing debilitating mobility and pain, RAD patients are also at a higher risk of suffering comorbidities such as cardiovascular or infectious events. Given the socioeconomic impact of RADs, broad research efforts have been dedicated to these diseases worldwide. In the present work, we applied literature mining platforms to identify "popular" proteins closely related to RADs. The platform is based on publicly available literature. The results not only will enable the systematic prioritization of candidates to perform targeted proteomics studies but also may lead to a greater insight into the key pathogenic processes of these disorders.