Small admixtures in water, e.g. of metal ions, often act as cell growth regulators. Here we report that enrichment of deuterium content in water, normally found at 8 mm concentration, two-three folds increases cell proliferation and lowers the oxidative stress level as well. Acting as an anti-oxidant, deuterium-enriched water prevents the toxic effect of such oxidative agents as hydrogen peroxide and auranofin. This action is opposite to that of deuterium depletion that is known to suppress cell growth and induce oxidative stress in mitochondria. We thus hypothesize that deuterium may be a natural cell growth regulator that controls mitochondrial oxidation-reduction balance. Because growth acceleration is reduced approximately by half by addition to water a minute amount (0.15%) of ¹⁸O isotope, at least part of the deuterium effect on cell growth can be explained by the isotopic resonance phenomenon. A slight (2-fold) enrichment of deuterium in water accelerates human cell growth. Quantitative MS based proteomics determined changes in protein abundances and redox states and found that deuterium-enriched water acts mainly through decreasing ROS production in mitochondria. This action is opposite to that of deuterium depletion that suppresses cell growth by inducing oxidative stress. Thus deuterium may be a natural cell growth regulator that controls mitochondrial oxidation-reduction balance. The role of isotopic resonance in this effect was validated by further experiments on bacteria.

Amino acid hydroxylation is a common post-translational modification, which generally regulates protein interactions or adds a functional group that can be further modified. Such hydroxylation is currently considered irreversible, necessitating the degradation and re-synthesis of the entire protein to reset the modification. Here we present evidence that the cellular machinery can reverse FIH-mediated asparagine hydroxylation on intact proteins. These data suggest that asparagine hydroxylation is a flexible and dynamic post-translational modification akin to modifications involved in regulating signaling networks, such as phosphorylation, methylation and ubiquitylation.

We performed an in-depth characterization and comparison of the pediatric and adult urinary glycomes using a nanoLC-MS/MS based glycomics method, which included normal healthy pediatric (1–10 years, n = 21) and adult (21–50 years, n = 22) individuals. A total of 116 N-glycan compositions were identified, and 46 of them could be reproducibly quantified. We performed quantitative comparisons of the 46 glycan compositions between different age and sex groups. The results showed significant quantitative changes between the pediatric and adult cohorts. The pediatric urinary N-glycome was found to contain a higher level of high-mannose (HM), asialylated/afucosylated glycans (excluding HM), neutral fucosylated and agalactosylated glycans, and a lower level of trisialylated glycans compared with the adult. We further analyzed gender-associated glycan changes in the pediatric and adult group, respectively. In the pediatric group, there was almost no difference of glycan levels between males and females. In adult, the majority of glycans were more abundant in males than females, except the high-mannose and tetrasialylated glycans. These findings highlight the importance to consider age-matching and adult sex-matching for urinary glycan studies. The identified normal pediatric and adult urinary glycomes can serve as a baseline reference for comparisons to other disease states affected by glycosylation.

Coronavirus disease 2019 (COVID-19) is a highly contagious infection and threating the human lives in the world. The elevation of cytokines in blood is crucial to induce cytokine storm and immunosuppression in the transition of severity in COVID-19 patients. However, the comprehensive changes of serum proteins in COVID-19 patients throughout the SARS-CoV-2 infection is unknown. In this work, we developed a high-density antibody microarray and performed an in-depth proteomics analysis of serum samples collected from early COVID-19 (n = 15) and influenza (n = 13) patients. We identified a large set of differentially expressed proteins (n = 132) that participate in a landscape of inflammation and immune signaling related to the SARS-CoV-2 infection. Furthermore, the significant correlations of neutrophil and lymphocyte with the CCL2 and CXCL10 mediated cytokine signaling pathways was identified. These information are valuable for the understanding of COVID-19 pathogenesis, identification of biomarkers and development of the optimal anti-inflammation therapy.

Pathway analyses are key methods to analyze 'omics experiments. Nevertheless, integrating data from different 'omics technologies and different species still requires considerable bioinformatics knowledge.

Here we present the novel ReactomeGSA resource for comparative pathway analyses of multi-omics datasets. ReactomeGSA can be used through Reactome's existing web interface and the novel ReactomeGSA R Bioconductor package with explicit support for scRNA-seq data. Data from different species is automatically mapped to a common pathway space. Public data from ExpressionAtlas and Single Cell ExpressionAtlas can be directly integrated in the analysis. ReactomeGSA greatly reduces the technical barrier for multi-omics, cross-species, comparative pathway analyses.

We used ReactomeGSA to characterize the role of B cells in anti-tumor immunity. We compared B cell rich and poor human cancer samples from five of the Cancer Genome Atlas (TCGA) transcriptomics and two of the Clinical Proteomic Tumor Analysis Consortium (CPTAC) proteomics studies. B cell-rich lung adenocarcinoma samples lacked the otherwise present activation through NFkappaB. This may be linked to the presence of a specific subset of tumor associated IgG+ plasma cells that lack NFkappaB activation in scRNA-seq data from human melanoma. This showcases how ReactomeGSA can derive novel biomedical insights by integrating large multi-omics datasets.

Despite the crucial function of the small intestine in nutrient uptake our understanding of the molecular events underlying the digestive function is still rudimentary. Recent studies demonstrated that enterocytes do not direct the entire dietary triacylglycerol toward immediate chylomicron synthesis. Especially after high-fat challenges, parts of the resynthesized triacylglycerol are packaged into cytosolic lipid droplets for transient storage in the endothelial layer of the small intestine. The reason for this temporary storage of triacylglycerol is not completely understood. To utilize lipids from cytosolic lipid droplets for chylomicron synthesis in the endoplasmic reticulum, stored triacylglycerol has to be hydrolyzed either by cytosolic lipolysis or lipophagy. Interestingly, triacylglycerol storage and chylomicron secretion rates are unevenly distributed along the small intestine, with the proximal jejunum exhibiting the highest intermittent storage capacity. We hypothesize that correlating hydrolytic enzyme activities with the reported distribution of triacylglycerol storage and chylomicron secretion in different sections of the small intestine is a promising strategy to determine key enzymes in triacylglycerol remobilization. We employed a serine hydrolase specific activity-based labeling approach in combination with quantitative proteomics to identify and rank hydrolases based on their relative activity in 11 sections of the small intestine. Moreover, we identified several clusters of enzymes showing similar activity distribution along the small intestine. Merging our activity-based results with substrate specificity and subcellular localization known from previous studies, carboxylesterase 2e and arylacetamide deacetylase emerge as promising candidates for triacylglycerol mobilization from cytosolic lipid droplets in enterocytes.

Extracellular vesicles (EVs) secreted by the epididymal epithelium transfer to spermatozoa key proteins that are essential in promoting motility and subsequent fertilization success. Using the domestic cat model, the objectives were to (1) characterize and compare protein content of EVs between segments of the epididymis, and (2) compare EV protein compositions between normo- and teratospermic individuals (producing >60% of abnormal spermatozoa). Epididymal EVs from adult cats were isolated and assessed via liquid chromatography tandem MS. Both male types shared 3008 proteins in total, with 98 and 20 EV proteins unique to normospermic and teratospermic males, respectively. Expression levels of several proteins changed between epididymal segments in both male types. Several proteins in both groups were related to sperm motility (e.g. hexokinase 1, adenylate kinase isoenzyme) and zona pellucida or oolemma binding (e.g. disintegrin and metalloproteinase domain proteins, zona binding proteins 1 and 2). Interestingly, seven cauda-derived EV proteins trended downward in teratospermic compared with normospermic males, which may relate to poor sperm quality. Collective results revealed, for the first time, EV proteins related to sequential sperm maturation with differences observed between normospermic and teratospermic individuals.

