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.

Splicing neoantigen discovery with SNAF reveals shared targets for cancer immunotherapy.

Immunotherapy has emerged as a crucial strategy to combat cancer by "reprogramming" a patient's own immune system. Although immunotherapy is typically reserved for patients with a high mutational burden, neoantigens produced from posttranscriptional regulation may provide an untapped reservoir of common immunogenic targets for new targeted therapies. To comprehensively define tumor-specific and likely immunogenic neoantigens from patient RNA-Seq, we developed Splicing Neo Antigen Finder (SNAF), an easy-to-use and open-source computational workflow to predict splicing-derived immunogenic MHC-bound peptides (T cell antigen) and unannotated transmembrane proteins with altered extracellular epitopes (B cell antigen). This workflow uses a highly accurate deep learning strategy for immunogenicity prediction (DeepImmuno) in conjunction with new algorithms to rank the tumor specificity of neoantigens (BayesTS) and to predict regulators of mis-splicing (RNA-SPRINT). T cell antigens from SNAF were frequently evidenced as HLA-presented peptides from mass spectrometry (MS) and predict response to immunotherapy in melanoma. Splicing neoantigen burden was attributed to coordinated splicing factor dysregulation. Shared splicing neoantigens were found in up to 90% of patients with melanoma, correlated to overall survival in multiple cancer cohorts, induced T cell reactivity, and were characterized by distinct cells of origin and amino acid preferences. In addition to T cell neoantigens, our B cell focused pipeline (SNAF-B) identified a new class of tumor-specific extracellular neoepitopes, which we termed ExNeoEpitopes. ExNeoEpitope full-length mRNA predictions were tumor specific and were validated using long-read isoform sequencing and in vitro transmembrane localization assays. Therefore, our systematic identification of splicing neoantigens revealed potential shared targets for therapy in heterogeneous cancers.

Connexin 37 sequestering of activated-ERK in the cytoplasm promotes p27-mediated endothelial cell cycle arrest.

Connexin37-mediated regulation of cell cycle modulators and, consequently, growth arrest lack mechanistic understanding. We previously showed that arterial shear stress up-regulates Cx37 in endothelial cells and activates a Notch/Cx37/p27 signaling axis to promote G1 cell cycle arrest, and this is required to enable arterial gene expression. However, how induced expression of a gap junction protein, Cx37, up-regulates cyclin-dependent kinase inhibitor p27 to enable endothelial growth suppression and arterial specification is unclear. Herein, we fill this knowledge gap by expressing wild-type and regulatory domain mutants of Cx37 in cultured endothelial cells expressing the Fucci cell cycle reporter. We determined that both the channel-forming and cytoplasmic tail domains of Cx37 are required for p27 up-regulation and late G1 arrest. Mechanistically, the cytoplasmic tail domain of Cx37 interacts with, and sequesters, activated ERK in the cytoplasm. This then stabilizes pERK nuclear target Foxo3a, which up-regulates p27 transcription. Consistent with previous studies, we found this Cx37/pERK/Foxo3a/p27 signaling axis functions downstream of arterial shear stress to promote endothelial late G1 state and enable up-regulation of arterial genes.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease.

A major fraction of loci identified by genome-wide association studies (GWASs) lead to alterations in alternative splicing, but interpretation of how such alterations impact proteins is hindered by the technical limitations of short-read RNA-seq, which cannot directly link splicing events to full-length transcript or protein isoforms. Long-read RNA-seq represents a powerful tool to define and quantify transcript isoforms, and recently, infer protein isoform existence. Here we present a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H 4 PP ≥ 0.75). We generated deep coverage PacBio long-read RNA-seq data (N=∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were novel. By casting the colocalized sQTLs directly onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Using these data, we created one of the first proteome-scale resources defining full-length isoforms impacted by colocalized sQTLs. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense mediated decay (NMD) and 190 that potentially resulted in the expression of new protein isoforms. Finally, we identified colocalizing sQTLs in TPM2 for splice junctions between two mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization. We expect our approach to be widely generalizable across diverse clinical traits and accelerate system-scale analyses of protein isoform activities modulated by GWAS loci.

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].

Endothelial cell cycle state determines propensity for arterial-venous fate.

