Workflow enabling deepscale immunopeptidome, proteome, ubiquitylome, phosphoproteome, and acetylome analyses of sample-limited tissues.

Serial multi-omic analysis of proteome, phosphoproteome, and acetylome provides insights into changes in protein expression, cell signaling, cross-talk and epigenetic pathways involved in disease pathology and treatment. However, ubiquitylome and HLA peptidome data collection used to understand protein degradation and antigen presentation have not together been serialized, and instead require separate samples for parallel processing using distinct protocols. Here we present MONTE, a highly sensitive multi-omic native tissue enrichment workflow, that enables serial, deep-scale analysis of HLA-I and HLA-II immunopeptidome, ubiquitylome, proteome, phosphoproteome, and acetylome from the same tissue sample. We demonstrate that the depth of coverage and quantitative precision of each 'ome is not compromised by serialization, and the addition of HLA immunopeptidomics enables the identification of peptides derived from cancer/testis antigens and patient specific neoantigens. We evaluate the technical feasibility of the MONTE workflow using a small cohort of patient lung adenocarcinoma tumors.

Automating UbiFast for High-throughput and Multiplexed Ubiquitin Enrichment.

Robust methods for deep-scale enrichment and site-specific identification of ubiquitylation sites are necessary for characterizing the myriad roles of protein ubiquitylation. To this end we previously developed UbiFast, a sensitive method for highly multiplexed ubiquitylation profiling where K-ϵ-GG peptides are enriched with anti-K-ε-GG antibody and labeled on-antibody with isobaric labeling reagents for sample multiplexing. Here, we present robotic automation of the UbiFast method using a magnetic bead-conjugated K-ε-GG antibody (mK-ε-GG) and a magnetic particle processor. We report the identification of ∼20,000 ubiquitylation sites from a TMT10-plex with 500 µg input per sample processed in ∼2 h. Automation of the UbiFast method greatly increased the number of identified and quantified ubiquitylation sites, improved reproducibility, and significantly reduced processing time. The automated method also significantly reduced variability across process replicates compared with the manual method. The workflow enables processing of up to 96 samples in a single day making it suitable to study ubiquitylation in large sample sets. Here we demonstrate the applicability of the method to profile small amounts of tissue using breast cancer patient-derived xenograft (PDX) tissue samples.

Profiling SARS-CoV-2 HLA-I peptidome reveals T cell epitopes from out-of-frame ORFs.

T cell-mediated immunity plays an important role in controlling SARS-CoV-2 infection, but the repertoire of naturally processed and presented viral epitopes on class I human leukocyte antigen (HLA-I) remains uncharacterized. Here, we report the first HLA-I immunopeptidome of SARS-CoV-2 in two cell lines at different times post infection using mass spectrometry. We found HLA-I peptides derived not only from canonical open reading frames (ORFs) but also from internal out-of-frame ORFs in spike and nucleocapsid not captured by current vaccines. Some peptides from out-of-frame ORFs elicited T cell responses in a humanized mouse model and individuals with COVID-19 that exceeded responses to canonical peptides, including some of the strongest epitopes reported to date. Whole-proteome analysis of infected cells revealed that early expressed viral proteins contribute more to HLA-I presentation and immunogenicity. These biological insights, as well as the discovery of out-of-frame ORF epitopes, will facilitate selection of peptides for immune monitoring and vaccine development.

Identification of protein interactions of grapevine fanleaf virus RNA-dependent RNA polymerase during infection of Nicotiana benthamiana by affinity purification and tandem mass spectrometry.

The RNA-dependent RNA polymerase (1EPol) is involved in replication of grapevine fanleaf virus (GFLV, Nepovirus, Secoviridae) and causes vein clearing symptoms in Nicotiana benthamiana. Information on protein 1EPol interaction with other viral and host proteins is scarce. To study protein 1EPol biology, three GFLV infectious clones, i.e. GHu (a symptomatic wild-type strain), GHu-1EK802G (an asymptomatic GHu mutant) and F13 (an asymptomatic wild-type strain), were engineered with protein 1EPol fused to a V5 epitope tag at the C-terminus. Following Agrobacterium tumefaciens-mediated delivery of GFLV clones in N. benthamiana and protein extraction at seven dpi, when optimal 1EPol:V5 accumulation was detected, two viral and six plant putative interaction partners of V5-tagged protein 1EPol were identified for the three GFLV clones by affinity purification and tandem mass spectrometry. This study provides insights into the protein interactome of 1EPol during GFLV systemic infection in N. benthamiana and lays the foundation for validation work.