Endometrial carcinoma (EC) is the most common gynecologic malignancy in the United States, with limited effective targeted therapies. Endometrial tumors exhibit frequent alterations in protein kinases, yet only a small fraction of the kinome has been therapeutically explored. To identify kinase therapeutic avenues for EC, we profiled the kinome of endometrial tumors and normal endometrial tissues using Multiplexed Inhibitor Beads and Mass Spectrometry (MIB-MS). Our proteomics analysis identified a network of kinases overexpressed in tumors, including Serine/Arginine-Rich Splicing Factor Kinase 1 (SRPK1). Immunohistochemical (IHC) analysis of endometrial tumors confirmed MIB-MS findings and showed SRPK1 protein levels were highly expressed in endometrioid and uterine serous cancer (USC) histological subtypes. Moreover, querying large-scale genomics studies of EC tumors revealed high expression of SRPK1 correlated with poor survival. Loss-of-function studies targeting SRPK1 in an established USC cell line demonstrated SRPK1 was integral for RNA splicing, as well as cell cycle progression and survival under nutrient deficient conditions. Profiling of USC cells identified a compensatory response to SRPK1 inhibition that involved EGFR and the up-regulation of IGF1R and downstream AKT signaling. Co-targeting SRPK1 and EGFR or IGF1R synergistically enhanced growth inhibition in serous and endometrioid cell lines, representing a promising combination therapy for EC.

The absence of the dystrophin protein in Duchenne muscular dystrophy (DMD) results in myofiber fragility and a plethora of downstream secondary pathologies. Although a variety of experimental therapies are in development, achieving effective treatments for DMD remains exceptionally challenging, not least because the pathological consequences of dystrophin loss are incompletely understood. Here we have performed proteome profiling in tibialis anterior muscles from two murine DMD models (mdx and mdx52) at three ages (8, 16, and 80 weeks of age), all n = 3. High-resolution isoelectric focusing liquid chromatography-tandem MS (HiRIEF-LC–MS/MS) was used to quantify the expression of 4974 proteins across all 27 samples. The two dystrophic models were found to be highly similar, whereas multiple proteins were differentially expressed relative to WT (C57BL/6) controls at each age. Furthermore, 1795 proteins were differentially expressed when samples were pooled across ages and dystrophic strains. These included numerous proteins associated with the extracellular matrix and muscle function that have not been reported previously. Pathway analysis revealed multiple perturbed pathways and predicted upstream regulators, which together are indicative of cross-talk between inflammatory, metabolic, and muscle growth pathways (e.g. TNF, INF, NF-B, SIRT1, AMPK, PGC-1α, PPARs, ILK, and AKT/PI3K). Upregulation of CAV3, MVP and PAK1 protein expression was validated in dystrophic muscle by Western blot. Furthermore, MVP was upregulated during, but not required for, the differentiation of C2C12 myoblasts suggesting that this protein may affect muscle regeneration. This study provides novel insights into mutation-independent proteomic signatures characteristic of the dystrophic phenotype and its progression with aging.

Sepsis-induced acute kidney injury (S-AKI) is the most common complication in hospitalized and critically ill patients, highlighted by a rapid decline of kidney function occurring a few hours or days after sepsis onset. Systemic inflammation elicited by microbial infections is believed to lead to kidney damage under immunocompromised conditions. However, although AKI has been recognized as a disease with long-term sequelae, partly because of the associated higher risk of chronic kidney disease (CKD), the understanding of kidney pathophysiology at the molecular level and the global view of dynamic regulations in situ after S-AKI, including the transition to CKD, remains limited. Existing studies of S-AKI mainly focus on deriving sepsis biomarkers from body fluids. In the present study, we constructed a mid-severity septic murine model using cecal ligation and puncture (CLP), and examined the temporal changes to the kidney proteome and phosphoproteome at day 2 and day 7 after CLP surgery, corresponding to S-AKI and the transition to CKD, respectively, by employing an ultrafast and economical filter-based sample processing method combined with the label-free quantitation approach. Collectively, we identified 2,119 proteins and 2950 phosphosites through multi-proteomics analyses. Among them, we identified an array of highly promising candidate marker proteins indicative of disease onset and progression accompanied by immunoblot validations, and further denoted the pathways that are specifically responsive to S-AKI and its transition to CKD, which include regulation of cell metabolism regulation, oxidative stress, and energy consumption in the diseased kidneys. Our data can serve as an enriched resource for the identification of mechanisms and biomarkers for sepsis-induced kidney diseases.

Specific E3 ligases target tumor suppressors for degradation. Inhibition of such E3 ligases may be an important approach to cancer treatment. RNF146 is a RING domain and PARylation-dependent E3 ligase that functions as an activator of the β-catenin/Wnt and YAP/Hippo pathways by targeting the degradation of several tumor suppressors. Tankyrases 1 and 2 (TNKS1/2) are the only known poly-ADP-ribosyltransferases that require RNF146 to degrade their substrates. However, systematic identification of RNF146 substrates have not yet been performed. To uncover substrates of RNF146 that are targeted for degradation, we generated RNF146 knockout cells and TNKS1/2-double knockout cells and performed proteome profiling with label-free quantification as well as transcriptome analysis. We identified 160 potential substrates of RNF146, which included many known substrates of RNF146 and TNKS1/2 and 122 potential TNKS-independent substrates of RNF146. In addition, we validated OTU domain-containing protein 5 and Protein mono-ADP-ribosyltransferase PARP10 as TNKS1/2-independent substrates of RNF146 and SARDH as a novel substrate of TNKS1/2 and RNF146. Our study is the first proteome-wide analysis of potential RNF146 substrates. Together, these findings not only demonstrate that proteome profiling can be a useful general approach for the systemic identification of substrates of E3 ligases but also reveal new substrates of RNF146, which provides a resource for further functional studies.

AAA+ ATPases constitute a large family of proteins that are involved in a plethora of cellular processes including DNA disassembly, protein degradation and protein complex disassembly. They typically form a hexametric ring-shaped structure with six subunits in a (pseudo) 6-fold symmetry. In a subset of AAA+ ATPases that facilitate protein unfolding and degradation, six subunits cooperate to translocate protein substrates through a central pore in the ring. The number and type of nucleotides in an AAA+ ATPase hexamer is inherently linked to the mechanism that underlies cooperation among subunits and couples ATP hydrolysis with substrate translocation. We conducted a native MS study of a monodispersed form of PAN, an archaeal proteasome AAA+ ATPase, to determine the number of nucleotides bound to each hexamer of the WT protein. We utilized ADP and its analogs (TNP-ADP and mant-ADP), and a nonhydrolyzable ATP analog (AMP-PNP) to study nucleotide site occupancy within the PAN hexamer in ADP- and ATP-binding states, respectively. Throughout all experiments we used a Walker A mutant (PAN^K217A) that is impaired in nucleotide binding as an internal standard to mitigate the effects of residual solvation on mass measurement accuracy and to serve as a reference protein to control for nonspecific nucleotide binding. This approach led to the unambiguous finding that a WT PAN hexamer carried – from expression host – six tightly bound ADP molecules that could be exchanged for ADP and ATP analogs. Although the Walker A mutant did not bind ADP analogs, it did bind AMP-PNP, albeit at multiple stoichiometries. We observed variable levels of hexamer dissociation and an appearance of multimeric species with the over-charged molecular ion distributions across repeated experiments. We posit that these phenomena originated during ESI process at the final stages of ESI droplet evolution.

Plasmodium, the malaria parasite, undergoes a complex life cycle alternating between a vertebrate host and a mosquito vector of the genus Anopheles. In red blood cells of the vertebrate host, Plasmodium multiplies asexually or differentiates into gamete precursors, the male and female gametocytes, responsible for parasite transmission. Sexual stage maturation occurs in the midgut of the mosquito vector, where male and female gametes egress from the host erythrocytes to fuse and form a zygote. Gamete egress entails the successive rupture of two membranes surrounding the parasite, the parasitophorous vacuole membrane and the erythrocyte plasma membrane. In this study, we used the rodent model parasite Plasmodium berghei to design a label-free quantitative proteomic approach aimed at identifying gender-related proteins differentially released/secreted by purified mature gametocytes when activated to form gametes. We compared the abundance of molecules secreted by wild type gametocytes of both genders with that of a transgenic line defective in male gamete maturation and egress. This enabled us to provide a comprehensive data set of egress-related molecules and their gender specificity. Using specific antibodies, we validated eleven candidate molecules, predicted as either gender-specific or common to both male and female gametocytes. All of them localize to punctuate, vesicle-like structures that relocate to cell periphery upon activation, but only three of them localize to the gametocyte-specific secretory vesicles named osmiophilic bodies. Our results confirm that the egress process involves a tightly coordinated secretory apparatus that includes different types of vesicles and may put the basis for functional studies aimed at designing novel transmission-blocking molecules.