During blood vessel development, endothelial cells become specified toward arterial or venous fates to generate a circulatory network that provides nutrients and oxygen to, and removes metabolic waste from, all tissues. Arterial-venous specification occurs in conjunction with suppression of endothelial cell cycle progression; however, the mechanistic role of cell cycle state is unknown. Herein, using Cdh5-CreERT2;R26FUCCI2aR reporter mice, we find that venous endothelial cells are enriched for the FUCCI-Negative state (early G1) and BMP signaling, while arterial endothelial cells are enriched for the FUCCI-Red state (late G1) and TGF-ß signaling. Furthermore, early G1 state is essential for BMP4-induced venous gene expression, whereas late G1 state is essential for TGF-ß1-induced arterial gene expression. Pharmacologically induced cell cycle arrest prevents arterial-venous specification defects in mice with endothelial hyperproliferation. Collectively, our results show that distinct endothelial cell cycle states provide distinct windows of opportunity for the molecular induction of arterial vs. venous fate.

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.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing.

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

Characteristics of Immune Memory and Effector Activity to Cancer-Expressed MHC Class I Phosphopeptides Differ in Healthy Donors and Ovarian Cancer Patients.

Elevated immunity to cancer-expressed antigens can be detected in people with no history of cancer and may contribute to cancer prevention. We have previously reported that MHC-restricted phosphopeptides are cancer-expressed antigens and targets of immune recognition. However, the extent to which this immunity reflects prior or ongoing phosphopeptide exposures was not investigated. In this study, we found that preexisting immune memory to cancer-expressed phosphopeptides was evident in most healthy donors, but the breadth among donors was highly variable. Although three phosphopeptides were recognized by most donors, suggesting exposures to common microbial/infectious agents, most of the 205 tested phosphopeptides were not recognized by peripheral blood mononuclear cells (PBMC) from any donor and the remainder were recognized by only 1 to 3 donors. In longitudinal analyses of 2 donors, effector immune response profiles suggested active reexposures to a subset of phosphopeptides. These findings suggest that the immunogens generating most phosphopeptide-specific immune memory are rare infectious agents or incipient cancer cells with distinct phosphoproteome dysregulations, and that repetitive immunogenic exposures occur in individual donors. Phosphopeptide-specific immunity in PBMCs and tumor-infiltrating lymphocytes from ovarian cancer patients was limited, regardless of whether the phosphopeptide was expressed on the tumor. However, 4 of 10 patients responded to 1 to 2 immunodominant phosphopeptides, and 1 showed an elevated effector response to a tumor-expressed phosphopeptide. As the tumors from these patients displayed many phosphopeptides, these data are consistent with lack of prior exposure or impaired ability to respond to some phosphopeptides and suggest that enhancing phosphopeptide-specific T-cell responses could be a useful approach to improve tumor immunotherapy.

The proteome of mouse vestibular hair bundles over development.

Development of the vertebrate hair bundle is a precisely orchestrated event that culminates in production of a tightly ordered arrangement of actin-rich stereocilia and a single axonemal kinocilium. To understand how the protein composition of the bundle changes during development, we isolated bundles from young (postnatal days P4-P6) and mature (P21-P25) mouse utricles using the twist-off method, then characterized their constituent proteins using liquid-chromatography tandem mass spectrometry with data-dependent acquisition. Using MaxQuant and label-free quantitation, we measured relative abundances of proteins in both bundles and in the whole utricle; comparison of protein abundance between the two fractions allows calculation of enrichment in bundles. These data, which are available via ProteomeXchange with identifier PXD002167, will be useful for examining the proteins present in mammalian vestibular bundles and how their concentrations change over development.

Accurate label-free protein quantitation with high- and low-resolution mass spectrometers.

Label-free quantitation of proteins analyzed by tandem mass spectrometry uses either integrated peak intensity from the parent-ion mass analysis (MS1) or features from fragment-ion analysis (MS2), such as spectral counts or summed fragment-ion intensity. We directly compared MS1 and MS2 quantitation by analyzing human protein standards diluted into Escherichia coli extracts on an Orbitrap mass spectrometer. We found that summed MS2 intensities were nearly as accurate as integrated MS1 intensities, and both outperformed MS2 spectral counting in accuracy and linearity. We compared these results to those obtained from two low-resolution ion-trap mass spectrometers; summed MS2 intensities from LTQ and LTQ Velos instruments were similar in accuracy to those from the Orbitrap. Data from all three instruments are available via ProteomeXchange with identifier PXD000602. Abundance measurements using MS1 or MS2 intensities had limitations, however. While measured protein concentration was on average well-correlated with the known concentration, there was considerable protein-to-protein variation. Moreover, not all human proteins diluted to a mole fraction of 10(-3) or lower were detected, with a strong falloff below 10(-4) mole fraction. These results show that MS1 and MS2 intensities are simple measures of protein abundance that are on average accurate but should be limited to quantitation of proteins of intermediate to higher fractional abundance.