Oncogenic KRAS engages an RSK1/NF1 pathway to inhibit wild-type RAS signaling in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy with limited treatment options. Although activating mutations of the KRAS GTPase are the predominant dependency present in >90% of PDAC patients, targeting KRAS mutants directly has been challenging in PDAC. Similarly, strategies targeting known KRAS downstream effectors have had limited clinical success due to feedback mechanisms, alternate pathways, and dose-limiting toxicities in normal tissues. Therefore, identifying additional functionally relevant KRAS interactions in PDAC may allow for a better understanding of feedback mechanisms and unveil potential therapeutic targets. Here, we used proximity labeling to identify protein interactors of active KRAS in PDAC cells. We expressed fusions of wild-type (WT) (BirA-KRAS4B), mutant (BirA-KRAS4BG12D), and nontransforming cytosolic double mutant (BirA-KRAS4BG12D/C185S) KRAS with the BirA biotin ligase in murine PDAC cells. Mass spectrometry analysis revealed that RSK1 selectively interacts with membrane-bound KRASG12D, and we demonstrate that this interaction requires NF1 and SPRED2. We find that membrane RSK1 mediates negative feedback on WT RAS signaling and impedes the proliferation of pancreatic cancer cells upon the ablation of mutant KRAS. Our findings link NF1 to the membrane-localized functions of RSK1 and highlight a role for WT RAS signaling in promoting adaptive resistance to mutant KRAS-specific inhibitors in PDAC.

SARS-CoV-2 infected cells present HLA-I peptides from canonical and out-of-frame ORFs.

T cell-mediated immunity may play a critical role in controlling and establishing protective immunity against SARS-CoV-2 infection; yet the repertoire of viral epitopes responsible for T cell response activation remains mostly unknown. Identification of viral peptides presented on class I human leukocyte antigen (HLA-I) can reveal epitopes for recognition by cytotoxic T cells and potential incorporation into vaccines. Here, we report the first HLA-I immunopeptidome of SARS-CoV-2 in two human cell lines at different times post-infection using mass spectrometry. We found HLA-I peptides derived not only from canonical ORFs, but also from internal out-of-frame ORFs in Spike and Nucleoprotein not captured by current vaccines. Proteomics analyses of infected cells revealed that SARS-CoV-2 may interfere with antigen processing and immune signaling pathways. Based on the endogenously processed and presented viral peptides that we identified, we estimate that a pool of 24 peptides would provide one or more peptides for presentation by at least one HLA allele in 99% of the human population. These biological insights and the list of naturally presented SARS-CoV-2 peptides will facilitate data-driven selection of peptides for immune monitoring and vaccine development.

Mito-oncology agent: fermented extract suppresses the Warburg effect, restores oxidative mitochondrial activity, and inhibits in vivo tumor growth.

Mitochondrial dysfunction and significant changes in metabolic pathways accompany cancer development and are responsible for maintaining the tumor microenvironment. Normal mitochondria can trigger intrinsic apoptosis by releasing cytochrome c into the cytosol. The survival of malignant cells highly depends on the suppression of this function. We validated that A250, a highly purified fraction of fermented wheat germ extract (FWGE), increases the carbon flux into the mitochondria, the expression of key elements of the Krebs cycle and oxidative phosphorylation (OXPHOS). The increased respiratory chain activity is related to the mitochondria's ability to release cytochrome c into the cytosol, which triggers the apoptotic cascade. The 68% tumor growth inhibitory effect observed in the murine melanoma study is related to this effect, as proteomic analysis validated similar changes in mitochondrial protein levels in the isolated tumor tissue samples. Blood count data indicated that this effect was not accompanied by general toxicity. This study is significant, as it shows that a highly concentrated form of FWGE is an effective agent that increases normal mitochondrial functionality. The lack of hepatotoxic and general toxic effects makes A250 an excellent candidate targeting mitochondria function in cancer therapy.

Tryp-N: A Thermostable Protease for the Production of N-terminal Argininyl and Lysinyl Peptides.

Bottom-up proteomics is a mainstay in protein identification and analysis. These studies typically employ proteolytic treatment of biological samples to generate suitably sized peptides for tandem mass spectrometric (MS) analysis. In MS, fragmentation of peptides is largely driven by charge localization. Consequently, peptides with basic centers exclusively on their N-termini produce mainly b-ions. Thus, it was long ago realized that proteases that yield such peptides would be valuable proteomic tools for achieving simplified peptide fragmentation patterns and peptide assignment. Work by several groups has identified such proteases, however, structural analysis of these suggested that enzymatic optimization was possible. We therefore endeavored to find enzymes that could provide enhanced activity and versatility while maintaining specificity. Using these previously described proteases as informatic search templates, we discovered and then characterized a thermophilic metalloprotease with N-terminal specificity for arginine and lysine. This enzyme, dubbed Tryp-N, affords many advantages including improved thermostability, solvent and detergent tolerance, and rapid digestion time.