Mallory-Denk-bodies (MDBs) are hepatic protein aggregates associated with inflammation both clinically and in MDB-inducing models. Similar protein aggregation in neurodegenerative diseases also triggers inflammation and NF-B activation. However, the precise mechanism that links protein aggregation to NF-B-activation and inflammatory response remains unclear. Herein we find that treating primary hepatocytes with MDB-inducing agents (N-methylprotoporphyrin (NMPP), protoporphyrin IX (PPIX), or Zinc-protoporphyrin IX (ZnPP)) elicited an IBα-loss with consequent NF-B activation. Four known mechanisms of IBα-loss i.e. the canonical ubiquitin-dependent proteasomal degradation (UPD), autophagic-lysosomal degradation, calpain degradation and translational inhibition, were all probed and excluded. Immunofluorescence analyses of ZnPP-treated cells coupled with 8 M urea/CHAPS-extraction revealed that this IBα-loss was due to its sequestration along with IBβ into insoluble aggregates, thereby releasing NF-B. Through affinity pulldown, proximity biotinylation by antibody recognition, and other proteomic analyses, we verified that NF-B subunit p65, which stably interacts with IBα under normal conditions, no longer binds to it upon ZnPP-treatment. Additionally, we identified 10 proteins that interact with IBα under baseline conditions, aggregate upon ZnPP-treatment, and maintain the interaction with IBα after ZnPP-treatment, either by cosequestering into insoluble aggregates or through a different mechanism. Of these 10 proteins, the nucleoporins Nup153 and Nup358/RanBP2 were identified through RNA-interference, as mediators of IBα-nuclear import. The concurrent aggregation of IBα, NUP153, and RanBP2 upon ZnPP-treatment, synergistically precluded the nuclear entry of IBα and its consequent binding and termination of NF-B activation. This novel mechanism may account for the protein aggregate-induced inflammation observed in liver diseases, thus identifying novel targets for therapeutic intervention. Because of inherent commonalities this MDB cell model is a bona fide protoporphyric model, making these findings equally relevant to the liver inflammation associated with clinical protoporphyria.

At neuronal synapses, activation of group I metabotropic glutamate receptors (mGluR1/5) triggers a form of long-term depression (mGluR-LTD) that relies on new protein synthesis and the internalization of AMPA-type glutamate receptors. Dysregulation of these processes has been implicated in the development of mental disorders such as autism spectrum disorders and therefore merit a better understanding on a molecular level. Here, to study mGluR-induced signaling pathways, we integrated quantitative phosphoproteomics with the analyses of newly synthesized proteins via bio-orthogonal amino acids (azidohomoalanine) in a pulsed labeling strategy in cultured hippocampal neurons stimulated with DHPG, a specific agonist for group I mGluRs. We identified several kinases with important roles in DHPG-induced mGluR activation, which we confirmed using small molecule kinase inhibitors. Furthermore, changes in the AMPA receptor endocytosis pathway in both protein synthesis and protein phosphorylation were identified, whereby Intersectin-1 was validated as a novel player in this pathway. This study revealed several new insights into the molecular pathways downstream of group I mGluR activation in hippocampal neurons, and provides a rich resource for further analyses.

Huanglongbing (HLB) is the most devastating and widespread citrus disease. All commercial citrus varieties are susceptible to the HLB-associated bacterium, Candidatus Liberibacter asiaticus (CLas), which resides in the phloem. The phloem is part of the plant vascular system and is involved in sugar transport. To investigate the plant response to CLas, we enriched for proteins surrounding the phloem in an HLB susceptible sweet orange variety, Washington navel (Citrus sinensis (L) Osbeck). Quantitative proteomics revealed global changes in the citrus proteome after CLas inoculation. Plant metabolism and translation were suppressed, whereas defense-related proteins such as peroxidases, proteases and protease inhibitors were induced in the vasculature. Transcript accumulation and enzymatic activity of plant peroxidases in CLas infected sweet orange varieties under greenhouse and field conditions were assessed. Although peroxidase transcript accumulation was induced in CLas infected sweet orange varieties, peroxidase enzymatic activity varied. Specific serine proteases were up-regulated in Washington navel in the presence of CLas based on quantitative proteomics. Subsequent activity-based protein profiling revealed increased activity of two serine proteases, and reduced activity of one protease in two C. sinensis sweet orange varieties under greenhouse and field conditions. The observations in the current study highlight global reprogramming of the citrus vascular proteome and differential regulation of enzyme classes in response to CLas infection. These results open an avenue for further investigation of diverse responses to HLB across different environmental conditions and citrus genotypes.

Stroke remains a leading cause of death and disability worldwide. Despite continuous advances, the identification of key molecular signatures in the hyper-acute phase of ischemic stroke is still a primary interest for translational research on stroke diagnosis, prognosis, and treatment. Data integration from high-throughput -omics techniques has become crucial to unraveling key interactions among different molecular elements in complex biological contexts, such as ischemic stroke. Thus, we used advanced data integration methods for a multi-level joint analysis of transcriptomics and proteomics data sets obtained from mouse brains at 2 h after cerebral ischemia. By modeling net-like correlation structures, we identified an integrated network of genes and proteins that are differentially expressed at a very early stage after stroke. We validated 10 of these deregulated elements in acute stroke, and changes in their expression pattern over time after cerebral ischemia were described. Of these, CLDN20, GADD45G, RGS2, BAG5, and CTNND2 were next evaluated as blood biomarkers of cerebral ischemia in mice and human blood samples, which were obtained from stroke patients and patients presenting stroke-mimicking conditions. Our findings indicate that CTNND2 levels in blood might potentially be useful for distinguishing ischemic strokes from stroke-mimicking conditions in the hyper-acute phase of the disease. Furthermore, circulating GADD45G content within the first 6 h after stroke could also play a key role in predicting poor outcomes in stroke patients. For the first time, we have used an integrative biostatistical approach to elucidate key molecules in the initial stages of stroke pathophysiology and highlight new notable molecules that might be further considered as blood biomarkers of ischemic stroke.

Middle East and North Africa

Research on the Middle East and North Africa region focuses on changes to politics and society, economics, and security issues.

This is a turbulent period for the region following the Arab Spring, with conflict in Syria continuing to impact its neighbours, governance in Libya under scrutiny, and increasing pressures on the Gulf monarchies, especially around human rights.   

Key research areas include the Gulf States and the Gulf Cooperation Council (GCC), the future of the state, mapping the region’s war economies, the Yemen conflict, Iraq’s reconstruction, and the influence of Saudi Arabia and Iran.

Colorectal cancer (CRC) arises as the consequence of progressive changes from normal epithelial cells through polyp to tumor, and thus is an useful model for studying metabolic shift. In the present study, we studied the metabolomic profiles using high analyte specific gas chromatography/mass spectrometry (GC/MS) and liquid chromatography tandem mass spectrometry (LC/MS/MS) to attain a systems-level view of the shift in metabolism in cells progressing along the path to CRC. Colonic tissues including tumor, polyps and adjacent matched normal mucosa from 26 patients with sporadic CRC from freshly isolated resections were used for this study. The metabolic profiles were obtained using GC/MS and LC/MS/MS. Our data suggest there was a distinct profile change of a wide range of metabolites from mucosa to tumor tissues. Various amino acids and lipids in the polyps and tumors were elevated, suggesting higher energy needs for increased cellular proliferation. In contrast, significant depletion of glucose and inositol in polyps revealed that glycolysis may be critical in early tumorigenesis. In addition, the accumulation of hypoxanthine and xanthine, and the decrease of uric acid concentration, suggest that the purine biosynthesis pathway could have been substituted by the salvage pathway in CRC. Further, there was a step-wise reduction of deoxycholic acid concentration from mucosa to tumors. It appears that to gain a growth advantage, cancer cells may adopt alternate metabolic pathways in tumorigenesis and this flexibility allows them to adapt and thrive in harsh environment.