Molecular architecture of the chick vestibular hair bundle.

Hair bundles of the inner ear have a specialized structure and protein composition that underlies their sensitivity to mechanical stimulation. Using mass spectrometry, we identified and quantified >1,100 proteins, present from a few to 400,000 copies per stereocilium, from purified chick bundles; 336 of these were significantly enriched in bundles. Bundle proteins that we detected have been shown to regulate cytoskeleton structure and dynamics, energy metabolism, phospholipid synthesis and cell signaling. Three-dimensional imaging using electron tomography allowed us to count the number of actin-actin cross-linkers and actin-membrane connectors; these values compared well to those obtained from mass spectrometry. Network analysis revealed several hub proteins, including RDX (radixin) and SLC9A3R2 (NHERF2), which interact with many bundle proteins and may perform functions essential for bundle structure and function. The quantitative mass spectrometry of bundle proteins reported here establishes a framework for future characterization of dynamic processes that shape bundle structure and function.

Trk activation of the ERK1/2 kinase pathway stimulates intermediate chain phosphorylation and recruits cytoplasmic dynein to signaling endosomes for retrograde axonal transport.

The retrograde transport of Trk-containing endosomes from the axon to the cell body by cytoplasmic dynein is necessary for axonal and neuronal survival. We investigated the recruitment of dynein to signaling endosomes in rat embryonic neurons and PC12 cells. We identified a novel phosphoserine on the dynein intermediate chains (ICs), and we observed a time-dependent neurotrophin-stimulated increase in intermediate chain phosphorylation on this site in both cell types. Pharmacological studies, overexpression of constitutively active MAP kinase kinase, and an in vitro assay with recombinant proteins demonstrated that the intermediate chains are phosphorylated by the MAP kinase ERK1/2, extracellular signal-regulated kinase, a major downstream effector of Trk. Live cell imaging with fluorescently tagged IC mutants demonstrated that the dephosphomimic mutants had significantly reduced colocalization with Trk and Rab7, but not a mitochondrial marker. The phosphorylated intermediate chains were enriched on immunoaffinity-purified Trk-containing organelles. Inhibition of ERK reduced the amount of phospho-IC and the total amount of dynein that copurified with the signaling endosomes. In addition, inhibition of ERK1/2 reduced the motility of Rab7- and TrkB-containing endosomes and the extent of their colocalization with dynein in axons. NGF-dependent survival of sympathetic neurons was significantly reduced by the overexpression of the dephosphomimic mutant IC-1B-S80A, but not WT IC-1B, further demonstrating the functional significance of phosphorylation on this site. These results demonstrate that neurotrophin binding to Trk initiates the recruitment of cytoplasmic dynein to signaling endosomes through ERK1/2 phosphorylation of intermediate chains for their subsequent retrograde transport in axons.

Matrix rigidity regulates cancer cell growth by modulating cellular metabolism and protein synthesis.

BACKGROUND: Tumor cells in vivo encounter diverse types of microenvironments both at the site of the primary tumor and at sites of distant metastases. Understanding how the various mechanical properties of these microenvironments affect the biology of tumor cells during disease progression is critical in identifying molecular targets for cancer therapy. METHODOLOGY/PRINCIPAL FINDINGS: This study uses flexible polyacrylamide gels as substrates for cell growth in conjunction with a novel proteomic approach to identify the properties of rigidity-dependent cancer cell lines that contribute to their differential growth on soft and rigid substrates. Compared to cells growing on more rigid/stiff substrates (>10,000 Pa), cells on soft substrates (150-300 Pa) exhibited a longer cell cycle, due predominantly to an extension of the G1 phase of the cell cycle, and were metabolically less active, showing decreased levels of intracellular ATP and a marked reduction in protein synthesis. Using stable isotope labeling of amino acids in culture (SILAC) and mass spectrometry, we measured the rates of protein synthesis of over 1200 cellular proteins under growth conditions on soft and rigid/stiff substrates. We identified cellular proteins whose syntheses were either preferentially inhibited or preserved on soft matrices. The former category included proteins that regulate cytoskeletal structures (e.g., tubulins) and glycolysis (e.g., phosphofructokinase-1), whereas the latter category included proteins that regulate key metabolic pathways required for survival, e.g., nicotinamide phosphoribosyltransferase, a regulator of the NAD salvage pathway. CONCLUSIONS/SIGNIFICANCE: The cellular properties of rigidity-dependent cancer cells growing on soft matrices are reminiscent of the properties of dormant cancer cells, e.g., slow growth rate and reduced metabolism. We suggest that the use of relatively soft gels as cell culture substrates would allow molecular pathways to be studied under conditions that reflect the different mechanical environments encountered by cancer cells upon metastasis to distant sites.