Regulation of PTP1B activation through disruption of redox-complex formation.

We have identified a molecular interaction between the reversibly oxidized form of protein tyrosine phosphatase 1B (PTP1B) and 14-3-3ζ that regulates PTP1B activity. Destabilizing the transient interaction between 14-3-3ζ and PTP1B prevented PTP1B inactivation by reactive oxygen species and decreased epidermal growth factor receptor phosphorylation. Our data suggest that destabilizing the interaction between 14-3-3ζ and the reversibly oxidized and inactive form of PTP1B may establish a path to PTP1B activation in cells.

The glycan CA19-9 promotes pancreatitis and pancreatic cancer in mice.

Glycosylation alterations are indicative of tissue inflammation and neoplasia, but whether these alterations contribute to disease pathogenesis is largely unknown. To study the role of glycan changes in pancreatic disease, we inducibly expressed human fucosyltransferase 3 and ß1,3-galactosyltransferase 5 in mice, reconstituting the glycan sialyl-Lewisa, also known as carbohydrate antigen 19-9 (CA19-9). Notably, CA19-9 expression in mice resulted in rapid and severe pancreatitis with hyperactivation of epidermal growth factor receptor (EGFR) signaling. Mechanistically, CA19-9 modification of the matricellular protein fibulin-3 increased its interaction with EGFR, and blockade of fibulin-3, EGFR ligands, or CA19-9 prevented EGFR hyperactivation in organoids. CA19-9-mediated pancreatitis was reversible and could be suppressed with CA19-9 antibodies. CA19-9 also cooperated with the KrasG12D oncogene to produce aggressive pancreatic cancer. These findings implicate CA19-9 in the etiology of pancreatitis and pancreatic cancer and nominate CA19-9 as a therapeutic target.

Muller's Ratchet and Ribosome Degeneration in the Obligate Intracellular Parasites Microsporidia.

Microsporidia are fungi-like parasites that have the smallest known eukaryotic genome, and for that reason they are used as a model to study the phenomenon of genome decay in parasitic forms of life. Similar to other intracellular parasites that reproduce asexually in an environment with alleviated natural selection, Microsporidia experience continuous genome decay that is driven by Muller's ratchet-an evolutionary process of irreversible accumulation of deleterious mutations that lead to gene loss and the miniaturization of cellular components. Particularly, Microsporidia have remarkably small ribosomes in which the rRNA is reduced to the minimal enzymatic core. In this study, we analyzed microsporidian ribosomes to study an apparent impact of Muller's ratchet on structure of RNA and protein molecules in parasitic forms of life. Through mass spectrometry of microsporidian proteome and analysis of microsporidian genomes, we found that massive rRNA reduction in microsporidian ribosomes appears to annihilate the binding sites for ribosomal proteins eL8, eL27, and eS31, suggesting that these proteins are no longer bound to the ribosome in microsporidian species. We then provided an evidence that protein eS31 is retained in Microsporidia due to its non-ribosomal function in ubiquitin biogenesis. Our study illustrates that, while Microsporidia carry the same set of ribosomal proteins as non-parasitic eukaryotes, some ribosomal proteins are no longer participating in protein synthesis in Microsporidia and they are preserved from genome decay by having extra-ribosomal functions. More generally, our study shows that many components of parasitic cells, which are identified by automated annotation of pathogenic genomes, may lack part of their biological functions due to continuous genome decay.

Error-prone protein synthesis in parasites with the smallest eukaryotic genome.

Microsporidia are parasitic fungi-like organisms that invade the interior of living cells and cause chronic disorders in a broad range of animals, including humans. These pathogens have the tiniest known genomes among eukaryotic species, for which they serve as a model for exploring the phenomenon of genome reduction in obligate intracellular parasites. Here we report a case study to show an apparent effect of overall genome reduction on the primary structure and activity of aminoacyl-tRNA synthetases, indispensable cellular proteins required for protein synthesis. We find that most microsporidian synthetases lack regulatory and eukaryote-specific appended domains and have a high degree of sequence variability in tRNA-binding and catalytic domains. In one synthetase, LeuRS, an apparent sequence degeneration annihilates the editing domain, a catalytic center responsible for the accurate selection of leucine for protein synthesis. Unlike accurate LeuRS synthetases from other eukaryotic species, microsporidian LeuRS is error-prone: apart from leucine, it occasionally uses its near-cognate substrates, such as norvaline, isoleucine, valine, and methionine. Mass spectrometry analysis of the microsporidium Vavraia culicis proteome reveals that nearly 6% of leucine residues are erroneously replaced by other amino acids. This remarkably high frequency of mistranslation is not limited to leucine codons and appears to be a general property of protein synthesis in microsporidian parasites. Taken together, our findings reveal that the microsporidian protein synthesis machinery is editing-deficient, and that the proteome of microsporidian parasites is more diverse than would be anticipated based on their genome sequences.