Oxidative stress has been implicated in aging and many human diseases, notably neurodegenerative disorders and various cancers. The reactive oxygen species that are generated by aerobic metabolism and environmental stressors can chemically modify proteins and alter their biological functions. Cells possess protein repair pathways to rescue oxidized proteins and restore their functions. If these repair processes fail, oxidized proteins may become cytotoxic. Cell homeostasis and viability are therefore dependent on the removal of oxidatively damaged proteins. Numerous studies have demonstrated that the proteasome plays a pivotal role in the selective recognition and degradation of oxidized proteins. Despite extensive research, oxidative stress-triggered regulation of proteasome complexes remains poorly defined. Better understanding of molecular mechanisms underlying proteasome function in response to oxidative stress will provide a basis for developing new strategies aimed at improving cell viability and recovery as well as attenuating oxidation-induced cytotoxicity associated with aging and disease. Here we highlight recent advances in the understanding of proteasome structure and function during oxidative stress and describe how cells cope with oxidative stress through proteasome-dependent degradation pathways.

Cover art by Julie Newdoll for MCP April issue.

As the capability of mass spectrometry-based proteomics has matured, tens of thousands of peptides can be measured simultaneously, which has the benefit of offering a systems view of protein expression. However, a major challenge is that with an increase in throughput, protein quantification estimation from the native measured peptides has become a computational task. A limitation to existing computationally-driven protein quantification methods is that most ignore protein variation, such as alternate splicing of the RNA transcript and post-translational modifications or other possible proteoforms, which will affect a significant fraction of the proteome. The consequence of this assumption is that statistical inference at the protein level, and consequently downstream analyses, such as network and pathway modeling, have only limited power for biomarker discovery. Here, we describe a Bayesian model (BP-Quant) that uses statistically derived peptides signatures to identify peptides that are outside the dominant pattern, or the existence of multiple over-expressed patterns to improve relative protein abundance estimates. It is a research-driven approach that utilizes the objectives of the experiment, defined in the context of a standard statistical hypothesis, to identify a set of peptides exhibiting similar statistical behavior relating to a protein. This approach infers that changes in relative protein abundance can be used as a surrogate for changes in function, without necessarily taking into account the effect of differential post-translational modifications, processing, or splicing in altering protein function. We verify the approach using a dilution study from mouse plasma samples and demonstrate that BP-Quant achieves similar accuracy as the current state-of-the-art methods at proteoform identification with significantly better specificity. BP-Quant is available as a MatLab ® and R packages at https://github.com/PNNL-Comp-Mass-Spec/BP-Quant.

Protein tyrosine kinases (PTKs) play key roles in cellular signal transduction, cell cycle regulation, cell division, and cell differentiation. Dysregulation of PTK-activated pathways, often by receptor overexpression, gene amplification, or genetic mutation, is a causal factor underlying numerous cancers. In this study, we have developed a parallel reaction monitoring (PRM)-based assay for quantitative profiling of 83 PTKs. The assay detects 308 proteotypic peptides from 54 receptor tyrosine kinases and 29 nonreceptor tyrosine kinases in a single run. Quantitative comparisons were based on the labeled reference peptide method. We implemented the assay in four cell models: 1) a comparison of proliferating versus epidermal growth factor (EGF)-stimulated A431 cells, 2) a comparison of SW480Null (mutant APC) and SW480APC (APC restored) colon tumor cell lines, and 3) a comparison of 10 colorectal cancer cell lines with different genomic abnormalities, and 4) lung cancer cell lines with either susceptibility (11-18) or acquired resistance (11-18R) to the epidermal growth factor receptor tyrosine kinase inhibitor erlotinib. We observed distinct PTK expression changes that were induced by stimuli, genomic features or drug resistance, which were consistent with previous reports. However, most of the measured expression differences were novel observations. For example, acquired resistance to erlotinib in the 11-18 cell model was associated not only with previously reported upregulation of MET, but also with upregulation of FLK2 and downregulation of LYN and PTK7. Immunoblot analyses and shotgun proteomics data were highly consistent with PRM data. Multiplexed PRM assays provide a targeted, systems-level profiling approach to evaluate cancer-related proteotypes and adaptations. Data are available through Proteome eXchange Accession PXD002706.

This article has been withdrawn by the authors. We discovered an error after this manuscript was published as a Paper in Press. Specifically, we learned that the structures of glycans presented for the PD-L1 peptide were drawn and labeled incorrectly. We wish to withdraw this article and submit a corrected version for review.

The histidine kinase Hik33 plays important roles in mediating cyanobacterial response to divergent types of abiotic stresses including cold, salt, high light (HL), and osmotic stresses. However, how these functions are regulated by Hik33 remains to be addressed. Using a hik33-deficient strain (hik33) of Synechocystis sp. PCC 6803 (Synechocystis) and quantitative proteomics, we found that Hik33 depletion induces differential protein expression highly similar to that induced by divergent types of stresses. This typically includes downregulation of proteins in photosynthesis and carbon assimilation that are necessary for cell propagation, and upregulation of heat shock proteins, chaperons, and proteases that are important for cell survival. This observation indicates that depletion of Hik33 alone mimics divergent types of abiotic stresses, and that Hik33 could be important for preventing abnormal stress response in the normal condition. Moreover, we found the majority of proteins of plasmid origin were significantly upregulated in hik33, though their biological significance remains to be addressed. Together, the systematically characterized Hik33-regulated cyanobacterial proteome, which is largely involved in stress responses, builds the molecular basis for Hik33 as a general regulator of stress responses.

This article has been withdrawn by the authors. This article did not comply with the editorial guidelines of MCP. Specifically, single peptide based protein identifications of 9-19% were included in the analysis and discussed in the results and conclusions. We wish to withdraw this article and resubmit a clarified, corrected manuscript for review.

Glycopeptides in peptide or digested protein samples pose a number of analytical and bioinformatics challenges beyond those posed by unmodified peptides or peptides with smaller posttranslational modifications. Exact structural elucidation of glycans is generally beyond the capability of a single mass spectrometry experiment, so a reasonable level of identification for tandem mass spectrometry, taken by several glycopeptide software tools, is that of peptide sequence and glycan composition, meaning the number of monosaccharides of each distinct mass, for example HexNAc(2)Hex(5) rather than man5. Even at this level, however, glycopeptide analysis poses challenges:  finding glycopeptide spectra when they are a tiny fraction of the total spectra; assigning spectra with unanticipated glycans, not in the initial glycan database; and finding, scoring, and labeling diagnostic peaks in tandem mass spectra.  Here we discuss recent improvements to Byonic, a glycoproteomics search program, that address these three issues. Byonic now supports filtering spectra by m/z peaks, so that the user can limit attention to spectra with diagnostic peaks, for example, at least two out of three of 204.087 for HexNAc, 274.092 for NeuAc (with water loss), and 366.139 for HexNAc-Hex, all within a set mass tolerance, for example, ± 0.01 Daltons. Also new is glycan "wildcard" search, which allows an unspecified mass within a user-set mass range to be applied to N- or O-linked glycans and enables assignment of spectra with unanticipated glycans. Finally the next release of Byonic supports user-specified peak annotations from user-defined posttranslational modifications. We demonstrate the utility of these new software features by finding previously unrecognized glycopeptides in publicly available data, including glycosylated neuropeptides from rat brain.