An experimentally derived database of candidate Ras-interacting proteins.

We used a TAP-tag approach to identify candidate binding proteins for the related Ras family GTPases: H-Ras, R-Ras, and Rap1A. Protein complexes were isolated from mouse fibroblasts, and component proteins were identified by a combination of nanoflow HPLC and tandem mass spectrometry. H-Ras was found to associate with numerous cytoskeletal proteins including talin-1. R-Ras and Rap1A each associated with various signaling molecules, many of which are membrane-associated. Thus, we have established the first database of potential Ras interactors in mammalian cells.

RLIP76 (RalBP1) is an R-Ras effector that mediates adhesion-dependent Rac activation and cell migration.

The Ras family of small GTPases regulates cell proliferation, spreading, migration and apoptosis, and malignant transformation by binding to several protein effectors. One such GTPase, R-Ras, plays distinct roles in each of these processes, but to date, identified R-Ras effectors were shared with other Ras family members (e.g., H-Ras). We utilized a new database of Ras-interacting proteins to identify RLIP76 (RalBP1) as a novel R-Ras effector. RLIP76 binds directly to R-Ras in a GTP-dependent manner, but does not physically associate with the closely related paralogues H-Ras and Rap1A. RLIP76 is required for adhesion-induced Rac activation and the resulting cell spreading and migration, as well as for the ability of R-Ras to enhance these functions. RLIP76 regulates Rac activity through the adhesion-induced activation of Arf6 GTPase and activation of Arf6 bypasses the requirement for RLIP76 in Rac activation and cell spreading. Thus, we identify a novel R-Ras effector, RLIP76, which links R-Ras to adhesion-induced Rac activation through a GTPase cascade that mediates cell spreading and migration.

The characterization of amphibian nucleoplasmins yields new insight into their role in sperm chromatin remodeling.

BACKGROUND: Nucleoplasmin is a nuclear chaperone protein that has been shown to participate in the remodeling of sperm chromatin immediately after fertilization by displacing highly specialized sperm nuclear basic proteins (SNBPs), such as protamine (P type) and protamine-like (PL type) proteins, from the sperm chromatin and by the transfer of histone H2A-H2B. The presence of SNBPs of the histone type (H type) in some organisms (very similar to the histones found in somatic tissues) raises uncertainty about the need for a nucleoplasmin-mediated removal process in such cases and poses a very interesting question regarding the appearance and further differentiation of the sperm chromatin remodeling function of nucleoplasmin and the implicit relationship with SNBP diversity The amphibians represent an unique opportunity to address this issue as they contain genera with SNBPs representative of each of the three main types: Rana (H type); Xenopus (PL type) and Bufo (P type). RESULTS: In this work, the presence of nucleoplasmin in oocyte extracts from these three organisms has been assessed using Western Blotting. We have used mass spectrometry and cloning techniques to characterize the full-length cDNA sequences of Rana catesbeiana and Bufo marinus nucleoplasmin. Northern dot blot analysis shows that nucleoplasmin is mainly transcribed in the egg of the former species. Phylogenetic analysis of nucleoplasmin family members from various metazoans suggests that amphibian nucleoplasmins group closely with mammalian NPM2 proteins. CONCLUSION: We have shown that these organisms, in striking contrast to their SNBPs, all contain nucleoplasmins with very similar primary structures. This result has important implications as it suggests that nucleoplasmin's role in chromatin assembly during early zygote development could have been complemented by the acquisition of a new function of non-specifically removing SNBPs in sperm chromatin remodeling. This acquired function would have been strongly determined by the constraints imposed by the appearance and differentiation of SNBPs in the sperm.