Organoid models of human and mouse ductal pancreatic cancer.

Pancreatic cancer is one of the most lethal malignancies due to its late diagnosis and limited response to treatment. Tractable methods to identify and interrogate pathways involved in pancreatic tumorigenesis are urgently needed. We established organoid models from normal and neoplastic murine and human pancreas tissues. Pancreatic organoids can be rapidly generated from resected tumors and biopsies, survive cryopreservation, and exhibit ductal- and disease-stage-specific characteristics. Orthotopically transplanted neoplastic organoids recapitulate the full spectrum of tumor development by forming early-grade neoplasms that progress to locally invasive and metastatic carcinomas. Due to their ability to be genetically manipulated, organoids are a platform to probe genetic cooperation. Comprehensive transcriptional and proteomic analyses of murine pancreatic organoids revealed genes and pathways altered during disease progression. The confirmation of many of these protein changes in human tissues demonstrates that organoids are a facile model system to discover characteristics of this deadly malignancy.

Dephosphorylation of tyrosine 393 in argonaute 2 by protein tyrosine phosphatase 1B regulates gene silencing in oncogenic RAS-induced senescence.

Oncogenic RAS (H-RAS(V12)) induces premature senescence in primary cells by triggering production of reactive oxygen species (ROS), but the molecular role of ROS in senescence remains elusive. We investigated whether inhibition of protein tyrosine phosphatases by ROS contributed to H-RAS(V12)-induced senescence. We identified protein tyrosine phosphatase 1B (PTP1B) as a major target of H-RAS(V12)-induced ROS. Inactivation of PTP1B was necessary and sufficient to induce premature senescence in H-RAS(V12)-expressing IMR90 fibroblasts. We identified phospho-Tyr 393 of argonaute 2 (AGO2) as a direct substrate of PTP1B. Phosphorylation of AGO2 at Tyr 393 inhibited loading with microRNAs (miRNAs) and thus miRNA-mediated gene silencing, which counteracted the function of H-RAS(V12)-induced oncogenic miRNAs. Overall, our data illustrate that premature senescence in H-RAS(V12)-transformed primary cells is a consequence of oxidative inactivation of PTP1B and inhibition of miRNA-mediated gene silencing.

Brainstem deficiency of the 14-3-3 regulator of serotonin synthesis: a proteomics analysis in the sudden infant death syndrome.

Impaired brainstem responses to homeostatic challenges during sleep may result in the sudden infant death syndrome (SIDS). Previously we reported a deficiency of serotonin (5-HT) and its key biosynthetic enzyme, tryptophan hydroxylase (TPH2), in SIDS infants in the medullary 5-HT system that modulates homeostatic responses during sleep. Yet, the underlying basis of the TPH2 and 5-HT deficiency is unknown. In this study, we tested the hypothesis that proteomics would uncover previously unrecognized abnormal levels of proteins related to TPH2 and 5-HT regulation in SIDS cases compared with controls, which could provide novel insight into the basis of their deficiency. We first performed a discovery proteomic analysis of the gigantocellularis of the medullary 5-HT system in the same data set with deficiencies of TPH2 and 5-HT levels. Analysis in 6 SIDS cases and 4 controls revealed a 42-75% reduction in abundance in 5 of the 6 isoforms identified of the 14-3-3 signal transduction family, which is known to influence TPH2 activity (p < 0.07). These findings were corroborated in an additional SIDS and control sample using an orthogonal MS(E)-based quantitative proteomic strategy. To confirm these proteomics results in a larger data set (38 SIDS, 11 controls), we applied Western blot analysis in the gigantocellularis and found that 4/7 14-3-3 isoforms identified were significantly reduced in SIDS cases (p ≤ 0.02), with a 43% reduction in all 14-3-3 isoforms combined (p < 0.001). Abnormalities in 5-HT and TPH2 levels and 5-HT(1A) receptor binding were associated with the 14-3-3 deficits in the same SIDS cases. These data suggest a potential molecular defect in SIDS related to TPH2 regulation, as 14-3-3 is critical in this process.