Sparkling wine is an alcoholic beverage enjoyed around the world. The sensory properties of sparkling wine depend on a complex interplay between the chemical and biochemical components in the final product. Glycoproteins have been linked to positive and negative qualities in sparkling wine, but the glycosylation profiles of sparkling wine have not been previously investigated in detail. We analysed the glyco/proteome of sparkling wines using protein- and glycopeptide-centric approaches. We developed an automated workflow that created ion libraries to analyse Sequential Window Acquisition of all THeoretical mass spectra (SWATH) Data Independent Acquisition (DIA) mass spectrometry data based on glycopeptides identified by Byonic. We applied our workflow to three pairs of experimental sparkling wines to assess the effects of aging on lees and of different yeast strains used in the Liqueur de Tirage for secondary fermentation. We found that aging a cuvée on lees for 24 months compared to 8 months led to a dramatic decrease in overall protein abundance and an enrichment in large glycans at specific sites in some proteins. Secondary fermentation of a Riesling wine with Saccharomyces cerevisiae yeast strain Siha4 produced more yeast proteins and glycoproteins than with S. cerevisiae yeast strain DV10. The abundance and glycosylation profiles of grape glycoproteins were also different between grape varieties. This work represents the first in-depth study into protein- and peptide-specific glycosylation in sparkling wines and describes a quantitative glycoproteomic SWATH/DIA workflow that is broadly applicable to other sample types.

Many cell surface and secreted proteins are modified by the covalent addition of glycans that play an important role in the development of multicellular organisms. These glycan modifications enable communication between cells and the extracellular matrix via interactions with specific glycan-binding lectins and the regulation of receptor-mediated signaling. Aberrant protein glycosylation has been associated with the development of several muscular diseases suggesting essential glycan- and lectin-mediated functions in myogenesis and muscle development but our molecular understanding of the precise glycans, catalytic enzymes and lectins involved remain only partially understood. Here, we quantified dynamic remodeling of the membrane-associated proteome during a time-course of myogenesis in cell culture. We observed wide-spread changes in the abundance of several important lectins and enzymes facilitating glycan biosynthesis. Glycomics-based quantification of released N-linked glycans confirmed remodeling of the glycome consistent with the regulation of glycosyltransferases and glycosidases responsible for their formation including a previously unknown di-galactose-to-sialic acid switch supporting a functional role of these glycoepitopes in myogenesis. Furthermore, dynamic quantitative glycoproteomic analysis with multiplexed stable isotope labelling and analysis of enriched glycopeptides with multiple fragmentation approaches identified glycoproteins modified by these regulated glycans including several integrins and growth factor receptors. Myogenesis was also associated with the regulation of several lectins most notably the up-regulation of galectin-1 (LGALS1). CRISPR/Cas9-mediated deletion of Lgals1 inhibited differentiation and myotube formation suggesting an early functional role of galectin-1 in the myogenic program. Importantly, similar changes in N-glycosylation and the up-regulation of galectin-1 during postnatal skeletal muscle development were observed in mice. Treatment of new-born mice with recombinant adeno-associated viruses to overexpress galectin-1 in the musculature resulted in enhanced muscle mass. Our data form a valuable resource to further understand the glycobiology of myogenesis and will aid the development of intervention strategies to promote healthy muscle development or regeneration.

CD4+ T cell responses are crucial for inducing and maintaining effective anti-cancer immunity, and the identification of human leukocyte antigen class II (HLA-II) cancer-specific epitopes is key to the development of potent cancer immunotherapies. In many tumor types, and especially in glioblastoma (GBM), HLA-II complexes are hardly ever naturally expressed. Hence, little is known about immunogenic HLA-II epitopes in GBM. With stable expression of the class II major histocompatibility complex transactivator (CIITA) coupled to a detailed and sensitive mass spectrometry based immunopeptidomics analysis, we here uncovered a remarkable breadth of the HLA-ligandome in HROG02, HROG17 and RA GBM cell lines. The effect of CIITA expression on the induction of the HLA-II presentation machinery was striking in each of the three cell lines, and it was significantly higher compared to interferon gamma (IFN) treatment. In total, we identified 16,123 unique HLA-I peptides and 32,690 unique HLA-II peptides. In order to genuinely define the identified peptides as true HLA ligands, we carefully characterized their association with the different HLA allotypes. In addition, we identified 138 and 279 HLA-I and HLA-II ligands, respectively, most of which are novel in GBM, derived from known GBM-associated tumor-antigens that have been used as source proteins for a variety of GBM vaccines. Our data further indicate that CIITA-expressing GBM cells acquired an antigen presenting cell-like phenotype as we found that they directly present external proteins as HLA-II ligands. Not only that CIITA-expressing GBM cells are attractive models for antigen discovery endeavors, but also such engineered cells have great therapeutic potential through massive presentation of a diverse antigenic repertoire.

The glycoprotein spike (S) on the surface of SARS-CoV-2 is a determinant for viral invasion and host immune response. Herein, we characterized the site-specific N-glycosylation of S protein at the level of intact glycopeptides. All 22 potential N-glycosites were identified in the S-protein protomer and were found to be preserved among the 753 SARS-CoV-2 genome sequences. The glycosites exhibited glycoform heterogeneity as expected for a human cell-expressed protein subunit. We identified masses that correspond to 157 N-glycans, primarily of the complex type. In contrast, the insect cell-expressed S protein contained 38 N-glycans, completely of the high-mannose type. Our results revealed that the glycan types were highly determined by the differential processing of N-glycans among human and insect cells, regardless of the glycosites’ location. Moreover, the N-glycan compositions were conserved among different sizes of subunits. Our study indicate that the S protein N-glycosylation occurs regularly at each site, albeit the occupied N-glycans were diverse and heterogenous. This N-glycosylation landscape and the differential N-glycan patterns among distinct host cells are expected to shed light on the infection mechanism and present a positive view for the development of vaccines and targeted drugs.

Esophageal squamous cell cancer (ESCC) is an aggressive malignancy with poor therapeutic outcomes. However, the alterations in proteins and post-translational modifications (PTMs) leading to the pathogenesis of ESCC remains unclear. Here, we provide the comprehensive characterization of the proteome, phosphorylome, lysine acetylome and succinylome for ESCC and matched control cells using quantitative proteomic approach. We identify abnormal protein and post-translational modification (PTM) pathways, including significantly downregulated lysine succinylation sites in cancer cells. Focusing on hyposuccinylation, we reveal that this altered PTM was enriched on enzymes of metabolic pathways inextricably linked with cancer metabolism. Importantly, ESCC malignant behaviors such as cell migration are inhibited once the level of succinylation was restored in vitro or in vivo. This effect was further verified by mutations to disrupt succinylation sites in candidate proteins. Meanwhile, we found that succinylation has a negative regulatory effect on histone methylation to promote cancer migration. Finally, hyposuccinylation is confirmed in primary ESCC specimens. Our findings together demonstrate that lysine succinylation may alter ESCC metabolism and migration, providing new insights into the functional significance of PTM in cancer biology.

The choice for adjuvant chemotherapy in stage II colorectal cancer (CRC) is controversial as many patients are cured by surgery alone and it is difficult to identify patients with high-risk of recurrence of the disease. There is a need for better stratification of this group of patients. Mass spectrometry imaging could identify patients at risk. We report here the N-glycosylation signatures of the different cell populations in a group of stage II CRC tissue samples. The cancer cells, compared to normal epithelial cells, have increased levels of sialylation and high-mannose glycans, as well as decreased levels of fucosylation and highly branched N-glycans. When looking at the interface between cancer and its microenvironment, it seems that the cancer N-glycosylation signature spreads into the surrounding stroma at the invasive front of the tumor. This finding was more outspoken in patients with a worse outcome within this sample group.

Paramphistomosis, caused by the rumen fluke, Calicophoron daubneyi, is a parasitic infection of ruminant livestock which has seen a rapid rise in prevalence throughout Western Europe in recent years. Following ingestion of metacercariae (parasite cysts) by the mammalian host, newly-excysted juveniles (NEJs) emerge and invade the duodenal submucosa which causes significant pathology in heavy infections. The immature larvae then migrate upwards, along the gastrointestinal tract, and enter the rumen where they mature and begin to produce eggs. Despite their emergence, and sporadic outbreaks of acute disease, we know little about the molecular mechanisms used by C. daubneyi to establish infection, acquire nutrients and to avoid the host immune response. Here, transcriptome analysis of four intra-mammalian life-cycle stages, integrated with secretome analysis of the NEJ and adult parasites (responsible for acute and chronic disease respectively), revealed how the expression and secretion of selected families of virulence factors and immunomodulators are regulated in accordance with fluke development and migration. Our data show that whilst a family of cathepsins B with varying S2 sub-site residues (indicating distinct substrate specificities) are differentially secreted by NEJs and adult flukes, cathepsins L and F are secreted in low abundance by NEJs only. We found that C. daubneyi has an expanded family of aspartic peptidases, which is up-regulated in adult worms, although they are underrepresented in the secretome. The most abundant proteins in adult fluke secretions were helminth defence molecules (HDMs) that likely establish an immune environment permissive to fluke survival and/or neutralise pathogen-associated molecular patterns (PAMPs) such as bacterial lipopolysaccharide in the microbiome-rich rumen. The distinct collection of molecules secreted by C. daubneyi allowed the development of the first coproantigen-based ELISA for paramphistomosis which, importantly, did not recognise antigens from other helminths commonly found as co-infections with rumen fluke.

Antibodies play essential roles in both diagnostics and therapeutics. Epitope mapping is essential to understand how an antibody works and to protect intellectual property. Given the millions of antibodies for which epitope information is lacking, there is a need for high-throughput epitope mapping. To address this, we developed a strategy, Antibody binding epitope Mapping (AbMap), by combining a phage displayed peptide library with next generation sequencing. Using AbMap, profiles of the peptides bound by 202 antibodies were determined in a single test, and linear epitopes were identified for >50% of the antibodies. Using spike protein (S1 and S2)-enriched antibodies from the convalescent serum of one COVID-19 patient as the input, both linear and conformational epitopes of spike protein specific antibodies were identified. We defined peptide-binding profile of an antibody as the Binding Capacity (BiC). Conceptually, the BiC could serve as a systematic and functional descriptor of any antibody. Requiring at least one order of magnitude less time and money to map linear epitopes than traditional technologies, AbMap allows for high-throughput epitope mapping and creates many possibilities.

Thermal proteome profiling (TPP) allows for the unbiased detection of drug – target protein engagements in vivo. Traditionally, one cell type is used for TPP studies, with the risk of missing important differentially expressed target proteins. The use of whole organisms would circumvent this problem. Zebrafish embryos are amenable to such an approach. Here, we used TPP on whole zebrafish embryo lysate to identify protein targets of napabucasin, a compound that may affect Signal transducer and activator of transcription 3 (Stat3) signaling through an ill-understood mechanism. In zebrafish embryos, napabucasin induced developmental defects consistent with inhibition of Stat3 signaling. TPP profiling showed no distinct shift in Stat3 upon napabucasin treatment, but effects were detected on the oxidoreductase, Pora, which might explain effects on Stat3 signaling. Interestingly, thermal stability of several aldehyde dehydrogenases (Aldhs) was affected. Moreover, napabucasin activated ALDH enzymatic activity in vitro. Aldhs have crucial roles in retinoic acid metabolism and functionally we validated napabucasin-mediated activation of the retinoic acid pathway in zebrafish in vivo. We conclude that TPP profiling in whole zebrafish embryo lysate is feasible and facilitates direct correlation of in vivo effects of small molecule drugs with their protein targets.

The increasing consumption of high-fat foods combined with a lack of exercise is a major contributor to the burden of obesity in humans. Aerobic exercise such as running is known to provide metabolic benefits, but how the over-consumption of a high fat diet (HFD) and exercise interact is not well characterized at the molecular level. Here, we examined the plasma proteome in mice for the effects of aerobic exercise as both a treatment and as a preventative regime for animals on either HFD or a healthy control diet. This analysis detected large changes in the plasma proteome induced by the HFD, such as increased abundance of SERPINA7, ALDOB, and down-regulation of SERPINA1E, CFD (adipsin). Some of these changes were significantly reverted using exercise as a preventative measure, but not as a treatment regime. To determine if either the intensity, or duration, of exercise influenced the outcome, we compared high-intensity interval training (HIIT) and endurance running. Endurance running slightly out-performed HIIT exercise, but overall, both provided similar reversion in abundance of plasma proteins modulated by the high-fat diet including SERPINA7, APOE, SERPINA1E, and CFD. Finally, we compared the changes induced by over-consumption of HFD to previous data from mice fed an isocaloric high saturated fat (SFA) or polyunsaturated fat (PUFA) diet. This identified several common changes including increased APOC2 and APOE, but also highlighted changes specific for either over-consumption of HFD (ALDOB, SERPINA7, CFD), SFA-based diets (SERPINA1E), or PUFA-based diets (Haptoglobin - Hp). Together, these data highlight the importance of early intervention with exercise to revert HFD-induced phenotypes and suggest some of the molecular mechanisms leading to the changes in the plasma proteome generated by high fat diet consumption. Web-based interactive visualizations are provided for this dataset (larancelab.com/hfd-exercise), which give insight into diet and exercise phenotypic interactions on the plasma proteome.

The complexity and dynamics of the immensely heterogeneous glycoproteome of the prostate cancer (PCa) tumour micro-environment remain incompletely mapped, a knowledge gap that impedes our molecular-level understanding of the disease. To this end, we have used sensitive glycomics and glycoproteomics to map the protein-, cell- and tumour grade-specific N- and O-glycosylation in surgically-removed PCa tissues spanning five histological grades (n = 10/grade) and tissues from patients with benign prostatic hyperplasia (n = 5). Quantitative glycomics revealed PCa grade-specific alterations of the oligomannosidic-, paucimannosidic- and branched sialylated complex-type N-glycans, and dynamic remodelling of the sialylated core 1- and core 2-type O-glycome. Deep quantitative glycoproteomics identified ~7,400 unique N-glycopeptides from 500 N-glycoproteins and ~500 unique O-glycopeptides from nearly 200 O-glycoproteins. With reference to a recent Tissue and Blood Atlas, our data indicate that paucimannosidic glycans of the PCa tissues arise mainly from immune cell-derived glycoproteins. Further, the grade-specific PCa glycosylation arises primarily from dynamics in the cellular makeup of the PCa tumour microenvironment across grades involving increased oligomannosylation of prostate-derived glycoproteins and decreased bisecting GlcNAcylation of N-glycans carried by the extracellular matrix proteins. Further, elevated expression of several oligosaccharyltransferase subunits and enhanced N-glycoprotein site occupancy were observed associated with PCa progression. Finally, correlations between the protein-specific glycosylation and PCa progression were observed including increased site-specific core 2-type O-glycosylation of collagen VI. In conclusion, integrated glycomics and glycoproteomics have enabled new insight into the complexity and dynamics of the tissue glycoproteome associated with PCa progression generating an important resource to explore the underpinning disease mechanisms.

Malaria elimination is still pending on the development of novel tools that rely on a deep understanding of parasite biology. Proteins of all living cells undergo a myriad number of posttranslational modifications (PTMs) that are critical to multifarious life processes. An extensive proteome-wide dissection revealed a fine PTM map of most proteins in both Plasmodium falciparum, the causative agent of severe malaria, and the infected red blood cells. More than two-thirds of proteins of the parasite and its host cell underwent extensive and dynamic modification throughout the erythrocytic developmental stage. PTMs critically modulate the virulence factors involved in the host-parasite interaction and pathogenesis. Furthermore, P. falciparum stabilized the supporting proteins of erythrocyte origin by selective de-modification. Collectively, our multiple omic analyses, apart from having furthered a deep understanding of the systems biology of P. falciparum and malaria pathogenesis, provide a valuable resource for mining new antimalarial targets.