Interlaboratory evaluation of automated, multiplexed peptide immunoaffinity enrichment coupled to multiple reaction monitoring mass spectrometry for quantifying proteins in plasma.

The inability to quantify large numbers of proteins in tissues and biofluids with high precision, sensitivity, and throughput is a major bottleneck in biomarker studies. We previously demonstrated that coupling immunoaffinity enrichment using anti-peptide antibodies (SISCAPA) to multiple reaction monitoring mass spectrometry (MRM-MS) produces Immunoprecipitation MRM-MS (immuno-MRM-MS) assays that can be multiplexed to quantify proteins in plasma with high sensitivity, specificity, and precision. Here we report the first systematic evaluation of the interlaboratory performance of multiplexed (8-plex) immuno-MRM-MS in three independent labs. A staged study was carried out in which the effect of each processing and analysis step on assay coefficient of variance, limit of detection, limit of quantification, and recovery was evaluated. Limits of detection were at or below 1 ng/ml for the assayed proteins in 30 µl of plasma. Assay reproducibility was acceptable for verification studies, with median intra- and interlaboratory coefficients of variance above the limit of quantification of 11% and <14%, respectively, for the entire immuno-MRM-MS assay process, including enzymatic digestion of plasma. Trypsin digestion and its requisite sample handling contributed the most to assay variability and reduced the recovery of target peptides from digested proteins. Using a stable isotope-labeled protein as an internal standard instead of stable isotope-labeled peptides to account for losses in the digestion process nearly doubled assay accuracy for this while improving assay precision 5%. Our results demonstrate that multiplexed immuno-MRM-MS can be made reproducible across independent laboratories and has the potential to be adopted widely for assaying proteins in matrices as complex as plasma.

Lack of association of the serotonin transporter polymorphism with the sudden infant death syndrome in the San Diego Dataset.

Dysfunction of medullary serotonin (5-HT)-mediated respiratory and autonomic function is postulated to underlie the pathogenesis of the majority of sudden infant death syndrome (SIDS) cases. Several studies have reported an increased frequency of the LL genotype and L allele of the 5-HT transporter (5-HTT) gene promoter polymorphism (5-HTTLPR), which is associated with increased transcriptional activity and 5-HT transport in vitro, in SIDS cases compared with controls. These findings raise the possibility that this polymorphism contributes to or exacerbates existing medullary 5-HT dysfunction in SIDS. In this study, we tested the hypothesis that the frequency of LL genotype and L allele are higher in 179 SIDS cases compared with 139 controls of multiple ethnicities in the San Diego SIDS Dataset. We observed no significant association of genotype or allele with SIDS cases either in the total cohort or on stratification for ethnicity. These observations do not support previous findings that the L allele and/or LL genotype of the 5-HTTLPR are associated with SIDS.

Neuroanatomic relationships between the GABAergic and serotonergic systems in the developing human medulla.

gamma-Amino butyric (GABA) critically influences serotonergic (5-HT) neurons in the raphé and extra-raphé of the medulla oblongata. In this study we hypothesize that there are marked changes in the developmental profile of markers of the human medullary GABAergic system relative to the 5-HT system in early life. We used single- and double-label immunocytochemistry and tissue receptor autoradiography in 15 human medullae from fetal and infant cases ranging from 15 gestational weeks to 10 postnatal months, and compared our findings with an extensive 5-HT-related database in our laboratory. In the raphé obscurus, we identified two subsets of GABAergic neurons using glutamic acid decarboxylase (GAD65/67) immunostaining: one comprised of small, round neurons; the other, medium, spindle-shaped neurons. In three term medullae cases, positive immunofluorescent neurons for both tryptophan hydroxylase and GAD65/67 were counted within the raphé obscurus. This revealed that approximately 6% of the total neurons counted in this nucleus expressed both GAD65/67 and TPOH suggesting co-production of GABA by a subset of 5-HT neurons. The distribution of GABA(A) binding was ubiquitous across medullary nuclei, with highest binding in the raphé obscurus. GABA(A) receptor subtypes alpha1 and alpha3 were expressed by 5-HT neurons, indicating the site of interaction of GABA with 5-HT neurons. These receptor subtypes and KCC2, a major chloride transporter, were differentially expressed across early development, from midgestation (20 weeks) and thereafter. The developmental profile of GABAergic markers changed dramatically relative to the 5-HT markers. These data provide baseline information for medullary studies of human pediatric disorders, such as sudden infant death syndrome.