Biological functions emerge from complex and dynamic networks of protein-protein interactions. Because these protein-protein interaction networks, or interactomes, represent pairwise connections within a hierarchically organized system, it is often useful to identify higher-order associations embedded within them, such as multi-member protein complexes. Graph-based clustering techniques are widely used to accomplish this goal, and dozens of field-specific and general clustering algorithms exist. However, interactomes can be prone to errors, especially when inferred from high-throughput biochemical assays. Therefore, robustness to network-level noise is an important criterion for any clustering algorithm that aims to generate robust, reproducible clusters. Here, we tested the robustness of a range of graph-based clustering algorithms in the presence of noise, including algorithms common across domains and those specific to protein networks. Strikingly, we found that all of the clustering algorithms tested here markedly amplified noise within the underlying protein interaction network. Randomly rewiring only 1% of network edges yielded more than a 50% change in clustering results, indicating that clustering markedly amplified network-level noise. Moreover, we found the impact of network noise on individual clusters was not uniform: some clusters were consistently robust to injected noise while others were not. To assist in assessing this, we developed the clust.perturb R package and Shiny web application to measure the reproducibility of clusters by randomly perturbing the network. We show that clust.perturb results are predictive of real-world cluster stability: poorly reproducible clusters as identified by clust.perturb are significantly less likely to be reclustered across experiments. We conclude that graph-based clustering amplifies noise in protein interaction networks, but quantifying the robustness of a cluster to network noise can separate stable protein complexes from spurious associations.

Protease activity has been associated with pathological processes that can lead to cancer development and progression. However, understanding the pathological unbalance in proteolysis is challenging since changes can occur simultaneously at protease, their inhibitor and substrate levels. Here, we present a pipeline that combines peptidomics, proteomics and peptidase predictions for studying proteolytic events in the saliva of seventy-nine patients and their association with oral squamous cell carcinoma (OSCC) prognosis. Our findings revealed differences in the saliva peptidome of patients with (pN+) or without (pN0) lymph node metastasis and delivered a panel of ten endogenous peptides correlated with poor prognostic factors plus five molecules able to classify pN0 and pN+ patients (ROC-AUC>0.85). In addition, endo- and exopeptidases putatively implicated in the processing of differential peptides were investigated using cancer tissue gene expression data from publicly repositories reinforcing their association with poorer survival rates and prognosis in oral cancer. The dynamics of the OSCC-related proteolysis was further explored via the proteomic profiling of saliva. This revealed that peptidase/endopeptidase inhibitors exhibited reduced levels in the saliva of pN+ patients, as confirmed by SRM-MS, whilst minor changes were detected in the level of saliva proteases. Taken together, our results indicated that proteolytic activity is accentuated in the saliva of OSCC patients with lymph node metastasis and, at least in part, this is modulated by reduced levels of salivary peptidase inhibitors. Therefore, this integrated pipeline provided better comprehension and discovery of molecular features with implications in the oral cancer metastasis prognosis.

Modulation of the host cell is integral to the survival and replication of microbial pathogens. Several intracellular bacterial pathogens deliver bacterial proteins, termed ‘effector proteins’ into the host cell during infection by sophisticated protein translocation systems, which manipulate cellular processes and functions. The functional contribution of individual effectors is poorly characterised, particularly in intracellular bacterial pathogens with large effector protein repertoires. Technical caveats have limited the capacity to study these proteins during a native infection, with many effector proteins having only been demonstrated to be translocated during over-expression of tagged versions. Here we developed a novel strategy to examine effector proteins in the context of infection. We coupled a broad, unbiased proteomics-based screen with organelle purification to study the host-pathogen interactions occurring between the host cell mitochondrion and the Gram-negative, Q fever pathogen Coxiella burnetii. We identify 4 novel mitochondrially-targeted C. burnetii effector proteins, renamed Mitochondrial Coxiella effector protein (Mce) B to E. Examination of the subcellular localisation of ectopically expressed proteins confirmed their mitochondrial localisation, demonstrating the robustness of our approach. Subsequent biochemical analysis and affinity enrichment proteomics of one of these effector proteins, MceC, revealed the protein localises to the inner membrane and can interact with components of the mitochondrial quality control machinery. Our study adapts high-sensitivity proteomics to study intracellular host-pathogen interactions, providing a robust strategy to examine the sub-cellular localisation of effector proteins during native infection. This approach could be applied to a range of pathogens and host cell compartments to provide a rich map of effector dynamics throughout infection.

Hirschsprung disease (HSCR) is a heterogeneous group of neurocristopathy characterized by the absence of the enteric ganglia along a variable length of the intestine. Genetic defects play a major role in the pathogenesis of HSCR while family studies of pathogenic variants in all the known genes (loci) only demonstrate incomplete penetrance and variable expressivity for unknown reasons. Here, we applied large-scale, quantitative proteomics of human colon tissues from 21 patients using iTRAQ method followed by bioinformatics analysis. Selected findings were confirmed by parallel reaction monitoring (PRM) verification. At last the interesting differentially expressed proteins were confirmed by western blot. A total of 5341 proteins in human colon tissues were identified. Among them, 664 proteins with >1.2-fold difference were identified in 6 groups: groups A1 and A2 pooled protein from the ganglionic and aganglionic colon of male, long-segment HSCR patients (L-HSCR, n=7); groups B1 and B2 pooled protein from the ganglionic and aganglionic colon of male, short-segment HSCR patients (S-HSCR, n=7); and groups C1 and C2 pooled protein from the ganglionic and aganglionic colon of female, S-HSCR patients (n=7). Based on these analyses, 49 proteins from 5 pathways were selected for PRM verification, including ribosome, endocytosis, spliceosome, oxidative phosphorylation and cell adhesion. The downregulation of three neuron projection development genes ARF4, KIF5B and RAB8A in the aganglionic part of the colon were verified in 15 paired colon samples using WB. The findings of this study will shed new light on the pathogenesis of HSCR and facilitate the development of therapeutic targets.

Thyroglobulin (Tg) is a secreted iodoglycoprotein serving as the precursor for T3 and T4 hormones. Many characterized Tg gene mutations produce secretion-defective variants resulting in congenital hypothyroidism (CH). Tg processing and secretion is controlled by extensive interactions with chaperone, trafficking, and degradation factors comprising the secretory proteostasis network. While dependencies on individual proteostasis network components are known, the integration of proteostasis pathways mediating Tg protein quality control and the molecular basis of mutant Tg misprocessing remain poorly understood. We employ a multiplexed quantitative affinity purification–mass spectrometry approach to define the Tg proteostasis interactome and changes between WT and several CH-variants. Mutant Tg processing is associated with common imbalances in proteostasis engagement including increased chaperoning, oxidative folding, and engagement by targeting factors for ER-associated degradation (ERAD). Furthermore, we reveal mutation-specific changes in engagement with N-glycosylation components, suggesting distinct requirements for one Tg variant on dual engagement of both oligosaccharyltransferase complex isoforms for degradation. Modulating dysregulated proteostasis components and pathways may serve as a therapeutic strategy to restore Tg secretion and thyroid hormone biosynthesis.

The Rhizobium-legume symbiosis is a beneficial interaction in which the bacterium converts atmospheric nitrogen into ammonia and delivers it to the plant in exchange for carbon compounds. This symbiosis implies the adaptation of bacteria to live inside host plant cells. In this work we apply RP-LC-MS/MS and  iTRAQ techniques to study the proteomic profile of endosymbiotic cells (bacteroids) induced by Rhizobium leguminosarum bv viciae strain UPM791 in legume nodules. Nitrogenase subunits, tricarboxylic acid cycle enzymes, and stress response proteins are amongst the most abundant from over one thousand rhizobial proteins identified in pea (Pisum sativum) bacteroids. Comparative analysis of bacteroids induced in pea and in lentil (Lens culinaris)nodules revealed the existence of a significant host-specific differential response affecting dozens of bacterial proteins, including stress-related proteins, transcriptional regulators, and proteins involved in the carbon and nitrogen metabolisms. A mutant affected in one of these proteins, homologous to a GntR-like transcriptional regulator, showed a symbiotic performance significantly  impaired in symbiosis with pea, but not with lentil plants. Analysis of the proteomes of bacteroids isolated from both hosts also revealed the presence of different sets of plant-derived nodule-specific cysteine rich (NCR) peptides, indicating that the endosymbiotic bacteria find a host-specific cocktail of chemical stressors inside the nodule. By studying variations of the bacterial response to different plant cell environments we will be able to identify specific limitations imposed by the host that might give us clues for the improvement of rhizobial performance.

The early detection of pancreatic ductal adenocarcinoma is a complex clinical obstacle yet is key to improving the overall likelihood of patient survival. Current and prospective carbohydrate biomarkers CA19-9 and sTRA are sufficient for surveilling disease progression yet are not approved for delineating PDAC from other abdominal cancers and non-cancerous pancreatic pathologies. To further understand these glycan epitopes, an imaging mass spectrometry approach was utilized to assess the N-glycome of the human pancreas and pancreatic cancer in a cohort of PDAC patients represented by tissue microarrays and whole tissue sections. Orthogonally, these same tissues were characterized by multi-round immunofluorescence which defined expression of CA19-9 and sTRA as well as other lectins towards carbohydrate epitopes with the potential to improve PDAC diagnosis. These analyses revealed distinct differences not only in N-glycan spatial localization across both healthy and diseased tissues but importantly between different biomarker-categorized tissue samples. Unique sulfated bi-antennary N-glycans were detected specifically in normal pancreatic islets. N-glycans from CA19-9 expressing tissues tended to be bi-, tri- and tetra-antennary structures with both core and terminal fucose residues and bisecting N-acetylglucosamines. These N-glycans were detected in less abundance in sTRA-expressing tumor tissues, which favored tri- and tetra-antennary structures with polylactosamine extensions. Increased sialylation of N-glycans was detected in all tumor tissues. A candidate new biomarker derived from IMS was further explored by fluorescence staining with selected lectins on the same tissues. The lectins confirmed the expression of the epitopes in cancer cells and revealed different tumor-associated staining patterns between glycans with bisecting GlcNAc and those with terminal GlcNAc. Thus, the combination of lectin-IHC and IMS techniques produces more complete information for tumor classification than the individual analyses alone. These findings potentiate the development of early assessment technologies to rapidly and specifically identify PDAC in the clinic that may directly impact patient outcomes.

Aspergillus flavus (A. flavus), a pathogenic fungus, can produce carcinogenic and toxic aflatoxins that are a serious agricultural and medical threat worldwide. Attempts to decipher the aflatoxin biosynthetic pathway have been hampered by the lack of a high-quality genome annotation for A. flavus. To address this gap, we performed a comprehensive proteogenomic analysis using high-accuracy mass spectrometry data for this pathogen. The resulting high-quality dataset confirmed the translation of 8,724 previously-predicted genes, and identified 732 novel proteins, 269 splice variants, 447 single amino acid variants, 188 revised genes. A subset of novel proteins was experimentally validated by RT-PCR and synthetic peptides. Further functional annotation suggested that a number of the identified novel proteins may play roles in aflatoxin biosynthesis and stress responses in A. flavus. This comprehensive strategy also identified a wide range of post-translational modifications (PTMs), including 3,461 modification sites from 1,765 proteins. Functional analysis suggested the involvement of these modified proteins in the regulation of cellular metabolic and aflatoxin biosynthetic pathways. Together, we provided a high quality annotation of A. flavus genome and revealed novel insights into the mechanisms of aflatoxin production and pathogenicity in this pathogen.

The molecular mechanism associated with mammalian meiosis has yet to be fully explored, and one of the main reasons for this lack of exploration is that some meiosis-essential genes are still unknown. The profiling of gene expression during spermatogenesis has been performed in previous studies, yet few studies have aimed to find new functional genes. Since there is a huge gap between the number of genes that are able to be quantified and the number of genes that can be characterized by phenotype screening in one assay, an efficient method to rank quantified genes according to phenotypic relevance is of great importance. We proposed to rank genes by the probability of their function in mammalian meiosis based on global protein abundance using machine learning. Here, nine types of germ cells focusing on continual substages of meiosis prophase I were isolated, and the corresponding proteomes were quantified by high-resolution mass spectrometry. By combining meiotic labels annotated from the MGI mouse knockout database and the spermatogenesis proteomics dataset, a supervised machine learning package, FuncProFinder, was developed to rank meiosis-essential candidates. Of the candidates whose functions were unannotated, four of ten genes with the top prediction scores, Zcwpw1, Tesmin, 1700102P08Rik and Kctd19, were validated as meiosis-essential genes by knockout mouse models. Therefore,  mammalian meiosis-essential genes could be efficiently predicted based on the protein abundance dataset, which provides a paradigm for other functional gene mining from a related abundance dataset.

The histopathological subtype of lung adenocarcinoma (LUAD) is closely associated with prognosis. Micropapillary or solid predominant LUAD tends to relapse after surgery at an early stage, whereas lepidic pattern shows a favorable outcome. However, the molecular mechanism underlying this phenomenon remains unknown. Here, we recruited 31 lepidic predominant LUADs (LR: low-risk subtype group) and 28 micropapillary or solid predominant LUADs (HR: high-risk subtype group). Tissues of these cases were obtained and label-free quantitative proteomic and bioinformatic analyses were performed. Additionally, prognostic impact of targeted proteins was validated using The Cancer Genome Atlas databases (n=492) and tissue microarrays composed of early-stage LUADs (n=228). A total of 192 differentially expressed proteins were identified between tumor tissues of LR and HR and three clusters were identified via hierarchical clustering excluding eight proteins. Cluster 1 (65 proteins) showed a sequential decrease in expression from normal tissues to tumor tissues of LR and then to HR and was predominantly enriched in pathways such as tyrosine metabolism and ECM-receptor interaction, and increased matched mRNA expression of 18 proteins from this cluster predicted favorable prognosis. Cluster 2 (70 proteins) demonstrated a sequential increase in expression from normal tissues to tumor tissues of LR and then to HR and was mainly enriched in pathways such as extracellular organization, DNA replication and cell cycle, and high matched mRNA expression of 25 proteins indicated poor prognosis. Cluster 3 (49 proteins) showed high expression only in LR, with high matched mRNA expression of 20 proteins in this cluster indicating favorable prognosis. Furthermore, high expression of ERO1A and FEN1 at protein level predicted poor prognosis in early-stage LUAD, supporting the mRNA results. In conclusion, we discovered key differentially expressed proteins and pathways between low-risk and high-risk subtypes of early-stage LUAD. Some of these proteins could serve as potential biomarkers in prognostic evaluation.

Extracellular vesicle (EV) proteins from acute myeloid leukemia (AML) cell lines were analyzed using mass spectrometry. The analyses identified 2450 proteins, including 461 differentially expressed proteins (290 upregulated and 171 downregulated). CD53 and CD47 were upregulated and were selected as candidate biomarkers. The association between survival of patients with AML and the expression levels of CD53 and CD47 at diagnosis was analyzed using mRNA expression data from The Cancer Genome Atlas database. Patients with higher expression levels showed significantly inferior survival than those with lower expression levels. Enzyme-linked immunosorbent assay results of the expression levels of CD53 and CD47 from EVs in the bone marrow of patients with AML at diagnosis and at the time of complete remission with induction chemotherapy revealed that patients with downregulated CD53 and CD47 expression appeared to relapse less frequently. Network model analysis of EV proteins revealed several upregulated kinases, including LYN, CSNK2A1, SYK, CSK, and PTK2B. The potential cytotoxicity of several clinically applicable drugs that inhibit these kinases was tested in AML cell lines. The drugs lowered the viability of AML cells. The collective data suggest that AML-derived EVs could reflect essential leukemia biology.